数学模型结业课程设计求解钢管订购和运输问题

《数学模型》课程结业论文

题 目 院 系 专 业 学 号 学生姓名 任课教师

钢管订购与运输

理学院 信息与计算科学

单锋

沈阳航空航天大学

2013年4月

任务及要求

任 务 书

[要求]

1、将所给的问题翻译成汉语；

2、给论文起个题目（名字或标题） 3、根据任务来完成数学模型论文；

4、论文书写格式要求按给定要求书写；

5、态度要认真，要独立思考，独立完成任务；

6、论文上交时间：5月30日前（要求交纸质论文和电子文档）。 7、严禁抄袭行为，若发现抄袭，则成绩记为“不及格”。

[任务]

钢管订购和运输

要铺设一条A1?A2???A15的输送天然气的主管道, 如图一所示(见下页)。经筛选后可以生产这种主管道钢管的钢厂有S1,S2,?S7。图中粗线表示铁路，单细线表示公路，双细线表示要铺设的管道(假设沿管道或者原来有公路，或者建有施工公路)，圆圈表示火车站，每段铁路、公路和管道旁的阿拉伯数字表示里程(单位km)。

为方便计，1km主管道钢管称为1单位钢管。 一个钢厂如果承担制造这种钢管，至少需要生产500个单位。钢厂Si在指定期限内能生产该钢管的最大数量为si个单位，钢管出厂销价1单位钢管为pi万元，如下表：

i 1 800 160 2 800 155 3 1000 155 4 2000 160 5 2000 155 6 2000 150 7 3000 160 si pi

1单位钢管的铁路运价如下表：

里程(km)

≤300 301～350 351～400 401～450 451～500 I

任务及要求

运价(万元)

20 23 26 29 32 里程(km) 501～600 601～700 701～800 801～900 运价(万元)

1000km以上每增加1至100km运价增加5万元。

37 44 50 55 901～1000 60 公路运输费用为1单位钢管每公里0.1万元（不足整公里部分按整公里计算）。

钢管可由铁路、公路运往铺设地点（不只是运到点A1,A2,?,A15，而是管道全线）。

（1）请制定一个主管道钢管的订购和运输计划，使总费用最小（给出总费用)。

（2）请就（1）的模型分析：哪个钢厂钢管的销价的变化对购运计划和总费用影响最大，哪个钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大，并给出相应的数字结果。

（3）如果要铺设的管道不是一条线，而是一个树形图，铁路、公路和管道构成网络，请就这种更一般的情形给出一种解决办法，并对图二按（1）的要求给出模型和结果。

II

任务及要求

290 S3 S2 690 1200 720 202 1100 20 195 1150 600 3060 5 10 194 10 31 201 A8A1A711S1 12 42 70 10 170 520 88 462 S5 10 220 A11 S4 320 160 70 30 70 62 S6 110 420 A15 500 A14 20 30 S7 20 690 160 A13 210 A12 480 680 11A9300 A10 450 80 A5 606 2 750 A4 3 A3 301 104 A2 A1 A6A11 205 图一 290 S3 S2 690 1200 720 A16 202 1100 20 195 1150 600 3060 5 10 194 A6 A5 606 750 2 A4 A3 301 A2 10 31 S1 12 42 170 520 88 S4 A18 160 130 160 70 A20 100 30 260 S6 70 (A21) 320 A19 62 110 420 30 20 S7 20 690 A15

190 462 500 A14 A17 70 10 S5 10 220 A13 210 A12 201 A8A1205 A7 300 A10 A9 680 480 A11 450 80 图二 3 104 A1

III

任务及要求

评语：

成绩

任课教师签字 年 月 日

成 绩 评 定 单

IV

摘要

摘 要

本文讨论了在铺设天然气管道的过程中如何合理订购与运输钢管以使总费用最小的优化问题。

问题一是在一定约束条件下以钢管订购和运输的总费用为目标函数的非线性规划问题。总费用由订购钢管的总费用、从钢厂到站点运输钢管的总费用及从站点开始铺设钢管的总费用三部分组成。订购钢管的总费用和从钢厂到各站点运输钢管的总费用分别通过在各厂购买量与各厂出厂销价和各厂购买量与从各钢厂到各站点运输单位钢管的最小费用的线性规划运算得到。从站点开始铺设钢管的总费用通过等差数列求和得到。在求从钢厂到站点的运输钢管的总费用时，关键是采用弗洛伊德算法，用MATLAB软件编程求出单位钢管从各钢厂运往各站点最小运输费用。利用LINGO软件求解此模型，得到钢管订购与运输的最小费用。 问题二是对问题一模型的灵敏度分析，通过控制变量法的方法即每次只让一家钢厂的销价或生产线发生变化并且每次的变化是相同，分别得出各变量对购运计划的影响。 问题三是对问题一的推广，要铺设的管道不是一条线，而是一个树形图，铁路、公路和管道构成网络，在问题一的模型中又增加了一些约束条件和变量，同时在目标函数中增加相应的铺设费用。利用LINGO软件编程求解新的模型。

关键词：非线性规划；弗洛伊德算法；灵敏度分析；

I

目录

目 录

钢管订购与运输 ................................................................................................................... 1 1.1 问题提出 .................................................................................................................... 1 1.2 模型假设 .................................................................................................................... 3 1.3 符号说明 .................................................................................................................... 4 1.4 问题一的模型建立：求钢管订购和运输最小运费 ................................................ 4 1.5 问题一的求解 ............................................................................................................ 5 2.1问题二的模型建立：钢管销价变化对购运计划的影响。 ..................................... 6 2.2问题二的求解 ............................................................................................................. 6 3.1问题三的模型建立：直线管道向管道网变化时的购运计划 ................................. 7 3.2题三的求解 ................................................................................................................. 9 4．优缺点改进 ................................................................................................................. 9 5．参考文献 ................................................................................................................... 10 6．附录............................................................................................................................ 10

I

数学模型课程结业论文

钢管订购与运输

1.1 问题提出

要铺设一条A1?A2???A15的输送天然气的主管道, 如图一所示(见下页)。经筛选后可以生产这种主管道钢管的钢厂有S1,S2,?S7。图中粗线表示铁路，单细线表示公路，双细线表示要铺设的管道(假设沿管道或者原来有公路，或者建有施工公路)，圆圈表示火车站，每段铁路、公路和管道旁的阿拉伯数字表示里程(单位km)。

为方便计，1km主管道钢管称为1单位钢管。 一个钢厂如果承担制造这种钢管，至少需要生产500个单位。钢厂Si在指定期限内能生产该钢管的最大数量为si个单位，钢管出厂销价1单位钢管为pi万元，如下表：

i 1 800 160 2 800 155 3 1000 155 4 2000 160 5 2000 155 6 2000 150 7 3000 160 si pi

1单位钢管的铁路运价如下表：

里程(km) 运价(万元)

里程(km) 501～600 601～700 701～800 801～900 运价(万元)

1000km以上每增加1至100km运价增加5万元。

1

≤300 20 301～350 351～400 401～450 451～500 23 26 29 32 901～1000 60 37 44 50 55

数学模型课程结业论文

公路运输费用为1单位钢管每公里0.1万元（不足整公里部分按整公里计算）。

钢管可由铁路、公路运往铺设地点（不只是运到点A1,A2,?,A15，而是管道全线）。

（1）请制定一个主管道钢管的订购和运输计划，使总费用最小（给出总费用)。

（2）请就（1）的模型分析：哪个钢厂钢管的销价的变化对购运计划和总费用影响最大，哪个钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大，并给出相应的数字结果。

（3）如果要铺设的管道不是一条线，而是一个树形图，铁路、公路和管道构成网络，请就这种更一般的情形给出一种解决办法，并对图二按（1）的要求给出模型和结果。

290 S3 S2 690 1200 720 202 1100 20 195 1150 600 3060 5 10 194 10 31 201 A8A1A711S1 12 42 70 10 170 520 88 462 S5 10 220 A11 S4 320 160 70 30 70 62 S6 110 420 A15 500 A14 20 30 S7 20 690 160 A13 210 A12 480 680 11A9300 A10 450 8A5 606 2 750 A4 3 A3 301 104 A2 0 A6A11 205 图一 A1

2

数学模型课程结业论文

290 S3 S2 690 1200 720 A16 202 1100 20 195 1150 600 3060 5 10 194 A6 A5 606 10 12 31 300 A10 A9 680 480 S1 42 170 520 88 S4 A18 160 130 160 70 A20 100 30 260 S6 70 (A21) 320 A19 110 62 420 30 20 S7 20 690 A15

190 462 500 A14 A17 70 10 S5 10 220 A13 210 A12 A11 201 A8A1205 A7 450 80 2 750 A4 A3 301 A2 图二 3 104 A1

1.2 模型假设

1.模型只考虑钢管销价费用和钢管从钢管厂运送到铺设点的钢管运费，而不考虑其它费用，

如不计换车、转站的时间和费用，不计装卸费用等。

2.要铺设的管道侧有公路，可运输所需钢管。

2.钢管单价与订购量、订购次数、订购日期无关,即在钢管订购与运输过程中，钢管的单价保持不变。

3.将每一单位的管道所在地看成一个需求点，向以单位管道的所在地运输钢管即向一个点运输钢管。 4.钢管在运送和使用中没有损耗。

5.不计运输时由于运输工具出现故障等意外事故引起工期延误造成损失。

3

数学模型课程结业论文

ti?0或1(i=1,..,7)

3.2题三的求解

S1 800 S2 800 S3 1000 S4 0 S5 1303 S6 2000 S7 0 费用 1406330 得到最优最小费用为W?1406330万元。

4．优缺点改进

由于总费用由订购费用和运输费用部分组成，运输费又由一般线路上的运输费和铺设管道上的运输费组成. 利用求网络中最短路径的弗洛伊德方法得到新的算法，可对含多种权重计算方式的网络进行搜索，得出最小费用路径（最短路径），算出两点之间的最优路径，进而根据非线性规划，借助于Lingo软件求解即可求出相应的结果.

1．优点

1）本问题中运用了求网络中最短路径的弗洛伊德思想，改进和修改得到新的算法，可对含多种权重计算方式的网络进行搜索，算出两点之间的最优路径，计算结果准确，从而得出相应的购运单价的矩阵.

2）本问题构造出的模型算法较简单，也可以运用相应的其他编程软件来得到比较满意的结果.

3）本模型计算步骤清晰，借助于Lingo软件求解，可靠性较高. 2．缺点

1）由于题意中不考虑铁路公路间转运的中转费用，也不限制转运次数，因此在算法设计中存在着考虑不周全的缺限，如我们考虑是先通过铁路再通过公路到铺设点，但这不一定是最小费用路径，有可能先通过公路，然后经铁路再经公路运到铺设点，费用更少，这里没有理论证明.

2) 问题二要求根据问题一的分析,指出哪家钢厂销价的变化对购运计划和总费用影响最大，哪家钢厂钢管产量的上限的变化对购运计划和总费用的影响最大，并给出相应的数字结果. 这个问题属于规划问题的灵敏度分析，在此模型中，只是通过说明销价增加一万元，减少一万元来说明，并没有给出一般的理论说明。

3．模型改进

这个数学模型可以应用于西部开发中“西气东送”问题，当然，西部开发中“西气东送”问题远比我们的假设还要复杂的多，但无论如何，他们的本质一样，我们可将本问题运用于时间的变化等范围的推广。文还可以把问题1归结为网络最小费用流问题，建立了线性和非线性最小费用流模型，并运用相应的解法和分支定界法求解，简洁，层次分明。

9

数学模型课程结业论文

5．参考文献

1．《数学模型》 单峰 朱丽梅 田贺民 国防工业出版社 2.《运筹学教程》 胡运权 清华大学出版社

3.《MATLAB程序设计与应用》 刘卫国 中国水利水电出版社

6．附录

附录一：

Floyd算法函数在matlab下的M函数文件如下： function [D,path]=floyd(a) n=size(a,1);

D=a;path=zeros(n,n); for i=1:n

for j=1:n

if D(i,j)~=inf path(i,j)=j; end end end

for k=1:n

for i=1:n

for j=1:n

if D(i,k)+D(k,j)

附录二：

ab=[1 1 2 3 4 5 6 7 8 9 10 11 12 13 15 16 17 18 19 20 20 22 23];

bb=[ 14 15 15 16 19 18 23 24 10 10 11 15 13 14 16 17 19 19 20 21 22 23 24]; w=[20 202 1200 690 690 462 70 30 450 80 1150 1100 306 195 720 520 170 88 160 70 320 160 290];

ab1=[1 2 4 5 6 7 8 9 10 11 14 15 16 17 18 33 34 35];

bb1=[19 20 21 22 23 24 25 26 27 28 29 30 31 31 19 24 31 32]; w1=[3 2 600 10 5 10 12 42 70 10 10 62 30 20 104 31 110 20]; a=sparse(ab,bb,w); a(24,24)=0; a=a+a'; a=full(a);

10

数学模型课程结业论文

for i=1:24

for j=1:24

if(a(i,j)==0&i~=j) a(i,j)=inf; end end end

[D,path]=floyd(a);

a1=sparse(ab1,bb1,w1); a1(35,35)=0; a1=a1+(a1)'; a1=full(a1); for i=1:35

for j=1:35

if(a1(i,j)==0&i~=j) a1(i,j)=inf; end end end

[D1,path1]=floyd(a1);

%距离转换为费用的程序

D1=D1*0.1; %把公路最短距离换算成公路最少费用 for k=1:300 m1(k)={k}; end

for k=1:50

m2(k)={300+k}; m3(k)={350+k}; m4(k)={400+k}; m5(k)={450+k}; end

for k=1:100

m6(k)={500+k}; m7(k)={600+k}; m8(k)={700+k}; m9(k)={800+k}; m0(k)={900+k}; end

for i=1:24

for j=1:24 %把铁路最短距离换算成铁路最少费用 switch D(i,j)

11

数学模型课程结业论文

case 0

D(i,j)=0; case m1

D(i,j)=20; case m2

D(i,j)=23; case m3

D(i,j)=26; case m4

D(i,j)=29; case m5

D(i,j)=32; case m6

D(i,j)=37; case m7

D(i,j)=44; case m8

D(i,j)=50; case m9

D(i,j)=55; case m0

D(i,j)=60; otherwise

D(i,j)=ceil((D(i,j)-1000)/100)*5+60; end end end

%c矩阵表示七个钢管生产厂到十五个铺设节点之间的距离，先把它们都设成20000（任意一个钢管厂到任意一个铺设节点之间的距离不会超过20000），然后用 for 循环求出最小值

c=20000*ones(7,15);

for i=1:7 %7个钢管生产厂 for k=18:32 个铺设节点 for j=8:24 %7个钢管生产厂和17个中转点，i=1，表示第一个钢管生产厂，j=8，表示第一个中转点

if c(i,k-17)>D(i,j)+D1(j-7,k) c(i,k-17)=D(i,j)+D1(j-7,k); %对于所有中转点，在铁路网和公路网上的下标相差8 end end end

12

数学模型课程结业论文

end

for i=1:7 for k=18:32

if c(i,k-17)>D(i,1)+D1(33,k) c(i,k-17)=D(i,1)+D1(33,k); 3代表第一个钢管生产厂S1点 end

if c(i,k-17)>D(i,6)+D1(34,k) c(i,k-17)=D(i,6)+D1(34,k); 4代表第六个钢管生产厂S6点 end

if c(i,k-17)>D(i,7)+D1(35,k) c(i,k-17)=D(i,7)+D1(35,k); 5代表第七个钢管生产厂S7点 end end

%因为S1,S6,S7这三个钢管厂有公路直接连接到铺设节点，所以把这三个点单独处理 end

13

本文来源：https://www.bwwdw.com/article/bzpg.html