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

1、精品文档functionf=result(t)% 求解问题 1 tic;x0=zeros(8,15);vlb=zeros(8,15);m=zeros(1,7);s=800 800 1000 2000 2000 2000 3000;s(t)=s(t)-50;N=1110110;% 每公里钢管从 Si到达 Ai 站点的最小费用C=330.7320.3000300.2000258.6000198.0000180.5000163.1000181.2000224.2000252.0000256.0000 266.0000281.2000288.0000302.0000;370.7360.3000345.

2、2000326.6000266.0000249.6000241.0000226.2000269.2000297.0000301.0000 311.0000326.2000333.0000347.0000;385.7375.3000355.2000336.6000276.0000260.5000251.0000241.2000203.2000237.0000241.0000 251.0000266.2000273.0000287.0000;420.7410.3000395.2000376.6000316.0000299.6000291.0000276.2000244.2000222.000021

3、1.0000 221.0000236.2000243.0000257.0000;410.7400.3000380.2000361.6000301.0000285.5000276.0000.精品文档266.2000234.2000212.0000188.0000 206.0000226.2000228.0000242.0000;415.7405.3000385.2000366.6000306.0000290.5000281.0000271.2000234.2000212.0000201.0000 195.0000176.2000161.0000178.0000;435.7425.3000405.

4、2000386.6000326.0000310.5000301.0000291.2000259.2000236.0000226.0000 216.0000198.2000186.0000162.0000;options=optimset('LargeScale' ,'off' ,'Algorithm','active-set' ,'MaxFunEvals',50000); %,'Tolx',1.0000e-032);x,f=fmincon('myfun',x0,vlb,'my

5、con',options,C,N,s);for i=1:7for j=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=104 301 750 606 194 205 201 680 480 300 220 210 420 500;f=0;for i=1:7for j=1:15f=f+N(i)*x(i,j)*C(i,j);% 运输费和成本费end

6、endfor i=1:14f=f+(rl(i)*(rl(i)+1)/2+(L(i)-rl(i)*(L(i)-rl(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=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:7for j=1:15m(i)=m(i)+N(i)*x(i,j);endc(i)=m(i)

7、-s(i);cc=cc+m(i);endfor i=1:14c(i+7)=rl(i)-L(i);endfor i=2:14for j=1:7a(i)=a(i)+N(j)*x(j,i);endceq(i-1)=a(i)-rl(i)+rl(i-1)-L(i-1);endt1=0;t2=0;for i=1:7.精品文档t1=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;结果：表一Si 到 Ai 的最小费用（单位：万元 / 单位）S1S2S3S4S5S6S7170.7215.72

8、30.7260.7255.7265.7275.7A1160.3205.3220.3250.3245.3255.3265.3A2140.2190.2200.2235.2225.2235.2245.2A398.6171.6181.6216.6206.6216.6226.6A438111121156146156166A5.精品文档20.564.6105.5139.6130.5140.5150.5A63.18696131121131141A721.271.286.2116.2111.2121.2131.2A864.2114.248.284.279.284.299.2A9921428262576276A

9、10961468651335166A111061569661514556A12121.2171.2111.276.271.226.238.2A1312817811883731126A141421921329787282A15表二各厂的生产量及总费用（生产量可小于500）（单位：单位 、万元）S1S2S3S4S5S6总费用S7800800100001190.51135.52451061.2754=0>500表三各厂的生产量及总费用（单位：单位、万元）S1S2S3S4S5S6总费用S7800800100001190.51180.501061.2786800800100001185.5885.55001061.2797表四问题一订购和运输方案（不足1km的按整数计）（单位：单位、万元）S1S2S3S4S5S6S7订购量80080010000119011810A10000000.精品文档017900000A20137141023000A31497479016600A41861101