黄世杰,豆彩霞,李金玲指导老师,数模组

1、 降落伞的选择摘 要本文针对降落伞的选择问题建立相应优化模型。首先，对降落伞进行受力分析，根据有关数据进行拟合，得出空气阻力系数。其次，题中已给出限制降落伞的最大落地速度为 ，所以当速度为20m/s时，伞的承载量最大。建立高度与时间，速度与时间，最大速度，伞面面积（题中已给出可能选购的每种伞的半径）的方程组，分别算出每种伞的最大承载量。22.533.54150.6787235.4355339.0272461.4536602.7150再次，运用软件进行线性规划求解得：即购买半径为3的降落伞6个时，降落伞总费用最少为4932元。最后，对所建模型进行评价，并对模型的改进方向提出建议。关键词：线性规划

2、、空气阻力系数、拟合、最大载重量 1 问题重述1.1 背景资料为向灾区空投救灾物资共，需选购一些降落伞。1.2 已知信息已知空投高度为，要求降落伞落地时的速度不能超过。降落伞面为半径的半球面，用每根长， 共16根绳索连接的载重的物体位于球心正下方球面处，每个降落伞的价格由三部分组成。伞面费用由伞的半径决定，见附录【1】；绳索费用由绳索总长度及单价4元/米决定；固定费用为200元。降落伞在降落过程中受到的空气阻力，可认为与降落速度和伞的受力面积的乘积成正比。为了确定阻力系数，用半径、载重的降落伞从高度作降落试验，测得各时刻的高度。见附录【2】1.3 提出问题空气阻力系数，最大载重量，每个降落伞的

3、单价，试根据以上条件确定降落伞的选购方案，即共需多少个，每个伞的半径多大（在附录【1】选择），在满足空投要求的条件下，使总费用最低。2 模型的假设1、 假设空投物资的瞬时伞已打开。2、 空投物资的总数可以任意分割。3、 空气的阻力系数与除空气外的其它因素无关。4、 降落伞和绳的质量可以忽略不计。5、 假设降落伞只受到竖直方向上的空气阻力作用。3 符号的说明1、:重力加速度2、:降落伞的载重量3、:空气阻力系数4、:降落伞的表面积5、:降落伞的加速度6、 :降落伞的速度7、 :降落伞的位移8、 :降落伞离地高度9、分别为每种伞的个数4 问题的分析由题可知，每个伞的价格由三部分组成。伞面费用，由其

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

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

6、的二次微分等于加速度建立方程得: 用解微分方程得：(见附录【3】)则题目已经给t-h数据为：时刻036912151821242730高度500470425372317264215160108551对给定的数据以为拟合函数进行拟合，,得出。（见附录【4】）5.2求解最大承载量用速度对时间的微分等于加速度，且建立方程组得：用解得(见附录【5】)由前面的和函数建立方程组得：最大载重量取伞在下降到地面的瞬间达到最大速度，此时，由方程组调用分别解得半径为r的降落伞在满足空投条件下的最大载重量如下表：（见附录【6】）表1 最大载重量22.533.54150.6787235.4355339.0272461.

7、4536602.71505.3线性规划求解数量和费用由分析可知每种伞的单价：由题可知为:22.533.54费用（元）651703506601000为： 为固定值即：由以上数据求得，如：表2 每种伞的单价22.533.54单价446.1596.3821.51176.81562取整44659682211771562我们设每种伞分别取个,则其目标函数为：对其进行优化求解z的最小值，就是所需的最小费用。由分析可知其限制条件如下：用LINGO求解得（程序见附件【7】）。最少总费用为4932元。6 模型的评价 改进 推广6.1 模型的评价优点 （1）本模型的求解过程大量的运用了电脑软件，使得计算更加精确。

8、（2）模型原理简单明了，容易理解与灵活运用（3）模型的建立中有成熟的理论基础和利用专业的软件进行求解，可信度较高。缺点 (1)本模型未考虑降落伞打开的时间，将其假设成在下降时伞就已经打开。(2)由于在实际生活中降落伞还受到风向的影响，本模型假设的是理想的状态下（无风）(3)在实际生活中载重物也受空气阻力的影响，本模型假设的是理想状态中（载重物随意分割）。6.2模型的改进由于本模型假设的是在物资抛落的瞬时伞已打开，而在实际情况中物资抛落后应有一段自由落体运动，载重物还受空气阻力影响，降落伞降落过程中受风向的影响。在模型的改进时应考虑到这三点，以便让模型更切合实际。6.3 模型的推广本案例建立的优

9、化模型解决了降落伞费用最少的问题。因此，本模型还可应用与其他类似，如：生产成本最低，利润最高等方面的问题,只需稍微改动模型即可。7 参考文献1、数学实验萧树铁主编 高等教育出版社 1999.72、孙祥 徐流美 吴清编著MATLAB7.0基础教程清华大学出版社 2005.53、 姜启源 谢金星 叶俊编数学模型高等教育出版社 2003.84、于义良总主编，运筹学，中国人民大学出版社，2004.38 附录附录【1】伞面费用22.533.54651703506601000附录【2】各时刻与高度的关系时刻036912151821242730高度500470425372317264215160108551

10、附录【3】求解位移的程序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附录【4】拟合k程序建立一个名为myfun的m文件function F=myfun(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

11、 window中输入下列命令：xdata=0 3 6 9 12 15 18 21 24 27 30;ydata=500470425372317264215160108551 ;x0=1;x=lsqcurvefit(myfun,x0,xdata,ydata)附录【5】求解速度程序dsolve('m*Dv+k*S*v-m*g=0','v(0)=0','t')解得：V(t)=g/k/s*m-exp(-k*s/m*t)*g/k/s*m 附录【6】 求最大载重量在matlab中建立一个名为myfun的m文件，如下：function F=myfun(x)r=

12、2;g=9.8;k=2.9377;s=2*pi*r2;的最F=x(1)2*g/(k2*s2)*exp(-k*s*x(2)/x(1)+x(1)*g*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中输入以下命令：量M(r)r（m）22.533.54x0 = 1; 1; % 初始点options=optimset('Display','iter'); % 显示输出信息x = fsolve(myfun,x0,opti

13、ons)在m文件中更改r的值，然后在命令窗口重复输入以上命令就可分别求出不同半径的降落伞大载重量。如下：最大载重m150.6787235.4355339.0272461.4536602.7150附录【7】优化求解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. Obje