数学建模降落伞的选购问题.docx

1、数学建模降落伞的选购问题数学建模降落伞的选购问题本模型研究的是降落伞的选购方案问题，目的是在满足空投要求的条件下，使费用最少。为了方便对降落伞进行受力分析，我们把降落伞和其负载的物资看做一个整体，忽略了伞和绳子的质量，并假设降落伞只受到竖直方向上空气阻力和重力的作用。通过对降落伞在空中的受力情况的分析建立起了高度与时间的方程，然后以高度与时间的方程作为拟合曲线与题中给出的时间与高度的数据进行拟合，得出阻力系数k的值。我们建立了速度与质量的方程，并证明其为严格增函数（证明过程见建模与求解）。由于题中已限制降落伞的最大落地速度为20m/s，所以当速度为20m/s时，伞的承载量最大。建立高度与时间，

2、速度与时间的方程组，代入最大速度20m/s，高度500m,伞的半径（题中已给出可能选购的每种伞的半径），分别计算出每种伞的最大承载量。最后运用LINGO软件进行线性规划求解得：x1=0,x2=0,x3=6,x4=0,x5=0.即购买半径为3m的降落伞6个时总费用最少为4932元。 为向灾区空投救灾物资共2000kg，需选购一些降落伞。已知空投高度为500m，要求降落伞落地时的速度不能超过20m/s。降落伞面为半径r的半球面，用每根长 L, 共16根绳索连接的载重m的物体位于球心正下方球面处，每个降落伞的价格由三部分组成。伞面费用C由伞的半径r决定，见表1；绳索费用1C由绳索总长度及单价4元/米

3、决定；固定费用C为200元。 23表1 r（m） 2 2.5 3 3.5 4 费用（元） 65 170 350 660 1000 降落伞在降落过程中受到重力作用外还受到的空气阻力，可以认为与降落速度和伞的受力面积的乘积成正比。为了确定阻力系数，用半径r=3m、载重m=300kg的降落伞从500m高度作降落试验，测得各时刻的高度 ，见表2。 表2 时刻t（s） 0 3 6 9 12 15 18 21 24 27 30 高度h（m） 500 470 425 372 317 264 215 160 108 55 1 试根据以上条件确定降落伞的选购方案，即共需多少个，每个伞的半径多大（在表1中选择），

4、在满足空投要求的条件下，使费用最低。1、假设空投物资的瞬时伞已打开。 2、空投物资的总数2000kg可以任意分割。 3、空气的阻力系数与除空气外的其它因素无关。 4、降落伞和绳的质量可以忽略不计。 5、假设降落伞只受到竖直方向上的空气阻力作用。 1、m 降落伞的负载重量 2、g 重力加速度 3、a 降落伞的加速度 4、k 空气的阻力系数 5、S 降落伞的伞面面积 6、v 降落伞的速度 7、H 降落伞的位移 8、h 降落伞离地高度 9、x1,x2,x3,x4,x5 分别为每种伞的个数 由题意可知每个伞的价格由三部分组成：三面费用C、绳索费用C、固定12费用C。伞面费用由伞的半径r决定；绳索费用C

5、由绳索的长度及单价决定，32由图一可知绳索的长度又由降落伞的半径决定即；固定费用为定值200。L,2r因为题中已给出每种伞面的半径，所以每种伞的价格为定值。要想确定选购方案，即共需半径（在题中给出的半径中选择）为多大的伞的数量，在满足空投物资要求的条件下使总费用最少。因此，我们需要确定每种伞的最大承载量。然后进行线性规划，确定总费用和每种伞的个数。 要确定最大载重量，我们需对降落伞进行受力分析（如图二）。降落伞在降落过程中除受到竖直向下的重力作用外还受到竖直向上的空气阻力的作用，而由题可知空气阻力又与阻力系数、运动速度、伞的受力面积有关。运动速度和受力面积是已知的，所以要想确定每种伞的最大承载

6、量，就必须先要确定空气的阻力系数。 图一 图二 对图二的分析可知降落伞的运动状态是做加速度趋近于0的加速运动。因此，我们可以建立一个位移与时间的函数关系式，在根据题中所给的数据拟合出阻力系数k的值。然后再建立一个速度与时间的函数关系式，两个关系式联立求解出最大载重量（其中高度和速度由题目已经给出）。最后用LINGO软件进行线性规划算出问题要的结果。 （1）首先确定阻力系数K 为了方便对物资进行受力分析，我们把降落伞和物资看作一个整体如图二。由假设5可知物体A只受到竖直向上的空气阻力和竖直向下的重力作用。又由题可知空气阻力与降落速度v和伞的受力面积S的乘积成正比。则物体A在竖直方向上受到的合外力

7、为： F,mg,kSv 合由运动学方程： F,ma合得 Fmg,kSv合 a,mm由物体位移H和时间的二次微分等于加速度建立方程得: 2dHmg,kSv, 2mdt用MATLAB解微分方程得：(程序见附录【1】) ,kS22tmgmgtmgmH(t),e，, 2222kSkSkS则 ,kS22tmgmgtmgmh(t),500,e,， 2222kSkSkS题目已经给t-h数据为： 时刻t（s） 0 3 6 9 12 15 18 21 24 27 30 高度h（m） 500 470 425 372 317 264 215 160 108 55 1 2对给定的数据以h(t)为拟合函数进行拟合，r=

8、3m,m=300kg,g=9.8,得出 S,2,rk=2.9377 。（程序见附录【2】） （2）求解最大承载量 用速度对时间的微分等于加速度，且v=0建立方程组得： 0dvmg,kSv ,dtmv,0 0用MATLAB解得(程序见附录【3】) ,kStmgmgmev(t), kSkSH(t)v(t)由前面的和函数建立方程组得： ,kst,mgmgm,v(t)e,ksks,kst22mmgtmgmg, ,，,H(t)e,2222ksksks,2,S2r,H500h,k=2.9458,g=9.8,r=2 2.5 3 3.5 4 关系式中t看做一个因为降落伞在下落过程中其质量是不变的，所以我们把v

9、(t)定值，则关于m的方程为 ,kStmgmgme(), vmkSkS从上式我们可以知道是关于m的单调递增函数（证明见附件【7】），v(m)mg并且如果存在平衡状态则必须满足，那么 而又通过对mg,kvsv,ks,kstmgmgmgmv(t),e 分析，只有在，这与实际矛盾，t,，,时，才有v(t),ksksks故降落伞是一直做加速度减小的加速运动，不存在平衡状态。因此，求最大载重量取伞在下降到地面的瞬间达到最大速度v(t),20m/s,此时H(t),500，由方程组调用MATLAB分别解得半径为r的降落伞在满足空投条件下的最大载重量M(r)如下表：（程序见附录【5】） r（m） 2 2.5

10、3 3.5 4 最大承载150.6787 235.4355 339.0272 461.4536 602.7150 （kg） 取整(kg) 150 235 339 461 602 （3）线性规划求解数量和费用 由分析可知每种伞的单价： C,C，C，C 123C由题可知为: 1r（m） 2 2.5 3 3.5 4 费用（元） 65 170 350 660 1000 为： C2C,16，2r，4 2C 为固定值即： 3C,200 3由以上数据求得每种伞的单价见下表： r（m） 2 2.5 3 3.5 4 单价C 446 596.3 821.5 1176.8 1562 取整 446 596 822 1

11、177 1562 ,x,x,x,x个,则其目标函数为： 我们设每种伞分别取x12345z=446x+596x+822x+1177x+1562x 12345对其进行优化求解z的最小值，就是所需的最小费用。由分析可知其限制条件如下： s.t. 150x+235x+339x+461x+602x=2000; 12345(x，x，xx，x); ,N123，54用LINGO求解得（程序见附件【6】） x=0,x=0,x=6,x=0,x=0。 12345最少总费用为4932元。 优点： 1、本模型的求解过程大量的运用了电脑软件，使得计算更加精确。 缺点： 1、本模型未考虑降落伞打开的时间，将其假设成在下降时

12、伞就已经打开。 2、由于在实际生活中降落伞还受到风向的影响，本模型假设的是理想的状态下（无风） 改进：由于本模型假设的是在物资抛落的瞬时伞已打开，而在实际情况中物资抛落后应有一段自由落体运动。在模型的改进时应考虑到这一点，以便让模型更切合实际。 1、数学实验萧树铁主编 高等教育出版社 1999 7 1 附录【1】求解位移的程序 H=dsolve(m*D2H+k*S*DH=m*g,H(0)=0,DH(0)=0,t) 解得：g/k2/S2*m2*exp(-k*S/m*t)+g/k/S*m*t-1/k2/S2*m2*g 附录【2】拟合k程序 建立一个名为myfun的m文件 function F=my

13、fun(x,xdata) s=2*pi*32; m=300; g=9.8; F=500-m2*g/(x(1)2*s2)*exp(-x(1)*s*xdata/m)-m*g*xdata/(x(1)*s)+m2*g/(x(1)2*s2); 在matlab command window中输入下列命令： xdata=0 3 6 9 12 15 18 21 24 27 30; ydata=500 470 425 372 317 264 215 160 108 55 1 ; x0=1; x=lsqcurvefit(myfun,x0,xdata,ydata) 附录【3】求解速度程序 v=dsolve(m*Dv

14、+k*S*v-m*g=0,v(0)=0,t) ,kStgmgmm解得：v(t),e kSkS附录【4】在v-t,m函数中对m求二阶导数 syms m t g S k f=g*m/(k*S)-g*m/(k*S)*exp(-k*S*t/m); diff(f,m2) 求得： -g/m3*t2*k*s*exp(-k*s/m*t) 附录【5】 求最大载重量 在matlab中建立一个名为myfun的m文件，如下： function F=myfun(x) r=2.5; g=9.8;k=2.9458; s=2*pi*r2; F=x(1)2*g/(k2*s2)*exp(-k*s*x(2)/x(1)+x(1)*g

15、*x(2)/(k*s)-x(1)2*g/(k2*s2)-500; g*x(1)/(k*s)-g*x(1)/(k*s)*exp(-k*s*x(2)/x(1)-20; 在matlab中command window中输入以下命令： x0 = 1; 1; % 初始点 options=optimset(Display,iter); % 显示输出信息 x = fsolve(myfun,x0,options) 在m文件中更改r的值，然后在命令窗口重复输入以上命令就可分别求出不同半径的降落伞的最大载重量。 分别求解可得最大载重量如下表： r（m） 2 2.5 3 3.5 4 m 150.6787 235.43

16、55 339.0272 461.4536 602.7150 附录【6】优化求解 min=446*x1+596*x2+822*x3+1177*x4+1562*x5; 150*x1+235*x2+339*x3+461*x4+602*x5=2000; x1=0; x2=0; x3=0; x4=0; x5=0; gin(x1);gin(x2);gin(x3);gin(x4);gin(x5); 求解得： Global optimal solution found. Objective value: 4932.000 Extended solver steps: 0 Total solver iterat

17、ions: 0 Variable Value Reduced Cost X1 0.000000 446.0000 X2 0.000000 596.0000 X3 6.000000 822.0000 X4 0.000000 1177.000 X5 0.000000 1562.000 附件【7】 证明速度与质量m成正比关系 v(m)由高数定理可知：函数的一阶导数大于零，则原函数是单调递增的。一阶导数小于零，则原函数是单调递减的。 ,kStmgmgme(), vmkSkS对求一阶导数得： v(m),kSt,kStmmgegteg,v(m),， kSmSk由上式分析可知无法确定其是否大于零，在对其求二阶导数为： ,kSt2gtkSmv(m),e,0 3m则一阶导数为单调递减函数，当m趋近于无穷大时对一阶导数求极限可知 ,kSt,kStmmgegteggg,lim(,，),，,0 m,kSmSkkSkS由此可得： v(m),0 则原函数是单调递增函数，即速度v和m是成正比关系的。