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

钢管订购和运输问题一
代码和结果
公司内部档案编码：[OPPTR-OPPT28-OPPTL98-OPPNN08]
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=[1 1 1 0 1 1 0];
%每公里钢管从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;
结果：
表一 i S 到i A 的最小费用（单位：万元/单位）
1S 2S 3S 4S 5S 6S 7S
1A
2A
3A
4A
5A 38
111
121
156
146
156
166
6A
7A 86
96
131
121
131
141
8A 9A
10A
92 142 82 62 57 62 76 11A 96 146 86 51 33 51 66 12A
106
156
96
61
51
45
56
13A
14A
128 178 118 83 73 11 26 15A
142
192
132
97
87
28
2
表二各厂的生产量及总费用（生产量可小于500）（单位：单位、万
元）
表三各厂的生产量及总费用（单位：单位、万元）
表四问题一订购和运输方案（不足1km的按整数计）（单位：单
位、万元）。