钢管订购和运输数学建模毕业论文.doc

钢管订购和运输数学建模论文钢管订购和运输摘 要本文建立了一个运输问题的最优化模型。通过对图(一)的分析，我们首先直观地将路线分成两段，将图分为两个子图建立了模型一, 利用分支定界法求得总费用最优解为1279496万元。然后对模型一进行优化，得到全线的最优模型二，求得总费用最优解为1278632万元。通过对最优模型二的分析，我们得出钢厂S1的上限产量和钢厂S6的销价的的变化对运购计划和总费用的影响最大，并给出了数据结果。我们利用截取和连接的方法将树形图转化成为对线性图进行分析，并给出了一般的解决方法。对图(二)给出的具体模型，类似与问题一,分别建立了模型三和模型四,求得最优解分别为1408859.4和1403948万元一、 问题的提出已知有7个钢厂，可生产输送天然气主管道的钢管，用Si表示(i=1,2,7)。现有个地点（A1,A2,A15），沿着这个地点铺设一条输送天然气的主管道。为方便计，1km主管道称为单位钢管。 一个钢厂如果承担制造这种钢管，至少需要生产500个单位。钢厂在指定期限内能生产该钢管的最大数量为个单位，钢管出厂销价1单位钢管为万元，如下表：1234567800800100020002000200030001601551551601551501601单位钢管的铁路运价如下表：里程(km)300301350351400401450451500运价(万元)2023262932里程(km)5016006017007018008019009011000运价(万元)37445055601000km以上每增加1至100km运价增加5万元。公路运输费用为1单位钢管每公里0.1万元（不足整公里部分按整公里计算）。钢管可由铁路、公路运往铺设地点（不只是运到点，而是管道全线）。（1）请制定一个主管道钢管的订购和运输计划，使总费用最小（给出总费用)。（2）请就（1）的模型分析：哪个钢厂钢管的销价的变化对购运计划和总费用影响最大，哪个钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大，并给出相应的数字结果。（3）如果要铺设的管道不是一条线，而是一个树形图，铁路、公路和管道构成网络，请就这种更一般的情形给出一种解决办法，并对图二按（1）的要求给出模型和结果。二、问题的分析 该问题是图论中运输问题的最优化问题。 经过分析，我们认为总费用可分为两种费用进行求解，分别为：（1）销价和运输钢管至管道结点Ai(i=1,2,15)的总费用（称为成本费用）（2）铺设过程中的运输费用； 要解决此问题，我们认为有两点关键：（1）如何求出Ai（i=1,2,15）至Si(i=1,2,7)的最小成本费用；（2）如何调整使得各路径满足题中的最优指标。 针对上述问题，我们分别运用了图上作业法、枚举法、逐次修正法、重绕最小生成树法等方法，在综合考虑算法的精度和算法的复杂度后，我们选择了图上作业法、枚举法、逐次修正法对模型进行逐次优化，直至求得最优解。三、模型的基本假设及符号说明（一） 基本假设：（1）运输方式的改变所花费用包含在运费中；（2）铁路线上任意两点可以直达，不需中途转车，即铁路线上两点间运费按线路总长计算；（3）假设一单位钢管可由任意长度钢管组成，购买钢管可以非整数单位购买;（4）不考虑其它外界因素对费用的影响；（5） 钢管在铺设时，先将钢管运到结点处，再由结点处向左右两方相邻结点铺设；（6） 在Si厂购买钢管要么为零,要么至少为500单位。（二）符号说明： xij：从钢厂Si运到结点Aj的单位钢管数； fij：单位钢管从钢厂Si运到结点Aj的的最少成本费用； tj：第Aj点与Aj+1点间的路线长度，j=1,2,14； yj：从结点Aj开始沿管道向右铺设的路线长度, j=1,2,3,15； w：钢管订购和运输的总费用； mi: 钢厂Si的最小产量，mi=500，i=1,2,3,4,5,6,7 ni: 钢厂Si的产量上限，ni=si，i=1,2, ，7 bi: 图一中各节点（见附录七） Vi: 图二中各节点（见附录八）四、模型的建立与求解（一）问题一及其求解：针对图（一），我们首先采用图上作业法对所给图进行分析，利用枚举法，我们求出一单位钢管由钢管厂Si运输至管道结点i的最小成本费用，具体数据如下表：表一 单位钢管由Si运输至i的最小成本费用 (单位：万元)S1S2S3S4S5S6S7A1330.7370.7385.7420.7410.7415.7435.7A2320.3360.3375.3410.3400.3405.3425.3A3300.2345.2355.2395.2380.2385.2405.2A4258.6326.6336.6376.6361.6366.6386.6A5198266276316301306326A6180.5250.5260.5300.5285.5290.5310.5A7163.1241251291276281301A8181.2226.2241.2276.2266.2271.2291.2A9224.2269.2203.2244.2234.2234.2259.2A10252297237222212212237A11256301241211188201226A12266311251221206195216A13281.2326.2266.2236.2226.2176.2198.2A14288333273243228161186A15302347287257242178162 这样，原问题即转化为通常的运输问题。我们进一步对上表数据进行分析，得到结论：结论1：将4、S5两列数据进行对比，发现同一行中S4列数据均大于S5列。由此得出，应优先考虑钢厂S5，只有当对S5的需求量超过其最大产量时，才需要考虑钢厂S4。 进一步对原图进行分析可得：结论2：由于A1点只有通过A2点才能与钢厂Si连接,故模型中不需考虑往A1运输钢管。模型(一)：对表一数据进一步分析。将各行数据由小到大排序，发现由A2至A9各行中费用最小的前四个数据均在S1,S2,S3,S5列中，由A10至A15各行中费用最小的前三个数据均在S5,S6,S7中，直观感觉到：A1至A9需要来自S6及S7等后一部分钢厂的钢管的可能性很小，而且A9至A15不太可能需要S1、S2、S3等前一部分钢厂生厂的钢管，故把原图分为两个子图考虑，分别为A1至A9和A9至A15，即分为两个子模型，分别进行最优化，得到两个子模型(非线性规划模型)如下： 子模型1： 子模型2： 由于这两个子模型属于非线性规划问题，且第一个约束条件很复杂，利用结论1、2，采用分支定界法用LINGO软件分别进行编程（见附录一、二）得到： w1=859629.3 w2=419866.7模型一的总费用为w=w1+w2=1279496（万元）。模型（二）：由于前面的模型直观将管道分为两端，即原图分为两个子图，可能存在一定的误差。下面对模型（一）进行优化，即对主管道建立统一模型，如下：同样利用结论1、2并采用分支定界法，利用LINGO软件编程（见附录三）运行得出:最优解W=1278632万元，各节点向右铺设的管道单位数：Y1=0 Y2=175 Y3=282 Y4=0 Y5=10 Y6=16 Y7=76 Y8=175 Y9=159 Y10=30 Y11=145 Y12=11 Y13=34 Y14=335 表二 各钢厂定购计划：S1S2S3S4S5S6S7数量800800100001237.51333.50表三 运输计划如下（bi含义参见附图七）:路 线 运量S1-b7-b6-b5-A5-A4; 335S1-b7-b6-A6; 200S1-A7; 265;S2-b8-b4-b2-b1-A2; 179S2-b8-b4-b2-b3-A3; 131.3S2-b8-S1-b7-b6-b5-A5-A4; 116S2-b8-S1-b7-b6-b5-A5; 73.7S2-b8-A8; 300 S3-b9-b8-b4-b2-b3-A3; 319S3-b9-b8-S1-b7-b6-b5-A5-A4; 11S3-b9-A9; 604S5-b12-b11-b10-b9-b8-b4-b2-b3-A3; 57.7S5-b12-b11-b10-b9-b8-S1-b7-b6-b5-A5; 542.3S5-b12-b11-b10-A10; 222.5S5-b12-A11; 415S6-b16-b15-b13-b11-b10-A10; 128.8S6-b16-b15-b13-b14-A12; 86S6-b16-b15-A13; 333S6-A14; 621S6-b16-b17-A15; 165（二）问题（2）的求解 通过模型（二）的求解，我们确定对图（一）仅需要S1、S2、S3、S5、S6承担生产任务即可取得最优解。对模型所用程序进行灵敏度分析，并具体考虑当钢厂钢管的销价增加1万元或产量上限增加1单位时，购运计划和总费用的变化情况，结果如下表：表四 钢厂Si的钢管产量上限增加1单位对总费用的影响S1S2S3S4S5S6S7总费用1278529127859712786071278632127863212786321278632减少量10335250000则可以发现：钢厂S1的钢管产量上限的变动对购运计划和总费用影响最大。表五 钢厂Si的钢管销价增加1万元对总费用的影响S1S2S3S4S5S6S7总费用1279432127943212796321278632127963913798341278632增加量80080010000100712020则可以发现：钢厂S6的钢管销价的微小变动对购运计划和总费用影响最大。（三） 问题（3）的求解 1、一般模型 经过对图（一）的分析求解可以看出，订购及运输钢管的总费用可由各段所需费用求和得到。若要铺设的管道是一个树形图， 则可以将其转化为线性管道进行分析。具体解题步骤如下：（1）运用图论的最小权匹配法（简单图可利用枚举法），求出从各钢厂定购并运输一单位钢管至主管道各结点的最小成本费用值。（2）将树型图转化为线性图： 在树形图中取其最长的线形段，称为主干线形段；将图中剩余分支截取，并通过一个虚拟段（长度t为零）联接至主干线形段上，组成一个新的线性管道；如，图二中A9至A16段，可将该段接到A15 点，再分别从原图A9、A16点连一条长度为0的公路，而在A15点之间连一条长度为0的公路并增加一新的结点A9，同时约定A15点不向右铺设管道。其它各段类似处理。（3）按照问题（一）的思路建立模型进行求解。建立如下非线性规划模型：其中K表示钢厂的总数，L+1表示管道节点总数。2、问题三图（二）的求解图（二）为树形图，将其转化为类似图（一）的线性图。利用枚举法，求出一单位钢管由钢管厂Si运输至管道结点i的最小成本费用，具体数据如下表：表六 单位钢管由Si运输至i的最小成本费用 (单位：万元)S1S2S3S5S6S7A1330.7370.7385.4410.7410.7435.7A2320.3360.3375400.3400.3425.3A3300.2345.2355.2380.2385.2405.2A4258.6326.6336.6361.6361.6386.6A5198216276301301326A6180.5250.5260.5285.5290.5310.5A7163.1241251276278.1301A8181.2226.2241.2266.2266.2291.2A9224.2269.2203.2234.2234.2259.2A10252297237212211237A11256301241188201224A12266311251206187216A13281.2326.2266.2226.2166.2198.2A14288333273228161186A15302347287242178162A16220265199230230255A17255300240187197223A18260305245200183210A19265310250205186215A20275320260220160192A21285330270230150186模型（三） 首先类似于模型一，按分段的思想建立模型。由A1A9、A16为一段可得如下模型：由A21-A20-A19-A17-A11-A12-A13-A14-A15-A9-A10-A11-A17-A18为一段可得如下模型：运用LINGO软件进行编程（见附录四、五）运行结果为： 费用w1=69127.2 w2=539732.2 则总费用为：w=w1+w2=1408859.4模型（四）： 对整个图进行考虑，可得如下模型：利用LINGO软件编程（见附录六）运行得出最优解w= 1403948万元，同时得出下表数据：表七 各钢厂定购计划：S1S2S3S4S5S6S7数量80080010000130320000表八 运输计划如下（Vi含义参见附录八）:路 线 运量S1-V7-V6-V5-A5; 335S1-V7-V6-A6; 200S1-A7; 265;S2-V8-V4-V2-V1-A2; 179S2-V8-V4-V2-V3-A3; 171S2-V8-S1-V7-V6-V5-A5; 150S2-V8-A8; 300 S3-A16-V8-V4-V2-V3-A3; 336S3-A16-A9; 664S5-A17-V10-V9-A16-V8-V4-V2-V3-A3; 1S5-A17-V10-V9-A16-V8-V4-A4; 468S5-A17-V10-V9-A16-V8-S1-V7-V6-V5-A5; 131S5-A17-V10-V9-A10; 218S5-A17-A11; 380S5-A17; 105S6-A20-A18-V10-V9-A10; 175S6-A20-A19-A17-A11; 111S6-A20-A13; 393S6-A14; 571S6-V11-V12-A15; 165S6-A20-A18; 130S6-A20-A19; 95S6-A20; 260S6-A21; 100（五）模型的评价1、本文从简单的角度入手建立模型，运用枚举法、图上作业法、图论等多种方法对模型进行逐步优化。过程严谨，理论性强，逻辑严密，而且易于理解。2、模型一、三为问题的近似解，但涉及变量少，易于求解，且精度较高；模型二、四为问题的精确解，变量多，求解稍难。3、一般模型（1）（2）式具有很强的通用性，网络、树型图转化为线性图的方法也具有一般性；4、本文大量运用了计算机程序，所有数据均由计算机处理，故误差由计算机精度产生，模型具有较好的稳定性。参考文献（1）李德、钱颂迪 运筹学 清华大学出版社 1982年；（2）许卓群等 数据结构 高等教育出版社 1987年；（3）滕传琳 管理运筹学 中国铁道出版社 1986年；（4）美E 米涅卡 网络和图的最优计算方法 中国道出铁版社 1984年；评注：本题是一个典型的离散优化问题，实际上最终归结为两个二次规划模型。求解的关键包括定购和运输单价的计算、二次规划模型的建立和求解。其中的难点是题中复杂图的恰当处理和模型中变量数特别是整型变量数过多的灵活解决。本优秀论文在对这两问题的处理上颇有特色：其一是对问题一的图，作者首先根据直观分析，将原图分成两个小图分别求解，再合二为一建立优化模型；而针对问题二的树形图，通过增加一个长度为零费用为零的线路转化为线性图后，问题二的求解方法就与问题一的求解方法相同了。这是“从简单到复杂，从局部到全局，化繁琐为简单”的建模思想的具体体现。其二是对变量数过多问题的解决，作者经过约束条件的松弛化处理后，充分利用非线性规划软件Lingo的强大功能，主要是循环语句的灵活使用，得出了最优解。本篇论文获得2000年数学建模竞赛的全国一等奖。1. 基于C8051F单片机直流电动机反馈控制系统的设计与研究2. 基于单片机的嵌入式Web服务器的研究 3. MOTOROLA单片机MC68HC（8）05PV8/A内嵌EEPROM的工艺和制程方法及对良率的影响研究 4. 基于模糊控制的电阻钎焊单片机温度控制系统的研制 5. 基于MCS-51系列单片机的通用控制模块的研究 6. 基于单片机实现的供暖系统最佳启停自校正（STR）调节器7. 单片机控制的二级倒立摆系统的研究8. 基于增强型51系列单片机的TCP/IP协议栈的实现 9. 基于单片机的蓄电池自动监测系统 10. 基于32位嵌入式单片机系统的图像采集与处理技术的研究11. 基于单片机的作物营养诊断专家系统的研究 12. 基于单片机的交流伺服电机运动控制系统研究与开发 13. 基于单片机的泵管内壁硬度测试仪的研制 14. 基于单片机的自动找平控制系统研究 15. 基于C8051F040单片机的嵌入式系统开发 16. 基于单片机的液压动力系统状态监测仪开发 17. 模糊Smith智能控制方法的研究及其单片机实现 18. 一种基于单片机的轴快流CO,2激光器的手持控制面板的研制 19. 基于双单片机冲床数控系统的研究 20. 基于CYGNAL单片机的在线间歇式浊度仪的研制 21. 基于单片机的喷油泵试验台控制器的研制 22. 基于单片机的软起动器的研究和设计 23. 基于单片机控制的高速快走丝电火花线切割机床短循环走丝方式研究 24. 基于单片机的机电产品控制系统开发 25. 基于PIC单片机的智能手机充电器 26. 基于单片机的实时内核设计及其应用研究 27. 基于单片机的远程抄表系统的设计与研究 28. 基于单片机的烟气二氧化硫浓度检测仪的研制 29. 基于微型光谱仪的单片机系统 30. 单片机系统软件构件开发的技术研究 31. 基于单片机的液体点滴速度自动检测仪的研制32. 基于单片机系统的多功能温度测量仪的研制 33. 基于PIC单片机的电能采集终端的设计和应用 34. 基于单片机的光纤光栅解调仪的研制 35. 气压式线性摩擦焊机单片机控制系统的研制 36. 基于单片机的数字磁通门传感器 37. 基于单片机的旋转变压器-数字转换器的研究 38. 基于单片机的光纤Bragg光栅解调系统的研究 39. 单片机控制的便携式多功能乳腺治疗仪的研制 40. 基于C8051F020单片机的多生理信号检测仪 41. 基于单片机的电机运动控制系统设计 42. Pico专用单片机核的可测性设计研究 43. 基于MCS-51单片机的热量计 44. 基于双单片机的智能遥测微型气象站 45. MCS-51单片机构建机器人的实践研究 46. 基于单片机的轮轨力检测 47. 基于单片机的GPS定位仪的研究与实现 48. 基于单片机的电液伺服控制系统 49. 用于单片机系统的MMC卡文件系统研制 50. 基于单片机的时控和计数系统性能优化的研究 51. 基于单片机和CPLD的粗光栅位移测量系统研究 52. 单片机控制的后备式方波UPS 53. 提升高职学生单片机应用能力的探究 54. 基于单片机控制的自动低频减载装置研究 55. 基于单片机控制的水下焊接电源的研究 56. 基于单片机的多通道数据采集系统 57. 基于uPSD3234单片机的氚表面污染测量仪的研制 58. 基于单片机的红外测油仪的研究 59. 96系列单片机仿真器研究与设计 60. 基于单片机的单晶金刚石刀具刃磨设备的数控改造 61. 基于单片机的温度智能控制系统的设计与实现 62. 基于MSP430单片机的电梯门机控制器的研制 63. 基于单片机的气体测漏仪的研究 64. 基于三菱M16C/6N系列单片机的CAN/USB协议转换器 65. 基于单片机和DSP的变压器油色谱在线监测技术研究 66. 基于单片机的膛壁温度报警系统设计 67. 基于AVR单片机的低压无功补偿控制器的设计 68. 基于单片机船舶电力推进电机监测系统 69. 基于单片机网络的振动信号的采集系统 70. 基于单片机的大容量数据存储技术的应用研究 71. 基于单片机的叠图机研究与教学方法实践 72. 基于单片机嵌入式Web服务器技术的研究及实现 73. 基于AT89S52单片机的通用数据采集系统 74. 基于单片机的多道脉冲幅度分析仪研究 75. 机器人旋转电弧传感角焊缝跟踪单片机控制系统 76. 基于单片机的控制系统在PLC虚拟教学实验中的应用研究77. 基于单片机系统的网络通信研究与应用 78. 基于PIC16F877单片机的莫尔斯码自动译码系统设计与研究79. 基于单片机的模糊控制器在工业电阻炉上的应用研究 80. 基于双单片机冲床数控系统的研究与开发 81. 基于Cygnal单片机的C/OS-的研究82. 基于单片机的一体化智能差示扫描量热仪系统研究 83. 基于TCP/IP协议的单片机与Internet互联的研究与实现 84. 变频调速液压电梯单片机控制器的研究 85. 基于单片机-免疫计数器自动换样功能的研究与实现 86. 基于单片机的倒立摆控制系统设计与实现 87. 单片机嵌入式以太网防盗报警系统 88. 基于51单片机的嵌入式Internet系统的设计与实现 89. 单片机监测系统在挤压机上的应用 90. MSP430单片机在智能水表系统上的研究与应用 91. 基于单片机的嵌入式系统中TCP/IP协议栈的实现与应用92. 单片机在高楼恒压供水系统中的应用 93. 基于ATmega16单片机的流量控制器的开发 94. 基于MSP430单片机的远程抄表系统及智能网络水表的设计95. 基于MSP430单片机具有数据存储与回放功能的嵌入式电子血压计的设计 96. 基于单片机的氨分解率检测系统的研究与开发 97. 锅炉的单片机控制系统 98. 基于单片机控制的电磁振动式播种控制系统的设计 99. 基于单片机技术的WDR-01型聚氨酯导热系数测试仪的研制 100. 一种RISC结构8位单片机的设计与实现 101. 基于单片机的公寓用电智能管理系统设计 102. 基于单片机的温度测控系统在温室大棚中的设计与实现103. 基于MSP430单片机的数字化超声电源的研制 104. 基于ADC841单片机的防爆软起动综合控制器的研究105. 基于单片机控制的井下低爆综合保护系统的设计 106. 基于单片机的空调器故障诊断系统的设计研究 107. 单片机实现的寻呼机编码器 108. 单片机实现的鲁棒MRACS及其在液压系统中的应用研究 109. 自适应控制的单片机实现方法及基上隅角瓦斯积聚处理中的应用研究110. 基于单片机的锅炉智能控制器的设计与研究 111. 超精密机床床身隔振的单片机主动控制 112. PIC单片机在空调中的应用 113. 单片机控制力矩加载控制系统的研究 项目论证，项目可行性研究报告，可行性研究报告，项目推广，项目研究报告，项目设计，项目建议书，项目可研报告，本文档支持完整下载，支持任意编辑！选择我们，选择成功！项目论证，项目可行性研究报告，可行性研究报告，项目推广，项目研究报告，项目设计，项目建议书，项目可研报告，本文档支持完整下载，支持任意编辑！选择我们，选择成功！单片机论文，毕业设计，毕业论文，单片机设计，硕士论文，研究生论文，单片机研究论文，单片机设计论文，优秀毕业论文，毕业论文设计，毕业过关论文，毕业设计，毕业设计说明，毕业论文，单片机论文，基于单片机论文，毕业论文终稿，毕业论文初稿，本文档支持完整下载，支持任意编辑！本文档全网独一无二，放心使用，下载这篇文档，定会成功！