。
=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号井只能在三年中的某一年打造或者不打造，故应有
?Xij <=1;
j?132) 每年的费用不能超过计划即
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〈=60
3) 每年的水量应满足要求，水量有三部分构成：现有水井的产水量，新建水井的产水量， 管 道铺好后的管道水量。
现有水井产水量可根据2001——2009数据拟合出2010——2014年的，程序编码及拟合图 见附录1，拟合结果如下图所示：
编号 产水量 年份 1号井 2号井 21.4472 1.069 20.2456 0 19.0439 0 17.8422 0 16.6406 0 2010 2011 2012 2013 2014 -可编辑修改-
。
3号井 4号井 总水量 9.1306 11.8286 43.4754 7.0456 4.5071 31.7983 4.9606 0 24.0045 2.8756 0 20.7178 0.7906 0 17.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+N1 W2=31.7983+N2 W3=24.0045+N3
-可编辑修改-
。
W4=20.7178+N4+Q W5=17.4312+N5+Q
满足，W1>=150, W2>=160, W3>=170, W4>=180, W5>=190 4)每年的铺管道费取整且总管道不小20公里即
Pj=0.66Q^0.51*Li Pj取整
L1+L2+L3 >= 20
6 模型求解
将上述模型输入LINGO可得到【2】
Local optimal solution found.
Objective value: 172.0000 Extended solver steps: 308 Total solver iterations: 10226
Variable Value Reduced Cost Z1 55.00000 0.000000 Z2 60.00000 0.000000 Z3 57.00000 0.000000 X11 1.000000 5.000000 X21 0.000000 7.000000
-可编辑修改-
。
X31 1.000000 5.000000 X41 0.000000 4.000000 X51 0.000000 6.000000 X61 1.000000 5.000000 X71 1.000000 5.000000 X81 0.000000 P1 35.00000 X12 0.000000 X22 1.000000 X32 0.000000 X42 0.000000 X52 0.000000 X62 0.000000 X72 0.000000 X82 0.000000 P2 53.00000 X13 0.000000 X23 0.000000 X33 0.000000 X43 0.000000 X53 1.000000 X63 0.000000 -可编辑修改-
3.000000 1.000000 5.000000 7.000000 5.000000 4.000000 6.000000 5.000000 5.000000 3.000000 1.000000 5.000000 7.000000 5.000000 4.000000 6.000000 5.000000