降落伞选择的数学模型.docVIP

降落伞选择的数学模型 摘要 本文建立了一个关于对降落伞选择使得费用最节省的模型。根据降落伞伞面的大小、材质、所悬挂重物的质量与降落伞下落速度之间的关系，建立在空气阻力作用下物体降落的模型，运用最小二乘法拟合出阻力系数。利用高度与时间的关系式，计算出不同半径的降落伞的最大载重量。 关键词 空气阻力系数；线性规划；matlab 数学建模是将数学知识、实际问题与计算机应用有机地结合起来，旨在提高学生的综合素质与分析问题、解决问题的能力。对于现实中的原型，为了某个特定目的，作出一些必要的简化和假设，运用适当的数学工具得到一个数学结构。数学建模是利用数学语言（符号、式子与图象）模拟现实的模型。把现实模型抽象、简化为某种数学结构是数学模型的基本特征。 本文建立了一个关于对降落伞选择使得费用最节省，最优化的模型。根据降落伞伞面的大小、材质、所悬挂重物的质量与降落伞下落速度之间的关系，建立在空气阻力作用下物体降落的物理模型，并对其进行有效分析。运用线性最小二乘法，拟合出空气的阻力系数。利用高度h与时间t的关系式，计算出不同半径的降落伞的最大载重量。 1 问题引出 为向灾区空投一批共2000 kg的救灾物资，需选购一些降落伞，空投高度为500 m，其落地时的速度不能超过20 m每秒，伞面为半径为r的半球面，用每根长l共16根绳索连接的重物m位于球心正下方球面处，如图： 每个降落伞的价格由三部分组成：伞面费用由伞的半径r决定，绳索费用由绳索总长度及单价4元每米决定，固定费用为200元；降落伞在降落过程中除受到重力外，受到空气的阻力，可认为与降落的速度和伞的面积的乘积成正比。为了确定阻力系数，用半径r=3 m，载重m=300 kg的降落伞以500 m高度作试验，测得各时刻t的高度x，试确定共需多少个伞，每个伞的半径多大（在给定半径的伞中选），在满足空投要求的条件下，使总费用最低。 我们要考虑降落伞伞面的大小、材质，所悬挂重物的质量与降落伞下落速度之间的关系，因为降落伞下降过程是一个物理模型，根据物理理论，系统在下降过程中做加速度减小的加速运动，直到所受阻力等于自身重力时，加速度为零，速度达到最大。由已知，降落伞在降落过程中受到的空气阻力与降落速度以及伞面积成正比，所以我们要先确定它们的比例系数k，在求k之前必须求出时间t与高度h的关系式。根据kvs=mg及20 m/s的最大速度求出不同半径的降落伞的最大载重量m，最后通过lingo软件得出结果。 2 问题求解 降落伞在下降过程中受到的力符合牛顿定律，即am=mg-f，a=（mg-kvs）/m，加速度、速度、位移之间通过微积分知识可知： ，。且时间与高度之间的关系如表3。 根据物理公式可以得到高度h(t)的表达式： ，。 再根据表3数据，对阻力系数k利用非线性最小二乘法拟合求解。其中重力加速度g=9.8 m/s，m=300 kg，在matlab中运行程序如下： fun=inline(‘9.8*300*t)/(56.52*k)-(90000*9.8)/(k*56.56)*exp (-56.52*t*k/300)-(90000*9.8)/(k*56.56)’,’k’,’t’);t=[0 3 6 9 12 15 18 21 24 27 30]; y=[0 30 75 128 183 236 285 340 392 445 499];k=lsqcurvefit(fun,2,t,y) 运行结果为k=2.9720 由如上分析，继而根据表1可建立线性规划模型如下： min=x1*(65+16*1.414*2*4+200)+x2*(170+16*1.414*2.5*4+ 200)+x3*(350+16*1.414*3*4+200)+x4*(660+16*1.414*3.5*4+200)+x5*(1000+16*1.414*4*4+200)； x1*151+x2*236+x3*340+x4*463+x5*605=2000;x1=0;x2=0;x3=0; x4=0;x5=0； @gin(x1);@gin(x2);@gin(x3);@gin(x4);@gin(x5)； 在lingo运行结果为：采用半径为2 m的降落伞1个，半径为2.5 m的降落伞2个，半径为3 m的降落伞4个，这时2000 kg物资能被完全投放，而且使得总费用最低为4924.7元。 3 结论 由上述结论可以看出，一个实际问题在经过数学方法进行分析之后，可以在操作之前经过简单的运算就进行很好的预测，这样可以节省大量的人力物力，并且杜绝浪费行为，利用最少的条件得到最好的结果，这就是数学模型的作用与魅力，希望本文方法与结果可以为读者提供帮助。 参考文献 [1]么焕民,孙秀梅,等.数学建模[m].哈尔滨:哈尔滨工业大学出版社,2003. [2]杨