抗旱最优解.docVIP

抗旱方案最优解 黄新颖 黄剑锋 杨锦华 （020900305）（020900310）（020900327） 福州大学机械学院09级机制3班 摘要 我们根据题目分析出该优化问题中水井产水量随时间变化，属于多阶段决策，故我们建立动态规划模型进行求解。求解这个动态优化问题分为两个步骤: 第一是动态规划基本方程，逆序地求出条件最优目标函数和条件最优决策。在此过程中，我们应用了线性回归、指数回归、逆回归等函数来求解出相应的数学模型，从而预测原有水井未来5年的产水量并以最小开支为目标按未来五年供水的要求求出目标函数。 第二是顺序地求出最优决策序列。在求解最优解的过程中，由题目要求每一年打井的决策按顺序排列得允许决策集合。由动态规划的基本方程验证得到其中达到最有效果的最优决策即打井和铺设管道计划。 问题重述 政府计划为某村制定抗旱方案分打新井与铺设管道两方面。打的新井的费用和预计的年产水量不同而且预计每口水井的年产水量会以平均每年10%左右的速率减少。铺设管道20公里外把河水引入该村。铺设管道的费用为 （万元）， 表示每年的可供水量（万吨/年）， 表示管道长度（公里）。铺设管道从开工到完成需要三年时间，且每年投资铺设管道的费用为万元的整数倍。要求完成之后，每年能够通过管道至少提供100万吨水。?????????从2010年开始，连续三年，每年最多可提供60万元用于该村打井和铺设管道，为了保证该村从2010至2014年这五年间每年分别能至少获得150、160、170、180、190万吨水，作出一个从2010年起三年的打井和铺设管道计划，以使整个计划的总开支尽量节省（不考虑小蓄水池的作用和利息的因素在内）。 从2010年算第一年 四、模型建立与求解 动态规划模型 由原有井产水量随年份的变化关系（如下图大体分析其函数类型）分析，通过对线性回归方程、指数回归、乘方回归、对数回归、逆回归等回归方程的比较： 原有水井在近几年的产水量折线图（万吨） 得出函数模型： 四组函数模型相关系数的平方均接近1，故所建立的函数模型与真实数据相关。 由函数模型得出未来五年的预测值，如表所示： 表一 2010~2014原有水井的预测产水量（单位：万吨） 2010 2011 2012 2013 2014 1号井 21.57 20.37 19.17 17.96 16.76 2号井 1.43 1.06 0.78 0.58 0.43 3号井 9.13 7.04 4.96 2.87 0.79 4号井 15.34 14.64 14.07 13.58 13.16 总计 47.47 43.11 38.98 34.99 31.14 要求供水量 150 160 170 180 190 尚需水量（） 102.53 116.89 131.02 145.01 158.86 表格中尚需水量一行在前三年内由于管道铺设尚未完成，所以需由新井产水提供，而在第四、第五年由新井与管道输水共同提供。所以总开支应有打井与铺设管道的两部分， 目标函数 min =P+W （Z=180） 其中与n的约束条件为： (=100) 依题可知，我们前三年只能从旧的井和刚打的井里获得水，由于旧井中的水会逐年减少，上面已经得出未来5年相互相应的预测值了，所以现在主要要解决的问题就是要找出一个最节省钱而又能满足要求的打井方案。至于水管的可供水量，由于它只能在第四年供水，所以只要第五年的时候井里的水与水管供水的总和大于190万吨就可以了。 因为在2014年的时候，旧井可提供的最少水量为31.14万吨，又2012新井的供水量至少要为131.02万吨，由此可以推算出2014年的时候新井可供的总水量为106.12万吨。由此可知，在2014年的时候，仅由井提供的水就可以达到137.26万吨，所以水管的供水量只要达到52.74万吨就足够了，但是题目要求管道的供水量至少为100万吨，又因为P要取万元的整数倍，因此我们取=101.1万吨，则此时的管道铺设费用为139万元。 由此，问题便转化为打井的最优决策了。 允许决策集合的确定 表二 8个位置打井费用（万元）和当年产水量（万吨） 编号 1 2 3 4 5 6 7 8 打井费用 5 7 5 4 6 5 5 3 当年产水 25 36 32 15 31 28 22 12 费用/产水量 0.200 0.194 0.156 0.267 0.194 0.179 0.227 0.25 第一年