抗旱方案的制定建模.docxVIP

----宋停云与您分享---- ----宋停云与您分享---- 福州大学第六届数学建模竞赛 题目：抗旱方案的制定 摘要 本文就抗旱方案的制定，从整个计划的总开支情况考虑，利用 MATLAB 拟合以及单纯形法定理、C 语言编程手段，设计了一个良好的打井与管道铺设计划，即铺设管道费用最小为 139 万元，供水量 Q 为 100 万吨，打井费用为 29 万元，即最小支出费用为 168 万元。同时结合具体实际，兼顾长远的经济效益与实用性，切实解决居民们的生活与农业用水问题。 一、模型假设： 1、忽略该地的年平均降水量。 2、打好的井在当年可立即投入使用。 3、铺设的管道在竣工两年时间内供水量稳定，不受其它因素影响。 4、不考虑小蓄水池的作用和利息的因素在内。 5、现有的四口井在未来的五年内供水量不会发生很大的突变。 6、每年铺设管道的费用可以不相等，若少，则铺设管道长度就短，反之则长。但总的铺设管道时间仍为三年不变。 二、问题分析： 1.若只采用专门打井而不铺设管道的方案，由于每口井年产水量逐年递减的因素，从长远供水效应出发，是不切实际的。所以模型的建立必须同时采用打井和铺设管道方案。 2.利用 MATLAB 对现有四口水井进行拟合，结果显示现有的 2 号井仅在 2010 年提供 1.1 万吨水，2010 年之后 2 号井的供水量可以忽略不计。同时，通过观察发现，现有 4 号井在 2002 至 2005 年供水量跳变范围太大，而在 2006 至 2009 年的供水量减小的速率较为稳定， 则使用其在 2006 至 2009 年的供水量来进行 MATLAB 拟合。 ----宋停云与您分享---- ----宋停云与您分享---- 年份产水量编号号井号井号井号井 年份 产水量编号 号井 号井 号井 号井 2010 2011 2012 2013 2014 21.2 1.1 8.8 15.1 19.4 0 6.2 13.8 17.2 0 3 12.6 14.3 0 0 11.3 10.5 0 0 10.1 3.现有四口水井在 2012 年的供水为 32.8 万吨,由于铺设管道需要三年，为了保证该村在 2012 年能至少获得 170 万吨水，则 2012 年通过 2010 至 2012 年打井所供应的水至少应为 170-32.8=137.2 万吨。同理，经计算可得 2013 年通过前三年打井所供应的水至少应为 0.9 137.2=123.48 万吨；2014 年通过前三年打井所供应的水应为 111.132 万吨。因为管道每年至少提供 100 万吨水，则2013 年总的供水量至少为 123.48+100=223.48 万吨；2014 年总的供水量至少为 111.132+100=211.132 万吨。因此满足 2013 和 2014 年的供水需求。 综上所述：只需考虑 2010 年至 2012 年的打井费用。为使开支最少，必须使管道费用 最少，即管道每年提供 100 万吨水。 三、符号说明： Y：计划总开支（万元） Ci：前三年(代表 2010 年、2011 年和 2012 年）每年打井费用（i=1,2,3）（万元） Q：管道每年的可供水量（万吨/年） L:管道长度（公里） P：铺管道费用 Si：所打的井当年的产水量（i=1,2,3）（万元） a ,a ,a ,a ,a ,a ,a ,a : a :1 表示 2010 年；i 表示按表 2 各井的编号顺序， 11 12 13 14 15 16 17 18 1i a 只能取 0 或 1,；0 代表该井该年没打，1 代表该井该年有打 1i a ,a ,a ,a ,a ,a ,a ,a : a :2 表示 2011 年；i 表示按表 2 各井的编号顺序， 21 22 23 24 25 26 27 28 2i a 只能取 0 或 1,；0 代表该井该年没打，1 代表该井该年有打 2i a ,a ,a ,a ,a ,a ,a ,a : a :3 表示 2012 年；i 表示按表 2 各井的编号顺序， 31 32 33 34 35 36 37 38 3i a 只能取 0 或 1,；0 代表该井该年没打，1 代表该井该年有打 3i 四、模型建立与求解 由于 Q 为 100 万吨，L=20 公里，得 P=[138.2]=139 万元，即前两年 60 万元必须全部支出， 用于管道铺设和打井费用。 ----宋停云与您分享---- ----宋停云与您分享---- ----宋停云与您分享---- ----宋停云与您分享---- 总开支 Y ? 0.66Q0.51 ? L ? ?3 i?1 Ci （Q=100） ----宋停云与您分享---- ----宋停云与您分享---- 2010 年: S1=25 ? a +36 ?