钢管订购和运输问题一代码和结果

1、钢管订购和运输问题一代码和结果function f=result(t) %求解问题1tic;x0=zeros(8,15);vlb=zeros(8/15);m=zeros(L7);s=800 800 1000 2000 2000 2000 3000;s(t)=s(t)-50;N = 1110 110;%每公里钢管从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.0

2、000;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;385.7 375.3000 355.2000 336.6000 276.0000 260.5000 251.0000 241.2000 203.2000 237.0000 241.0000 251.0000 266.2000 273.0000 287.0000;420.7 410.3000 395.2000 376.6000

3、 316.0000 299.6000 291.0000 276.2000 244.2000 222.0000 211.0000 221.0000 236.2000 243.0000 257.0000;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

4、 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(,LargeScale7off,/Algorithm, /active-sef/MaxFunEvals1,50000);%/ TolxM.OOOOe-032);xrfl=fm

5、incon(,myfunrxOffJr/vlbr/,mycon,roptionsfC/N/s);for i=l:7for j = l:15m(i)=m(i)+N(i)*x(i,j);end endb=(f-1278600)/1278600*(s(t)+50)/50 toe function f=myfun(XXJC/N,s)%问题1的目标函数x=XX(l:7fl:15);rl=XX(8,l:15);L=104 301 750 606 194 205 201 680 480 300 220 210 420 500;f=0;for i = l:7for j = l:15f=f+N(i)*x(i,j

6、)*C(i,j);% 运输费和成本费end endfor i=l:14f=f+(rl(i)*M(i)+l)/2+(L(i)川(i)*(L(i)-rl(i)+l)/2)*0.1;% 铺设时的运输费 endffunction c/ceq = mycon(XX/C/Nfs)%问题1的约束条件x=XX(l:7/l:15);rl=XX(8,l:15);L=104 301 750 606 194 205 201 680 480 300 220 210 420 500;m=zeros(l,7);a=zeos(L15);cc=O;for i=l:7for j = l:15m(i)=m(i)+N(i)*x(ij

7、);endc(i)=m(i)-s(i);cc=cc+m(i);endfor i=l:14c(i+7)=rl(i)-L(i);endfor i=2:14for j = 1:7a(i)=a(i)+N(j)*x(j,i);endceq(i-l)=a(i)-rl(i) + rl(i-l)-L(i-l);endtl=0;t2=0;for i=l:7tl=tl+N(i)*x(ifl);t2=t2+N(i)*x(ir15);endceq(14)=tl-rl(l);ceq(15)=rl(15);ceq(16)=cc-5171;Command Windowfigure 1iVTUsrlAiRK：A4A5MATA

8、3出AlOAllAt 2A13AUHKfl K2St810 00闭62Mi010a0052BIO】而2；60JW010q00幻nw 055E106(0c“410Q00Ma o00ci01010a0a%iW 0l?T153rott0IT441500054DW 00o(a(Q1730并加42i5700000d4taxinc inHMible dimg to rithlfi the 4e,wli viluc ” fsrHcc -lerss, idar* satisfied e.Mn tb- def cult sis of，be snmiP, “leg： 3】加，za、so5*2ef/freoB-廿卬

9、“2。IS -IM -513功衿*12：-11 ANtt*S do、iJbdotL _iibi劭，的，板工dml33 31313, dsvl4A A2Ui?=oi*-l ” x,J.fzncMt S：I0JW010000幻 UW 0104195(0(如10000M 00Q00I0I010a0a$S 13 02l015：幻3i0(0l：8410000X 313 0000(0(017302：93t$7 510 0000(0(01036490r218 R l n-MM.950A 1212091eT4)R1”30140IC)34umR3结果:表一 s,到4的最小费用（单位：万元/单位）Sis?S3S4

10、工r邑a170. 7215. 7230.7260.7255. 7265. 7275.7a2160. 3205. 3220.3250.3245. 3255. 3265.3A3140. 2190.2200.2235.2225.2235.2245.2A498.6171.6181.6216.6206.6216.6226.6438111121156146156166A620.564.6105.5139.6130.5140.5150.543. 18696131121131141421.271.286.2116.2111.2121.2131.2464.2114.248.284.279.284.299.2Ao

11、921428262576276小961468651335166A 21061569661514556小121.2171.2111.276.271.226.238.2A,12817811883731126小1421921329787282表二 各厂的生产量及总费用（生产量可小于500）（单位：单位、万元）4&$336S7总费用800800100001190. 51135. 52451.2754xl06表三 各厂的生产量及总费用（单位：单位、万元）Ss2S3S*SsSb当总费用二0800800100001190.51180. 501.2786xl06500800800100001185. 5885.55001.2797 xlO6表四 问题一订购和运输方案（不足1km的按整数计）（单位：单位、万元）Sis,3S5Sb订购量80080010000119011810a0000000A?01790