降落伞的采购问题.docVIP

降落伞的采购问题 摘要 本模型研究的是降落伞的选购方案问题，目的是在满足空投要求的条件下，使费用最少。为了方便对降落伞进行受力分析，我们把降落伞和其负载的物资看做一个整体，忽略了伞和绳子的质量，并假设降落伞只受到竖直方向上空气阻力和重力的作用。通过对降落伞在空中的受力情况的分析建立起了高度与时间的方程，然后以高度与时间的方程作为拟合曲线与题中给出的时间与高度的数据进行拟合，得出阻力系数k的值。我们建立了速度与质量的方程，并证明其为严格增函数（证明过程见附录）。由于题中已限制降落伞的最大落地速度为20m/s，所以当速度为20m/s时，伞的承载量最大。建立高度与时间，速度与时间的方程组，代入最大速度20m/s，高度500m,伞的半径（题中已给出可能选购的每种伞的半径），分别计算出每种伞的最大承载量。最后运用LINGO软件进行线性规划求解得：x1=0,x2=0,x3=6,x4=0,x5=0.即购买半径为3m的降落伞6个时总费用最少为4932元。 关键字：线性规划、空气阻力系数、拟合、Lingo软件。 一.问题重述 为空投救灾物资2000kg，需选购一些降落伞。 空投高度为500m，降落伞落地时的速度不能超过20m/s。降落伞伞面是半径为m的半球面，通过每根长m的16根绳索连接的载重位于球心正下方的球面处，如下图一所示： 图一 每个降落伞的价格由三部分组成：伞面费用由伞的半径决定，见下面附表1；绳索费用由绳索总长度及单价4$/m决定；固定费用为200$。降落伞在降落过程中受到的空气阻力可以认为与降落速度和伞面积的乘积成正比， 为了确定比例系数，用半径为m，载重为kg的降落伞从500m高度处做降落试验，测得各时刻的高度见下面附表2. 要求我们确定降落伞的最佳选购方案，即选购多少，每个降落伞的半径多大（在附表1中选取）. 附表1： （m） 2 2.5 3 3.5 4 （$） 65 170 350 660 1000 附表2： （s） 0 3 6 9 12 15 18 21 24 27 30 （m） 500 470 425 372 317 264 215 160 108 55 1 二.模型假设 (1)救灾物资总量2000可以任意分割. (2)降落伞的落地速度不超过20. (3)降落伞与绳索的质量可以忽略. (4)降落伞在降落过程中只受到重力和空气阻力的作用. (5)空气阻力系数是定值. 三.符号说明 ： 时刻时降落伞的速度 ： 时刻时降落伞的位移 ： 时刻时降落伞距离地面的高度 ：空气阻力系数 ：降落伞的面积 ：重力加速度 四、问题的分析 由题意可知每个伞的价格由三部分组成：伞面费用C1、绳索费用C2、固定费用C3。伞面费用C1由伞的半径r决定；绳索费用C2由绳索的长度及单价决定，由图一可知绳索的长度又由降落伞的半径决定即；固定费用为定值C200。因为题中已给出每种伞面的半径，所以每种伞的价格为定值。要想确定选购方案，即共需半径（在题中给出的半径中选择）为多大的伞的数量，在满足空投物资要求的条件下使总费用最少。因此，我们需要确定每种伞的最大承载量。然后进行线性规划，确定总费用和每种伞的个数。 要确定最大载重量，我们需对降落伞进行受力分析（如图二）。降落伞在降落过程中除受到竖直向下的重力作用外还受到竖直向上的空气阻力的作用，而由题可知空气阻力又与阻力系数、运动速度、伞的受力面积有关。运动速度和受力面积是已知的，所以要想确定每种伞的最大承载量，就必须先要确定空气的阻力系数。 图一 图二 对图二的分析可知降落伞的运动状态是做加速度趋近于0的加速运动。因此，我们可以建立一个位移与时间的函数关系式，在根据题中所给的数据拟合出阻力系数k的值。然后再建立一个速度与时间的函数关系式，两个关系式联立求解出最大载重量（其中高度和速度由题目已经给出）。最后用LINGO软件进行线性规划算出问题要的结果。 五、模型的建立 （1）求各种不同半径的降落伞的价格,见下表： 2 2.5 3 3.5 4 65 170 350 660 1000 181.0193 226.2742 271.5290 316.7838 362.0387 200 200 200 200 200 446.0 596.3 821.5 1176.8 1562.0 （取整） 446 596 822 1177 1562 此处各个数据的计算可利用MatLab7.1来做，如的计算： r=[2,2.5,3,3.5,4]; 64*sqrt(2)*r ans = 181.0193 226.2742 271.5290 316.7838