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

钢管订购和运输问题一代码和结果 钢管订购和运输问题一代码和结果 PAGE / NUMPAGES 钢管订购和运输问题一代码和结果 function f=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=[ ; ; ; ; ; ; ]; options=optimset( LargeScale,off ,Algorithm ,active-set ,MaxFunEvals ,50000);%,Tolx,; [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 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)*; %铺设时的运输费 end f function [c,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: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; 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; 结果： A1 A2 A3 A4 A5 A6 A7 A8 A9 A10 A11  表一 Si 到 Ai 的最小费用（单位：万元 / 单位） S1 S2 S3 S4 S5 S6 S7 38 111 121 156 146 156 166 86 96 131 121 131 141 92 142 82 62 57 62 76 96 146 86 51 33 51 66 A12 106 156 96 61 51 45 56 A13 A14 128 178 118 83 73 11 26 A15 142 192 132 97 87 28 2 表二 各厂的生产量及总费用（生产量可小于 500）（单位：单位 、万元） S1 S2 S3 S4 S5 S6 S7 总费用 800 800 1000 0 245 1.2754 106 =0 500  表三 各厂的生产量及总费用（单位：单位 、万元） S1 S2 S3 S4S5 S6 S7 总费用 800 800 1000 0 0 1.2786 106 800 800 1000 0 500 1.2797 106 表四 问题一订购和运输方案（不足 1km 的按整数计）（单位：单位 、万元） S1 S2 S3 S4 S5 S6 S7 订购量 800 800 1000 0 1190 1181 0 A1 0 0 0 0 0 0 0 A2 0 179 0 0 0 0 0 A3 0 137 141 0 230 0 0 A4 149 74 79 0 166 0 0 A5 186 110 116 0 203 0 0 A6 200 0 0 0 0 0 0 A7 265 0 0 0 0 0 0 A8 A9 A10 A1