泄洪设施修建具体计划

1、泄洪设施修建计划摘要 针对本题提出的如何修建泄洪河道使总费用最省以及维护人员在各村留宿 的概率的问题，分别建立了非线性规划模型、 马氏链模型，并运用 matlab 和 lingo 数学软件，对模型进行求解， 得出修建河道的最省方案和维护人员在各村留宿的 概率。最后还对原来建立的模型进行了评价，并加以推广。 在考虑修建的泄洪道路径和泄洪量的情况下， 得到修建泄洪道的最省花费的 0-1 规划模型。并通过 lingo 求解得到最优泄洪道网络连接图（见图 1 ）和修建 新泄洪河道最省总花费资金： 547.0804 万元。维护人员是在问题一中解得的新泄洪河道上移动的， 从一个村移动到与之相 连的一个村，

2、符合马氏链，所以建立了马氏链模型。通过分析得出，该马氏链是 正则链。根据正则链的性质可知， 正则链存在唯一的极限状态概率， 所以维护人 员在各村留宿的概率分布是稳定的。并运用 matlab 软件编程求解出维护人员在 各村留宿的稳态概率（见表 3）。由于前面的模型仅是从修建泄洪道花费最省来建立的。 而没有考虑建设后的 维护成本。 且由于上游地势高的村庄的水流要汇入下游地势低的村庄， 从而会使 得下游泄洪道的泄洪压力增大。 有可能洪涝来临时对下游村庄带来危险。 所以从 安全和维护等因素来看， 综合该乡地势由西向东逐渐降低的地势特点。 可以考虑 在该乡中间人口相对较少地区修建一由西向东的逐渐加宽的主

3、渠道。 再由各村庄 各自修建泄洪渠道与主渠道相连，最终将洪水排出。关键词 ：lingo0-1变量线性规划模型马氏链模型 matlab、问题重述位于我国南方的某个偏远贫困乡， 地处山区，一旦遇到暴雨， 经常发生洪涝 灾害，以往下雨时，完全是依靠天然河流进行泄洪。 2010 年入夏以来，由于史 无前例的连日大雨侵袭， 加上这些天然河流泄洪不畅， 造成大面积水灾， 不仅夏 粮无收，而且严重危重到当地群众的生命财产安全。为此，乡政府打算立即着手解决防汛水利设施建设问题。 从长远考虑， 可以 通过修建新泄洪河道的办法把洪水引出到主干河流。 经测算，修建新泄洪河道的 费用为2p 3L Q （万元）其中 Q

4、表示新泄洪河道的可泄洪量（万立方米 / 小时），L 表示新泄洪河道的 长度（公里）。该乡共有 10个村，分别标记为 - ，下图给出了它们大致的相对地理位置，海拔高度总体上呈自西向东逐渐降低的态势。其中村距离主干河流最近， 且海拔高度最低。 乡政府打算拟定一个修建在 各村之间互通的新泄洪河道网络计划， 将洪水先通过新泄洪河道引入村后， 再 经村引出到主干河流。 要求完成之后， 每个村通过新泄洪渠道能够达到可泄洪 量 100 万立方米 / 小时以上的泄洪能力。表 1 各村之间修建新泄洪河道的距离（单位：公里）2 3 4 5 6 7 8 9 1017 4 8 11 13 12 16 17 2229

5、14 16 8 11 18 14 2337 9 11 7 12 12 1744 17 10 7 15 1858 10 6 15 1569 16 8 1578 6 11813 11912请通过数学建模的方法，解决以下问题：问题 1：根据表 1 数据，为该乡提供一个各村之间修建新泄洪河道网络的合 理方案，使得总费用尽量最省。 （提示：从村 A 村 B 的新泄洪河道，一般要求 能够承载村 A 及上游新泄洪河道的泄洪量） 。问题 2：新泄洪河道网络铺设完成后，打算安排一位维护人员，每天可以从 一个村到与之直接有新泄洪河道连接的相邻村进行设施维护工作， 并在到达的村1 / 11 留宿，次日再随机地选择一

6、个与该村直接有新泄洪河道连接的相邻村进行维护工作。试分析长此以往，他在各村留宿的概率分布是否稳定？问题 3：是否能够为该乡提出一个更加合理的修建新泄洪河道的办法？二、问题分析针对问题一，要求使得 总费用尽量最省，而修建新泄洪河道的费 用2p 2 L Q(万元)，由此可知费用与新泄洪河道的可泄洪量 Q 和泄洪河道的长3 度 L 有关，要使费用 P最小，即泄洪量 Q尽量的小，且泄洪河道的长度 L 尽量的 短。对于此，可以运用 0-1 变量建立规划模型， 进而运用 lingo 软件编程求出最针对问题二， 维护人员是在问题一中解得的新泄洪河道网络上移动的， 从一 个村移动到与之相连的一个村， 符合马氏

7、链， 所以建立了马氏链模型。 通过分析 得出，该马氏链是正则链。 根据正则链的性质可知， 正则链存在唯一的极限状态 概率，所以维护人员在各村留宿的概率分布是稳定的。再用 matlab 软件编程求 解出维护人员在各村留宿的稳态概率。针对问题三， 由于上游地势高的村庄的洪水要汇入下游地势低的村庄， 从而 会使得下游泄洪道的泄洪压力增大。 有可能洪涝来临时对下游村庄带来危险。 且 由于渠道较多后期维护较难。 所以从安全和维护等因素来看， 综合该乡地势由西 向东逐渐降低的地势特点。 可以考虑修建一由西向东的逐渐加宽的主渠道。 最终 将洪水排出。三、问题假设1、村子 1-10 的海拔高度自西向东递减。2

8、、A-B 的泄洪河道，其中 B村的泄洪河道能够承载村 A泄洪量及上游所有 流入 A 的泄洪量。3、若泄洪河道相交，假设互不影响各自泄洪量。4、假设维修人员选择第一个村庄是随机的，且概率是相同的。四、符号说明Qij：第 i 村庄到第 j 村庄泄洪河道的流量。Lij ：i 村到 j 村的距离0,表示 i j不修河道Xij 1,表示i j 修建河道p ：修建河道总费用Pij ：维护员从村子 i 到相邻村子 j 的概率( i, j 1,2,3, ,8,9,10,且i j )；2 / 11wi ：在 i 村留宿的概率( i 1,2, ,10 )X n ：表示维护人员所处的状态( Xn可以取 10 个离散

9、值 Xn 1,2, ,10 ) ai (n) P(X n i) ：状态概率，维护人员处在村 i 的概率( i 1,2, ,10 )Pij P(Xn 1 j| Xn i) ：转移概率，维护人员从村 i转移到村 j 的概率 (i, j 1,2,3, ,8,9,10,且i j)P ：转移概率矩阵Qi ：各段主泄洪道的泄洪量( 100-900)。Li ：各段主泄洪道的长度。Qj ：各村的泄洪量。Lj ：各村到主泄洪道的泄洪道长度。五、模型的建立与求解5.1 ：根据表 1 数据，为该乡提供一个各村之间修建新泄洪河道网络的合理方案， 使得总费用尽量最省。根据表 1 数据，为该乡提供一个各村之间修建新泄洪河

10、道网络的合理方案， 使得 总费用尽量最省。其中为了表述的方便将村庄做如下编号即：3215796 1048于是，表一数据整理为：3 / 11表22 3 4 5 6 7 8 9 10 18 13 8 11 9 17 15 8 1627 14 9 11 14 23 16 18317 4 12 8 22 11 16412 6 15 12 15 1357 7 17 9 12610 11 10 8718 4 7815 1196为使费用最小，依据题意引入 0-1 变量，其中 0 表述不修河道， 1 表示要修建河道。故建立如下优化模型：minp23XijLij Qiji=1 j=i+1 3Qij (0 1 i

11、 j 10)i-110s.t.Qki +100Q(ij i=2 9)k=1j=i+110Q1j 100j=1由 lingo 软件可以得到（程序见附录 1 ）：X23 1,X35 1,X56 1,X46 1,X19 1,X79 1, X 9,10 1, X6,10 1, X8,10 1 则：Q21 100,Q13 200,Q37 300,Q97 100,Q78 500,Q58 300,Q10,8 100,Q45 100,Q65 100L21 7,L13 4,L37 7,L97 6,L78 8,L58 6, L10,8 11, L45 4,L65 8由此可绘出下图：4 / 11图 1. 各村之间互

12、通的新泄洪河道网络 可得出最小费用为：min p 547.0804万元5.2.1 ：维护人员的转移路线就是问题一中建立的新泄洪河道网络。为求维护人 员在各村留宿的概率分布以及是否稳定，建立马氏链模型。由问题一得出的新泄洪河道修建方案可知转移概率矩阵为：01/21/2000000010000000001/ 2000001/20000000100000P0001/301/ 301/3000000100000001/300001/31/3000001/ 301/3001/300000010000000000100由马氏链的性质可知： Xn 1的取值只取决于 Xn 的取值及转移概率，而与Xn 1,Xn

13、 2 , 的取值无关。由状态转移的无后效性和全概率公式可以写出马氏 链的基本方程为10ai(n 1)a j (n)Pji , i 1,2, ,10 j1并且5 / 1110ai (n) 1, n 0,1,2, i1s.t. Pij 0, i, j 1,2, ,1010Pij 1, j 1,2, ,10则状态概率向量(行向量)和转移概率矩阵a n a1 n ,a2 n , ,a10 n01/ 21/ 21001/ 200000P000000001/ 3000000000则基本方程( 1)可以表示为000000000000000001/ 200001000001/ 301/ 301/ 300010

14、000000001/ 31/ 3001/ 301/ 3001/ 300010000000100a(n 1) a n P由该递推关系式还可以得到a(n) a0 Pn5.2 3判断该马氏链是否是正则链正则链的定义为：一个有 k个状态的马氏链如果存在正整数 N ，使从任意 状态i经N 次转移，都以大于零的概率到达状态 j(i, j 1,2, , k ) ，则这样的马 氏链称为正则链。因为修建的新泄洪河道网络连接着这十个村，当维护人员沿着新泄洪河道 网络转移时，每个村都有可能到达。即假设维护人员在村 i ，一定可以经过正整 数N 次转移到任意村 j(i j),即都以大于零的概率到达状态 j。由正则链的

15、定义可以知道，本问题中建立的马氏链模型是正则链。5.2.4 求解极限状态概率由定理可知，正则链存在唯一的极限状态频率 w w1 , w2 ,.,wk ，使得当n 时状态概率 a n w ， w与初始状态概率 a 0 无关。 w 满足6 / 11wP wk wi 1 i1所以，长此以往，维护人员在各村留宿的概率就是极限状态频率w w1,w2,.,w9 ,w10 ，则由正则链的性质可知，维护人员在各村留宿的概率是 稳定的。由式可得：w1 , w2 ,.,w9 , w10 P w1 , w2 ,., w9 , w10 10wi 1i1联立可解得极限状态概率为（ matlab 程序见附录 2 ）：w

16、0.1111 0.0556 0.1111 0.0556 0.1667 0.0556 0.1667 0.1667 0.0556 0.0556表 3 即可得维护人员随机在各村留宿的稳态概率w1w2w3w4w51/91/181/91/183/18w6w7w8w9w101/183/183 /181/181/185.3 ：由于该乡地势由西向东逐渐降低。且各村庄由西向东分布。故可以在其中 间人口分布密度较少区建一由西向东的主泄洪道。泄洪道泄洪量依次递增。从 100-900。最终经村庄 8 注入主干河流。 各村再各自修建泄洪道与主泄洪道相连 各村最大泄洪量用 Q 表示。 L 为主泄洪道的长度。由几段构成。泄

17、洪量 100-900。目标函数：min p 2/3* Qi *Li2/3* Qj *Lj六、模型的评价和优化模型的优点：1、用 0-1 变量刻画村庄之间的河道修理与否，结合线性规划模型简单易懂用 lingo 软件求解也较方便。2、问题二建立的马氏链模型， 即符合题意也具有理论依据用 matlab 求解其程序 也比较简便。3、该模型可运用性强，不仅仅运用于河道的修理也可以用于电缆的铺设、道路7 / 11的修建等等类似的模型。 模型的缺点： 1、该模型模型考虑的因素还不多，如有些地方可能因为地貌不能挖泄洪河道、有的地方河道没法加宽、以及修建排洪沟这样会使建立的模型对实际产 生较大的误差。七、参考文

18、献1 孙祥 徐流美 吴清，matlab7.0 基础教程，清华大学出版社， 2005 年.2 刘卫国， MATLAB程序设计与应用（第二版） ，北京：高等教育出版社， 2006 年。3 姜启源 谢金星 叶俊，数学模型（第三版） ，北京：高等教育出版社， 2006 年。4 戴明强 李卫军 杨鹏飞，数学模型及其应用，北京：科学出版社， 2007 年.5 谢金星、薛毅编著， 优化建模与 lingo/lindo 软件，北京：清华大学出版社， 2005.7.8 / 11附件：附录 1model:sets :cun/1.10/;link(cun,cun) |&2#gt#&1:q,x,d;endsetsdata:d =8 13 8 11 9 17 15 8 167 14 9 11 14 23 16 1817 4 12 8 22 11 1612 6 15 12 15 137 7 17 9 1210 11 10 818 4 715 116;enddatamin=sum(link(k,i)|k#lt#i:x(k,i)*(2/3)* sqrt(q(k,i)