降落伞选择的优化模型

《降落伞选择的优化模型》由会员分享，可在线阅读，更多相关《降落伞选择的优化模型（8页珍藏版）》请在金锄头文库上搜索。

1、降落伞选择优化模型摘 要 本文针对灾区空投物资时选取降落伞的问题，建立了降落伞选择优化模型。在保证 空投任务完成的前提下， 选择不同大小的降落伞个数需满足费用最小。 （结尾这句不好） 首先，确定伞降落时的阻力系数k。根据降落伞降落的情形进行受力分析，由牛顿定 律得出关系， 建立微分方程， 利用拟合，运用数学软件 Matlab 求出阻力系数 k=3.0035。 其次，计算出不同半径的最大载重量分别为150.0726、235.8947、235.8947、462.3535、603.8903及每把伞的总费用分别为446.0193、596.2742 、821.529、1176.7838 、1562.08

2、37 。 然后，满足使用费最小的情况下，采用线性规划建立优化模型，并运用Lingo 软件 求解出只需用半径为3 的降落伞。 最后，运通 limgo22.533.540,0,6,0,0xxxxx，得出需要半径为3 的降落伞为6 个，满足使用费最小，且最小费用为4929.2。1、利用拟合，运用数学软件Matlab 求出阻力系数 k=3.0035。这句话不好！2 为150.0726、 235.8947、 235.8947、 462.3535、 603.8903及每把伞的总费用分别为446.0193、 596.2742、821.529、1176.7838、1562.0837最好用表格描述关键词： 阻力

3、系数拟合优化模型 MATLAB LINGO 1 1 问题重述 （OK） 为向灾区空投救灾物资，需选购一些降落伞。1.1 题目给出的条件1 空投物资重2000kg，空投高度500m。2 降落伞落地速度不能超过20/m s。 3 降落伞面为半径 r 的半球面，用每根长 L 共 16 根绳索连接的载重m位于球心正下 方球面处。 4 降落伞价格 （下面这个有问题！）=+cc123降落伞价格伞面费（,由r 决定）伞绳费（,由绳索总长度及单价4m/元决定） +固定费用 (=200元)c5 降落伞在降落过程中受到空气阻力与降落速度和伞面积的乘积成正比。然后用 半径 r=3m,载重 m=300kg的降落伞从

4、500m高度作降落试验，测得各时刻t 的高度 x（见 附录 A） ，从而确定阻力系数。 1.2 所求问题 试确定降落伞的选购方案，即共需多少个降落伞，每个伞的半径多大（在表1 中 选择） ，在满足空投要求的条件下，使费用最低。2 模型假设及符号说明 （OK） 2.1 模型假设 1 空投物资的 2000kg 可以任意分割； 2 假设空投物资的瞬间伞打开； 3 降落伞和绳的质量可以忽略不计； 4 降落伞的落地速度不会超过20m/s； 5 空气的阻力系数只与空气有关与其它因素无关； 6 每个降落伞载的物重都不会超过降落过程中的最大载重。 7 偶然情况忽略不计，空投情况正常。 8 g=10 2.2 符

5、号说明f空气阻力k阻力系数Mr半径为 r 的降落伞的最大载重rs半径为 r 的降落伞的伞面积( )x tt时刻降落伞的下降高度rx径为 r 的降落伞数目1C伞面费2C绳索费3C固定费用L 根绳索的长度 a降落伞的加速度g 重力加速度F合物体受的合外力3 问题分析 根据题目所给出的条件进行分析 （下面有些字母没有用公式编辑器，而且字线并未做符号说明，其他还好） 首先，由条件 4 分析得出降落伞所受的阻力与降落速度和伞面积的关系，伞在降落2 过程中做减速运动。同时利用二阶微分方程求出 x 与 t 的关系。再对 k 进行拟合，求出 k 值。然后用一阶微分方程求出v 与 t 之间的关系。由前面的两个关

6、系运用Matlab 解出 不同半径降落伞的最大载重量。 其次，由条件 3 分析得出每条伞绳的长度L 与三面半径 r 之间的关系，可以确定好 伞绳索的费用，再分别计算出不同半径降落伞的总费用，根据费用最小化，利用目标规 划求出不同半径的降落伞所需的个数。 最后，求所需的降落伞的总数量，并且确定了每个伞的半径，从而满足了空投条件 下，使用费用最低。4 建立模型及求解 4.1 求解阻力系数k 因为降落伞在降落的过程中所受的阻力与降落伞下降的速度和伞的面积的乘积成 正比，得 fkvs(1) 由受力分析可知货物受竖直空气阻力和竖直向下的重力作用，则货物在竖直方向的 合力为=Fmgf合(2)由牛顿第二定律

7、得=Fma合（3）由（1） （2） （3）联立得出mgkvsam（4）由物体位移 x和时间 t 的二次微分等于加速度建立方程得: 22(0)0(0)0d Xmgf dtm vx（5）用 Matlab 解微分方程得2 2 2 22 2kst mggmt m gXmek sksk s（t） (6) 则222222( )500kst mm gmgtm gx tek sksk s(7）(具体过程见附录 B) 题目已经 给 t-x 数据（见附录 A）对给定的数据以)(th为拟合函数进行拟合，23 ,300,10,2rm mkg gsr得出（详细过程见附录C） 3.0035k3 4.2求降落伞最大载重量

8、由速度对时间的微分等于加速度得：(0)0dvmgkvs dtm v（8）用 Matlab 解得（过程见附录D）kstmgmgmveksks（9）由（6） （7）联立建立方程组2222222( )( )2500kst mkst mmgmgv teksksmgm gm gX teksk sk ssrXx（10）k=3.0035 , g=10 , r=2 2.5 3 3.5 4因为降落伞在下落过程中其质量是不变的，所以我们把v(t)关系式中t看做一值，则关于 m的方程为v(m)=kstmgmgmeksks（11）从上式我们可以知道)(mv是关于 m的单调递增函数 : 又如果存在平衡状态则必须满足kv

9、smg， 那么ksmgv而又通过对mkst eksmgksmgtv )(分析，只有在 ksmgtvt)(时，才有，这与实际矛盾， 故降落伞是一直做加速度减小的加速运动，不存在平衡状态。因此，求最大载重量取伞在下降到地面的瞬间达到最大速度smtv/20)(, 此时( )500X tm，由方程组调用MATLAB 分别解得半径为r 的降落伞在满足空投条件下的最大载重量)(rM如下表（程序见附录E），（下面表格有问题）r （m ）2 2.5 3 3.5 4 最大承载)(rM（kg）150.9726 235.8947 339.6883 462.3535 603.8903 表 4-1 4 4.3 线性规划

10、求解降落伞数量和费用 由分析可知每种伞的单价321CCCC由分析可得2Lr，2C为: 42162rC则购买每把不同半径的降落伞所需总费用C如下: ()r m2. 2.5 3 3.5 4 1c65 170 350 660 1000 2c181.0193 226.2742 271.5290 316.7838 362.0837 3c200 200 200 200 200 总费用 （元）446.0193 596.2742 821.529 1176.7838 1562.0837 表 4-2 由上表利用目标规划求出最优解，22.533.5422.533.55min446.0193596.2742821.5

11、2911761562.0387150.072235.894339.6833462.3535603.89032000.rcxxxxxxxxxxstx目标函数：为自然数用 Lingo 软件解出22.533.540060=0xxxxx最小费用 min4929.2c（程序见附录 F ） 结果说明 共需 6 个降落伞，每个伞的半径3m,最小费用费用为4929.2 元。5 5 模型检验 1当36x时， 由表 4-1 可算出降落伞的总载重量为339.68836=2038.1kg2000kg件成立。 2 当m=339.6883, t=30s,r=3时，带入（ 7）式得v=20m/s，此时降落伞还没有落地。所以

12、条件成立。 6 模型评价及推广 本模型根据题目所给的条件综合考虑了不同半径的伞的载重量以及所需的费用，运 用 Matlab 和 Lingo 两个软件进行科学求解出模型的最优解。 此模型还可以运用到码头集装箱运送货物，以及海上救援等情况。 7 参考文献 1 吴建国等 . 数学建模案例精编 . 北京，中国水利水电出版社，2005,5. 2 孙祥 徐流美 吴清编著 MATLAB7.0 基础教程清华大学出版社2005.5 3 袁新生等 . Lingo 和Excel 在数学建模中的应用 . 北京，科学出版社， 2007. 6 8 附录 附录 A t(s) 0 3 6 9 12 15 18 21 24 2

13、7 30 x(m) 500 470 425 372 317 264 215 160 108 55 1 附录 B X=dsolve(m*D2X+k*S*DX=m*g,X(0)=0,DX(0)=0,t) 附录 C 拟合 k function F=myfun1(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); 输入下列命令：xdata=0 3 6 9 12 15 18 21 24 27 30; ydata=500 470 42

14、5 372 317 264 215 160 108 55 1 ; x0=1; x=lsqcurvefit(myfun1,x0,xdata,ydata) 附录 D 求解 V v=dsolve(m*Dv+k*S*v-m*g=0,v(0)=0,t) 附录 E 求最大承载量function F=myfun(x) r=2.5; %依次输入不同半径g=10;k=3.0035; 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

15、(-k*s*x(2)/x(1)-20; 输入以下命令：x0 = 1; 1; % 初始点options=optimset(Display,iter); % 显示输出信息x = fsolve(myfun,x0,option) %在 m文件中更改r 的值，然后在命令窗口重复输入以上命令就可分别求出不同半径的降落伞的最大载重量。附录 F 程序：7 min=446.0193*x1+596.2742*x2+821.529*x3+1176.7838*x4+1562.0837*x5; 150.9726*x1+235.8947*x2+339.6883*x3+462.3535*x4+603.8903*x5=2000; x1=0; x2=0; x3=0; x4=0; x5=0; gin(x1); gin(x2); gin(x3); gin(x4); gin (x5);