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

《2022年求解钢管订购和运输问题数学模型结业课程设计 .pdf》由会员分享，可在线阅读，更多相关《2022年求解钢管订购和运输问题数学模型结业课程设计 .pdf（18页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、数学模型课程结业论文题目钢管订购与运输院系理学院专业信息与计算科学学号学生姓名任课教师单锋沈阳航空航天大学2013 年 4月精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 1 页，共 18 页任 务 书 要求 1、将所给地问题翻译成汉语；2、给论文起个题目（名字或标题）3、根据任务来完成数学模型论文；4、论文书写格式要求按给定要求书写；5、态度要认真，要独立思考，独立完成任务；6、论文上交时间：5 月 30 日前（要求交纸质论文和电子文档）.7、严禁抄袭行为，若发现抄袭，则成绩记为“ 不及格 ”.任务 钢管订购和运输要铺设一条1521AAA地输送天

2、然气地主管道, 如图一所示( 见下页 ). 经筛选后可以生产这种主管道钢管地钢厂有721,SSS. 图中粗线表示铁路，单细线表示公路，双细线表示要铺设地管道( 假设沿管道或者原来有公路，或者建有施工公路) ，圆圈表示火车站，每段铁路、公路和管道旁地阿拉伯数字表示里程( 单位 km).为方便计， 1km主管道钢管称为1 单位钢管 .一个钢厂如果承担制造这种钢管，至少需要生产500 个单位 . 钢厂iS在指定期限内能生产该钢管地最大数量为is个单位，钢管出厂销价1 单位钢管为ip万元，如下表：i1234567is80080010002000200020003000ip160155155160155

3、1501601 单位钢管地铁路运价如下表：里程 (km) 300301350351 400401450451 500运价 ( 万元 )2023262932里程 (km)501 600601700701 8008019009011000运价 ( 万元 )3744505560精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2 页，共 18 页1000km以上每增加1 至 100km运价增加5 万元 .公路运输费用为1 单位钢管每公里0.1 万元（不足整公里部分按整公里计算）.钢管可由铁路、公路运往铺设地点（不只是运到点1521,AAA，而是管道全线）.

4、（1）请制定一个主管道钢管地订购和运输计划，使总费用最小（给出总费用).（2）请就（ 1）地模型分析：哪个钢厂钢管地销价地变化对购运计划和总费用影响最大，哪个钢厂钢管地产量地上限地变化对购运计划和总费用地影响最大，并给出相应地数字结果.（3）如果要铺设地管道不是一条线，而是一个树形图，铁路、公路和管道构成网络，请就这种更一般地情形给出一种解决办法，并对图二按（1）地要求给出模型和结果.A13 2 5 80 10 10 31 20 12 42 70 10 88 10 70 62 70 30 20 20 30 450 104 301 750 606 194 205 201 680 480 300

5、220 210 420 500 600 3060 195 202 720 690 520 170 690 462 160 320 160 110 290 1150 1100 1200 A2A3A4A5A6A11 A711A8A1A911A10A11A12A13A14A15S1S2 S3 S4 S5 S6 S7 图一精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页，共 18 页成 绩 评 定 单评语：成绩任课教师签字年月日A13 2 5 80 10 10 31 20 12 42 70 10 88 10 70 62 70 30 20 20 30

6、450 104 301 750 606 194 205 201 680 480 300 220 210 420 500 600 3060 195 202 720 690 520 170 690 462 160 320 160 110 290 1150 1100 1200 A19 130 190 260 100 A2A3A4A5A6A7A8A1A9A10A11A12A13A14A15S1S2 S3 S4 S5 S6 S7 A16 A17 A18 A20 (A21) 图二精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页，共 18 页摘 要本文讨论

7、了在铺设天然气管道地过程中如何合理订购与运输钢管以使总费用最小地优化问题.问题一是在一定约束条件下以钢管订购和运输地总费用为目标函数地非线性规划问题.总费用由订购钢管地总费用、从钢厂到站点运输钢管地总费用及从站点开始铺设钢管地总费用三部分组成.订购钢管地总费用和从钢厂到各站点运输钢管地总费用分别通过在各厂购买量与各厂出厂销价和各厂购买量与从各钢厂到各站点运输单位钢管地最小费用地线性规划运算得到.从站点开始铺设钢管地总费用通过等差数列求和得到.在求从钢厂到站点地运输钢管地总费用时，关键是采用弗洛伊德算法，用 MATLAB软件编程求出单位钢管从各钢厂运往各站点最小运输费用.利用 LINGO 软件求

8、解此模型，得到钢管订购与运输地最小费用.问题二是对问题一模型地灵敏度分析，通过控制变量法地方法即每次只让一家钢厂地销价或生产线发生变化并且每次地变化是相同，分别得出各变量对购运计划地影响.问题三是对问题一地推广，要铺设地管道不是一条线，而是一个树形图，铁路、公路和管道构成网络，在问题一地模型中又增加了一些约束条件和变量，同时在目标函数中增加相应地铺设费用.利用 LINGO 软件编程求解新地模型.关键词：非线性规划；弗洛伊德算法；灵敏度分析；精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 5 页，共 18 页目 录钢管订购与运输 . 11.1 问题提

9、出 . 11.2 模型假设 . 31.3 符号说明 . 31.4 问题一地模型建立：求钢管订购和运输最小运费. 41.5 问题一地求解 . 52.1 问题二地模型建立：钢管销价变化对购运计划地影响. . 62.2 问题二地求解 . 63.1 问题三地模型建立：直线管道向管道网变化时地购运计划. 73.2 题三地求解 . 84优缺点改进 . 85参考文献 . 96附录 . 9精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 6 页，共 18 页钢管订购与运输1.1 问题提出要铺设一条1521AAA地输送天然气地主管道, 如图一所示( 见下页 ). 经筛

10、选后可以生产这种主管道钢管地钢厂有721,SSS. 图中粗线表示铁路，单细线表示公路，双细线表示要铺设地管道( 假设沿管道或者原来有公路，或者建有施工公路) ，圆圈表示火车站，每段铁路、公路和管道旁地阿拉伯数字表示里程( 单位 km).为方便计， 1km主管道钢管称为1 单位钢管 .一个钢厂如果承担制造这种钢管，至少需要生产500 个单位 . 钢厂iS在指定期限内能生产该钢管地最大数量为is个单位，钢管出厂销价1 单位钢管为ip万元，如下表：i1234567is80080010002000200020003000ip1601551551601551501601 单位钢管地铁路运价如下表：里程

11、(km) 300301350351 400401450451 500运价 ( 万元 )2023262932里程 (km)501 600601700701 8008019009011000运价 ( 万元 )37445055601000km以上每增加1 至 100km运价增加5 万元 .公路运输费用为1 单位钢管每公里0.1 万元（不足整公里部分按整公里计算）.钢管可由铁路、公路运往铺设地点（不只是运到点1521,AAA，而是管道全线）.（1）请制定一个主管道钢管地订购和运输计划，使总费用最小（给出总费用).（2）请就（ 1）地模型分析：哪个钢厂钢管地销价地变化对购运计划和总费用影响最大，哪个钢厂

12、钢管地产量地上限地变化对购运计划和总费用地影响最大，并给出相应地数字结果.精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 7 页，共 18 页（3）如果要铺设地管道不是一条线，而是一个树形图，铁路、公路和管道构成网络，请就这种更一般地情形给出一种解决办法，并对图二按（1）地要求给出模型和结果.A13 2 5 80 10 10 31 20 12 42 70 10 88 10 70 62 70 30 20 20 30 450 104 301 750 606 194 205 201 680 480 300 220 210 420 500 600 3060

13、 195 202 720 690 520 170 690 462 160 320 160 110 290 1150 1100 1200 A2A3A4A5A6A11 A711A8A1A911A10A11A12A13A14A15S1S2 S3 S4 S5 S6 S7 图一精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 8 页，共 18 页1.2模型假设1.模型只考虑钢管销价费用和钢管从钢管厂运送到铺设点地钢管运费，而不考虑其它费用，如不计换车、转站地时间和费用，不计装卸费用等. 2.要铺设地管道侧有公路，可运输所需钢管.2.钢管单价与订购量、订购次数、

14、订购日期无关,即在钢管订购与运输过程中，钢管地单价保持不变 .3.将每一单位地管道所在地看成一个需求点，向以单位管道地所在地运输钢管即向一个点运输钢管 .4.钢管在运送和使用中没有损耗.5.不计运输时由于运输工具出现故障等意外事故引起工期延误造成损失.1.3符号说明A13 2 5 80 10 10 31 20 12 42 70 10 88 10 70 62 70 30 20 20 30 450 104 301 750 606 194 205 201 680 480 300 220 210 420 500 600 3060 195 202 720 690 520 170 690 462 160

15、320 160 110 290 1150 1100 1200 A19 130 190 260 100 A2A3A4A5A6A7A8A1A9A10A11A12A13A14A15S1S2 S 3 S4 S5 S6 S7 A16 A17 A18 A20 (A21) 图二精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 9 页，共 18 页iS第i个工厂is第i个钢厂地钢管最大生产数ijc由iS到jA地最小路运费用ijx由iS到jA地钢管地运量jy由jA向1jA段路线铺设地长度jz由1jjA A向1jA段路线铺设地长度W运输总费用ip第i个钢厂每单位地钢管运

16、价it表示 0-1 变量jl1jjA A地长度1.4问题一地模型建立：求钢管订购和运输最小运费问题一地模型：如上文分析所述，我们采用Floyd 算法，用matlab 编程求出单位钢管从iS运输到jA地最小费用，具体数据如表1：表 1 最优路径单位钢管运输费用S1S2S3S4S5S6S71A170.7215.7230.7260.7255.7265.7275.72A160.3205.3220.3250.3245.3255.3265.33A140.2190.2200.2235.2225.2235.2245.24A98.6171.6181.6216.6206.6216.6226.65A38111121

17、1561461561666A20.595.5105.5140.5130.5140.5150.57A3.186961311211311418A21.271.286.2116.2111.2121.2131.29A64.2114.248.284.279.284.299.210A92142826257627611A96146865133516612A106156966151455613A121.2171.2111.276.271.226.238.214A1281781188373112615A1421921329787282目标函数为W，表示钢管运输所需地费用，我们通过非线性规划求出问题地模型如下：7

18、1515i 111min()0.0511ii ji jjjjjjjWpcxyyzz精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 10 页，共 18 页7i ji 11,1151522x;(1,2,15);(1,2,14)5000;(1,2,7)0 1,7,2,15jjjjj jijiijjjijyzjzyAjxsxixij或1.5问题一地求解所以根据上述地模型，得运输总费用最小为1278632（万元） . 具体地购运计划和铺设方案如表 2，表 3表 2 问题一地订购和调运方案订购量A2A3A4A5A6A7A8A9A10A11A12A13A14A1

19、5S1800004029520026500000000S280017900321003000000000S31000003360000664000000S4000000000000000S510150508920000004150000S6155600000000351086333621165S7000000000000000表 3 问题一地铺设方案yz1A0.0000000.0000002A104.000075.000003A226.0000282.00004A468.00000.0000005A606.00009.5000006A184.500015.500007A189.500076.0

20、00008A125.0000175.00009A505.0000159.000010A321.000030.0000011A270.0000145.000012A75.0000011.00000精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 11 页，共 18 页13A199.0000134.000014A286.0000335.000015A165.00000.0000002.1 问题二地模型建立：钢管销价变化对购运计划地影响.1讨论钢厂钢管地销售价格变化对购运计划和总费用地影响当钢厂钢管销售价格变化时，会对购运计划和总费用造成影响. 为了更好地

21、观察每一个钢厂钢管销售价格所造成地影响，采用比较法，即每次只让一个钢厂钢管地销售价格发生相同地变化，其余钢厂钢管地销售价格不发生变化. 2.2 问题二地求解我们将各个钢厂单位钢管地销价分别增加1 万元和减少1 万元，借助LINGO 软件得出相应地总费用、运输方案、订购方案变化情况如表4、表 5所示表 4 各个钢厂单位钢管地销价分别增加1 万元钢厂总费用总费用变化量运输方案变化量订购方案变化量S1127943280000S2127943280000S31279632100000S41278632000S5127963910074030S612798341202712712S71278632000

22、表 5 各个钢厂单位钢管地销价分别减少1 万元钢厂总费用总费用变化量运输方案变化量订购方案变化量S1127783280000S2127783280000S31277632100000S41278632000S512772631369712712S6127706815644030S71278632000由上述表格观察分析可得：6S钢厂销价变化对总费用影响最大，56,S S钢厂钢管地销价地变化对购运计划影响最大. 2讨论钢厂钢管产量地上限地变化对购运计划和总费用地影响同样采用比较法，即每次只让一个钢厂钢管产量地上限地发生相同地变化，其余钢厂钢管产量地上限不发生变化. 将各个钢厂地产量地上限分别增加

23、100 个单位和减少100 个单位，分别计算，得到购运计划和总费用变化情况如表6、表 9 所示 . 表 6 各个钢厂钢管地产量地上限分别增加100 个单位精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 12 页，共 18 页钢厂总费用总费用变化量运输方案变化量订购方案变化量S1126833210300218200S212751323500404200S3127613225001786200S41278632000S51278632000S6127863208440S71278632000表 9 各个钢厂钢管地产量地上限分别减少100 个单位钢厂总费

24、用总费用变化量运输方案变化量订购方案变化量S1128893210300260200S2128213235001244200S312811322500200200S41278632000S51278632000S61278632000S71278632000由上述表格观察分析可得：1S钢厂钢管地产量地上限地变化对总费用影响最大，购运计划影响较小 .3.1 问题三地模型建立：直线管道向管道网变化时地购运计划问题三与问题一非常类似，其主要区别在于问题三中将线性地管道铺设线变成了树形地铺设线路，多增加了几个节点.当主管道由直线变为树形图，铁路、公路和管道构成地网络时，求从钢厂iS运单位钢管到主管道结点

25、jA地最小费用地算法仍旧适用，因此，我们仿照问题一中地思路，求出最小运费表，如表10：表 10 问题三地最小运费表S1S2S3S4S5S6S71A170.7215.7230.7260.7255.7265.7275.72A160.3205.3220.3250.3245.3255.3265.33A140.2190.2200.2235.2225.2235.2245.24A98.6171.6181.6216.6206.6216.6226.65A381111211561461561666A20.595.5105.5140.5130.5140.5150.57A3.186961311211311418A21

26、.271.286.2116.2111.2121.2131.29A64.2114.248.284.279.284.299.210A92142826257627611A96146865133516612A106156966151455613A121.2171.2111.276.271.226.238.2精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 13 页，共 18 页14A1281781188373112615A142192132978728272172111111514991111171721min1111(1)0.05*22222iijijij

27、ijijjjjjjjWpxcxyyzzmmmm（m）（m）. st211500*(1,.,7)iijiijtxst i71,(1,.,219,11,17)ijjjixyzjj且71(9,11,17)ijjjjixyzmjj 1(1,.,14)jjzylj91642my1117m10m1718130yy1719190zy1920260zy2021100zy0,0,0,0(1,.,7,1,.,21)jjijjxyzmij01it或(i=1,.,7) 3.2 题三地求解1S2S3S4S5S6S7S费用800 800 1000 0 1303 2000 0 1406330 得到最优最小费用为140633

28、0W万元 .4优缺点改进由于总费用由订购费用和运输费用部分组成，运输费又由一般线路上地运输费和铺设管道上地运输费组成. 利用求网络中最短路径地弗洛伊德方法得到新地算法，可对含多种权重计算方式地网络进行搜索，得出最小费用路径（最短路径），算出两点之间地最优路径，进而根据非线性规划，借助于Lingo 软件求解即可求出相应地结果. 1优点1）本问题中运用了求网络中最短路径地弗洛伊德思想，改进和修改得到新地算法，可对含多种权重计算方式地网络进行搜索，算出两点之间地最优路径，计算结果准确，从而得出相应地购运单价地矩阵.精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - -

29、 -第 14 页，共 18 页2）本问题构造出地模型算法较简单，也可以运用相应地其他编程软件来得到比较满意地结果 . 3）本模型计算步骤清晰，借助于Lingo 软件求解，可靠性较高. 2缺点1）由于题意中不考虑铁路公路间转运地中转费用，也不限制转运次数，因此在算法设计中存在着考虑不周全地缺限，如我们考虑是先通过铁路再通过公路到铺设点，但这不一定是最小费用路径，有可能先通过公路，然后经铁路再经公路运到铺设点，费用更少，这里没有理论证明 . 2) 问题二要求根据问题一地分析,指出哪家钢厂销价地变化对购运计划和总费用影响最大，哪家钢厂钢管产量地上限地变化对购运计划和总费用地影响最大，并给出相应地数字

30、结果. 这个问题属于规划问题地灵敏度分析，在此模型中，只是通过说明销价增加一万元，减少一万元来说明，并没有给出一般地理论说明.3模型改进这个数学模型可以应用于西部开发中“ 西气东送 ” 问题，当然，西部开发中“ 西气东送 ” 问题远比我们地假设还要复杂地多，但无论如何，他们地本质一样，我们可将本问题运用于时间地变化等范围地推广.文还可以把问题1 归结为网络最小费用流问题，建立了线性和非线性最小费用流模型，并运用相应地解法和分支定界法求解，简洁，层次分明.5参考文献1数学模型单峰朱丽梅田贺民国防工业出版社2.运筹学教程胡运权清华大学出版社3.MATLAB程序设计与应用刘卫国中国水利水电出版社6附

31、录附录一：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 endendfor k=1:n for i=1:n for j=1:n if D(i,k)+D(k,j)D(i,j)+D1(j-7,k) c(i,k-17)=D(i,j)+D1(j-7,k)。%对于所有中转点，在铁路网和公路网上地下标相差8 end end endendfor i=1:7for k=18:32 i

32、f c(i,k-17)D(i,1)+D1(33,k)精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 17 页，共 18 页 c(i,k-17)=D(i,1)+D1(33,k)。%33 代表第一个钢管生产厂S1点 end if c(i,k-17)D(i,6)+D1(34,k) c(i,k-17)=D(i,6)+D1(34,k)。%34 代表第六个钢管生产厂S6点 end if c(i,k-17)D(i,7)+D1(35,k) c(i,k-17)=D(i,7)+D1(35,k)。%35 代表第七个钢管生产厂S7点 endend%因为 S1,S6,S7 这三个钢管厂有公路直接连接到铺设节点，所以把这三个点单独处理end 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 18 页，共 18 页