数学建模——网络优化 钢管运输.pptVIP

优化方法与建模 钢管订购与运输问题 问题： 要铺设一条 的输送天然气的主管道, 如图一所示(见下页)。经筛选后可以生产这种主管道钢管的 钢厂有 . 图中粗线表示铁路，单细线表示公路，双细线表示要铺设的管道(假设沿管道或者原来有公路，或者建有施工公路)，圆圈表示火车站，每段铁路、公路和管道旁的阿拉伯数字表示里程(单位km)。 为方便计，1km主管道钢管称为1单位钢管。 一个钢厂如果承担制造这种钢管，至少需要生产500个单位。钢厂 在指定期限内能生产该钢管的最大数量为 个单位，钢管出厂销价1单位钢管为 万元，如下表： 1单位钢管的铁路运价如下表： （1）请制定一个主管道钢管的订购和运输计划，使总费用最小(给出总费用)。 （2）请就（1）的模型分析：哪个钢厂钢管的销价的变化对购运计划和总费用影响最大，哪个钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大，并给出相应的数字结果。 （3）如果要铺设的管道不是一条线，而是一个树形图，铁路、公路和管道构成网络，请就这种更一般的情形给出一种解决办法，并对图二按（1）的要求给出模型和结果。 1.购买、运输钢管 都是整单位。 2. 沿铺设主管道已有公路或者有施工公路。 3.钢厂先将钢管运输到各个结点Aj,再由Aj向各个方向运输。 4. 在主管道上，每公里卸1单位的钢管。 5.在求解钢厂的价格对总价的影响时，认为钢管的单价只会在一个小范围内变化，在求解钢厂的生产上限对总价的影响时，亦是如此。 　　用Lingo软件编程求解二次规划问题，得出如下结果： 当产量限制条件为： 时，求解出最小花费为127.53亿， S4的产量为0，但此时，S7的产量为245,不符合大于500的条件，故我们在二次规划问题中再加上两个限制条件再次求解。 额外限制一： , 再次求解得出最小花费为: 127.86亿. 额外限制二： , 再次求解得出最小花费为: 127.97亿. 故：问题1的最小运费为fm=127.86 亿，此时S1到S7的产量＝（800，800，1000，0，1366，1205，0） xij 是变量, 当xij =1时，表示弧（i，j）位于s-t路上，当xij =0时，表示弧（i，j）不在s-t路上. wij 表示弧（i，j）上的权. 4.3 最短路问题的数学描述 关联矩阵是全么模矩阵，因此0-1变量可以松弛为区间[0,1]中的实数 不含负圈，变量直接松弛为所有非负实数 思考：为什么xij 可以不限定为{0，1}？ Bellman方程(最短路方程) 一般情况下直接求解最短路方程是相当困难的. 定理5.1 对于只含正有向圈的连通有向网络，从起点s到任一顶点j都存在最短路，它们构成以起点s为根的树形图（称为最短路树（Tree of Shortest Paths）或最短路树形图（Shortest Path Arborescence）长度可以由Bellman方程唯一确定. 最短路树（树形图） A F 5 6 6 7 4 4 5 1 3 B E D C 4.4最短路算法 算法通过不断修改这些标号，进行迭代计算. 当算法结束时，距离标号表示的是从起点到该顶点的最短路长度. (1) Dijkstra算法 (1959, 正费用网络） 基本思想： 对于V 中每一个顶点j，赋予两个数值（通常称为“标号”）： 1.距离标号u(j):记录的是从起点到该顶点的最短路长度的上界； 2.前趋标号p（j）:记录的是当起点s到该顶点j 的一条路长 取到该上界时，该条路中顶点j 前面的那个节点. Dijkstra算法 STEP1. 如果S=V, 则 为节点s到节点j的最短路路长(最短路可以通过数组p所记录的信息反向追踪获得), 结束. STEP0. (初始化) 令S= ，=V， ；对V 中的顶点j(j s)令初始距离标号 . STEP2. 从 中找到距离标号最小的节点i，把它从 删除，加入S. 对于所有从i出发的弧 , 若 ，则令