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

1、functionf=result(t)%求解问题1tic；x0=zeros(8,15);vlb=zeros(8,15);m=zeros(1,7);s=80080010002000200020003000;s(t尸s(t)-50;N=1110110;%每公里钢管从Si到达Ai站点的最小费用C=330.7320.3000300.2000258.6000198.0000180.5000163.1000181.2000224.2000252.0000256.0000266.0000281.2000288.0000302.0000;370.7360.3000345.2000326.6000266.000

2、0249.6000241.0000226.2000269.2000297.0000301.0000311.0000326.2000333.0000347.0000;385.7375.3000355.2000336.6000276.0000260.5000251.0000241.2000203.2000237.0000241.0000251.0000266.2000273.0000287.0000;420.7410.3000395.2000376.6000316.0000299.6000291.0000276.2000244.2000222.0000211.0000221.0000236.200

3、0243.0000257.0000;410.7400.3000380.2000361.6000301.0000285.5000276.0000266.2000234.2000212.0000188.0000206.0000226.2000228.0000242.0000;415.7405.3000385.2000366.6000306.0000290.5000281.0000271.2000234.2000212.0000201.0000195.0000176.2000161.0000178.0000;435.7425.3000405.2000386.6000326.0000310.50003

4、01.0000291.2000259.2000236.0000226.0000216.0000198.2000186.0000162.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);fori=1:7fo

5、rj=1:15m(i)=m(i)+N(i)*x(i,j);endendx,m,f;b=(f-1278600)/1278600*(s(t)+50)/50tocfunctionf=myfun(XX,C,N,s)%'可题1的目标函数x=XX(1:7,1:15);rl=XX(8,1:15);L=104301750606194205201680480300220210420500;f=0;fori=1:7forj=1:15f=f+N(i)*x(i,j)*C(i,j);%运输费和成本费endendfori=1:14f=f+(rl(i)*(rl(i)+1)/2+(L(i)-rl(i)*(L(i)-r

6、l(i)+1)/2)*0.1;%铺设时的运输费endffunctionc,ceq=mycon(XX,C,N,s)%问题1的约束条件x=XX(1:7,1:15);rl=XX(8,1:15);L=104301750606194205201680480300220210420500;m=zeros(1,7);a=zeros(1,15);cc=0;fori=1:7forj=1:15m(i)=m(i)+N(i)*x(i,j);endc(i)=m(i)-s(i);cc=cc+m(i);endfori=1:14c(i+7)=rl(i)-L(i);endfori=2:14forj=1:7a(i)=a(i)+N

7、(j)*x(j,i);endceq(i-1)=a(i)-rl(i)+rl(i-1)-L(i-1);endt1=0;t2=0;fori=1:7t1=t1+N(i)*x(i,1);t2=t2+N(i)*x(i,15);endceq(14)=t1-rl(1);ceq(15)=rl(15);ceq(16)=cc-5171;2c 必 m an 加 113cl 3Cbflha 府! WhmJjw盯K*屿MA7jj城uDXUJillqJAK4工巧SI BIDQ«mX加2dri&4aJiD4DQeE2 mI,H6In1Q1WQ1Q中S3 JigQFIKm1q4m1DJDq.M fiD.口二l

8、gVQi口口dwK 11610I7Tmrn1*4a,-4U»中*a541M<0.0 1."ad4二 FO14力11AMa4E幅 fl中*q1ar0104at：>hT13；w1xE1IIIClm33V.支整用二-ral .J" tT 窜fQirH TAipr +i!if：14nF4 T|产 “旧近丁： rr 小b«uw*nn rwM«< bwi th* *iFhF fwrtivn it (w4Mftww hi*3r5：_m lLL+rt“iM, TO WLTlUE file Ctf&lll 工,Ulf <jf th

9、* funKLijE Fl + L 1一九UdUf 461才把。，口 VllJUM tM drt'EuLC mK rf 1 l>f i lTUt 113»1七的口Li；!/MKUTtuRjnj tr=L+-ri!i d<*inln£ JIM一1比 口 出J1山-'山击出田田<3Jb -Air Ji AZ - Tr«i 1-T-al i.> = LAd-«- 3r 11 I. bfl3J5 -|星 JrLjJiuxKni ihE 叽 T, M. bL» )Hr*Bh«s 14,4 J fel wi

10、T fcP 11i2' -fpcTirt力£；主打斗于!"工11嗣*： n <* fflr ill .后AT阜 ”1注上血总4恭- 他J Efng】fk. uoIIIm is uhHitM生户，干却厢二m七 二 1u HIlfib4"2理IIQ|71Hi*53imQ174.*M 4t舞0力。157XBl?Q«0r910Q40Tf MM c-nL -1kM.则Mi结果：表一Si到A的最小费用（单位：万元/单位）SiS2S3S4S5S6S7A170.7215.7230.7260.7255.7265.7275.7A2160.3205.3220.3

11、250.3245.3255.3265.3A3140.2190.2200.2235.2225.2235.2245.2A98.6171.6181.6216.6206.6216.6226.6A538111121156146156166A620.564.6105.5139.6130.5140.5150.5A3.18696131121131141A821.271.286.2116.2111.2121.2131.2A964.2114.248.284.279.284.299.2A10921428262576276A11961468651335166A121061569661514556AI3121.2171

12、.2111.276.271.226.238.2A|412817811883731126AI51421921329787282表二各厂的生产量及总费用（生产量可小于500）（单位：单位、万元）S1S2S3S4S5S6S7总费用800800100001190.51135.52451.2754106表三各厂的生产量及总费用（单位：单位、万元）S1S2S3S4S5S6S7总费用=0800800100001190.51180.501.2786106>500800800100001185.5885.55001.2797106表四问题一订购和运输方案（不足1km的按整数计）（单位：单位、万元）S1S2S3S4S5S6S7订购量80080010000119011810A0000000A017900000A30137141023000A14974790