钢管订购和运输问题一代码和结果_精品文档文档格式.docx

钢管订购和运输问题一代码和结果_精品文档文档格式.docx

《钢管订购和运输问题一代码和结果_精品文档文档格式.docx》由会员分享，可在线阅读，更多相关《钢管订购和运输问题一代码和结果_精品文档文档格式.docx（7页珍藏版）》请在冰豆网上搜索。

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:

L=[104301750606194205201680480300220210420500];

f=0;

f=f+N（i）*x（i,j）*C（i,j）;

%运输费和成本费

end

14

f=f+（rl（i）*（rl（i）+1）/2+（L（i）-rl（i））*（L（i）-rl（i）+1）/2）*0.1;

%铺设时的运输费

f

function[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）;

fori=2:

a（i）=a（i）+N（j）*x（j,i）;

ceq（i-1）=a（i）-rl（i）+rl（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.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

245.2

98.6

171.6

181.6

216.6

206.6

226.6

38

111

121

156

146

166

20.5

64.6

105.5

139.6

130.5

140.5

150.5

3.1

86

96

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

99.2

92

142

82

62

57

76

51

33

66

106

61

45

56

171.2

76.2

26.2

38.2

128

178

118

83

73

11

26

192

132

97

87

28

2

表二各厂的生产量及总费用（生产量可小于500）（单位：

单位、万元）

总费用

800

1000

1190.5

1135.5

245

表三各厂的生产量及总费用（单位：

=0

1180.5

>

500

1185.5

885.5

表四问题一订购和运输方案（不足1km的按整数计）（单位：

订购量

1190

1181

179

137

230

149

74

79

186

110

116

203

200

265

300

664

176

415

333

621

165

订购总量

5171