泄洪设施修建计划2文档格式.docx

泄洪设施修建计划2文档格式.docx

《泄洪设施修建计划2文档格式.docx》由会员分享，可在线阅读，更多相关《泄洪设施修建计划2文档格式.docx（42页珍藏版）》请在冰点文库上搜索。

分析长此以往，他在各村留宿的概率分布是否稳定？

问题4：

你们是否能够为该乡提出一个更加合理的解决泄洪的办法？

说明：

1、以上问题必须建立一般的数学模型，不能仅按照题目中提供的数据计算一个结果。

2、建模过程中，可自行提出合理的模型假设。

3

数学建模第三次论文

队号：

198 题目：

A

题 队员：

孙团伟

杨刚 张端正

2010年8月17日

摘要

修建泄洪设施时，如何在满足安全泄洪的情况下使总费用最少，是政府部门十分关注的问题。

针对本题提出的如何修建泄洪河道使总费用最省以及维护人员在各村留宿的概率的问题，分别建立了线性规划模型、最小生成树模型、马氏链模型、基于天然河道的最优泄洪基本模型,并运用MATLAB7.0和LINGO8.0数学软件，对模型进行求解，得出修建河道的最省方案和维护人员在各村留宿的概率。

最后还对原来建立的模型进行了评价，并加以推广。

问题一，开挖排洪沟和修建新泄洪河道计划按照以下五步骤进行求解（见图

1）：

1、拟合出四条天然河道的泄洪量与时间的关系式；

2、根据拟合的关系式，预测未来五年四条天然河道的泄洪量；

3、确定开挖排洪沟的费用和预计的可泄洪量以及要达到的泄洪要求；

4、建立线性规划模型，进行求解；

5、检验第

四、第五年是否满足泄洪要求。

计算得到整个方案的总开支最省为172万元。

问题二，在最短路径和泄洪量的要求下，建立了最小生成树模型。

先编写MATLAB程序代码（见附录17）,通过prim算法求出最小生成树，使得泄洪量达到最小。

然后根据最短路径的要求，再进行修改得到最优网络图（见图6）。

通过计算最终得出修建新泄洪河道网络最省总花费资金为593万元。

问题三，维护人员是在问题二中解得的新泄洪河道网络上移动的，从一个村移动到与之相连的一个村，符合马氏链，所以建立了马氏链模型。

通过分析得出，该马氏链是正则链。

根据正则链的性质可知，正则链存在唯一的极限状态概率，所以维护人员在各村留宿的概率分布是稳定的。

联立关系式（见P,求15③式）解出维护人员在各村留宿的稳态概率wi（见表3）。

问题四，考虑到天然河道有着很大的泄洪潜力，如果能够及时的对天然河道清淤以及在天然河道的基础上建立泄洪河道网络，可以显著地减少泄洪工程的花费，所以建立了基于天然河道的最优泄洪基本模型。

首先找出模型的影响因素：

天然河道的曲折；

不同地方土质的不同造成修建泄洪河道的成本不同；

有的村可能没有天然河道流过；

天然河道还要定期进行清淤；

根据地貌，修建人工湖或者水库。

这些因素都会对泄洪方案的花费造成影响。

因此，以这些影响因素建立约束函数，以工程的总花费建立目标函数（见P。

对这些影响因素进行17关系式④）

分别讨论，找出约束函数?

1、？

2、？

3、？

4最小和?

5最大的修建方案，此时目标函数z取最小值，则此时的方案是最优的，也是最合理的。

关键词：

线性规划模型最小生成树马氏链最优泄洪基本模型

1问题重述

位于我国南方的某个偏远贫困乡，地处山区，一旦遇到暴雨，经常发生洪涝5

灾害。

造成大面积水灾，不仅夏粮颗粒无收，而且严重危害到当地群众的生命财产安全。

要求通过数学建模方法，解决以下问题：

（1）该乡的某个村区域内原有四条天然河流，由于泥沙沉积，其泄洪能力逐年减弱。

在附录1中给出它们在近年来的可泄洪量（万立方米/小时）粗略统计数字。

由于它们的地质构造、长度不同，因而开挖的费用和预计的可泄洪量也不同（详见附录2）,而且预计每条排洪沟的可泄洪量还会以平均每年10%左右的速率减少。

同时开始修建一段20公里长的新泄洪河道。

修建工程从开工到完成需要三年时间，且每年投资修建的费用为万元的整数倍。

要求完成之后，通过新泄洪河道能够达到可泄洪量100万立方米/小时的泄洪能力。

乡政府从2010年开始，连续三年，每年最多可提供60万元用于该村开挖排洪沟和修建新泄洪河道，为了保证该村从2010至2014年这五年间每年分别能至少达到可泄洪量150、160、170、180、190万立方米/小时的泄洪能力，请作出一个从2010年起三年的开挖排洪沟和修建新泄洪河道计划，以使整个方案的总开支出尽量节省（不考虑利息的因素在内）。

（2）该乡共有10个村,分别标记为①一⑩，它们大致的相对地理位置（见附录3）,海拔高度总体上呈自西向东逐渐降低的态势。

根据附录4中的数据，为该乡提供一个各村之间修建新泄洪河道网络的合理方案，使得总费用尽量节省。

从村A一村B的新泄洪河道，一般要求能够承载村A及上游新泄洪河道的泄洪量。

（3）新泄洪河道网络铺设完成后，打算安排一位维护人员，每天可以从一个村到与之直接有新泄洪河道连接的相邻村进行设施维护工作，并在到达的村留宿，次日再随机地选择一个与该村直接有新泄洪河道连接的相邻村进行维护工作。

试分析长此以往，他在各村留宿的概率分布是否稳定？

（4）是否能够为该乡提出一个更加合理的解决泄洪的办法？

2基本假设

6

1、新泄洪河道在三年完工后，才能用于泄洪；

2、不受人力、物力等因素的影响，修建计划能够在指定的时间内完成；

3、在天然河道的泄洪能力非常小的情况下，认为其泄洪量为零。

3符号说明

xi：

第i条排洪沟的开挖情况（当xi=l时，表示该条排洪沟需要开挖；

当xi=0时，表示该条排洪沟不需要开挖）；

mi：

第i条排洪沟开挖费用;

qi：

第i条排洪沟当年泄洪量；

Rj：

第j年4条天然河道泄洪总泄洪量；

Pj：

第j年修建新泄洪河道的费用；

M：

每年的流动资金60万；

Pij：

维护员从村子i到相邻村子j的概率（叮?

1,23?

8940,且i?

j）；

wi：

在i村留宿的概率（i?

l,2,7,10）

4开挖排洪沟和修建新泄洪河道计划

4.1问题分析

整个开挖排洪沟和修建新泄洪河道计划应包括五个环节：

工、拟合出四条天然河道的泄洪量与时间的关系式；

2、根据拟合的关系式，预测未来五年四条天然河道的泄洪量；

3、确定开挖排洪沟的费用和预计的可泄洪量以及要达到的泄洪要求；

4、建立线性规划的数学模型，进行求解；

5、检验第四、第五年是否满足泄洪要求。

可用下面的流程图（图1）来表示，其中环节1、2、4、5是本文要做的工作，即应用拟合方法，预测出未来五年四条天然河道的泄洪量，建立数学模型进行求解出最佳方案，对第四、第五年进行检验是否满足泄洪要求。

7

图1开挖排洪沟和修建新泄洪河道计划的流程图

4.2拟合并预测未来五年四条天然河道的泄洪量

（1）由附录一中表格所给的数据对每条河道进行拟合，得到各条河道泄洪量随时间变化的关系式，然后对各河道未来五年的泄洪量进行预测。

当天然河道的泄洪量很低时，认为其泄洪量为零。

1号天然河道泄洪量随时间的变化近似为直线关系，所以使用MATLAB7.0对天然河道1号进行一次拟合（源程序见附录5）得到：

P=-1.2017 33.4639

Y=32.2622 31.0606 29.8589 28.6572 27.4556 26.2539

25.0522 23.850622.6489

拟合的图像为：

图21号天然河道泄洪量随时间变化的拟合图

则拟合的直线方程为：

y?

?

1.2017x?

33.4639

使用MATLAB7.0对未来五年该天然河道的泄洪量进行预测（源程序见附录

6）,得到：

y=21.4469 20.2452 19.0435 17.8418 16.6401

则未来五年1号天然河道的泄洪量分别为：

21.4469万立方米/小时、20.2452万立方米/小时、19.0435万立方米/小时、17.8418万立方米/小时、16.6401万立方米/小时。

（2）2号天然河道泄洪量随时间的变化近似为三次曲线关系，所以使用MATLAB7Q对天然河道2号进行三次拟合（源程序见附录7）得到：

P=-0.0311 0.7966 -7.5771 28.2357

9

Y=21.4242 16.0197 11.8357

8.6859

6.3840

4.7435

3.5781

2.7015

1.9273

图32号天然河道泄洪量随时间变化的拟合图

0.0311x3?

0.7966x2?

7.5771x?

28.2357

8）,得到：

y=1.0247 -0.1179 -1.7199 -3.9679 -7.048

从计算结果中可以看出，第二年以后结果为负数，则第二年以后2号天然河道的泄洪量已经变得很小，可以认为其泄洪量为0

则未来五年2号天然河道的泄洪量分别为：

1.0247万立方米/小时、0、

0、0、0o

（3）3号天然河道泄洪量随时间的变化近似为一次直线关系，所以使用

MATLAB7Q对天然河道2号进行一次拟合（源程序见附录9）得到：

P=-2.0850 29.9806

Y=27.8956 25.8106 23.7256 21.6406 19.5556 17.4706

15.3856 13.3006 11.2156

10

图43号天然河道泄洪量随时间变化的拟合图

2.0850x?

29.9806

10）,得到：

y=9.1306 7.0456 4.9606 2.8756 0.7906

则未来五年3号天然河道的泄洪量分别为：

9.1306万立方米/小时、

7.0456万立方米/小时、4.9606万立方米/小时、2.8756万立方米/小时、0.7906万立方米/小时。

（4）4号天然河道泄洪量随时间的变化近似为三次曲线关系，所以使用MATLAB7Q对天然河道4号进行三次拟合（源程序见附录11）得到：

P=-0.1136 2.5333-19.9078 73.0095

Y=42.4182 33.0180 26.6390 22.5994 20.2173 18.8110

17.6987 16.1985

11

图54号天然河道泄洪量随时间变化的拟合图

0.1136x3?

2.5333x2?

19.9078x?

73.0095

使用MATLAB7.0对未来五年该天然河道的泄洪量进行预测（源程序见

附录12）,得到：

y=13.6615 9.3514 2.6103 -7.2434-20.8913

由计算结果可知：

第四、第五年时，河道的泄洪量已经很小，可以认为其泄洪量为0。

则未来五年4号天然河道的泄洪量分别为：

13.6615万立方米/小时、

9.3514万立方米/小时、2.6103万立方米/小时、0、0。

（5）、根据以上的计算结果，可以列出未来五年四条天然河道的泄洪量情况：

表1未来五年四条天然河道的泄洪量（万立方米/小时）预测表

所以可以得到未来五年4条天然河道泄洪总泄洪量分别为

12

表2未来五年4条天然河道泄洪总泄洪量表

4.3线性规划模型的建立与求解

分析可知，本题研究的是，如何在资金使用最省的情况下，建立泄洪方案，达到每年的泄洪要求。

所以，应先建立线性规划模型，找出目标函数和约束函数，然后根据题目中的数据，进行求解。

再用第四、第五年的数据进行检验该方案是否满足该年的泄洪量。

影响最优泄洪方案的主要因素有两个，一是泄洪要求，一是资金花费。

泄洪的方式有：

天然河道泄洪、挖排洪沟泄洪、新泄洪河道泄洪（三年完工后，方可用于泄洪）。

费用被使用的途径：

挖排洪沟、建新泄洪河道。

则根据泄洪要求和资金花费可以建立约束函数。

而前两年是在修建的排洪沟满足当年的泄洪量要求时，剩下的资金全部用来修建新泄洪河道，只有在第三年时，修建的排洪沟满足当年的泄洪量要求时，剩下的资金部分用于修建新泄洪河道，因为新泄洪河道的泄洪量要求是确定的，所以修建新泄洪河道时，应该先按照最低的泄洪量要求进行修建，然后用第四、第五年的泄洪量要求进行检验。

若是满足，则可以按照最低的泄洪量要求进行修建新泄洪河道；

若不满足，则要增加排洪沟或者增大新泄洪河道的泄洪量。

所以建立线性规划模型中，新泄洪河道先按照最低的泄洪量来修建。

修建新泄洪河道的费用为P?

新泄洪河道能够达到可泄洪量100万立方米/小时的泄洪能力，长为20公里。

则求解可得PP138.2210万元，因为要求每年投入到修建新泄洪河道的资金为万元的整数倍，则修建新泄洪河道的钱最少为139万元。

该方案追求的是如何在达到泄洪量的要求下使资金费用最省，则如何选择开挖排洪沟是影响资金使用的关键因素。

首先建立开挖排洪沟资金使用的函数，该函数就是建立的线性规划模型的目标函数；

接着分阶段进行计算。

建立第一年的线性规划模型

目标函数：

设每年的开挖排洪沟的费用为乙则

minz?

xi?

mi?

i?

18

13

约束条件

泄洪要求：

x

18i?

qi?

150?

Rl

资金约束：

60

非约束：

O或1,当选择挖第i条排洪沟时，xi?

l；

否则xi?

O,i?

l,2,?

7,8综上可得

8?

Rl?

l?

60st?

i条升F洪沟时,xi?

否贝ljxi?

12?

78?

0或1,当选择挖第

?

第一年线性规划模型的求解：

带入数据，使用LINGO8.0求解（见附录13）可得

xl?

x3?

x6?

x7?

l

即在第一年，只需挖1、3、6、7条排洪沟。

挖排洪沟的花费为20万，即投入到建新泄洪河道的资金为40万。

建立第二年的线性规划模型

160?

R2?

60s.t.?

0?

0或1,当选择挖第i条排洪沟时，xi?

0,i?

8

第二年线性规划模型的求解：

14

带入数据，使用LINGO8.0求解（见附录14）可得，x5?

l,那么在第二年，只需挖第5条排洪沟。

挖排洪沟的花费为6万元，则有54万元被用来挖新泄洪河道。

建立第三年的线性规划模型

170?

R3?

s.t.?

x5?

x?

0 1,当选择挖第i条排洪沟时,x?

否

则x?

8ii?

i

第三年线性规划模型的求解：

带入数据，使用LINGO8.0求解（见附录15）可得，x2?

l,那么在第三年，只需挖第2条排洪沟。

挖排洪沟的花费为7万元，剩余的资金为53万元。

139-（40+54）=45

则第三年被用来修建新泄洪河道的资金为45万。

4.4模型的检验

使用LINGO8.0进行检验（见附录16）,结果输出为：

x4?

72.94?

31.96?

则能满足第四、第五年的泄洪要求，则第四、第五年不用再挖排洪沟。

所以该模型符合泄洪量的要求，该方案能达到最省钱的目的。

资金花费总额为：

139+5+7+5+6+5+5=172万元。

根据修建新泄洪河道的费用为P?

0.66Q0.51L（万元），要使费用P最小,即泄洪量Q尽量的小，且泄洪河道的长度L尽量的短。

根据prim算法得到最小生成树算法，根据图的基本定义，一个有n个

点的图，15

它的最小生成树必定含有n个点,

（n-1）条边,

运用MATLAB7.0软件

编程（源程序见附录17）

求得:

result=

根据以上的结果,

再结合从西向东流和最短路径的方法，可以得到如

下图所示的各村之间互通的新泄洪河道网络（图6）

图6.各村之间互通的新泄洪河道网络

泄洪量Q按照100万立方米/小时的倍数计算，得到:

P?

0.66Q0.51L?

0.66*[1000.51*（5?

9?

ll）?

2000.51*8?

3000.51*7?

4000.51*3?

6000.51*6]

7592.65万元

6在各村留宿的概率

6.1建立马氏链模型

维护人员从一个村到与之直接有新泄洪河道连接的相邻村进行设施维护工

作，并在到达的村留宿，次日再随机地选择一个与该村直接有新泄洪河道连接的相邻村进行维护工作，那么维护人员的转移路线就是问题二中建立的新泄洪河道网络。

为求维护人员在各村留宿的概率分布以及是否稳定，所以要建立马氏链模型。

16

则维护人员留宿的十个村就是建立的马氏链模型的十个状态，用随机变量

Xn表示维护人员所处的状态，则Xn可以取10个离散值Xn?

40,记

ai（n）?

P（Xn?

i）,即状态概率，维护人员处在村i的概率。

从Xn?

i到Xn?

j的概率记为Pij?

j|Xn?

i）,即转移概率，维护人员从村i转移

到村j的概率。

0000p=

0000

000000000

110

001/3

000

00001/30

000000

1/31/3

1/201/2001/301/300100000

01/300

01/21/201/3001/310

01

00

由马氏链的性质可知：

Xn?

l的取值只取决于Xn的取值及转移概率，而与

lzXn?

2,?

的取值无关。

由状态转移的无后效性和全概率公式可以写出马氏链的基本方程为

ai（n?

l）?

aj（n）Pji,

j?

并且ai（n）和Pij应满足

a（n）?

l,

ii?

n?

0,l,2,?

Pij?

0z

P

ij

l,j?

L2,?

」0

则状态概率向量（行向量）和转移概率矩阵

a?

al?

a2?

?

al0?

10?

则基本方程

（1）可以表示为

17

a（n?

由该递推关系式还可以得到

Pn

6.2判断该马氏链是否是正则链

正则链的定义为：

一个有k个状态的马氏链如果存在正整数N,使从任意状态i经N次转移，都以大于零的概率到达状态j（ij?

k），则这样的马氏链称为正则链。

因为修建的新泄洪河道网络连接着这十个村，当维护人员沿着新泄洪河道网络转移时，每个村都有可能到达。

即假设维护人员在村i,一定可以经过正整数N次转移到任意村j（i?

j）,即都以大于零的概率到达状态j。

由正则链的定义可以知道，本问题中建立的马氏链模型是正则链。

6.3求解极限状态概率

由定理可知，正则链存在唯一的极限状态频率w?

wl,w2,…,wk?

使得当

时状态概率a?

w,w与初始状态概率a?

无关。

w满足

wP?

w

wi?

lk ②

所以，长此以往，维护人员在各村留宿的概率就是极限状态频率

w?

wl,w2,...,w9,wl0?

则由正则链的性质可知，维护人员在各村留宿的概率是稳定的。

由②式可得：

wl,w27.../w9/wlO?

wl/w2,.../w9/wlO?

110 ③

联立可解得极限状态概率为（MATLAB源程序见附录18）：

x=0.0556 0.0556 0.1667 0.1111 0.1667 0.0556

0.1111 0.1667 0.0556