降落伞选购的优化模型精讲.doc

降落伞选购的优化模型 [摘要]: 本文对降落伞的选购问题建立了一个优化模型。对所给数据得载重为、半径为的降落伞从高空下降时的运动曲线，发现降落伞后期匀速直线运动。降落伞作匀速直线运动时，由求出空气阻力系数，落地速度为。空气阻力系数再通过，求出降落伞绳索长度，进而算出种半径的降落伞的绳索费。最后根据每种降落伞的总成本关系把问题转化为整数线性规划问题用MATLAB解购买半径的降落伞数量为6把总费用为元。 关键词: 数据拟合；整数线性规划 1．问题的提出 要向灾区空投救灾物资共，需要选购必要的降落伞。已知空投高度为，要求降落伞的落地速度不超过。降落伞的伞面为半径为(m)的半球面，用每根长(m)的共16根绳索连接的质量为(kg)的物件位于球心正下方球面处，如图所示。每个降落伞的费用由三部分组成：伞面费(元)由伞的半径决定，见表1；绳索费(元)由绳索总长度及单价4元/m来决定；固定费用为200元。降落伞下降时受到空气阻力的作用，阻力的大小可认为与降落伞速度和伞面的面积的乘积成正比。为测定阻力系数，用了半径，的降落伞从的高度作降落实验，测得各时刻(sec)的离地高度(m)，见表2。 试确定降落伞的选购方案，即共需多少个降落伞，每个降落伞的半径多大(在表1中选择)，在满足空投条件下，要使总的费用最低。 表1 2.5 3 3.5 4 (元) 65 170 350 660 1000 表2 0 3 6 9 12 15 18 21 24 27 30 500 470 425 372 317 264 215 160 108 55 1 2．模型的基本假设 1．降落伞的质量和绳索质量均忽略; 2．降落伞下降过程中只受重力及空气阻力的作用，其他因素忽略； 3．空气阻力与降落伞的伞面面积成正比，与下降速度也成正比，空气阻力系数是常数； 4．降落伞的落地速度不允许超过20m/sec； 5．每个降落伞载的物重都不会超过降落过程中的允许的最大载重； 6．救灾物资的总重量2000kg可以任意分拆处理。 3. 符号的约定 降落伞的载重的质量，单位。 重力加速度，单位。。 空气阻力 空气的阻力系数 降落伞伞面的面积，单位。 降落伞的伞面半径，单位。 降落伞的下降速度函数，单位。 降落伞的加速度函数，单位。 运动时刻，单位。为初始时刻。 降落伞的降落速度的上界，单位。 降落伞的高度。 空投高度，单位单位。。 降落伞的高度。。 降落伞落地时刻，单位。 降落伞落地速度，单位。 降落伞允许最大落地速度，单位。。 降落伞最大允许载重，单位。 一个降落伞的伞面费 一个降落伞的绳索费 一个降落伞的固定费 一套降落伞的总费用，单位元。 各种降落伞的个数。 空投的载重总量，单位。 降落伞购买总费用，单位元。 4．模型的建立与求解 4.1 降落伞降落的运动规律 降落伞降落运动时，注意到假设1，只考虑载重物的重量，可将问题简化为质点降落问题。 设质点（降落伞的载重）的质量为。则运动中质点所受的力为，其中为阻力，为重力加速度。又根据假设3，，为降落伞伞面的面积，为下降速度。若降落伞伞面半径为，则。从而质点所受的力为。 由牛顿第三定律，则有或，其中为质点的运动加速度。所以，得到速度函数需满足的运动方程或。 记 ， （1） 则降落速度函数满足方程 。 （2） 注意到初始（时）降落速度为零，即。所以满足微分方程初值问题： 。（3）是个，可解得降落伞的降落速度函数 。 （） 因（），故。 又因 ， 所以，当相当大时， ，且。 （） 于是 。 （） 这就是说，。 为时刻时降落伞降落高度。因，所以由公式（4）有 。 即 。 （7） 因（），故降落伞降落高度函数是单增函数。 又