数学建模论文范文[1]

11问题的提出位于我国西南地区的某个偏远贫困村，年平均降水量不足20mm，是典型的缺水地区。过去村民的日常生活和农业生产用水一方面靠的是每家每户自行建造的小蓄水池，用来屯积每逢下雨时获得的雨水，另一方面是利用村里现有的四口水井。由于近年来环境破坏，经常是一连数月滴雨不下，这些小蓄水池的功能完全丧失。而现有的四口水井经过多年使用后，年产水量也在逐渐减少，在表1中给出它们在近9年来的产水量粗略统计数字。2009年以来，由于水井的水远远不能满足需要，不仅各种农业生产全部停止，而且大量的村民每天要被迫翻山越岭到相隔十几里外去背水来维持日常生活。为此，今年政府打算着手帮助该村解决用水难的问题。从两方面考虑，一是地质专家经过勘察，在该村附近又找到了8个可供打井的位置，它们的地质构造不同，因而每个位置打井的费用和预计的年产水量也不同，详见表2，而且预计每口水井的年产水量还会以平均每年10%左右的速率减少。二是从长远考虑，可以通过铺设管道的办法从相隔20公里外的地方把河水引入该村。铺设管道的费用为L66Q.0P0.51（万元），其中Q表示每年的可供水量（万吨/年），L表示管道长度（公里）。铺设管道从开工到完成需要三年时间，且每年投资铺设管道的费用为万元的整数倍。要求完成之后，每年能够通过管道至少提供100万吨水。政府从2010年开始，连续三年，每年最多可提供60万元用于该村打井和铺设管道，为了保证该村从2010至2014年这五年间每年分别能至少获得150、160、170、180、190万吨水，请作出一个从2010年起三年的打井和铺设管道计划，以使整个计划的总开支尽量节省（不考虑小蓄水池的作用和利息的因素在内）。表1现有各水井在近几年的产水量（万吨）2表28个位置打井费用（万元）和当年产水量（万吨）编号12345678打井费用57546553当年产水25363215312822122问题的分析题中要求制定一个总费用（决策目标）最小的抗旱（打井，铺设管道）方案，属于优化问题，并且使得该村从2010至2014年这五年间每年分别能至少获得150、160、170、180、190万吨水，每年费用不超过60万元。（此两点为主要约束条件）其他的约束条件有：a.每口井只能在2010年开始，连续三年中的其中一年施工b.铺设管道费用为万元整数倍c.由于河位于与该村相隔20公里外的地方,所以管道总长度不小于20公里d.铺设管道需要3年时间,故前3年管道供水量为0,而第4,5年供水量不小于100万吨。故此模型即为基于以上约束条件的整数规划（最优决策目标）问题。3模型的假设a．忽略小蓄水池的作用和利息因素b．不考虑意外情况导致所需经费增加年份产水量编号2001200220032004200520062007200820091号井32.231.329.728.627.526.125.323.722.72号井21.515.911.88.76.54.83.52.62.03号井27.925.823.821.619.517.415.513.311.24号井46.232.626.723.020.018.917.516.33c．假设井在年初修建且时间很短，修完之后即可利用，管道铺好后即可用于供水d．假设这五年之内村民需水量基本稳定e．假设井供水量呈稳定规律变化，不考虑其他因素对产水量的影响f．从长远利益考虑，打井和铺设管道两个方案应同时协调进行4符号说明Xij0—1变量，表示第i号井在第j年的施工情况，Xij=1第i号井在第j年施工，Xij=0表示不施工Zj第j年的总费用Pj第j年的铺管道费用Lj第j年铺管道公里数Wj第j年的水量Q管道供水量Nj所有新建的水井在第j年的产水量5模型建立决策变量为三年间铺设管道和打井的总费用。0—1变量Xij表示i号井j年是否施工，为1则施工,产生费用，Pj表示第j年的铺路费用。所以第j年的总费用Zj=5*X1j+7*X2j+5*X3j+4*X4j+6*X5j+5*X6j+5*X7j+3*X8j+Pj三年费用minZ=Z1+Z2+Z3=5*X11+7*X21+5*X31+4*X41+6*X51+5*X61+5*X71+3*X81+P1+5*X12+7*X22+5*X32+4*X42+6*X52+5*X62+5*X72+3*X82+P2+5*X13+7*X23+5*X33+4*X43+6*X53+5*X63+5*X73+3*X83+P3约束条件：1）由于第i号井只能在三年中的某一年打造或者不打造，故应有431jXij=1;2)每年的费用不能超过计划即Z1=5*X11+7*X21+5*X31+4*X41+6*X51+5*X61+5*X71+3*X81+P1;Z2=5*X12+7*X22+5*X32+4*X42+6*X52+5*X62+5*X72+3*X82+P2;Z3=5*X13+7*X23+5*X33+4*X43+6*X53+5*X63+5*X73+3*X83+P3;Z1〈=60，Z2〈=60，Z3〈=603)每年的水量应满足要求，水量有三部分构成：现有水井的产水量，新建水井的产水量，管道铺好后的管道水量。现有水井产水量可根据2001——2009数据拟合出2010——2014年的，程序编码及拟合图见附录1，拟合结果如下图所示：编号产水量年份201020112012201320141号井21.447220.245619.043917.842216.64062号井1.06900003号井9.13067.04564.96062.87560.79064号井11.82864.5071000总水量43.475431.798324.004520.717817.4312新建水井产水量：第一年：N1=25*X11+36*X21+32*X31+15*X41+31*X51+28*X61+22*X71+12*X81;第二年：N2=25*X12+36*X22+32*X32+15*X42+31*X52+28*X62+22*X72+12*X82+25*X11*0.9+36*X21*0.9+32*X31*0.9+15*X41*0.9+31*X51*0.9+28*X61*0.9+22*X71*0.9+12*X81*0.9第三年：N3=25*X13+36*X23+32*X33+15*X43+31*X53+28*X63+22*X73+12*X83+25*X12*0.9+36*X22*0.9+32*X32*0.9+15*X42*0.9+31*X52*0.9+28*X62*0.9+22*X72*0.9+12*X82*0.9+25*X11*0.81+36*X21*0.81+32*X31*0.81+15*X41*0.81+31*X51*0.81+28*X61*0.81+22*X71*0.81+12*X81*0.81;第四年：N4=N3*0.9第五年：N5=N3*0.81管道水量：前三年为0，后两年为Q故每年的总水量W1=43.4754+N1W2=31.7983+N2W3=24.0045+N3W4=20.7178+N4+Q5W5=17.4312+N5+Q满足，W1=150,W2=160,W3=170,W4=180,W5=1904)每年的铺管道费取整且总管道不小20公里即Pj=0.66Q^0.51*LiPj取整L1+L2+L3=206模型求解将上述模型输入LINGO可得到【2】Localoptimalsolutionfound.Objectivevalue:172.0000Extendedsolversteps:308Totalsolveriterations:10226VariableValueReducedCostZ155.000000.000000Z260.000000.000000Z357.000000.000000X111.0000005.000000X210.0000007.000000X311.0000005.000000X410.0000004.000000X510.0000006.000000X611.0000005.000000X711.0000005.000000X810.0000003.000000P135.000001.000000X120.0000005.000000X221.0000007.000000X320.0000005.000000X420.0000004.000000X520.0000006.000000X620.0000005.000000X720.0000005.000000X820.0000003.000000P253.000001.000000X130.0000005.000000X230.0000007.000000X330.0000005.000000X430.0000004.000000X531.0000006.0000006X630.0000005.000000X730.0000005.000000X830.0000003.000000P351.000001.000000Q100.00000.000000L15.0643550.000000L27.6688800.000000L37.3794880.000000W1150.47540.000000W2164.09830.000000W3174.07450.000000W4255.78080.000000W5238.98780.000000RowSlackorSurplusDualPrice1172.0000-1.00000020.000000-1.00000030.000000-1.00000040.000000-1.00000050.0000000.00000060.0000000.00000070.0000000.00000081.0000000.00000090.0000000.000000100.0000000.000000110.0000000.000000121.0000000.00000013-0.4462500E-070.00000014-0.6757500E-070.00000015-0.6502500E-070.000000160.11272260.000000175.0000000.0000001855.000000.0000001960.000000.000000200.0000000.000000213.0000000.0000002257.000000.000000230.0000000.000000240.0000000.000000250.0000000.000000260.0000000.000000270.0000000.000000280.0000000.000000290.47540000.0000007304.0983000.000000314.0745000.0000003275.780800.0000003348.987820.000000即最小总费用为172万元第一年花费20万元打造1,3,6,7号井；花费35万元铺管道5.06公里，共计55万元；第二年花费7万元打造2号井，花费53万元铺管道7.669公里，共计60万元；第三年花费6万元打造5号井，花费51万元铺管道7.379公里，共计577结果分析由结果可知第一年打井1，3，6，7号。产生水量150.47万吨。由各井的产水量可知无论是减少井量，或是替换成其他的井，在保证费用不增加的情况下都会使产水量减小，所以第一年只能打井1，3，6，7号。第二年新增水井2号，总水量164.098，可供替换的井为4，5，7号，与2号水量之差分别为21，5，24皆大于4万吨，故也无法满足水量只能打2号井。同理第三年也只能打5号井。这样方案费用是最小的。8方案评价1）本文把所解决的问题归结为优化问题，建立的数学模型清晰合理。2）运用MATLAB和LINGO软件处理数据和进行运算，降低运算量，简单易行，有很大的可操作性。且所得数据较为合理可靠。3）运用0—1模型解题，全面可靠4）但在实