。

a=polyfit(x,y,1) t=1:1:14 s=polyval(a,t) hold on

plot(t,s,'r-','linewidth',2)

一号井水量模拟图线

grid

二号井： x=1:1:9

y=[21.5,15.9,11.8,8.7,6.5,4.8,3.5,2.6, 2.0] plot(x,y,'k.','markersize',25) a=polyfit(x,y,3)

-可编辑修改-

2000）（万吨）（年份减去

。

t=1:1:14 s=polyval(a,t) hold on

plot(t,s,'r-','linewidth',2) grid

二号井水量模拟图线

三号井： x=1:1:9

y=[ 27.9,25.8,23.8,21.6,19.5,17.4,15.5,13.3,11.2] plot(x,y,'k.','markersize',25) a=polyfit(x,y,1) t=1:1:14

-可编辑修改-

2000）（万吨）（年份减去

。

s=polyval(a,t) hold on

plot(t,s,'r-','linewidth',2) grid

三号井水量模拟图线

四号井： x=2:1:9

y=[ 46.2,32.6,26.7,23.0,20.0,18.9,17.5,16.3] plot(x,y,'k.','markersize',25)

-可编辑修改-

2000）（万吨）（年份减去

。

a=polyfit(x,y,3) t=2:1:14 s=polyval(a,t) hold on

plot(t,s,'r-','linewidth',2) grid

四号井水量模拟图线

附录二：

min=Z1+Z2+Z3;

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;

-可编辑修改-

2000）（万吨）（年份减去