1、武警部队大规模抗洪抢险中水的输送问题【摘要】 随着自然灾害的频繁发生，武警部队的职责使命越来越重，肩负着维护社会稳定，保障人民安居乐业和财产安全的任务。因此，灾区救援尤显重要。但是在救灾任务中如何确保救灾水的顺利输送，需要我们用科学的方式，合理的统筹安排，搞好水源输送问题。本论文就将输送水源这一问题进行研究，对抗洪抢险中水的输送这一问题进行合理的假设以及简化，建立相应模型。之后，我们使用 Matlab 对该典型线性规划进行了求解与结果分析。结论显示，引水管理费的差异是导致获利大小的关键因素。最后，本文对该模型还可引入的影响条件进行了改进讨论，并换用 LINGO 对结果进行了验证。关键词：自来水
2、输送问题 数学规划 线性规划 LP Matlab一、问题重述某市有甲、乙、丙、丁四个受灾区，由 A、B、C 三个分队对灾区输送水。四个灾区每天必须的基本生活用水分别为 30、70、10 、10 千吨，但三个送水分队每天最多只能分别送 50、 60、50 千吨水。由于地理位置的差别，往各灾区送水过程中所需要的兵力不同（如表，其中 C 水库与丁区间无输水管道） ，其它管理费均为 450 元/千吨。各区用户每千吨收费 900 元。此外，各区用户都向公司申请了额外用水量，分别为每天 50、70 、20、40 千吨。问三个分队应如何分配兵力输送水，才能用最少的兵力在最短时间内将水送到灾区？引水管理费(元
3、/千吨) 甲 乙 丙 丁A 160 130 220 170B 140 130 190 150C 190 230 200 /表 1.1二、问题假设（一）输送到各区的自来水只要在基本用水与额外用水量以内，各区即全额付费。三、符号说明1. x1,x2,x3,x4,y1,y2,y3,y4,z1,z2,z3：各水库向各居民区的供水量（详见表 1.2）2. u1,u2,u3：公司从 A、B、 C 的获利3. u：公司的总获利四、问题分析、模型的建立与求解1.问题的分析该问题为典型的数学规划问题，决策变量、目标函数都较为明显，求解过程较为简单。2.模型的建立设 A、B、C 各分队向甲、乙、丙、丁四个受灾区的
4、供水量如下，供水量（千吨） 甲 乙 丙 丁A x1 x2 x3 x4B y1 y2 y3 y4C z1 z2 z3 /表 1.2则公司从 A 水库的获利为：1 =900(1+2+3+4)(160+450)1(130+450)2(220+450)3（ 170+450）4公司从 B 水库的获利为：2 =900(1+2+3+4)(140+450)1(130+450)2(190+450)3（ 150+450）4公司从 C 水库的获利为：3=900(1+2+3)(190+450)1(200+450)2(230+450)3公司的总获利为：=1+2+3限定条件如下，各区每天的供水量：甲区： 503130zy
5、x乙区： 727z丙区： 01310zyx丁区： 44水库每天供水量的限定：A 水库：4=1=50B 水库：4=1=60C 水库： 3150iz3.模型的求解合并 u1,u2,u3 三式，得到总的目标函数：20z35260z1y463y8xx9uMax限定条件为： 3zyx70270131zyx40404=1=504=1=603150iz用 Matlab 写出线性规划程序求解（源程序详见附录） 。因 A 矩阵，b 矩阵的对应不等式为大于关系，为化为标准形式，故在 linprog 函数中 A,b 前加入负号。且 linprog 函数默认求解的是线性规划模型的标准形式，即最小量。故在取值范围允许的
6、情况下，在 f 矩阵前加负号，以求得负最小值。最终结果 fval 取相反数后即为所得结果。4.结果分析求解的结果如下：各输送管道的供水量：供水量（千吨） 甲 乙 丙 丁A 0 50 0 0B 0 50 0 10C 40 0 10 /表 1.3最大利润为 47600 元。对该供水量规划与表 1.1 的引水管理费数据作对比，不难发现获得供水量的规划方案都集中在引水管理费最低的运输途径上。由此可以得出初步结论，即引水管理费的差异是引起规划结果变化的关键因素，在满足供水需求的基础上合理规划各管道供水量，使管理费降到最低，即可获得最大利润。模型求解的结果亦支持了这一结论。五、模型的评价优点该模型较为简单
7、明了，易于理解，原理清晰。缺点程序的可重用性差。六、模型的推广与改进在实际的供水问题中，各居民区的基本用水量在人口一定的情况下是满足小范围浮动规律的，但在满足基本供水的基础之上，各居民区的额外用水量亦应获得一定比例的满足。额外用水量的调控可简单地以各居民区申请额外用水量的比例来加以满足，亦可引入新的额外需求量浮动数据。另外，该模型未考虑在受不同的日期影响下，基本用水量的起伏情况（如节假日用水量的激增，与离乡务工潮等周期性的大规模人口迁移来临后用水量的变化情况） ，以及各区域水价差异的调控政策等因素，故还有较大的改进空间。七、附录matlab 源程序：clear,clcf=290,320,230
8、,280,310,320,260,300,260,250,220A=1 0 0 0 1 0 0 0 1 0 0;0 1 0 0 0 1 0 0 0 1 0;0 0 1 0 0 0 1 0 0 0 1;0 0 0 1 0 0 0 1 0 0 0;1 0 0 0 1 0 0 0 1 0 0;0 1 0 0 0 1 0 0 0 1 0;0 0 1 0 0 0 1 0 0 0 1;0 0 0 1 0 0 0 1 0 0 0b=30;70;10;10;-80;-140;-30;-50Aeq=1 1 1 1 0 0 0 0 0 0 0;0 0 0 0 1 1 1 1 0 0 0;0 0 0 0 0 0 0
9、 0 1 1 1beq=50;60;50ub=50 50 50 50 60 60 60 60 50 50 50lb=zeros(11,1)x,fval,exitflag,output,lambda=linprog(-f,-A,-b,Aeq,beq,lb,ub) f =290 320 230 280 310 320 260 300 260 250 220A =1 0 0 0 1 0 0 0 1 0 00 1 0 0 0 1 0 0 0 1 00 0 1 0 0 0 1 0 0 0 10 0 0 1 0 0 0 1 0 0 01 0 0 0 1 0 0 0 1 0 00 1 0 0 0 1 0 0
10、 0 1 00 0 1 0 0 0 1 0 0 0 10 0 0 1 0 0 0 1 0 0 0b =30701010-80-140-30-50Aeq =1 1 1 1 0 0 0 0 0 0 00 0 0 0 1 1 1 1 0 0 00 0 0 0 0 0 0 0 1 1 1beq =506050ub =50 50 50 50 60 60 60 60 50 50 50lb =00000000000Optimization terminated.x =0.000050.00000.00000.00000.000050.00000.000010.000040.00000.000010.0000fval =-4.7600e+004exitflag =1output = iterations: 8algorithm: large-scale: interior pointcgiterations: 0message: Optimization terminated.lambda = ineqlin: 8x1 doubleeqlin: 3x1 doubleupper: 11x1 doublelower: 11x1 double 【参考文献】1 姜启源,谢金星,叶俊. 数学模型 M.北京: 高等教育出版社,2003.82 曾建军 MATLAB 语言与数学建模M
