1、LargeScale,offAlgorithm ,active-setMaxFunEvals ,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:L=104 301 750 606 194 205
2、 201 680 480 300 220 210 420 500;f=0; f=f+N(i)*x(i,j)*C(i,j);%运输费和成本费end14 f=f+(rl(i)*(rl(i)+1)/2+(L(i)-rl(i)*(L(i)-rl(i)+1)/2)*0.1;%铺设时的运输费ffunction c,ceq=mycon(XX,C,N,s)%问题1的约束条件a=zeros(1,15);cc=0; c(i)=m(i)-s(i); cc=cc+m(i); c(i+7)=rl(i)-L(i);for i=2: a(i)=a(i)+N(j)*x(j,i); ceq(i-1)=a(i)-rl(i)+rl
3、(i-1)-L(i-1);t1=0;t2=0; t1=t1+N(i)*x(i,1); t2=t2+N(i)*x(i,15);ceq(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.2245.298.6171.6181.6216.6206.6226.63811112115614616620.564.6105.5139.6130.
4、5140.5150.53.1869613114121.271.286.2116.2111.2121.2131.264.2114.248.284.279.299.29214282625776513366106614556171.276.226.238.2128178118837311261921329787282表二 各厂的生产量及总费用(生产量可小于500)(单位:单位 、万元)总费用80010001190.51135.5245表三 各厂的生产量及总费用(单位:=01180.55001185.5885.5表四 问题一订购和运输方案(不足1km的按整数计)(单位:订购量119011811791372301497479186110116203200265300664176415333621165订购总量5171