1、1运筹学模型与数学建模竞赛一、引言一般来说，大学生数学建模竞赛所涉及到的运筹学模型包括数学规划（线性规划和非线性规划） ，网络优化（含网络计划技术） ，排队模型，动态规划等，请看下表年份 题号 题名 模型分类1994 A 逢山开路 网络优化1995 A 一个飞行管理问题 数学规划1995 B 天车与冶炼炉的作业调度 网络计划技术1996 A 最优捕鱼策略 涉及数学规划1996 B 节水吸引机 涉及动态规划1997 A 零件的参数设计 数学规划1997 B 横断切割 涉及数学规划或网络优化1998 A 投资的收益和风险 数学规划1998 B 灾情巡视路线 网络优化1999 B 钻井布局 数学规划
2、2000 B 钢管的订购和运输 数学规划和网络优化2001 B 公交车调度 数学规划2002 A 车灯线光源的优化设计 涉及数学规划2003 B 露天矿生产的车辆安排 数学规划2004 A 奥运场馆的设计 涉及数学规划2004 B 电力市场的输电阻塞管理 涉及数学规划注：从 1999 年起，全国大学生数学建模竞赛开始设置专供大专院校学生做的 C，D 题。下面重点介绍运筹学模型的数学规划。二、数学规划的一般形式 )(max()inforxfubxljglihtsji ,21,0)(. 线性规划：整数规划：非线性规划：三、数学规划问题举例1 下料问题现要用 10050 厘米的板料裁剪出规格分别为
3、4040 厘米与 5020 厘米的零件，前者需要 25 件，后者需要 30 件。问如何裁剪，才能最省料？2解：先设计几个裁剪方案记 A-4040；B-5020A A/BA/B B BB B B B B注：还有别的方案吗？显然，若只用其中一个方案，都不是最省料的方法。最佳方法应是三个方案的优化组合。设方案 i 使用原材料 xi 件( i=1,2,3)。共用原材料 f 件。则根据题意，可用如下数学式子表示：)3,21(0532.min31321jxtsxj， 整 数这是一个整数线性规划模型。2 运输问题现要从两个仓库（发点）运送库存原棉来满足三个纺织厂（收点）的需要，数据如下表，试问在保证各纺织厂
4、的需求都得到满足的条件下应采取哪个运输方案，才能使总运费达到最小？（运价（元/吨）如下表）工厂 j仓库 i 1 号 2 号 3 号 库存量（吨）1 号2 号2212345030需求量（吨） 40 15 25解：题意即要确定从 i 号仓库运到 j 号工厂的原棉数量。故设 表示从 i 号仓运到ijx方案 1方案 2方案 33j 号工厂的原棉数量（吨）f 表示总运费.则运输模型为：运 输 量 非 负 约 束； 需 求 量 约 束运 出 量 受 存 量 约 束),j,i(xx.ts xxxfminij 3210543054222312211 2312一般地，对于有 m 个发点和 n 个收点的运输模型为
5、),.21;,.(0),.321(.i111njmixjbiaxtscfijmijnjijminjij其中 ai 为 i 号发点的运出量，b j 为 j 号收点的需求量，c ij 为从 i 号发点到 j 号收点的单位运价。特别当 时，存货必须全部运走，故上述约束条件中的 可改为等式：minj11 njiiax1),.21(1miaxnji3 选址问题某地区有 m 座煤矿，i # 矿每年产量为 ai 吨，现有火力发电厂一个，每年需用煤 b0 吨，每年运行的固定费用（包括折旧费，但不包括煤的运费）为 h0 元。现规划新建一个发电厂，m 座煤矿每年开采的原煤将全部供给这两个电厂发电用。现有 n 个备
6、选的厂址。若在 j#备选厂址建电厂，每年运行的固定费用为 hj 元，每吨原煤从 i# 矿运送到 j#备选厂址的运费为cij 元 (i=1,2,m, j=1,2n)。每吨原煤从 i# 矿运送到原有电厂的运费为 ci0 (i=1,2,m)。试问：1 应把新电厂厂址选在何处？2 m 座煤矿开采的原煤应如何分配给两个电厂？才能使每年的总费用（电厂运行的固定费用与原煤运费之和）为最小？4模型的建立(1) 变量的设置为了解决问题1 ，我们使用 0-1 变量 njjyj ,20#1否 则 备 选 厂 址选 中为了解决问题2，设从 i#煤矿运到 j#备选的厂址的运量为 xij 吨（i=1,2,m， j=0,1
7、,2,n）备 选 厂 址备选厂址 j煤矿 i现有电厂0 1 2 j n年产量1 1/xc1/jxc1a2 20 2j 2 m 0/mxc/mjxcma年需求量b 01iba1i总 供 应(2)目标函数的表达总运费： （对不被选中的备选厂址运费 xij,将由约束条件限制为 0）.ijminojxc1固定费用 h0+ njjy1每年总费用 z = 0110 hyxcnjjminjij (3)约束条件的表达（i）煤矿产量约束 m,iaxnji 210（ii）旧电厂用煤量约束 01bmi（iii ）新电厂用煤量约束 记 ，当 j#备选厂址被选中时 ，当01ai mijbx15j#备选厂址没被选中时 ，
8、综合表达为mijx10njybxjmij ,.211（iv）选址约束 由于只选一个厂址，所以 1njj(v)非负及整数约束 njynmixji ,20,20或综合得数学规划模型： 010100101min,.,1,.,;,.01;mnnijjijjijijmiiijjimiinjjijjzcxhyaxbynstbabyxijnj4 布点问题某市有 6 个区，每个区都可建消防站，为了节省开支，市政府希望设置的消防站最少，但必须保证在该市任何地区发生火警时，消防车能在 15 分钟内赶到现场。假定各区的消防站要建的话,就建在区的中心，根据实地测量,各区之间消防车行使的最长时间如下表:(单位:分钟)1
9、 区 2 区 3 区 4 区 5 区 6 区1 区 4 10 16 28 27 202 区 10 5 24 32 17 103 区 16 24 4 12 27 214 区 28 32 12 5 15 255 区 27 17 27 15 3 146 区 20 10 21 25 14 66请你为该市制定一个设置消防站的最节省的计划。建模并求解。解：本题实际上是要确定各个区是否要建立消防站，使其既满足要求，又最节省。这自然可引入 0-1 变量，故设 ）（区 建 消 防 站 时当 不 在 第区 建 消 防 站 时当 在 第 621 01 ,jj,xj 目标是 最少。以下考虑约束条件。61jf若 1 区
10、发生火警，按照“消防车要在 15 分钟内赶到现场”的要求，则 l 区和 2 区至少应有一个消防站，即 。同理得：121x 1 1 1 6526545434362 x,x,x,x从而得模型为： 62101652436261),j(,xx.tsfminjj再仔细观察知，若满足第一、三个约束条件，则必然满足第二、四个约束条件，故后者是多余的，可省略。从而化简得： 6210165243216),j(,xx.tsfminjj此模型当然可用软件求解，但由于比较简单,故可直接试算。若要求只有一个 ，则显1jx然不可行；若要求只有两个 ，则有唯一可行解 ，1jx 0 653142,x故这就是最优解，即只需在
11、2 区和 4 区建立消防站。附程序：c=1 1 1 1 1 17a=-1 1 0 0 0 0;0 0 1 1 0 0;0 0 0 1 1 1;0 1 0 0 1 1b=-1 1 1 1v=zeros(1,6)u=1 1 1 1 1 1x fval=linprog(c,a,b,v,u)Optimization terminated.x =0.00001.00000.00001.00000.00000.0000fval =2.00005 配套问题设有 n 个车间要生产 m 种产品，第 j 车间每天生产第 i 种产品至多 aij 件（即全天只安排生产产品 i 而不安排生产其他产品时的最大产量） ，假
12、设这 m 种产品每种一件配成一套，问如何安排生产任务才能使生产出的成套产品最多？（i =1,2,.,m； j=1,2,.,n）8（1） 建模方法（一）设 xij车间 j 安排用于生产产品 i 的数量，Z每天生产的成套产品数目，车间 j产品 i1 2 j n 全厂生产产品 I 的数量1 1x 1jx 1nx2 2 2j 2 i 1ix ijx inx全厂每天生产产品 I 的总量1ij m 1x mjx mnx车间 j每天生产产品I 的总量车间 j 每天1i生产产品 I 的总量原问题可以转化为以下数学模型： 1maxinijmfzx.0(1,;,)ijijfstxjn 模型改进为： mafz1(,
13、).0(,;1,)nijijijxmst jn 模型问题：没有将一天的生产时间约束考虑到！92、建模方法（二）设 x ij 车间 j 安排用于生产产品 i 的时间（占全天的比例）Z每天生产的成套产品数目则 aijxij车间 j 每天生产产品 i 的数目。例如：车间 2 每天至多生产某产品 6 件，若安排1/3 天时间去生产，则可产出 2 件） ，而 每天全厂产出产品 i 的总量。xaijnj1车间 j产品 i1 2 j n 全厂生产产品 I 的数量1 1xa 1jxa 1nxa2 2 2j 2 i 1ixa ijxa inxa全厂每天生产产品 I 的1nij总量 m 1xa mjxa mnxa
14、车间 j每天生产产品I 的总量车间 j1ij每天生产产品 I 的总量于是则有模型max Z ( )1minijzx(*) 整 数02121 1Z)n,.j;m,.i(x),j,.i(.tsijmijnjija其中常数 1 表示 1 天。10注： (1)此模型着重考虑安排生产的时间； (2)从实际情况考虑，安排生产的时间必须是每件产品耗用生产时间的整数倍才合适。例如，每 3 分钟生产一件产品，安排 5 分钟，也只能生产 1 件，不能认为生产了 1.67 件。模型（ *）没有考虑到这些因素，故是不合适的。（2）建模方法（二）改进（*）显然， 车间 j 生产每件产品 i 的耗用时间（天） 。从以上分
15、析应有ija(?是非负整数)axijij1?从而令 y ij=aijxij , yij 是非负整数 ,表示车间 j 每天生产产品 i 的数目，将它代入（*）后得max Z(*) 整 数整 数 021121Z)n,.j;m,.i(),.i(.tsyaijmiijijnjij这是一个整数线性规划模型。注：此模型着重考虑安排生产产品的数目。四、数学规划在数学建模中的应用举例1钻井布局勘探部门在某地区找矿。初步勘探时期已零散地在若干位置上钻井，取得了地质资料。进入系统勘探时期后，要在一个区域内按纵横等距的网格点来布置井位，进行“撒网式”全面钻探。由于钻一口井的费用很高，如果新设计的井位与原有井位重合（
16、或相当接近），便可利用旧井的地质资料，不必打这口新井。因此，应该尽量利用旧井，少打新井，以节约钻探费用。比如钻一口新井的费用为 500 万元，利用旧井资料的费用为 10 万元，则利用一口旧井就节约费用 490 万元.设平面上有 n 个点 Pi，其坐标为( ai,bi),i=1,2,n，表示已有的 n 个井位。新布置的井位是一个正方形网格 N 的所有结点（所谓“正方形网格”是指每个格子都是正方形的网格；结点是指纵线和横线的交叉点）。假定每个格子的边长（井位的纵横间距）都是 1单位（比如 100 米）。整个网格是可以在平面上任意移动的。若一个已知点 Pi与某个网格结点 Xi的距离不超过给定误差 （
17、=0.05 单位），则认为 Pi处的旧井资料可以利用，不必在结点 Xi处打新井。 为进行辅助决策，勘探部门要求我们研究如下问题： 1)假定网格的横向和纵向是固定的（比如东西向和南北向），并规定两点间的距离为其横向距离（横坐标之差绝对值）及纵向距离（纵坐标之差绝对值）的最大值。在平面上11平行移动网格 N，使可利用的旧井数尽可能大。试提供数值计算方法，并对下面的数值例子用计算机进行计算。 2)在欧氏距离的误差意义下，考虑网格的横向和纵向不固定（可以旋转）的情形，给出算法及计算结果。 3)如果有 n 口旧井，给出判定这些井均可利用的条件和算法（你可以任意选定一种距离）。 数值例子 n=12 个点的
18、坐标如下表所示：I 1 2 3 4 5 6 7 8 9 10 11 12ai 0.50 1.41 3.00 3.37 3.40 4.72 4.72 5.43 7.57 8.38 8.98 9.50bi 2.00 3.50 1.50 3.51 5.50 2.00 6.24 4.10 2.01 4.50 3.41 0.80问题重述（略）问题分析布局问题是在一定约束条件下的最优化问题，勘探部门要求尽可能多地利用旧井，因此我们必须围绕这个问题来移动网格。问题（1）：网格的横向和纵向是固定的，网格只能平移。如果考虑网格上所有的结点通过平移后与旧井位的距离，由于网格结点数较多，运算量较大，虽然可以解决问题
19、，但并非一个好的解法。因此我们从旧井位考虑，通过四舍五入法找到与其相邻的网格结点（相互关系如图所示），然后我们通过两个控制变量（角度与距离）对这些结点进行平移，使网格结点尽可能多的与旧井位重合。问题（2）：网格不固定移动，需要我们考虑平移和旋转，对于旋转我们同样从旧井位出发，通过旧井位与网格结点的相对关系，不难得出：旧井位的旋转是网格旋转的逆过程，因此我们利用旧井位旋转来替代网格的旋转问题，然后，根据问题（1）的方法，再解决网格平移问题，这样就将网格的移动问题简化成为旧井位的旋转问题和网格的平移问题。问题（3）：以上两个问题都是在正方形网格情况下进行的求解，在此问中我们不妨假设网格是长方形的，
20、通过调整长和宽的比例关系使 n 个旧井位都可以被利用。因此问题（3），就是在问题（2）的基础上，使 n 个旧井位都可以被利用的长方形的长、宽比例问题。具体实现方法如下：把目标函数值定为 12，通过对问题（2）的反向求解来确定长方形的长、宽比例。模型假设和符号说明模型假设：1 若一个旧井位与某个网格结点不超过 0.05 单位，则视为重合；2 网格是足够大且平铺的；3 模型采用直角坐标系，在网格未移动时，网格的横向和纵向就是直角坐标系的两坐标轴方向。4 逆时针旋转为正，顺时针方向为负；5 长方形网格的长是 1 单位。符号说明：ii 旧井位,i=1,2,12;Pi第 i 个旧井位;ai,bi第 i
21、个旧井位的横坐标、纵坐标;Ai,Bi第 i 个旧井位旋转的横坐标、纵坐标;Wi与第 i 个旧井位相邻的网格结点;12xi,yi与第 i 个旧井位相邻的网格结点的横坐标、纵坐标;xli,yli与第 i 个旧井位相邻的网格结点移动后的横坐标、纵坐标;di第 i 个旧井位与其相邻的网格结点距离;s网格平移长度;an网格平移角度;bn旧井位点旋转角度;pi第 i 个旧井位到坐标原点的距离;qi第 i 个旧井位与坐标原点连线的倾角;T可被利用的旧井位数;ti第 i 个旧井位数是否被利用的函数(1 表示可被利用;0 表示不可被利用);c长方形的宽与长的比例;um, tm分别为宽与长之比的上、下线.模型的建
22、立我们的目标是使旧井位利用数最多的 T, 它由 ti(d)决定，即 T= 。而由要求12)(idt知 ti(d)= (1)otherwisd05.)(01d(i)的表达式在下面给出。由于问题（1）是问题（2）的特殊情况，所以这里只给出问题（2）的模型。我们先将旧井位进行旋转。第 i 个旧井位到坐标原点的距离为p(i)= . (2)2)(iba第 i 个旧井位与坐标原点连线的倾角为q(i)=arctan(b(i)/a(i). (3)所以第 i 个旧井位旋转的横坐标、纵坐标分别为A(i)=p(i)cosq(i)+bn(i), B(i)=p(i)sinq(i)+bn(i) (4)(bn 旧井位点旋转
23、角度)。然后用四舍五入法求旧井位的相邻网格结点坐标：x(i)=A(i)+0.5, y(i)=B(i)+0.5 (5)由此得出网格结点平移后的坐标为xl(i)=x(i)+s*cos(an), yl(i)=y(i)+s*sin(an) (6)所以，网格结点与旧井位的误差距离为d(i)= (7)22)()(iByliAxl综上所述，本问题的数学模型可表示为：目标函数： maxT= 12)(idt约束条件是满足( 1)-(7)式。注：以上模型说明：当 bn=0, d(i)=max|xl(i)-A(i)|, |yl(i)-B(i)|时，目标函数的解即为问题（1）的解。13模型的求解我们的目标是使目标函数
24、最大化，这是一个有约束的优化问题，我们采用穷举法来求它的极值。所谓穷举法就是逐点确定寻优方向，然后利用固定的步长，进行搜索的方法。为使目标函数值最大，我们列出主要步骤如下：1定一个能包含所有解的初始范围与固定的搜索步；2据运行时间，再对搜索步长用穷举法进行精度调整。根据这种方法我们得到了问题答案：问题（1）：最多可利用的旧井位数 t=4, 网格的平移方向 an0=5.400000(弧度)，网格的平移距离 s0=0.583000(单位)；问题（2）：最多可利用的旧井位数 t=6, 网格的平移方向 an0=2.760000(弧度)，网格的平移距离 s0=0.230000(单位), 网格旋转角度：b
25、n0=-0.780000(弧度)（旧井点的相对旋转角度为：5.500000 弧度）。模型的改进降落伞的选择一、问题的重述：为向灾区空投一批救灾物资共 2000kg，需选购一些降落伞。已知空投高度为 500m，要求降落伞落地时的速度不能超过 20m/s，降落伞的伞面是半径为 r 的半球面，用每根长L 共 16 根绳索连接的重 m 位于球心正下方球面处，如下图：每个降落伞的价格由三部分组成，伞面费用 C1 由伞的半径 r 决定，见下表；绳索费用C2 由绳索总长度及单价 4 元每米决定，固定费用 C3 为 200 元。r 2 2.5 3 3.5 4C1 65 170 350 660 1000降落伞在
26、降落过程中除受到重力外，还受到空气的阻力，可以认为与降落的速度和伞的面积的乘积成正比。为了确定阻力系数，用半径 r=3m，载重 m=300kg 的降落伞从 500m高度作降落试验，测得各个时刻的高度 x 见下表：t(s) 0 3 6 9 12 15 18 21 24 27 30x(m) 500 470 425 372 317 264 215 160 108 55 1试确定降落伞的选购方案，即共需多少个伞，每个伞的半径多大（在给定半径的伞中选） ，在满足空投要求的条件下，使费用最低。二、问题的分析：根据题意，每种伞的价格是确定的。要确定降落伞的选购方案，即共需多少个伞，每个伞的半径多大（在给定半
27、径的伞中选） ，在满足空投要求的条件下，使费用最低。首先，必须知道每种伞在满足空投要求的条件下最大的载重量 M(r)，然后就是一个线性整数规划了。欲得到 M(r)，必须先求出空气阻力系数 k，然后根据运动方程得出 M(r)。最后运用分支定界法求解线性整数规划，得出问题要求的解。三、基本假设：1救灾物资 2000kg 可任意分割；2降落伞落地时的速度不能超过 20m/s；143降落伞与绳索的质量可以忽略；4降落伞下落过程中，只受到重力和空气阻力的作用；5空气阻力的阻力系数 k 是定值，与其他因素无关。四、符号说明：M(r)- 半径为 r 的伞在满足空投要求的条件下最大的载重量；t- 降落伞从开始
28、下降开始记时的时间；k-空气阻力系数H(t)-降落伞从降落位置到 t 时刻所下降的距离；m-降落伞负重重量；g-重力加速度；s-降落伞伞面面积；nr-选购的半径为 r 的降落伞的个数五、模型的建立与求解1、计算每种伞的单价如下：单位为元，C 2= 416r半径 r 2 2.5 3 3.5 4C1 65 170 350 660 1000C2 181.0193 226.2742 271.5290 316.7838 362.0387C3 200 200 200 200 200C 446 596.3 821.5 1176.8 1562C(取整) 446 596 822 1177 15622、 定空气阻
29、力系数 k。对给定的 r=3m，m=300kg ，取 g=9.8m/s2，s=2r 2，有数据t(s) 0 3 6 9 12 15 18 21 24 27 30x(m) 500 470 425 372 317 264 215 160 108 55 1处理得：t(s) 0 3 6 9 12 15 18 21 24 27 30H(m) 500 30 75 128 183 236 285 340 392 445 499作出 H(m)t(s)的关系图 1:15图 1从图 1 可看出：一方面，H(m)t(s)在后阶段基本是线性关系，即降落伞作匀速运动。从mg=kVs有又由于若在 500 米的物体做匀速运
30、动，则 ，5017/3svmstA代回 mg=kVs 中 ， 估算出 k=2.9。另一方面，根据题意，物体在降落过程中一直做匀速直线运动是不可能的。题中注意到：降落伞在下降过程中只受到重力和空气阻力的作用，而且初速度为 0，由牛顿第二定理知：，()dVtFmagkvs()dVtkvsgm（1）0)(kvst解得: （2）ksmgetVst)(由（2）式，并代入 k=2.9，r=3m，m=300kg，取 g=9.8m/s2，s=2r 2，作出下图 216图 2从图 2 可以看出当 时，V(t)迅速增长，但由于负项是成负指数衰减的，所以09t很快就接近极限值 mg/ks。当 时， 所以 9 秒以t
31、0,kstmkstge(),mgvtks后可以看作近似的匀速运动。下面就 9 秒以后的数据运用最小二乘法进行线性拟合。设 H(t)=Vt+b+ ，其中 服从正态分布。 Matlab 程序如下：x=9 12 15 18 21 24 27 30;H=128 183 236 285 340 392 445 499;P=polyfit(x,H,1)从而得到 p=17.5794 29.2976，所以 V=17.5794(m/s)，并由 mg=kVs 得到 k=2.9575。3、 瞬时速度与高度设降落伞从降落位置到 t 时刻所下降的距离为 H (t)，则有（3）tdVH0)(积分求得（4）22)(skgm
32、eksgttst4、求半径为 r 的降落伞在满足空投要求的条件下最大的载重量 M(r).V(t)V(t)S0 t0 tS=H(t0)= ()tv1722()kstmkstgeVttgHtks与 直接相关，也与 相关。 越大，则 越大。由（2）式知：V(m) 是关于 m 的增()vtmv函数。当 最大时， 也达到最大 ，即为最大的载重量。M特别的，在给定从 500m 的高空空投时，降落伞落地瞬间的速度在给定 g, k, s 后又有等式约束 H(t)=500，即 2250skgmeksgtst的情况下，V(m)是关于 m 的增函数；反之，其反函数也是关于 V 的增函数。所以要求半径为 r 的降落伞
33、在满足空投要求的条件下最大的载重量 M(r)，就是要在 V 取最大值时取得，即取 V=20m/s 时求出指定半径 r 的 M(r)，于是得到如下方程组： ksgetVmst)(（5）22)(sgttHt由此导出：（6）ksmVgksmt )/(1ln()( 22如前所述，取 H(t)=500，V=20m/s，得到方程：（7）ksks)/(l(5022将参数 g=9.8，k=2.9575 代入上式，有 smms 957.20)957./()8.91ln(8.922由 s=2r 2，分别代入 r=2，r=2.5；r=3；r=3.5；r=4，并调用 Matlab 的命令 solve，分别解得半径为
34、r 的降落伞在满足空投要求的条件下最大的载重量 M(r)如下：M(2)kg M(2.5)kg M(3)kg M(3.5)kg M(4)kg151.6947 237.0229 341.3130 464.5649 606.7787(取整)151 237 341 464 6063 4原问题就是如下的一个线性整数规划问题：设 nr 为选购的半径为 r 的降落伞的个数，则有 15627825964min 4.33.2 nn18s.t. （8）22.533.541741602,r rnnnN运用分支定界法可求得： ,345.3.2 总费用为 4932 元，即在需要空投 2000kg 的情况下，需要选购半径
35、为 3 米的降落伞 6 把。六、模型的检验及推广1在求解空气阻力系数 k 时，我们分析数据得出在运动后期降落伞作近似的匀速运动，并由此为前提对数据进行拟合求出了 k。下面在问题已经基本得到解决后，运用获得数据对降落伞的运动情况进行检验，如下表：r 2 2.5 3 3.5 4M(r) 151.6947 237.0229 341.3130 464.5649 606.7787s=2r 2 25.1327 39.2699 56.5487 76.9690 100.5310M/s 6.035740 6.035739 6.035739 6.035739 6.035739可以看出，M/s 几乎是常数，又有 m
36、g=kVs 是降落伞后期运动为匀速运动的充分必要条件，即 m/s=kV/g 为常数，而 k，g 为常数，所以降落伞在运动后期为近似匀速运动。因此，在求解空气阻力系数 k 时，假设后期运动为近似匀速运动是合理的。2从图 2 可以看出降落伞的大幅度加速过程很快就结束了，对给定的伞和给定的承载质量很快就进入近似匀速运动，而且速度与空投高度基本上是无关的，所以空投高度并不十分重要，只要能保证空投位置准确，高一些投放也可以，这样可降低空投高度。3题目中要求的是空投 2000kg 的救灾物资，若数据有变动，例如 3000kg，5000kg，则只需将（8）式的第一个约束右端改为 3000kg，5000kg，然后再运用分支定界法可求得结果。
