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