集大计院数学建模实验.docVIP

第二次实验内容 初等模型实验 课程：：计算1412： 实验目的： 掌握初等模型的建立的基本思路和方法，并了解其求解过程。对给定的初等模型问题能够借助Matlab工具进行求解。 实验内容： 实验1 用Matlab验证划艇比赛成绩模型的结果，通过数值结果来检验你所得到的模型正确性。(首先要阅读本目录中的Matlab数据拟合和matlab数据处理的相关材料) 模型检验：代码 Matlab 拟合函数图像： 结果分析： 代码： function E=fun(a,x,y) Y=a(1)*x.*x+0.75*x; E=y-Y; end x=[29.3 44 58.7 73.3 88 102.7 117.3]; y=[44 78 124 186 268 372 506]; a0=[0.5]; options=optimset(lsqnonlin); a=lsqnonlin(@fun,a0,[],[],options,x,y) ; %绘图 plot(x,y,o); hold on; x=[0:200]; y=a(1)*x.*x+0.75*x; plot(x,y,-); hold off Matlab拟合函数图像： 结果分析： 实验3 .从教材P56中的第7，13，14题，任选一题，建立相应的初等模型，并借助matlab进行求解，并给出合理的模型解释。 生物学家认为，对于休息状态的热血动物，消耗能量主要用于维持体温，能量与从心脏到全身的血流量成正比，而体温主要通过身体表面散失，现在要求建立一个动物体重与心率之间关系的模型，并用下面的数据加以检验 动物? 体重（g）? 心率（次/分）? 田鼠? 25? 670? 家鼠? 200 ?420? 兔? 2000 ?205? 小狗 ?5000 ?120? 大狗 ?30000 ?85? 羊? 50000? 70? 人? 70000? 72? 马 ?450000 ?38 模型假设： 1 4.?假定动物（包括人）体表的散热率完全恒等于心跳的产热率；? 5.?假定心脏体积与动物（包括人）的体积成正比，体积与重量成正比； 模型的建立： 由数学知识可知：体积V正比于长度的立方，表面积S正比于长度的平方，于是有 VL 3 S∝L2 则有S=K?*V^(2/3) 由体表散热率=心跳产热率，求得： 体表散热率?*S=K?*K?*V^(2/3) (1) 由心脏体积正比于动物的体积得，令心率为K?*V所以 心跳产热率=K?*V*P （） 则(1)(2)得： ?*S=K?*K?*V^(2/3)=K?*V*P 所以 P=(K?*K?*V^(2/3))/K? (3) 由体重和体积成正比都，令体重为 V=K?*G (4) 整理化简得 (K?*K?*G^-(1/3))/(K?*K?^(1/3)) 即 ^(-1/3) (5) 模型的求解： 由5)式用 G=[25 200 2000 5000 30000 50000 70000 450000]; P=[670 420 205 120 85 70 72 38]; plot(G,P,ro); title(已给数据的散点图); xlabel(生物的体重); ylabel(生物的心率); 为了得到更精确的信息，对P=K*G^(-1/3)进行修改得到lgP=lg(K*G^(-1/3))即lgP=lgK-1/lgG。令y=lgP，x=lgG，a=lgK，得y=a-1/G=[25 200 2000 5000 30000 50000 70000 450000]; P=[670 420 205 120 85 70 72 38]; x=log(G); y=log(P); plot(x,y,go); 再利用最小二乘法直线拟合，当所测各Yi值与拟合直线上的a+bxi之间的偏差的平方和最小，即yi-(a-1/3xi)]^2最小，系数a最好，拟合公式即为最佳经验公式。即a=(x*y*