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

钢管订购和运输数学建模论文PAGEPAGE１９钢管订购和运输摘要本文建立了一个运输问题的最优化模型。通过对图(一)的分析，我们首先直观地将路线分成两段，将图分为两个子图建立了模型一,利用分支定界法求得总费用最优解为1279496万元。然后对模型一进行优化，得到全线的最优模型二，求得总费用最优解为1278632万元。通过对最优模型二的分析，我们得出钢厂S1的上限产量和钢厂S6的销价的的变化对运购计划和总费用的影响最大，并给出了数据结果。我们利用截取和连接的方法将树形图转化成为对线性图进行分析，并给出了一般的解决方法。对图(二)给出的具体模型，类似与问题一,分别建立了模型三和模型四,求得最优解分别为1408859.4和1403948万元．

问题的提出已知有7个钢厂，可生产输送天然气主管道的钢管，用Si表示(i=1,2,…,7)。现有１５个地点（A1,A2,…,A15），沿着这１５个地点铺设一条输送天然气的主管道。为方便计，1km主管道称为１单位钢管。一个钢厂如果承担制造这种钢管，至少需要生产500个单位。钢厂在指定期限内能生产该钢管的最大数量为个单位，钢管出厂销价1单位钢管为万元，如下表：1234567800800100020002000200030001601551551601551501601单位钢管的铁路运价如下表：里程(km)≤300301～350351～400401～450451～500运价(万元)2023262932里程(km)501～600601～700701～800801～900901～1000运价(万元)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,7ni:钢厂Si的产量上限，ni=si，i=1,2,…，7bi:图一中各节点（见附录七）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.3w2=419866.7模型一的总费用为w=w1+w2=1279496（万元）。模型（二）：由于前面的模型直观将管道分为两端，即原图分为两个子图，可能存在一定的误差。下面对模型（一）进行优化，即对主管道建立统一模型，如下：同样利用结论1、2并采用分支定界法，利用LINGO软件编程（见附录三）运行得出:最优解W=1278632万元，各节点向右铺设的管道单位数：Y1=0Y2=175Y3=282Y4=0Y5=10Y6=16Y7=76Y8=175Y9=159Y10=30Y11=145Y12=11Y13=34Y14=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;300S3-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模型（三）首先类似于模型一，按分段的思想建立模型。由A1…A9、A16为一段可得如下模型：由A21-A20-A19-A17-A11-A12-A13-A14-A15-A9-A10-A11-A17-A18为一段可得如下模型：运用LINGO软件进行编程（见附录四、五）运行结果为：费用w1=８69127.2w2=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;300S3- 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年数学建模竞赛的全国一等奖。基于C8051F单片机直流电动机反馈控制系统的设计与研究基于单片机的嵌入式Web服务器的研究MOTOROLA单片机MC68HC（8）05PV8/A内嵌EEPROM的工艺和制程方法及对良率的影响研究基于模糊控制的电阻钎焊单片机温度控制系统的研制基于MCS-51系列单片机的通用控制模块的研究基于单片机实现的供暖系统最佳启停自校正（STR）调节器单片机控制的二级倒立摆系统的研究基于增强型51系列单片机的TCP/IP协议栈的实现基于单片机的蓄电池自动监测系统基于32位嵌入式单片机系统的图像采集与处理技术的研究基于单片机的作物营养诊断专家系统的研究基于单片机的交流伺服电机运动控制系统研究与开发基于单片机的泵管内壁硬度测试仪的研制基于单片机的自动找平控制系统研究基于C8051F040单片机的嵌入式系统开发基于单片机的液压动力系统状态监测仪开发模糊Smith智能控制方法的研究及其单片机实现一种基于单片机的轴快流CO〈,2〉激光器的手持控制面板的研制基于双单片机冲床数控系统的研究基于CYGNAL单片机的在线间歇式浊度仪的研制基于单片机的喷油泵试验台控制器的研制基于单片机的软起动器的研究和设计基于单片机控制的高速快走丝电火花线切割机床短循环走丝方式研究基于单片机的机电产品控制系统开发基于PIC单片机的智能手机充电器基于单片机的实时内核设计及其应用研究基于单片机的远程抄表系统的设计与研究基于单片机的烟气二氧化硫浓度检测仪的研制基于微型光谱仪的单片机系统单片机系统软件构件开发的技术研究基于单片机的液体点滴速度自动检测仪的研制基于单片机系统的多功能温度测量仪的研制基于PIC单片机的电能采集终端的设计和应用基于单片机的光纤光栅解调仪的研制气压式线性摩擦焊机单片机控制系统的研制基于单片机的数字磁通门传感器基于单片机的旋转变压器-数字转换器的研究基于单片机的光纤Bragg光栅解调系统的研究单片机控制的便携式多功能乳腺治疗仪的研制基于C8051F020单片机的多生理信号检测仪基于单片机的电机运动控制系统设计Pico专用单片机核的可测性设计研究基于MCS-51单片机的热量计基于双单片机的智能遥测微型气象站MCS-51单片机构建机器人的实践研究基于单片机的轮轨力检测基于单片机的GPS定位仪的研究与实现基于单片机的电液伺服控制系统用于单片机系统的MMC卡文件系统研制基于单片机的时控和计数系统性能优化的研究基于单片机和CPLD的粗光栅位移测量系统研究单片机控制的后备式方波UPS提升高职学生单片机应用能力的探究基于单片机控制的自动低频减载装置研究基于单片机控制的水下焊接电源的研究基于单片机的多通道数据采集系统基于uPSD3234单片机的氚表面污染测量仪的研制基于单片机的红外测油仪的研究96系列单片机仿真器研究与设计基于单片机的单晶金刚石刀具刃磨设备的数控改造基于单片机的温度智能控制系统的设计与实现基于MSP430单片机的电梯门机控制器的研制基于单片机的气体测漏仪的研究基于三菱M16C/6N系列单片机的CAN/USB协议转换器基于单片机和DSP的变压器油色谱在线监测技术研究基于单片机的膛壁温度报警系统设计基于AVR单片机的低压无功补偿控制器的设计基于单片机船舶电力推进电机监测系统基于单片机网络的振动信号的采集系统基于单片机的大容量数据存储技术的应用研究基于单片机的叠图机研究与教学方法实践基于单片机嵌入式Web服务器技术的研究及实现基于AT89S52单片机的通用数据采集系统基于单片机的多道脉冲幅度分析仪研究机器人旋转电弧传感角焊缝跟踪单片机控制系统基于单片机的控制系统在PLC虚拟教学实验中的应用研究基于单片机系统的网络通信研究与应用基于PIC16F877单片机的莫尔斯码自动译码系统设计与研究基于单片机的模糊控制器在工业电阻炉上的应用研究基于双单片机冲床数控系统的研究与开发基于Cygnal单片机的μC/OS-Ⅱ的研究基于单片机的一体化智能差示扫描量热仪系统研究基于TCP/IP协