410.7 400.3000 380.2000 361.6000 301.0000 285.5000 276.0000 266.2000 234.2000 212.0000 188.0000 206.0000 226.2000 228.0000 242.0000;
415.7 405.3000 385.2000 366.6000 306.0000 290.5000 281.0000 271.2000 234.2000 212.0000 201.0000 195.0000 176.2000 161.0000 178.0000;
435.7 425.3000 405.2000 386.6000 326.0000 310.5000 301.0000 291.2000 259.2000 236.0000 226.0000 216.0000 198.2000 186.0000 162.0000];
options=optimset('LargeScale','off','Algorithm' ,'active-set','MaxFunEvals' ,50000);%,'Tolx',1.0000e-032);
[x,f]=fmincon('myfun',x0,[],[],[],[],vlb,[],'mycon',options,C,N,s); for i=1:7 for j=1:15
m(i)=m(i)+N(i)*x(i,j); end; end; x,,m,f
b=(f-1278600)/1278600*(s(t)+50)/50 toc
16
function f=myfun(XX,C,N,s) %问题1的目标函数 x=XX(1:7,1:15); rl=XX(8,1:15);
L=[104 301 750 606 194 205 201 680 480 300 220 210 420 500]; f=0; for i=1:7 for j=1:15
f=f+N(i)*x(i,j)*C(i,j);%运输费和成本费 end; end; for i=1:14
f=f+(rl(i)*(rl(i)+1)/2+(L(i)-rl(i))*(L(i)-rl(i)+1)/2)*0.1;%铺设时的运输费 end; f
function [c,ceq]=mycon(XX,C,N,s) %问题1的约束条件
x=XX(1:7,1:15); rl=XX(8,1:15);
17
L=[104 301 750 606 194 205 201 680 480 300 220 210 420 500]; m=zeros(1,7);
a=zeros(1,15); cc=0; for i=1:7 for j=1:15
m(i)=m(i)+N(i)*x(i,j); end; c(i)=m(i)-s(i); cc=cc+m(i); end; for i=1:14 c(i+7)=rl(i)-L(i); end; for i=2:14 for j=1:7
a(i)=a(i)+N(j)*x(j,i); end;
ceq(i-1)=a(i)-rl(i)+rl(i-1)-L(i-1);end; t1=0;t2=0;
18
for i=1:7
t1=t1+N(i)*x(i,1); t2=t2+N(i)*x(i,15); end;
ceq(14)=t1-rl(1); ceq(15)=rl(15); ceq(16)=cc-5171; 附录2
function f=result2%求解问题3 x0=zeros(10,21);vlb=zeros(10,21); m=zeros(1,7); N=[1 1 1 0 1 1 0]; tic
%每公里钢管从Si到达Ai站点的最小费用
C=[330.7 320.3000 300.2000 258.6000 198.0000 180.5000 163.1000 181.2000 224.2000 252.0000 256.0000 266.0000 281.2000 288.0000 302.0000 220 255 260 265 275 290;
370.7 360.3000 345.2000 326.6000 266.0000 249.6000 241.0000 226.2000 269.2000 297.0000 301.0000 311.0000 326.2000 333.0000 347.0000 265 300 305 310 320 335;
385.7 375.3000 355.2000 336.6000 276.0000 260.5000 251.0000
19
241.2000 203.2000 237.0000 241.0000 251.0000 266.2000 273.0000 287.0000 199 240 245 250 260 270;
420.7 410.3000 395.2000 376.6000 316.0000 299.6000 291.0000 276.2000 244.2000 222.0000 211.0000 221.0000 236.2000 243.0000 257.0000 240 210 215 220 230 240;
410.7 400.3000 380.2000 361.6000 301.0000 285.5000 276.0000 266.2000 234.2000 212.0000 188.0000 206.0000 226.2000 228.0000 242.0000 230 187 205 205 220 230;
415.7 405.3000 385.2000 366.6000 306.0000 290.5000 281.0000 271.2000 234.2000 212.0000 201.0000 195.0000 176.2000 161.0000 178.0000 230 200 187 194 170 150;
435.7 425.3000 405.2000 386.6000 326.0000 310.5000 301.0000 291.2000 259.2000 236.0000 226.0000 216.0000 198.2000 186.0000 162.0000 255 225 210 215 192 186];
options=optimset('LargeScale','off','Algorithm' ,'active-set','MaxFunEvals' ,50000);%,'Tolx',1.0000e-032);
[x,f]=fmincon('myfun2',x0,[],[],[],[],vlb,[],'mycon2',options,C,N); for i=1:7 for j=1:21
20