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

function f=result(t)%求解问题1tic;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=1 1 1 0 1 1 0; %每公里钢管从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; 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 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 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 endx,m,f;b=(f-1278600)/1278600*(s(t)+50)/50tocfunction 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);%运输费和成本费 endend for i=1:14 f=f+(rl(i)*(rl(i)+1)/2+(L(i)-rl(i)*(L(i)-rl(i)+1)/2)*0.1;%铺设时的运输费endffunction 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);endfor i=1:14 c(i+7)=rl(i)-L(i);endfor 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);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;结果：表一 到的最小费用（单位：万元/单位）170.7215.7230.7260.7255.7265.7275.7160.3205.3220.3250.3245.3255.3265.3140.2190.2200.2235.2225.2235.2245.298.6171.6181.6216.6206.6216.6226.63811112115614615616620.564.6105.5139.6130.5140.5150.53.1869613112113114121.271.286.2116.2111.2121.2131.264.2114.248.284.279.284.299.29214282625762769614686513351661061569661514556121.2171.2111.276.271.226.238.2128178118837311261421921329787282表二 各厂的生产量及总费用（生产量可小于500）（单位：单位 、万元）总费用800800100001190.51135.5245表三 各厂的生产量及总费用（单位：单位 、万元）总费用=0800800100001190.51180.50500800800100001185.5885.5500表四 问题一订购和运输方案（不足1km的按整数计）（单位：单位 、万元）订购量80080010000119011810000000001790000001