《数学建模》实验指导书.doc

《《数学建模》实验指导书.doc》由会员分享，可在线阅读，更多相关《《数学建模》实验指导书.doc（16页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、Four short words sum up what has lifted most successful individuals above the crowd: a little bit more.-author-date数学建模实验指导书汇编语言程序设计实验指导书数学与计算机科学学院数学建模实验指导书2011年9月1日-目 录实验一 “商人们安全过河”的MATLAB程序1实验二 初等模型求解2实验三 数学规划模型求解3实验四 微分方程模型求解4实验五 离散模型求解6实验六 统计回归模型的求解7附件：数学建模实验报告9实验一 “商人们安全过河”的MATLAB程序一、实验目的复习Matl

2、ab编程；掌握编写简单的Matlab程序，掌握条件、循环和选择三种语句的用法。二、实验类型:设计三、实验环境计算机、软件Matlab7.0以上的环境四、实验内容1. 建立M-文件: 已知函数 计算 ，并作出该函数的曲线图。 2. 编写利用顺序Guass消去法求方程组解的M-函数文件,并计算方程组的解3. 编写“商人们安全过河”的Matlab程序五、实验总结根据实验操作和实验报告要求，完成实验报告;实验二 初等模型求解一、实验目的学会使用Matlab软件进行一维插值、二维插值运算，会进行多项式拟合、一般非线性拟合。二、实验类型:验证三、实验环境计算机、软件Matlab7.0以上的环境四、实验内容

3、1、 用生成一组数据，并用一维数据插值的方法（插值方法为：三次样条插值）对给出的数据进行曲线拟合，并在图像上显示出拟合效果。2、 假设已知的数据点来自函数，试根据生成的数据用5次多项式拟合的方法拟合函数曲线，并画出图形。3、 下表中给出的数据满足原型，试用最小二乘法求出， 的值，并用得出的函数将函数曲线绘制出来，观察拟合效果。（假设已知数据已读入，并存储在变量x和y中了）X-2-1.7-1.4-1.1-0.8-0.5-0.20.10.40.711.3Y0.102890.117410.131580.144830.156560.166220.173320.17750.178530.176350.1

4、71090.16302X1.61.92.22.52.83.13.43.744.34.64.9Y0.152550.14020.126550.112190.097680.083530.070150.057860.046870.037290.029140.02236五、实验总结根据实验操作和实验报告要求，做好实验报告。 实验三 数学规划模型求解一、实验目的学会根据实际问题建立线性规划模型，求解线性极值问题，掌握用Matlab、Lindo软件求解线性规划问题。二、实验类型:设计三、实验环境计算机、软件Matlab7.0、Lindo5.0以上的环境四、实验内容1、 求解线性规划问题：2、 某车间有甲、乙

5、两台机床，可用于加工三种工件。假定这两台车床的可用台时数分别为800和900，三种工件的数量分别为400、600和500，且已知用三种不同车床加工单位数量不同工件所需的台时数和加工费用如下表。问怎样分配车床的加工任务，才能既满足加工工件的要求，又使加工费用最低？车床类 型单位工件所需加工台时数单位工件的加工费用可用台时数工件1工件2工件3工件1工件2工件3甲0.41.11.013910800乙0.51.21.3111289003、 某工厂生产每件产品需经A，B，C三个车间，每个车间所需的工时数如下表所示，已知生产单位甲产品工厂可获利4万元，生产单位乙产品工厂可获利3万元，问该厂如何安排生产才能

6、使每周获得的利润最大？车间ABC生产单位甲产品需工时数210生产单位乙产品需工时数111一周可用工时数1087五、实验总结根据实验操作和实验报告要求，做好实验报告。 实验四 微分方程模型求解一、实验目的1掌握用Matlab解微分方程的方法；2理解微分方程的数值解原理。二、实验类型:验证三、实验环境多媒体计算机、WINDOWS XP系统、Matlab软件（7.0以上版本）四、实验内容（一）、示例一个慢跑者在平面上沿椭圆以恒定的速率v=1跑步,设椭圆方程为: x=10+20cos t, y=20+5sin t突然有一只狗攻击他 这只狗从原点出发,以恒定速率w跑向慢跑者,狗的运动方向始终指向慢跑者分

7、别求出w=20,w=5时狗的运动轨迹1 模型建立设t 时刻慢跑者的坐标为(X(t)，Y(t)，狗的坐标为(x(t)，y(t)则 X=10+20cos t, Y=20+15sin t. 狗从(0,0)出发, 与导弹追踪问题类似，狗的运动轨迹的参数方程为:2 模型求解(1) w=20时,建立文件eq3m如下: function dy=eq3(t,y) dy=zeros(2,1); dy(1)=20*(10+20*cos(t)-y(1)/sqrt (10+20*cos(t)-y(1)2+(20+15*sin(t)-y(2)2); dy(2)=20*(20+15*sin(t)-y(2)/sqrt (1

8、0+20*cos(t)-y(1)2+(20+15*sin(t)-y(2)2);取t0=0，tf=10，建立主程序chase3m如下： t0=0;tf=10; t,y=ode45(eq3,t0 tf,0 0); T=0:01:2*pi; X=10+20*cos(T); Y=20+15*sin(T); plot(X,Y,-) hold on plot(y(:,1),y(:,2),*)在chase3m中，不断修改tf的值,分别取tf=5, 25, 35,至315时,狗刚好追上慢跑者(2) w=5时建立M文件eq4m如下: function dy=eq4(t,y) dy=zeros(2,1); dy(

9、1)=5*(10+20*cos(t)-y(1)/sqrt (10+20*cos(t)-y(1)2+(20+15*sin(t)-y(2)2); dy(2)=5*(20+15*sin(t)-y(2)/sqrt (10+20*cos(t)- y(1)2+(20+15*sin(t)-y(2)2);取t0=0，tf=10，建立主程序chase4m如下： t0=0;tf=10; t,y=ode45(eq4,t0 tf,0 0); T=0:01:2*pi; X=10+20*cos(T); Y=20+15*sin(T); plot(X,Y,-) hold on plot(y(:,1),y(:,2),*)在ch

10、ase3m中，不断修改tf的值,分别取tf=20, 40, 80,可以看出,狗永远追不上慢跑者（二）、实验题1一个小孩借助长度为a的硬棒拉(或推)某玩具此小孩沿某曲线行走，计算并画出玩具的轨迹2讨论资金积累、国民收入与人口增长的关系（1）若国民平均收入x与人口平均资金积累y成正比，说明仅当总资金积累的相对增长率k大于人口的相对增长率r时，国民平均收入才是增长的（2）作出k(x)和r(x)的示意图，分析人口激增会导致什么后果五、实验总结写出其操作步骤及程序的运行结果。实验五 离散模型求解一、实验目的学会用Matlab软件对矩阵进行一些数值计算，学会用Matlab软件解线性方程组。二、实验类型:设

11、计三、实验环境计算机、软件Matlab7.0以上的环境四、实验内容1、 产生一个4阶的随机矩阵，执行下面的操作：(1) 求其行列式，检验其是否可逆；若可逆，求其逆矩阵。(2) 计算该矩阵的特征值、特征向量。(3) 将该矩阵化为行最简的阶梯形。(4) 验证矩阵的特征值之和等于矩阵主对角元之和，特征值之积等于矩阵的行列式。2、 判断下面的线性方程组是否有解，若有解求其通解。(1)(2) （3）3、 求矩阵的特征值和特征向量。五、实验总结根据实验操作和实验报告要求，完成实验报告。 实验六 统计回归模型的求解一、实验目的掌握Matlab求解统计回归模型的方法。二、实验类型:验证三、实验环境计算机、软件

12、Matlab7.0以上的环境四、实验内容(一)、示例用观测物体降落的距离s与时间t的关系，得到数据如下表，求s关于t的回归方程.t (s)1/302/303/304/305/306/307/30s (cm)11.8615.6720.6026.6933.7141.9351.13t (s)8/309/3010/3011/3012/3013/3014/30s (cm)61.4972.9085.4499.08113.77129.54146.48法一、直接作二次多项式回归：t=1/30:1/30:14/30;s=11.86 15.67 20.60 26.69 33.71 41.93 51.13 61.4

13、9 72.90 85.44 99.08 113.77 129.54 146.48; p,S=polyfit(t,s,2)得：法二、 化为多元线性回归： t=1/30:1/30:14/30;s=11.86 15.67 20.60 26.69 33.71 41.93 51.13 61.49 72.90 85.44 99.08 113.77 129.54 146.48;T=ones(14,1) t (t.2);b,bint,r,rint,stats=regress(s,T);b,stats得：(二)、实验题1、 设某商品的需求量与消费者的平均收入、商品价格的统计数据如下，建立回归模型，预测平均收入为

14、1000、价格为6时的商品需求量.需求量10075807050659010011060收入1000600 1200500300400130011001300300价格57668754392、 某零件上有一段曲线，为了在程序控制机床上加工这一零件，需要求这段曲线的解析表达式，在曲线横坐标xi处测得纵坐标yi共11对数据如下：xi02468101214161820yi0.62.04.47.511.817.123.331.239.649.761.7求这段曲线的纵坐标y关于横坐标x的二次多项式回归方程.3、根据课本（P294页）牙膏的销售量有关数据建立模型,并分析。五、实验总结根据实验操作和实验报告要求，做好实验报告。 附件：数学建模实验报告姓 名学 号系 别班级指导教师实验日期课程名称数学建模同组实验者一、实验名称： 二、实验目的：三、实验类型：验证、演示、设计、综合四、实验材料、工具、或软件多媒体计算机、Matlab软件五、实验要求六、试验内容与步骤：（附上必要的调试过程及运行结果截图，可加附页）七、实验总结八、教师评语（或成绩） 教师签字： 年 月 日