钢管订购和运输问题一代码和结果_精品文档.docx
- 文档编号:973102
- 上传时间:2022-10-14
- 格式:DOCX
- 页数:7
- 大小:120.44KB
钢管订购和运输问题一代码和结果_精品文档.docx
《钢管订购和运输问题一代码和结果_精品文档.docx》由会员分享,可在线阅读,更多相关《钢管订购和运输问题一代码和结果_精品文档.docx(7页珍藏版)》请在冰豆网上搜索。
functionf=result(t)
%求解问题1
tic;
x0=zeros(8,15);vlb=zeros(8,15);
m=zeros(1,7);
s=[80080010002000200020003000];
s(t)=s(t)-50;
N=[1110110];
%每公里钢管从Si到达Ai站点的最小费用
C=[330.7320.3000300.2000258.6000198.0000180.5000163.1000181.2000224.2000252.0000256.0000266.0000281.2000288.0000302.0000;
370.7360.3000345.2000326.6000266.0000249.6000241.0000226.2000269.2000297.0000301.0000311.0000326.2000333.0000347.0000;
385.7375.3000355.2000336.6000276.0000260.5000251.0000241.2000203.2000237.0000241.0000251.0000266.2000273.0000287.0000;
420.7410.3000395.2000376.6000316.0000299.6000291.0000276.2000244.2000222.0000211.0000221.0000236.2000243.0000257.0000;
410.7400.3000380.2000361.6000301.0000285.5000276.0000266.2000234.2000212.0000188.0000206.0000226.2000228.0000242.0000;
415.7405.3000385.2000366.6000306.0000290.5000281.0000271.2000234.2000212.0000201.0000195.0000176.2000161.0000178.0000;
435.7425.3000405.2000386.6000326.0000310.5000301.0000291.2000259.2000236.0000226.0000216.0000198.2000186.0000162.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);
fori=1:
7
forj=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
functionf=myfun(XX,C,N,s)
%问题1的目标函数
x=XX(1:
7,1:
15);
rl=XX(8,1:
15);
L=[104301750606194205201680480300220210420500];
f=0;
fori=1:
7
forj=1:
15
f=f+N(i)*x(i,j)*C(i,j);%运输费和成本费
end
end
fori=1:
14
f=f+(rl(i)*(rl(i)+1)/2+(L(i)-rl(i))*(L(i)-rl(i)+1)/2)*0.1;%铺设时的运输费
end
f
function[c,ceq]=mycon(XX,C,N,s)
%问题1的约束条件
x=XX(1:
7,1:
15);
rl=XX(8,1:
15);
L=[104301750606194205201680480300220210420500];
m=zeros(1,7);
a=zeros(1,15);
cc=0;
fori=1:
7
forj=1:
15
m(i)=m(i)+N(i)*x(i,j);
end
c(i)=m(i)-s(i);
cc=cc+m(i);
end
fori=1:
14
c(i+7)=rl(i)-L(i);
end
fori=2:
14
forj=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;
fori=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;
结果:
表一到的最小费用(单位:
万元/单位)
170.7
215.7
230.7
260.7
255.7
265.7
275.7
160.3
205.3
220.3
250.3
245.3
255.3
265.3
140.2
190.2
200.2
235.2
225.2
235.2
245.2
98.6
171.6
181.6
216.6
206.6
216.6
226.6
38
111
121
156
146
156
166
20.5
64.6
105.5
139.6
130.5
140.5
150.5
3.1
86
96
131
121
131
141
21.2
71.2
86.2
116.2
111.2
121.2
131.2
64.2
114.2
48.2
84.2
79.2
84.2
99.2
92
142
82
62
57
62
76
96
146
86
51
33
51
66
106
156
96
61
51
45
56
121.2
171.2
111.2
76.2
71.2
26.2
38.2
128
178
118
83
73
11
26
142
192
132
97
87
28
2
表二各厂的生产量及总费用(生产量可小于500)(单位:
单位、万元)
总费用
800
800
1000
0
1190.5
1135.5
245
表三各厂的生产量及总费用(单位:
单位、万元)
总费用
=0
800
800
1000
0
1190.5
1180.5
0
>500
800
800
1000
0
1185.5
885.5
500
表四问题一订购和运输方案(不足1km的按整数计)(单位:
单位、万元)
订购量
800
800
1000
0
1190
1181
0
0
0
0
0
0
0
0
0
179
0
0
0
0
0
0
137
141
0
230
0
0
149
74
79
0
166
0
0
186
110
116
0
203
0
0
200
0
0
0
0
0
0
265
0
0
0
0
0
0
0
300
0
0
0
0
0
0
0
664
0
0
0
0
0
0
0
0
176
176
0
0
0
0
0
415
0
0
0
0
0
0
0
86
0
0
0
0
0
0
333
0
0
0
0
0
0
621
0
0
0
0
0
0
165
0
订购总量
5171
总费用
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 钢管 订购 运输 问题 代码 结果 精品 文档