快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

留学>>

tag 标签: 解决方案

相关帖子

 版块 作者 回复/查看 最后发表
【数模招聘猎头】数据分析师招聘信息 人才招聘 厚积薄发 2011-7-15 3 3328 alair005 2012-2-7 16:55
DIY系列之——某校图书证条形码制作 attach_img 全国大学生数学建模竞赛(CUMCM) buguniaoo1991 2011-9-28 5 5619 书呆子 2011-10-2 22:53
2011年第一期数学中国SAS职业统计分析实训开班通知 attach_img MADIO数模培训 厚积薄发 2011-11-22 36 15032 流雨星月 2014-8-31 09:43
常见问题 美国大学生数学建模竞赛(MCM/ICM) 萧邦很忙 2012-2-1 0 1975 萧邦很忙 2012-2-1 15:59
为什么16号才收到电子解决方案！！！ 美国大学生数学建模竞赛(MCM/ICM) 550406088 2012-2-17 22 4866 garfieldme 2012-2-20 00:06
运行spss，弹出：试图连接至远程服务器失败。求解决方案？具体点 attach_img SPSS论坛 txwgzyzz 2012-4-2 1 4 txwgzyzz 2012-4-4 12:18
运行spss，弹出：试图连接至远程服务器失败。求解决方案？ attach_img SPSS论坛 txwgzyzz 2012-4-3 1 12143 txwgzyzz 2012-4-4 12:15
MATLAB安装过程中常见问题及解决方案 MATLAB论坛 hopeoflight 2012-5-29 5 24297 369254huang 2013-1-24 19:32
从数学建模到经理人的跨越归功于能力的提升 数模经验分享 飞翔的娇龙 2012-10-16 12 4022 becareful 2013-12-17 22:45
大规模工程计算解决方案-并行程序设计 培训班 算法类综合讨论版 塔木德2026 2012-11-13 0 8761 塔木德2026 2012-11-13 19:18
大规模工程计算解决方案-并行程序设计 培训班 attachment 并行算法 塔木德2026 2012-11-23 3 4511 空木葬花 2014-3-9 22:30
2013年“深圳杯”数学建模夏令营 B题：深圳关内外交通拥堵探究与治理 attachment 深圳杯数学建模夏令营 厚积薄发 2013-4-15 44 26634 wekun 2018-5-15 18:17
中考数学备考：数学问题的学习策略及解决方案 初中数学学习方法汇总 梦想在飞 2013-5-6 0 1984 梦想在飞 2013-5-6 17:48
matlab simulink 2009a和2011a的兼容问题 MATLAB论坛 王小丫丫 2013-7-31 0 748 王小丫丫 2013-7-31 16:49
关于选购问题的解决方案 层次分析法 不要回头不要悔 2013-8-2 0 1299 不要回头不要悔 2013-8-2 15:45
[书籍] 数据挖掘解决方案集锦 attachment 数据挖掘 挚爱/Sunsgne° 2013-8-13 3 1542 李钊2013 2013-9-28 07:12
投票 快来谈谈自己的国赛Idea吧，回帖就可能会有免费参加国际赛与美赛数学中国赛前培训~！ 全国大学生数学建模竞赛(CUMCM) madio 2013-9-16 0 11 madio 2013-9-16 11:36
关于美赛的常见问题解答 美国大学生数学建模竞赛(MCM/ICM) 百年孤独 2013-9-27 21 14296 zt1222 2014-2-10 10:06
LaTeX技巧81：WinEdt中文乱码问题解决方案 attachment LaTeX/ScienceWord论坛 挚爱/Sunsgne° 2013-11-13 0 827 挚爱/Sunsgne° 2013-11-13 14:57
LaTeX技巧344：一个warning里的简单解决方案 attachment LaTeX/ScienceWord论坛 梦里花111 2013-11-21 0 1129 梦里花111 2013-11-21 11:00
更多...

相关日志

分享 你好数模
热度 1 我是学渣 2016-9-6 21:18
数学建模竞赛对专业的限制较小，任何专业的学生都可以参与这项活动，而不需要有什么专业背景。如果非要说背景，当然理工科的会有一定优势，但是文科学生也可以在这项比赛中发挥他们的长处，比如在论文写作方面。 在数学建模竞赛中，参赛队伍被要求在三到四天时间内解决一个实际问题，在这个过程中，可以使用任何工具、软件，可以随意跑动而没有固定的比赛场所，可以查阅任何资料、文献，但是不能与本队伍之外的其他任何 有生命 的物体就赛题进行交流（可以上网搜索，但是不能通过网络向别人求助）。最终参赛队必须提交一篇科技论文，以展示他们对问题的理解和做出的解决方案。而队伍成绩评定的依据，也全部来自他们提交的这篇论文。 。
115 次阅读|0 个评论
分享 泄洪河道的解决方案
huitaolang 2015-9-2 11:33
泄洪设施修建计划 摘 要 位于我国南方的某个偏远贫困乡，地处山区，一旦遇到暴雨，经常发生洪涝灾害。以往下雨时，完全是依靠天然河流进行泄洪。 2010 年入夏以来，由于史无前例的连日大雨侵袭，加上这些天然河流泄洪不畅，造成大面积水灾，不仅夏粮颗粒无收，而且严重危害到当地群众的生命财产安全。本文利用数学建模知识并结合生活实际设计一系列最优的方案。 对于问题一，首先根据原天然河流近几年排洪能力的数据，利用灰度预测的 MATLAB 程序求出未来五年的排洪能力。将排洪沟的决策变量假设成 0-1 变量，建立 0-1 规划模型。根据投入费用及排洪量的约束条件，建立线性规划模型。利用 MATLAB 求出从 2010 年起，第一年开挖 2 、 3 、 4 、 5 四条排洪沟，其费用为 22 万元，修建泄洪河道的费用为 38 万元，第二年开挖第 6 条排洪沟，其费用为 5 万元，修建泄洪河道的费用为 55 万元，第三年开挖第 1 条排洪沟，其费用为 5 万元，修建泄洪河道的费用为 46 万元。在整个修建计划中所花费的金额为 171 万元。 对于问题二，首先根据地势的高低自西向东排列村庄，利用 Dijkstra 算法求出最小路径，将泄洪河道的决策变量假设成 0-1 变量，根据排洪能力以及费用最小的约束条件， 列出目标函数，并对最短路径进行人工干预求出的路径为 ⑥ - ⑤，② - ① - ③ - ④ - ⑤ - ⑧，⑨ - ⑦ - ⑧，⑩ - ⑧。计算所得的最少费用为 576.99 万元。 对于问题三，利用问题二中建立的泄洪河道网络考虑维修工路线进而考虑留在各村的概率，建立卡尔科夫链模型。利用 LINGO 求解出留在各村的概率，再利用 MATLAB 求出转移概率矩阵的特征值和特征向量。利用 MATLAB 分析稳定性。验证结果说明在各村留宿概率为上述 LINGO 求解结果。判定留宿结果稳定，当进行足够多次维修后在各村留宿概率稳定。 对于问题四，因为影响泄洪设施的因素有多种，因此找出影响泄洪设施的主要因素就行，主要因素有 1. 天然河道的曲折； 2. 天然河道需要定期进行清淤； 3. 建立一个主泄洪干道。这些因素都会影响修建泄洪河道的成本，对这些影响因素建立约束函数，找出约束函数，求出相应因素所需的最小费用 、 、 ，然后相加即为最少费用。 关键词 ： 0-1 规划 线性规划 最小生成树 马尔科夫链模型 一、 问题背景和重述 1.1 问题的背景 位于我国南方的某个偏远贫困乡，地处山区，一旦遇到暴雨，经常发生洪涝灾害。以往下雨时，完全是依靠天然河流进行泄洪。 2010 年入夏以来，由于史无前例的连日大雨侵袭，加上这些天然河流泄洪不畅，造成大面积水灾，不仅夏粮颗粒无收，而且严重危害到当地群众的生命财产安全。 1.2 问题的重述 为此，乡政府打算立即着手解决防汛水利设施建设问题。从两方面考虑，一是在各村开挖一些排洪沟 , 以满足近两三年的短期防汛需要；二是从长远考虑，可以通过修建新泄洪河道的办法把洪水引出到主干河流。经测算，修建新泄洪河道 的费用为 （万元），其中 表示 泄洪河道 的可 泄洪 量（万立方米/小时）， 表示 泄洪河道 的长度（公里）。 请你们通过数学建模方法，解决以下问题： 问题 1 ： 该乡的某个村区域内原有四条天然河流，由于泥沙沉积，其泄洪能力逐年减弱。在表 1 中给出它们在近年来的 可 泄洪 量（万立方米/小时） 粗略统计数字。水利专家经过勘察，在该村区域内规划了 8 条可供开挖排洪沟的路线。由于它们的地质构造、长度不同，因而开挖的费用和预计的 可 泄洪 量 也不同，详见表 2 ，而且预计每条排洪沟的 可 泄洪 量 还会以平均每年 10% 左右的速率减少。同时开始修建一段 20 公里长的新泄洪河道。修建工程 从开工到完成需要三年时间，且每年投资修建的费用为万元的整数倍。要求完成之后，通过 新泄洪河道 能够达到可 泄洪 量100万立方米/小时的泄洪能力。 乡政府从2010年开始，连续三年，每年最多可提供60万元用于该村 开挖排洪沟和修建新泄洪河道 ，为了保证该村从2010至2014年这五年间每年分别能至少达到可 泄洪 量150、160、170、180、190万立方米/小时的泄洪能力，请作出一个从2010年起三年的 开挖排洪沟和修建新泄洪河道 计划，以使整个方案的总开支尽量节省（不考虑利息的因素在内）。 表 1 现有四条天然河道在近几年的可泄洪量（万立方米 / 小时） 年份 编号 2001 2002 2003 2004 2005 2006 2007 2008 2009 1 号 32.2 31.3 29.7 28.6 27.5 26.1 25.3 23.7 22.7 2 号 21.5 15.9 11.8 8.7 6.5 4.8 3.5 2.6 2.0 3 号 27.9 25.8 23.8 21.6 19.5 17.4 15.5 13.3 11.2 4 号 46.2 32.6 26.7 23.0 20.0 18.9 17.5 16.3 表 2 开挖各条排洪沟费用（万元）和预计当年可泄洪量（万立方米 / 小时） 编号 1 2 3 4 5 6 7 8 开挖费用 5 7 5 4 6 5 5 3 当年泄洪量 25 36 32 15 31 28 22 12 问题 2 ： 该乡 共有10个村,分别标记为①— ⑩ ，下图给出了它们大致的相对地理位置，海拔高度总体上呈自西向东逐渐降低的态势。 ① ② ③ ④ ⑤ ⑥ ⑦ ⑧ ⑨ ⑩ 其中村 ⑧ 距离 主干河流 最近 ， 且海拔高度最低 。乡政府打算拟定一个 修建 在各村之间互通的 新泄洪河道网络 计划， 将洪水先通过 新泄洪河道 引入村 ⑧ 后，再经村 ⑧ 引出到主干河流。要求 完成之后，每个村通过 新泄洪河道 能够达到可 泄洪 量100万立方米/小时的泄洪能力。 表 3 各村之间修建新泄洪河道的距离（单位：公里） 2 3 4 5 6 7 8 9 10 1 2 3 4 5 6 7 8 9 8 5 9 12 14 12 16 17 22 9 15 17 8 11 18 14 22 7 9 11 7 12 12 17 3 17 10 7 15 18 8 10 6 15 15 9 14 8 16 8 6 11 11 11 10 请你们根据表 3 中的数据，为该乡提供一个各村之间修建新泄洪河道网络的合理方案，使得总费用尽量节省。（说明：从村 A→ 村 B 的新泄洪河道，一般要求能够承载村 A 及上游新泄洪河道的泄洪 量 。） 问题 3 ： 新泄洪河道网络铺设完成后，打算安排一位维护人员，每天可以从一个村到与之直接有新泄洪河道连接的相邻村进行设施维护工作，并在到达的村留宿，次日再随机地选择一个与该村直接有新泄洪河道连接的相邻村进行维护工作。试分析长此以往，他在各村留宿的概率分布是否稳定？ 问题 4 ： 你们是否能够为该乡提出一个更加合理的解决泄洪的办法？ 说明： 1 、以上问题必须建立一般的数学模型，不能仅按照题目中提供的数据计算一个结果。 2 、建模过程中，可自行提出合理的模型假设。 二、问题分析 针对问题一：该村区域内原有四条天然河流，但是由于泥沙沉积，其泄洪能力逐年减弱。根据表 1 中四条天然河流近几年来的 可 泄洪 量并利用灰度预测算法求出未来 5 年中每条河道的可泄洪量。假设修建的新泄洪河道在修建完成之前不能投入使用，在此期间修建 条排洪沟道，每条排洪沟道只能开挖一次。用 MALTAB 软件及 0-1 模型求出每年需要开挖的排洪沟道，按照每年的排洪量以及费用的约束条件，列出总费用与修建项目的目标函数，求解总费用最少的优化问题。 针对问题二：该村地势自西向东，泄洪河道只能自西向东的修建，建立 0-1 模型，利用 Dijkstra 算法的 MALTAB 程序计算最小生成树，再根据排洪能力和排洪量的约束条件进行人工干预，求解出泄洪河道的修建结果，已知河道的修建网络图，计算出各村排洪量的能力以及对应的距离，在根据 求解各村间修建泄洪河道的费用，最后将各村的费用相加即为最终结果。 针对问题三： 假设维护人员从某村出发，随机选择一个与该村直接有新泄洪河道连接的相邻进行维修工作。根据问题二中的泄洪干道分布图可以的得到村落之间的邻接矩阵，村落之间的邻接矩阵就是维修工留宿可转移矩阵。根据矩阵建立马尔科夫链模型，求出转接概率矩阵，进而求出在各村留宿概率。在检验概率稳定性。 针对问题四： 1. 考虑到天然河道有些地段可能过于曲折， 2. 由于泥沙的逐年沉积，天然河道还要定期进行清淤，以维持其泄洪量， 3. 建立一个主泄洪道来缓解其泄洪压力，三种主要因素的费用只要都达到最小即可。 三、模型假设 1. 新建泄洪河道在竣工前不投入使用，并且不停止投入工期费用。 2 ．不受人力、物力等因素的影响，修建计划按时完成。 3. 修完的排洪道排洪能力不随着年份衰减，即第 4 、 5 年排洪量为 100 万立方米/小时。 4. 修建的泄洪河道水流流向自西向东。 5 ．各村之间只能修建一条泄洪河道。 四、符号说明 符号 符号说明 第 年开挖第 条排洪沟道 第 年原四条天然河流排洪总量 第 年修建新泄洪河道的费用 第 年新建泄洪河道的长度 第 条排洪沟道的开挖费用 第 条排洪沟道的排洪量 开挖排洪沟和修建新泄洪道的总费用 第 年开挖排洪沟排洪能力 第 个村庄到第 个村庄的排洪量 第 个村庄到第 个村庄的距离 修建泄洪河道的总费用 \* MERGEFORMAT 维修人员转移状态概率矩阵 \* MERGEFORMAT 维修人员可转移矩阵 \* MERGEFORMAT 第 \* MERGEFORMAT 留宿稳定概率 \* MERGEFORMAT 第 \* MERGEFORMAT 天的状态量 五、模型建立与求解 5.1 问题一的模型与求解 整个排洪沟的开挖和新泄洪河道的修建计划可分为以下步骤：1.预测原四条天然河流在未来5年的排洪能力；2.求解排洪沟的开挖费用、修建新泄洪河道费用；3.预计开挖排洪沟及新泄洪河道的排洪能力；4.根据每年提供的费用及每年排洪能力的要求进行约束；5.建立线性规划模型求解目标函数；6.检验第4、5年的排洪能力是否达到要求。以上分析步骤可用流程图表示： 开始 预测原四条天然河流在未来5年的排洪能力 求解排洪沟的开挖费用、修建新泄洪河道费用 预计开挖排洪沟及新泄洪河道的排洪能力 根据每年提供的费用及每年排洪能力的要求进行约束 建立线性规划模型求解目标函数 检验第4、5年的排洪能力是否达到要求 结束 图5.1 求解问题一步骤的程序框图 5.1.1原天然河流泄洪能力的预测 该乡的某个村区域内原有四条天然河流，由于泥沙沉积，其泄洪能力逐年减弱。根据表 1 中每个天然河流前 9 年的排洪能力数据，写出灰度预测算法的程序（附录 1 ）并利用 MALTAB 软件求解。求解的结果如下 : 表 5.1.1 天然河流泄洪未来 5 年间泄洪能力的预测结果 年份 编号 2010 2011 2012 2013 2014 1 号 21.5 20.9 19.9 19.1 18.2 2 号 1.5 1.1 0.8 0.6 0.4 3 号 10.9 9.7 8.7 7.8 7.0 4 号 13.2 11.7 10.3 9.2 8.1 总量 47.0 43.4 39.8 36.6 33.8 利用 MALTAB 软件求解四条天然河流排洪能力变化的图形如下： 图 5.1.2 四条天然河流排洪能力 5.1.2 排洪沟的开挖费用、修建新泄洪河道费用 乡政府从 2010 年开始，连续三年，每年最多可提供 60 万元用于该村 开挖排洪沟和修建新泄洪河道。根据这一条件可列公式： 该乡从长远考虑，通过修建新泄洪河道的办法把洪水引出到主干河流。经测算，修建新泄洪河道 的费用为 （万元），其中 表示 泄洪河道 的可 泄洪 量（万立方米 / 小时）， 表示 泄洪河道 的长度（公里）。 修建一段 20 公里长的新泄洪河道。修建工程 从开工到完成需要三年时间，且每年投资修建的费用为万元的整数倍。根据该条件可列公式： 5.1.3 预计开挖排洪沟及新泄洪河道的排洪能力 水利专家经过勘察，在该村区域内规划了 8 条可供开挖排洪沟的路线。由于它们的地质构造、长度不同，因而开挖的费用和预计的 可 泄洪 量 也不同，而且预计每条排洪沟的 可 泄洪 量 还会以平均每年 10% 左右的速率减少。同时开始修建一段 20 公里长的新泄洪河道。新泄洪河道 能够达到可 泄洪 量 100 万立方米 / 小时的泄洪能力。根据该条件可列公式： 5.1.4 建立线性规划模型求解目标函数 乡政府从 2010 年开始，连续三年，每年最多可提供 60 万元用于该村 开挖排洪沟和修建新泄洪河道 ，为了保证该村从 2010 至 2014 年这五年间每年分别能至少达到可 泄洪 量 150 、 160 、 170 、 180 、 190 万立方米 / 小时的泄洪能力。新建泄洪河道完工后，通过 新泄洪河道 能够达到可 泄洪 量 100 万立方米 / 小时的泄洪能力。根据该条件可列公式： 约束条件 表示每条排洪沟只能开挖一次； 表示三年刚好修建 20 公里的泄洪河道。 利用 MATLAB 软件（附录 2 ）可以求出三年间需要开挖的排洪沟，其结果如下： 表 5.1.3 三年间需要开挖的排洪沟 年份 编号 1 2 3 4 5 6 7 8 2010 0 1 1 1 1 0 0 0 2011 0 0 0 0 0 1 0 0 2012 1 0 0 0 0 0 0 0 （注）： 0 表示不需要开挖的排洪沟， 1 表示需要开挖的排洪沟。 由表格的数据可知：第一年需要开挖第 2 、 3 、 4 、 5 四条排洪沟，其费用为 22 万元；第二年需要开挖第 6 条排洪沟，其费用为 5 万元；第三年需要开挖第 1 条排洪沟，其费用为 5 万元。 已知排洪沟的开挖时间和数目可以求出每一年中的排洪量，其结果如下： 表 5.1.4 每年的总排洪量 年份 2010 2011 2012 2013 2014 旧河道总排洪量 47.0 43.6 39.8 36.6 33.8 2 、 3 、 4 、 5 号排洪量 114 102.6 92.34 83.106 74.8 6 号排洪量 0 28 25.2 22.68 20.4 1 号排洪量 0 0 25 22.5 20.25 总排洪量 161.0 174.0 182.3 265.0 249.3 三年间开挖的排洪沟所用的总费用为 32 万元， 等于 100 ， 等于 20, 由 可求出 等于 138.221 万元。修建泄洪河道花费为万元整数倍， 等于 139 万元。因此修建泄洪河道及开挖排洪沟所需的总费用为 171 万元。新建的泄洪河道在未竣工之前不停工期，所以每年用于修建新泄洪河道的费用如表所示： 表 5.1.5 两种排洪方案各自所需费用 年份 2010 2011 2012 排洪沟费用 22 5 5 泄洪河道费用 38 55 46 5.2 问题二的模型与求解 对于问题二的求解可分为以下步骤： 1. 按自西向东的顺序排列各村间的距离； 2. 建立 0-1 模型并利用 Dijkstra 算法求出各村间是否修建泄洪河道； 3. 根据第二步的结果求解各村的排洪量； 4. 根据 求解各村间修建费用，其结果相加得到最终费用。 5.2.1自西向东排列各村间距离 由图可知，⑥的地势最高，⑧的地势最低，各个村之间地势顺序为：⑥②①=⑨③⑦④=⑩⑤⑧,各村间的距离如下表： 表5.2.1 自西向东各村间的距离 6 2 9 1 3 7 10 4 5 8 6 2 9 1 3 7 10 4 5 8 0 8 14 15 11 9 17 16 8 14 0 8 14 9 11 15 22 17 18 0 17 5 12 9 22 12 16 0 12 6 15 10 15 11 0 7 7 17 9 12 0 10 11 10 8 0 18 3 7 0 15 11 0 6 0 5.2.2 各个村庄是否修建泄洪河道的求解 根据第一步的距离数据可以利用 Dijkstra 算法（附录 3 ）的 MALTAB 软件求出最短路径，其结果如下图所示： 图 5.2.2 最小生成树 该村 ⑧ 距离 主干河流 最近 ， 且海拔高度最低 。乡政府打算拟定一个 修建 在各村之间互通的 新泄洪河道网络 计划， 将洪水先通过 新泄洪河道 引入村 ⑧ 后，再经村 ⑧ 引出到主干河流。根据该条件并建立 0-1 模型可得以下的数学模型。 目标函数： 约束条件 表示 6 号村庄排洪能力为 100 万立方米 /小时。 要求 完成之后，每个村通过 新泄洪河道 能够达到可 泄洪 量 100 万立方米 / 小时的泄洪能力 （说明：从村 A → 村 B 的 新泄洪河道 ，一般要求能够承载村 A 及上游 新泄洪河道 的 泄洪 量 。）根据该条件可列公式： 假设各个村庄间只能修建一条泄洪河道，对于图 5.2.2 的结果，根据约束条件进行人工干预，得到最终的结果如下图： ① ② ③ ④ ⑤ ⑥ ⑦ ⑧ ⑨ ⑩ 图 5.2.3 修建泄洪河道网络图 由图可知其结果为：⑥ - ⑤，② - ① - ③ - ④ - ⑤ - ⑧，⑨ - ⑦ - ⑧，⑩ - ⑧。又可以表示出相连村庄间的距离。 表 5.2.4 相连村庄间的距离 村庄 距离 村庄 距离 （ ⑥ - ⑤ ） 8 （ ⑤ - ⑧ ） 6 （ ② - ① ） 8 （ ⑨ - ⑦ ） 6 （ 道 4 angANG ① - ③ ） 5 （ ⑦ - ⑧ ） 8 （③ - ④ ） 7 （ ⑩ - ⑧ ） 11 （ ④ - ⑤ ） 3 5.2.3 各个村庄排洪量的求解 泄洪河道 完成之后，每个村通过 新泄洪河道 能够达到可 泄洪 量100万立方米/小时的泄洪能力 （说明：从村 A→ 村 B 的新泄洪河道，一般要求能够承载村 A 及上游新泄洪河道的泄洪 量 。）再结合步骤 1 和步骤 2 可以得到各村的排洪能力如下表： 表 5.2.5 各村排洪能力 村庄 排洪能力 村庄 排洪能力 （ ⑥ - ⑤ ） 100 （ ⑤ - ⑧ ） 500 （ ② - ① ） 100 （ ⑨ - ⑦ ） 100 （① - ③ ） 200 （ ⑦ - ⑧ ） 200 （ ③ - ④ ） 300 （ ⑩ - ⑧ ） 100 （ ④ - ⑤ ） 400 5.2.4 各村修建泄洪河道的费用及总费用 根据以上步骤先进行各村费用的求解，最后各村间费用的和即为最少的费用，其结果如下： 表 5.2.6 修建泄洪河道各村的费用及总的费用 村庄 费用 村庄 费用 （ ⑥ - ⑤ ） 55.29 （ ⑤ - ⑧ ） 94.23 （ ② - ① ） 55.29 （ ⑨ - ⑦ ） 41.47 （ 道 4 angANG ① - ③ ） 49.21 （ ⑦ - ⑧ ） 78.73 （② - ④ ） 84.72 （ ⑩ - ⑧ ） 76.02 （ ④ - ⑤ ） 42.05 总费用 576.99 5.3 模型建立求解 针对问题三 , 可分为三步， 1. 分析状态量变化过程建立马尔科夫链。 2.LINGO 求解维修人员留宿各村概率。 3. 分析 LINGO 求得的解的稳定性。 5.3.1 建立建立马尔科夫链模型 在题中维修工是一天为单位进行维修和留宿。可以将时间以天为单位离散化 \* MERGEFORMAT 单位是天，对 \* MERGEFORMAT 系统即维护人员第 \* MERGEFORMAT 天所在村庄，用随机变量 \* MERGEFORMAT 表示 \* MERGEFORMAT 表示第 \* MERGEFORMAT 天维护人员赶到一村， \* MERGEFORMAT 表示维修人员赶到二号村，于是 \* MERGEFORMAT 可取离散值为 \* MERGEFORMAT 。用 \* MERGEFORMAT 表示第 \* MERGEFORMAT 个村庄留宿的概率。记为 \* MERGEFORMAT ，表示第 \* MERGEFORMAT 天在第 \* MERGEFORMAT 个村庄留宿的概率。用 \* MERGEFORMAT 表示已知今天在第 \* MERGEFORMAT 个村庄留宿明天到第 \* MERGEFORMAT 个村庄留宿的概率。且 \* MERGEFORMAT , 记 \* MERGEFORMAT 。有题意可知 , 第 \* MERGEFORMAT 天的留宿位置 , 即状态量 \* MERGEFORMAT 仅仅取决于第 \* MERGEFORMAT 天的状态量 \* MERGEFORMAT 和状态转移矩阵。即维修工留宿转移概率矩阵。而与以前的状态量 \* MERGEFORMAT 无关，即状态转移具有无后效性。这种离散状态按照离散时间的随机转移过程满足马尔科夫链条件，有状态转移无后效性和全概率性。 马尔科夫链模型： \* MERGEFORMAT \* MERGEFORMAT 并且 \* MERGEFORMAT 和 \* MERGEFORMAT 满足： \* MERGEFORMAT \* MERGEFORMAT \* MERGEFORMAT \* MERGEFORMAT \* MERGEFORMAT \* MERGEFORMAT 5.3.2 求解维修人员留宿各村概率 引入状态概率向量和转移概率矩阵： \* MERGEFORMAT \* MERGEFORMAT 则基本方程表示为 \* MERGEFORMAT 有此可以得到 \* MERGEFORMAT 。 显然第 \* MERGEFORMAT 状态量 \* MERGEFORMAT 取决于 \* MERGEFORMAT 和 \* MERGEFORMAT 而且与前状态量无关具有无后效性，用 \* MERGEFORMAT 表示第 \* MERGEFORMAT 天维护人员留宿状态概率矩阵，即 \* MERGEFORMAT ，用 \* MERGEFORMAT 表示矩阵 \* MERGEFORMAT 表示转移概率矩阵。根据泄洪道网络分布图用 \* MERGEFORMAT 表示表示村落之间是否修建泄洪道，可以得到留宿人员转移矩阵 \* MERGEFORMAT 。其中： \* MERGEFORMAT 根据马尔科夫链模型存在状态概率 \* MERGEFORMAT 当 \* MERGEFORMAT 时状态量 \* MERGEFORMAT ， \* MERGEFORMAT 与初始状态无关。满足： \* MERGEFORMAT \* MERGEFORMAT 计算出 \* MERGEFORMAT 再对转移概率矩阵进行稳定性分析。 模型求解： 由问题二得到的可转移矩阵： 表 5.3.1 可转移矩阵 0 1 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 1 0 0 0 0 0 0 0 0 1 0 1 0 0 0 0 0 0 0 0 1 0 1 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 1 0 1 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 维修工转移概率矩阵： 表 5.3.2 维修工转移概率矩阵 0 0.5 0.5 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0.5 0 0 0.5 0 0 0 0 0 0 0 0 0.5 0 0.5 0 0 0 0 0 0 0 0 0.33333 0 0.33333 0 0.33334 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0.5 0.5 0 0 0 0 0 0.33333 0 0.33333 0 0 0.33334 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 用 LINGO 计算的各村留宿概率 \* MERGEFORMAT 的结果为： 5.3.3 留宿各村概率 1 2 3 4 5 6 7 8 9 10 0.1052 0.0526 0.1052 0.1052 0.1578 0.0526 0.1052 0.1578 0.0526 0.1052 （源代码见附录 4 ） 5.3.3 稳定性分析 转移矩阵为 \* MERGEFORMAT ，且 \* MERGEFORMAT 。利用 MATLAB 求出特征值 \* MERGEFORMAT 特征向量 \* MERGEFORMAT 。如果特征向量可逆，则可以求出 \* MERGEFORMAT 的无穷次幂。求得 \* MERGEFORMAT 可逆，有 \* MERGEFORMAT 可以看出 \* MERGEFORMAT 的无穷次幂存在。从而得到的 \* MERGEFORMAT 稳定存在（程序源代码和 \* MERGEFORMAT 值矩阵见附录 5 ）。 5.4 问题四的模型与求解 5.4.1 模型影响因素 根据问题分析可知，影响天然河道的最优泄模型的主要因素为： 1. 天然河道的曲折； 2. 天然河道需定期清淤； 3. 修建主泄洪干道； 5.4.2 天然河道最优泄洪模型及求解 综合以上主要影响因素，建立天然河道最优泄洪模型， 表示拓宽天然河道费用； 表示天然河道清淤的费用； 表示修建主干泄洪河道的费用。建立修建泄洪河道网络的资金花费的目标函数，即： 对于因素 1 ，当天然河道过于曲折时，会造成拓深拓宽天然河道的成本增加，若在此段距离修建新的河道费用小于拓深拓宽天然河道的花费，则可以修建新的河道绕过这段天然河道， 此时的费用 最小。 对于因素 2 ，由于泥沙的沉积，天然河道还要定期进行清淤，以维持其泄洪量。则可以根据河流的泥沙年沉积速率和泄洪量的要求，通过计算找出最优的泥沙清理时间间隔，这样能很好的减少平均的年清淤花费，此时的费用 最小。 对于因素 3 ，考虑到该乡处在山区，且村落分步基本成南北两向分散状，可以在中间人口分布密度较少的地带修建一条主泄洪干道，各村分别修建泄洪河道将引入主泄洪干道，主泄洪道自西向东汇入村落排洪量逐渐增大，求其最少的费用 。 六、模型的评价和推广 6.1 模型的优点 建立的 0-1 规划一般模型大大减低了建模的难度，简化了计算编程过程，并使论文叙述更加清晰易懂。 马氏链模型很接近于实际情况，模型与实际情况误差很小。而且建立的模型简单易行，不仅可以解决十个村庄的随机留宿的概率分布问题，也可以解决更多的村庄留宿概率分布问题，体现了此模型具有通用性。 问题 4 的模型的建立使泄洪方案更具可行性，设计更具可行性。 6.2 模型的缺点 问题一：模型是在假设洪水来之前当年的排洪沟的开挖工程已经结束，但是在现实中并不排除洪水来之前，排洪沟还没有开挖或尚未完成的情况。 问题二：村庄之间并不是一定能在直线上对新泄洪河道进行修建，可能会有一些无法解决的问题使新泄洪河道沿曲线修建。 6.3 模型的推广 1. 可以推广到 n 个自高到低的村庄或山坡等 。 2. 此模型也能应用于救灾紧急物资购置，如灾区暂住房与长期住房建设， 救灾应急线路修建等。 3. 此模型还能应用于考虑到线路承载能力的网状 流运输模型，如互联网流量，物流运输等。 4. 问题三为马氏链模型，具有无后效性，能广泛的应用于经济、社会、生态、 遗传等众多领域。 5. 此类模型能应用于大范围的水利工程修 建，从整体上起到调控作用。还能应用于长期投资与短期投资。 七、参考文献 . 姜启源 . 数学模型 . 北京，高等教育出版社， 2003 。 . 胡量剑，丁晓东，孙晓君 . 数学实验，上海科学技术出版社， 2003 。 . 杨超，运筹学，北京：科学出版社， 2004 。 . 谢金星、薛毅，优化建模与 LINDO 软件，北京：清华大学出版社。 . 刘克，《实用马尔科夫接侧过程》，清华大学出版社， 2004 。 . 韩中庚 . 数学建模方法及其应用， 高等教育出版社 2008 。 附 录 附录 1 syms a b; c= '; A= ; B=cumsum(A); n=length(A); for i=1:n-1; C(i)=(B(i)+B(i+1))/2; end D=A;D(1)= ; c=inv(E*E')*E*D; c=c'; a=c(1);b=c(2); F= ;G(1)=A(1); for i=2:(n+5) G(i)=F(i)-F(i-1); end t1=2001:2009; t2=2010:2014; G plot(t1,A,'o'); xlabel(' 年份 '),ylabel(' 排洪能力 '); title(' 第 4 条天然河流 '); hold on plot(t2,G); 表 1 现有四条天然河道在近几年的可泄洪量（万立方米 / 小时） 1 号 32.2 31.3 29.7 28.6 27.5 26.1 25.3 23.7 22.7 2 号 21.5 15.9 11.8 8.7 6.5 4.8 3.5 2.6 2.0 3 号 27.9 25.8 23.8 21.6 19.5 17.4 15.5 13.3 11.2 4 号 46.2 32.6 26.7 23.0 20.0 18.9 17.5 16.3 附录 2 n=10; w=inf*ones(10); w(6, )= ; w(2, )= ; w(1, )= ; w(9, )= ; w(3, )= ; w(7, )= ; w(4, )= ; w(10, )= ; w(5,8)= ; =mintreek(n,w) p=@(q,l) 0.66*q^0.51*l; price= ; A= ; b= ; A= ; b= ; LL= ; A1=- ; b1=- '+ '; A= ; b= ; temp=eye(8); temp= ; A= ; b= ; =bintprog( ,A,b); answer=reshape(x,8,3)'; Optimization terminated. 附录 3 function =mintreek(n,W) tmpa=find(W~=inf); =find(W~=inf); w=W(tmpa); e= ; =sort(w); E= ; =size(E); temp=find(E(:,1)-E(:,2)); E=E(temp,:); P=E(1,:); k=length(E(:,1)); while(rank(E)0) temp1=max(E(1,2),E(1,1)); temp2=min(E(1,2),E(1,1)); for i=1:k; if(E(i,1)==temp1),E(i,1)=temp2;end; if(E(i,2)==temp1),E(i,2)=temp2;end end; a=find(E(:,1)-E(:,2)); E=E(a,:); if(rank(E)0),P= ;k=length(E(:,1));end; end; Wt=sum(P(:,3)); Pp= ; for i=1:length(P(:,3)); disp( ); end; % 画图程序 axis equal;hold on =cylinder(1,n);xm=min(x(1,:));ym=min(y(1,:)); xx=max(x(1,:));yy=max(y(1,:)); axis( );plot(x(1,:),y(1,:),'ko') for i=1:n;temp= ; text(x(1,i),y(1,i),temp);end; for i=1:length(P(:,4)), plot(x(1,Pp(i,1:2)),y(1,Pp(i,1:2)),'k');end; text(-0.35,-1.2, ); title(' 最终的 Steiner 树 ');axis('off');hold off 附录 4 model : sets : cun/A..J/:p; ! 用大写字母表示村庄 ， 定义名称为拼音村 “cun”, 用 p 表示留宿在各村的概率 ; link(cun,cun):t; ! 用 t 表示转移的概率 ; endsets data : t=0 0.5 0.5 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0.5 0 0 0.5 0 0 0 0 0 0 0 0 0.5 0 0.5 0 0 0 0 0 0 0 0 0.333333 0 0.333333 0 0.333334 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0.5 0.5 0 0 0 0 0.333333 0 0.333333 0 0 0.333334 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0; enddata @for (cun(i):p(i)= @sum (link(j,i):p(j)*t(j,i))); @sum (cun:p)=1; @for (cun(i): @warn (' 矩阵每行之和是 1', @abs (1- @sum (link(i,j):t(i,j)))#GT#0.0000001)); End 附录 5 D= 1 0 0 0 0 0 0 0 0 0 0 0.911 0 0 0 0 0 0 0 0 0 0 0.707 0 0 0 0 0 0 0 0 0 0 0.442 0 0 0 0 0 0 0 0 0 0 0.292 0 0 0 0 0 0 0 0 0 0 -1 0 0 0 0 0 0 0 0 0 0 -0.91 0 0 0 0 0 0 0 0 0 0 -0.71 0 0 0 0 0 0 0 0 0 0 -0.29 0 0 0 0 0 0 0 0 0 0 -0.44 V= -0.32 -0.46 -0.25 0.11 0.14 -0.32 -0.46 0.25 -0.14 -0.11 -0.32 -0.50 -0.35 0.24 0.47 0.32 0.50 -0.35 0.47 0.24 -0.32 -0.33 0.00 -0.15 -0.39 0.32 0.33 0.00 -0.39 -0.15 -0.32 -0.15 0.25 -0.24 -0.36 -0.32 -0.15 -0.25 0.36 0.24 -0.32 0.06 0.35 -0.06 0.17 0.32 -0.06 0.35 0.17 -0.06 -0.32 0.07 0.50 -0.14 0.60 -0.32 0.07 -0.50 -0.60 0.14 -0.32 0.35 -0.35 -0.22 0.03 0.32 -0.35 -0.35 0.03 -0.22 -0.32 0.25 0.00 0.30 -0.08 -0.32 0.25 0.00 0.08 -0.30 -0.32 0.38 -0.50 -0.49 0.10 -0.32 0.38 0.50 -0.10 0.49 -0.32 0.28 0.00 0.67 -0.27 0.32 -0.28 0.00 -0.27 0.67 det （ v ） =-0.165 clear =eig(A) det(v)
258 次阅读|0 个评论
分享 京东方科技集团股份有限公司
李崇森 2013-12-10 22:38
公司网址：http://www.boe.com.cn/ 一、集团简介（上市国企） 京东方科技集团股份有限公司（京东方A：000725；京东方B：200725）创立于1993年4月，是一家光电显示技术、产品与解决方案的提供商。经过多年的技术积累与创新，京东方可使用专利数量达7500余项，2011年度新增专利申请数量突破1000项，居国内同业前列，成为中国大陆显示领域最具综合实力的高科技企业。目前，京东方拥有一条月产能为4.5万片玻璃基板的第4.5代TFT-LCD生产线（成都）、一条月产能为10万片玻璃基板的第5代TFT-LCD生产线（北京）、一条月产能为10万片玻璃基板的第6代TFT-LCD生产线（合肥）、一条月产能为9万片玻璃基板的第8.5代TFT-LCD生产线（合肥）、一条月产能为9万片玻璃基板的第8.5代TFT-LCD生产线（北京）、一条月产能为5.4万片玻璃基板的第5.5代AMOLED生产线（鄂尔多斯）、一条月产能为9万片玻璃基板的第8.5代TFT-LCD生产线（重庆）。 截至2013年9月30日，京东方注册资本135.2亿元，净资产272.5亿元，总资产906.1亿元，员工人数22980人。在北京、四川成都、安徽合肥、内蒙古鄂尔多斯、重庆、河北固安、江苏苏州、福建厦门等地拥有多个制造基地，营销和服务体系覆盖欧、美、亚等全球主要地区。网址：http://www.boe.com.cn/ 二、招聘分公司简介： 重庆京东方光电科技有限公司是中国显示产业龙头的京东方携手重庆合作兴建的第8.5代新型半导体显示器件及系统项目，是京东方完善西南地区产业布局、占据新型显示技术战略制高点的重要一步。项目选址两江新区，总投资328亿元人民币，建设用地1035亩，形成集制造与研发于一体、面板与整机于一体的综合性产业基地。项目设计月产能为9万片玻璃基板，计划于2015年二季度正式投产。（2014年首批骨干员工培训实习） 三、招聘岗位及薪资 岗位类别：生产技术，大专以上学历 岗位名称： 1、设备运行助工 2、管理储备干部 要求专业：电子类、机械类、自动化等 应聘条件： 1.专业基础扎实，获得奖学金者优先。 2.具有良好团队意识，学生干部优先。 3.身体健康（能适应倒班）无色盲、色弱。 薪资及福利： 1.工资：2600-3200元，年底发放年终奖。 2.根据国家规定为员工缴纳五险一金。 3.免费工作餐，公寓式宿舍及生日、节日贺金等 四、联系方式 公司地址：重庆市北碚区水土高新技术产业园云汉大道5号附12号 集团招聘电话：010-57676219 简历投递邮箱：cqzp@boe.com.cn
60 次阅读|0 个评论
更多...
qq
收缩

  • 电话咨询

  • 04714969085

关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

手机版|Archiver| |繁體中文 手机客户端  

蒙公网安备 15010502000194号

Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问：王兆丰

GMT+8, 2024-4-26 06:17 , Processed in 0.273098 second(s), 26 queries .

回顶部