第2讲第3讲作业答案.pptVIP

Matlab与数学建模 计算机科学与信息工程系 软件工程教研室 许晓华 xuxiaohua@tjau.edu.cn 作业1 问题: 在医院的外科手术室,往往需要将病人安置到活动病床上,沿走廊推到手术室或送回病房.然而有的医院走廊较窄,病床必须沿过道推过直角拐角(如图所示). 设标准病床长2米,宽1米，拐弯前的过道宽1.5米，拐弯后的过道宽1.2米，问标准的病床能否安适的推过拐角? 首先，将问题简化。想象床缩为一根杆子，问题转化为要想让杆子绕过拐角，杆子的长度要受到多大的限制。如图所示，过道的宽度分别是a和b。PQ线段的长度即固定? 角时能安放的最大杆长。 此问题即梯子长度问题 Matlab求得最大杆长3.8105米 现在，回到病床问题。 当病床的宽度一定时可以很方便地求出病床的长度ＡＢ 但是，当病床的长度一定时，要计算病床的宽度就比较麻烦 因此从定宽变长的角度来考虑问题 根据图中的三角关系可以得到： AB = AN + NB syms x; f=(1.2-sin(x))*sec(x)+(1.5-cos(x))*csc(x); %直接求函数的极小值 [xmin,fmin]=fminbnd(f,0,pi/2) xmin = 0.9137 fmin = 1.7914 即，当床宽度固定为1米时，长度大于1.7914的床都是过不去的 看成直线 (1460^2+380^2)^(1/2) ans = 1.5086e+003 9个点拟和成一条抛物线 x=[0,220,420,560,680,780,900,1040,1240,1460]; y=[0,50,95,135,173,207,245,285,330,380]; plot(x,y,o) hold on s=polyfit(x,y,2) syms x; f=s(1)*x^2+s(2)*x+s(3) g=diff(f,x) h=(1+g^2)^(1/2) L=int(h,x,0,1460) vpa(L) 加上每段的中点 x=[0,110,220,320,420,490,560,620,680,730,780,840,900,970,1040,1140,1240,1350,1460]; y=[0,-1,50,49,95,94,135,134,173,172,207,206,245,244,285,284,330,329,380]; plot(x,y,o) hold on s=polyfit(x,y,2) syms x; f=s(1)*x^2+s(2)*x+s(3) g=diff(f,x) h=(1+g^2)^(1/2) L=int(h,x,0,1460) vpa(L) 循环拟和成9段抛物线 d=[220,200,140,120,100,120,140,200,220]; h=[50,45,40,38,34,38,40,45,50]; sum=0 for i=1:9 a=[0,d(i)]; b=[0,h(i)]; s=polyfit(a,b,2); syms x; f=s(1)*x^2+s(2)*x+s(3); g=diff(f,x); j=(1+g^2)^(1/2); L=int(j,x,0,d(i)); sum=sum+L; end vpa(sum) 加上每段的中点（循环） d=[220,200,140,120,100,120,140,200,220]; h=[50,45,40,38,34,38,40,45,50]; sum=0 for i=1:9 a=[0,d(i)/2,d(i)]; b=[0,-1,h(i)]; s=polyfit(a,b,2); syms x; f=s(1)*x^2+s(2)*x+s(3); g=diff(f,x); j=(1+g^2)^(1/2); L=int(j,x,0,d(i)); sum=sum+L; end vpa(sum) 人口问题 首先建立一个M文件，保存fun函数 function x=fun(c,t) x=c(1)./(1+(c(1)./c(2)-1)*exp(-c(3)*t)); 再建立一个M文件，调用fun函数 clear all; close all; t=0:1:20; x=[3.9 5.3 7.2 9.6 12.9 17.1 23.2 31.4 38.6 50.2 62.9 76.0 92.0 106.5 123.2 131.7 150.7 179.3 204.0 226.5 251.4]; plot(t,x,o); hold on; c(1)=1000; c(2)=3.9; c(3)=mean(diff(x)./diff(t)./x(1: