运筹学第8讲 运输问题.pptVIP

表上作业法 解 引入虚拟的销地B4,其中从各产地到虚拟销地B4的运 费单价为0元。 根据运输表写出初始基本可行解(最小元素法)。 求解过程如下表所示： 表上作业法 继续迭代： 表中各检验数均非负，故当前解为最优解。 即A1运送10吨到B1；A1运送15吨到B3；A2运送25吨到 B2，A2就地储存10吨。 Lingo软件求解 例8.5 有四个煤产地A1,A2,A3,A4，今年产量分别为35， 45，55，65万吨，另外有六个煤销地B1,B2,B3,B4,B5, B6，今年销量分别为40，20，30，40，30，40。如果 从产地到销地的单位运输成本如下表，问如何运输可使 总运费最小？ B1 B2 B3 B4 B5 B6 A1 5 7 2 4 1 8 A2 9 1 3 5 6 7 A3 2 4 8 1 3 5 A4 7 6 1 2 4 9 设Ai产地的产量为ai，Bj销地的销量为bj，从Ai到Bj的 单位运价为cij，xij为从产地Ai运到Bj的数量。 Lingo软件求解 一般运输问题的数学模型： 定义集合产地（Production）； 定义集合销地（Distribution）； 由产地和销地派生出来的集合（LinkPD）； 在每个集合上可以定义变量（包括已知变量和未知 的决策变量）; 在Lingo程序中必须在集合段内定义集合。 Lingo软件求解 用Lingo软件求解： Model: Sets: Production/A1,A2,A3,A4/: Pvolumn; Distribution/B1,B2,B3,B4,B5,B6/:Dvolumn; LinkPD(Production,Distribution): Cost,Tvolumn; ！Tvolumn为决策变量 Endsets Data： Pvolumn ＝ 35，45，55，65； Dvolumn ＝40，20，30，40，30，40； Cost ＝ 5 7 2 4 1 8 9 1 3 5 6 7 2 4 8 1 3 5 7 6 1 2 4 9； Enddata Lingo软件求解 运输问题Lingo程序： ！对每个产地运出的货物量不超过其产量：@for(Production(i): @sum(Distribution(j):Tvolumn(i,j))=Pvolumn(i)); ！对每个销地运到的货物量应超过其需求量：@for(Distribution(j): @sum(Production(i):Tvolumn(i,j)=Dvolumn(j)); ！目标函数是单位运价乘以运量。 min = @sum(LinkPD: Cost*TVOLUMN); end Lingo软件求解 运输问题Lingo程序： 程序见运输问题1.lg4。 Lingo软件求解 问题的推广： 如果明年各销地的销量都增加5万吨，总销量增加了30万吨。为了产销平衡，打算增加一套年产30万吨的采煤设备，如果这套采煤设备放在A1，A2，A3，A4的生产成本分别为20，30，30，25万元。问题：这套采煤设备应当放在何处，可使总费用最小？ Lingo软件求解 建模分析： 总费用包括生产成本费用和运输成本费用。 设定变量yi：yi=1，表示采煤设备放在Ai，否则yi=0 yi称为0-1变量。 Lingo软件求解 建模分析： 总费用包括生产成本费用和运输成本费用。 设定变量yi：yi=1，表示采煤设备放在Ai，否则yi=0 yi称为0-1变量。 Sets: Production/A1,A2,A3,A4/: Pvolumn,y,Pcost; Distribution/B1,B2,B3,B4,B5,B6/:Dvolumn; LinkPD(Production,Distribution): Cost,Tvolumn; Endsets ！产地集合中，定义一个决策变量y和已知变量Pcost。 ！对每个产地运出的货物量不超过其产量：@for(Production(i): @sum(Distribution(j):Tvolumn(i,j))=Pvolumn(i)＋30*y(i);@bin(y(i))); ！@bin(y)是限定变量y为0－1变量。 ！对每个销地运到的货物量应超过其需求量：@for(Distribution(j): @sum(Production(i):Tvolumn(i,j)=Dvolumn(j)); ！目标函数运输成本加生产成本 Lingo软件求解 推广的运输问题的Lingo程序： 程序见运输问题2.lg4。 运输问题 小结： 运输问题模型的构建； 初始基本可行解的三种确定