降落伞的选择问题.doc

降落伞的选择问题一， 问题的提出与重述1.1问题提出在物资救援中,空投已经成为一种十分重要且便利的方式,由于降落伞难以多次利用,所以如何减少空投的成本,让人们有更多的资金购买救援物资已经成为了一个不可忽视的课题。1.2问题重述为向灾区空投救灾物资共2000kg，需选购一些降落伞。已知空投高度过500m，要求降落伞落地时的速度不能超过20m/s。降落伞面为半径r的半球面，用每根长共16根绳索连接的载重m位于球心正下方球面处，如下图： 每个降落伞的价格由三部分组成。伞面费用C1由伞的半径r决定，；绳索费用C2由绳索总长度及单价6元/米决定；固定费用C3为400元。降落伞在降落过程中受到的空气阻力，可以认为与降落速度和伞面积的乘积成正比。为了确定阻力系数，用半径r=3m、载重m=300kg的降落伞从500m高度作降落试验，测得各时刻t的高度。试确定降落伞的选购方案，即共需多少个，每个伞的半径多大（在表1中选择），在满足空投要求的条件下，使费用最低。二，问题分析本文主要解决的是在满足空投要求下的降落伞的选择问题，是典型的优化问题，通过对题目的分析可以进一步确定是整数线性规划问题。本题所建的模型的目标函数比较简单，主要是约束条件，而在约束条件中每种降落伞的最大载重质量又与空气阻力系数是有一定的量化关系的，因此此模型的关键在于求空气阻力系数。三，模型假设1 降落伞和绳索的质量均不计;2 救灾物资的大小不计,可以看作质点处理;3 降落伞下落的初速度为0;4 救灾物资可以任意分割.四，变量及符号说明第i种降落伞:半径:,伞面费用:;所需绳索长:;绳索费用:6;最大载重质量:;费用:;选用的个数:.总的费用:Z.空气阻力系数:k.重力加速度:g(取).五，模型建立与求解由载重m位于球心正下方球面处可知：绳索与竖直方向的夹角为45度。每种降落伞的费用由三部分组成，所以第i种降落伞的费用为：，又每种降落伞选用的个数为，目标函数为,(i=1,2,3,4,5).=约束条件为：此问题的关键在于求每种降落伞的最大载重质量。而最大载重质量与空气阻力系数k有关，归根结底，想要求得目标函数的最优解必须先得求出空气阻力系数k的值。一 求解空气阻力系数k。对物体做受力分析,物体受重力mg和空气阻力f，0fmg 物体在这两个力的作用下以初速度m/s，（1）向下运动，由牛顿第二定律知：，（2）又题中知，降落伞在降落过程中受到的空气阻力，可以认为与降落速度和伞面积的乘积成正比，所以，（3）又由牛顿第二定律的微分形式得：（4）由（1）（2）（3）（4）得：应用MATLAB可求得此微分方程：syms g m k pi r v;v=dsolve(Dv=g-2*k*pi*r2*v/m,v(0)=0,t) 结果：v =1/2*g/k/pi/r2*m-1/2*exp(-2*k*pi*r2/m*t)*g/k/pi/r2*m即。由速度的微分定义知：，所以利用MATLAB可求得syms t;z=int(1/2*g/k/pi/r2*m-1/2*exp(-2*k*pi*r2/m*t)*g/k/pi/r2*m);N,D=numden(z)N =g*m*(2*k*pi*r2*t+m*exp(-2*k*pi*r2/m*t);D =4*k2*pi2*r4;H=N/D所以将m=300,g=9.8,r=3代入得：利用Origin进行H与t的非线性拟合，可求得：k=2.945.k求出了可进一步通过分析求得降落伞的最大载重质量。当降落伞的半径为r时，最大载重质量为m。由可得v随着m的增大而增大，由反函数的性质可知若m是v的函数，则m随着v的增大而增大。而在此题中最大的落地速度为20。因此v=20所对应的质量就是降落伞的最大载重质量。因为v中含有t所以无法解出m。此时将v和H联立消去t之后，可得：syms g m k pi r v;v=dsolve(Dv=g-2*k*pi*r2*v/m,v(0)=0,t);a=solve(v=1/2*g/k/pi/r2*m-1/2*exp(-2*k*pi*r2/m*t)*g/k/pi/r2*m,t);syms t;z=int(1/2*g/k/pi/r2*m-1/2*exp(-2*k*pi*r2/m*t)*g/k/pi/r2*mV);N,D=numden(z);N =g*m*(2*k*pi*r2*t+m*exp(-2*k*pi*r2/m*t);D =4*k2*pi2*r4;H=N/D;subs替换。(H,t,a)利用MATLAB可解得：将H=500，v=20，k=2.945代入可得solve(500=1/4*9.8*m*(-log(-(2*20*2.945*3.14*r2-9.8*m)/9.8/m)*m-(2*20*2.945*3.14*r2-9.8*m)/9.8)/2.9452/3.142/r4,m)solve命令主要是用来求解代数方程（即多项式）的解。ans =37.744262265421074325678782244657*r2所以。由此代入不同的即可得与之相对应的最大载重质量利用MATLAB可求得:或者现在空气阻力系数k，和各种降落伞的最大载重质量均已求得。接下来的工作就是求线性规划的最优解问题了。目标函数：min z=s.t.利用Lingo编程min=736.529*x1+909.411*x2+1157.294*x3+1535.176*x4+1943.058*x5;150.976*x1+235.900*x2+339.696*x3+462.364*x4+603.904*x5=2000;gin边界限定函数：bin(x) 限制x为0或1 bnd(L,x,U) 限制LxUfree(x) 取消对变量x的默认下界为0的限制，即x可以取任意实数gin(x) 限制x为整数(x1);gin(x2);gin(x3);gin(x4);gin(x5);Global optimal solution found. Objective value: 6578.586 Extended solver steps: 0 Total solver iterations: 0 Variable Value Reduced Cost X1 0.000000 736.5290 X2 0.000000 909.4110 X3 1.000000 1157.294 X4 1.000000 1535.176 X5 2.000000 1943.058 Row Slack or Surplus Dual Price 1 6578.586 -1.000000 2 9.868000 0.000000可求得：=0，=1，=2。此时求得最优解：z=6578.586元。即用半径为3m，3.5m的降落伞各一个，用半径为4m的降落伞2个，既能满足空投需要，又能使总费用最小。六，模型的应用与推广在实际生活中，空投已经成为一种比较重要的物资救援方式了，因此降落伞的选择问题具有比较重要的实际意义，此模型能够解决在不同的高度投下，具有不同的落地速度要求的降落伞的选择问题，只需要改变H和v的值即可。进一步，当降落伞的半径连续时，此问题就要先对降落伞的半径和费用进行数据拟合，根据已知数据拟合出降落伞的半径和费用的函数关系则此问题又变为了非线性优化问题，利用MATLAB一样可以求得目标函数的最优解。七，模型的评价与改进此模型比较简单，求解也