降落伞选购问题.docVIP

下面是几道Lingo软件的练习题，请同学们练习。 要求：1、编写lingo程序并能正确运行，特别要注意第3题模型及结果是否正确； 2、将问题、解题思路、lingo程序和运算结果（注意结果的正确性）写成word文档2003版或pdf文件于8月15日发送到dinggenhong@126.com, 文件名同前几次的论文要求。 3，向灾区空投救灾物资共2000kg，需选购一些降落伞。已知空投高度为500m，要求降落伞落地时的速度不能超过20m/s。降落伞面是半径为r的半球面，用每根长l共16根绳索连接的载重m位于球心正下方球面处，如图所示。 每个降落伞的价格由三部分组成。伞面费用C1由伞的半径r决定，见表1；绳索费用C2由绳索总长度及单价4元/m决定；固定费用C3为200元。 降落伞在坠落过程中受到的空气阻力，可以认为与坠落速度和伞面积的乘积成正比。为了确定阻力系数。用半径r=3m、载重m=300kg的降落伞从500m高度做降落试验，测得各时刻t的高度x，见表1。 表1 时刻t的高度x t/s 0 3 6 9 12 15 18 21 24 27 30 x/m 500 470 425 372 317 264 215 160 108 55 1 试确定降落伞的选购方案，即共需多少个，每个伞的半径多大（在表2中选择），在满足空投要求的条件下，使费用最低。 表2 降落伞的选购方案 r/m 2 2.5 3 3.5 4 C1/元 65 170 350 660 1000 问题分析 由题意可知每个伞的价格由三部分组成：三面费用C1、绳索费用C2、固定费用C3。伞面费用由伞的半径r决定；绳索费用C2由绳索的长度及单价决定，绳索的长度又由降落伞的半径决定即；固定费用为定值200。因为题中已给出每种伞面的半径，所以每种伞的价格为定值。要想确定选购方案，即共需半径（在题中给出的半径中选择）为多大的伞的数量，在满足空投物资要求的条件下使总费用最少。因此，我们需要确定每种伞的最大承载量。然后进行线性规划，确定总费用和每种伞的个数。 要确定最大载重量，我们需对降落伞进行受力分析（如图二）。降落伞在降落过程中除受到竖直向下的重力作用外还受到竖直向上的空气阻力的作用，而由题可知空气阻力又与阻力系数、运动速度、伞的受力面积有关。运动速度和受力面积是已知的，所以要想确定每种伞的最大承载量，就必须先要确定空气的阻力系数。 模型建立与求解 （1）确定阻力系数k：降落伞在坠落过程中受到的空气阻力，可以认为与坠落速度和伞面积的乘积成正比，根据牛顿第二定律：物体位移H和时间的二次微分等于加速度建立方程得,用matlab进行常微分方程的求解，得到用MATLAB解微分方程得：(程序见附录【1】) 则 题目已经给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 为拟合函数用lingo进行求解，r=3m,m=300kg,g=9.8,,得出 k=2.941020。(程序见附录【4】)。 （2）确定降落伞的载重； 已知k，用速度对时间的微分等于加速度，且v（0）=0建立方程组得： 用MATLAB解得: (程序见附录【2】) 由前面的和函数建立方程组得： k=2.941020,g=9.8,r=[2 2.5 3 3.5 4] 因为降落伞在下落过程中其质量是不变的，所以我们把关系式中t看做一个定值，则关于m的方程为用mathmatica对v(m)进行求导之后得（程序见附录【3】），当,所以可以知道是关于m的单调递增函数，由于限定降落伞落地时的速度不能超过20m/s。因此，求最大载重量取伞在下降到地面的瞬间达到最大速度,此时，由方程组调用lingo分别解得半径为r的降落伞在满足空投条件下的最大载重量如下表：（程序见附录【5】） r（m） 2 2.5 3 3.5 4 ）461.4536 602.7150 取整(kg) 150 235 339 461 602 （3）分配方案： 由分析可知每种伞的单价： 由题可知为: r（m） 2 2.5 3 3.5 4 65 170 350 660 1000 ; 由以上数据求得每种伞的单价见下表： r（m） 2 2.5 3 3.5 4 取整 446 596 822 1177 1562 我们设每种伞分别取x1,x2,x3