1、2019/1/3,1,运输问题(Transportation Problem),运输规划的数学模型 表上作业法 产销不平衡的运输问题 有转运的运输问题,2019/1/3,2,一、运输问题的数学模型,例：某运输问题的资料如下：,2019/1/3,3,2019/1/3,4,运输问题数学模型的一般形式,已知某产品有m个生产地点Ai（i=1,2,，m），其供应量(产量)分别为ai（i=1，2，m），有n个销地Bj（j=1,2,n），其需要量分别为bj（j=1,2,n），从Ai到Bj运输单位物资的运价(单价)为cij，这些数据可汇总于产销平衡表和单位运价表中，见表4-1，表4-2。有时可把这两表合二为一
2、。,2019/1/3,5,表4-1 产销平衡表,2019/1/3,6,表4-2 单位运价表,2019/1/3,7,运输问题数学模型的一般形式,若用xij表示从Ai到Bj的运量，那么在产销平衡的条件下，要求得总运费最小的调运方案，其数学模型为,2019/1/3,8,运输问题数学模型的特点,1.运输问题有有限个最优解 2.运输问题约束条件的特点：运输问题的数学模型包含mn个变量，(m+n)个约束方程，其系数矩阵的结构比较松散且特殊。,2019/1/3,9,运输问题数学模型的特点,该系数矩阵中对应于某一变量的系数向量， 其分量中除第i个和第m+j个为1以外，其余的都为零， 在约束条件中，前m个约束条
3、件的含义是：由某一产地运往各销地的产品数量之和等于该产地的产量；后n个约束条件的含义是：由各产地运往某一销地的产品数量之和等于该销地的销量。,2019/1/3,10,运输问题的对偶问题,运输问题的对偶问题可按照前面写线性规划问题的对偶问题的方法给出（略）,2019/1/3,11,运输问题的解,运输问题也是一个线性规划问题，其求解时仍然可以先找一个基可行解，进行解的最优性检验，若不是最优，就进行迭代，继续检验和调整直到最优，因此要求每步得到的解都是基可行解，需满足(1)满足所有约束条件；(2)基变量对应的系数列向量线性无关；(3)解中非零变量的个数不能大于m+n-1个，原因是运输问题中虽有m+n
4、个约束条件，但由于总产量等于总销量，故只有m+n-1个约束条件是线性独立的；(4)保持基变量的个数在迭代过程中为m+n-1个。,2019/1/3,12,第2节 表上作业法,表上作业法是单纯形法在求解运输问题时的一种简化方法，其实质是单纯形法,但具体计算和术语有所不同，其步骤可归纳为： (1) 找出初始基可行解：即在(mn)产销平衡表上用西北角法或最小元素法、Vogel法给出m+n-1个数字，称为数字格，它们就是初始基变量的取值。 (2) 求各非基变量的检验数，即在表上计算空格的检验数，判别是否达到最优解。如已是最优解，则停止计算，否则转到下一步。 (3) 确定换入变量和换出变量，找出新的基可行
5、解。在表上用闭回路法调整。 (4) 重复(2)，(3)直到得到最优解为止。,2019/1/3,13,例题,例：某食品公司下属的A1、A2、A3 ，3个厂生产方便食品，要运输到B1、B2、B3、B4 ，4个销售点，数据如下： 求最优运输方案。,2019/1/3,14,一 初始基可行解的确定,确定初始基可行解的方法很多，有西北角法，最小元素法和沃格尔(Vogel)法等，一般希望的方法是既简便，又尽可能接近最优解。下面分别予以介绍： 1. 最小元素法此方法的基本思想是就近供应，即从单位运价表中最小的运价开始确定供销关系，然后次小，一直到给出初始基可行解为止。以上例进行讨论得下表：,2019/1/3,
6、15,3,11,3,10,1,9,2,7,4,10,5,8,3,4,1,6,3,3,总的运输费用（31）（64）（43）（12）（310）（35）86元,2019/1/3,16,2.西北角法（或左上角法）：此法是纯粹的人为规定，没有理论依据和实际背景，但它易操作，特别适合在计算机上编程计算，因而受欢迎。方法如下：,3 6 5 6,7 4 9,3,4 4 9,0 6 5 6,4,0 4 9,0 2 5 6,2,0 2 9,0 0 5 6,2,0 0 9,0 0 3 6,3 6,0 0 0 0,0 0 0,3 4 0 0 0 2 2 0 0 0 3 6,总的运费(33)(411)(29)(22)(
7、310)(65)135元,2019/1/3,17,3. 沃格尔法,最小元素法的缺点是：为了节省一处的费用，有时造成在其他处要多花几倍的运费。沃格尔法考虑到，一产地的产品假如不能按最小运费就近供应，就考虑次小运费，这就有一个差额。差额越大，说明不能按最小运费调运时，运费增加越多。因而对差额最大处，就应当采用最小运费调运。 沃格尔法的步骤是：第一步：在原表中分别计算出各行和各列的最小运费和次小运费的差额，并填入该表的最右列和最下行，见表4-1,2019/1/3,18,表4-1,第二步：从行或列差额中选出最大者，选择它所在行或列中的最小元素。在表4-1中B2列是最大差额所在列。B2列中最小元素为4，
8、可确定A3的产品先供应B2的需要。得表4-2：,2019/1/3,19,同时将运价表中的B2列数字划去。如表4-3所示（表4-3）：,表4-2,2019/1/3,20,第三步：对表4-3中未划去的元素再分别计算出各行、各列的最小运费和次小运费的差额，并填入该表的最右列和最下行，重复第一、二步，直到给出初始解为止。用此法给出上例的初始解列于表4-4（下表）。,2019/1/3,21,说明,由以上可见：沃格尔法同最小元素法除在确定供求关系的原则上不同外，其余步骤相同。沃格尔法给出的初始解比用最小元素法给出的初始解更接近最优解。本例用沃格尔法给出的初始解就是最优解。,2019/1/3,22,二 最优
9、解的判别,最优性检验就是检查所得到的方案是不是最优方案。检查的方法与单纯形方法中的原理相同，即计算检验数。由于目标要求极小，因此，当所有的检验数都大于或等于零时该调运方案就是最优方案；否则就不是最优，需要进行调整。下面介绍两种求检验数的方法：1.闭回路法；2.位势法,2019/1/3,23,1 闭回路法,在给出调运方案的计算表上，如上例表，从每一空格出发找一条闭回路。它是以某空格为起点，用水平或垂直线向前划，当碰到一数字格时可以转90后，继续前进，直到回到起始空格为止。闭回路如图(a),(b),(c)等所示。,2019/1/3,24,闭回路法计算检验数的经济解释为：在已给出初始解的表中， 可从
10、任一空格出发，如(A1，B1)，若让A1的产品调运1吨给B1。为了保持产销平衡，就要依次作调整： 在(A1，B3)处减少1吨，(A2，B3)处增加1吨，(A2，B1)处减少1吨，即构成了以(A1，B1)空格为起点，其他为数字格的闭回路，如下表虚线所示。在这表中闭回路各顶点所在格的右上角数字是单位运价,2019/1/3,25,可见这调整的方案使运费增加 (+1)3+(-1)3+(+1)2+(-1)=1(元) 将“1”这个数填入(A1，B1)格，这就是检验数。按以上所述，可找出所有空格的检验数，见下表。,2019/1/3,26,0,0,0,0,0,1,2,1,-1,12,10,0,当检验数还存在负
11、数时，说明原方案不是最优解，要继续改进。 闭回路法的主要缺点是：当变量个数较多时，寻找闭回路以及计算两方面都会产生困难。,运输问题的约束条件共有m+n个，其中：m是产地产量的限制；n是销地销量的限制。其对偶问题也应有m+n个变量，据此： ij=cij(ui+vj) ,其中前m个计为ui（i=1.2m）,前n个计为vj （j=1.2n）由单纯形法可知，基变量的ij 0 cij(ui+vj) 0 因此ui ,vj可以求出。,.位势法,2019/1/3,28,仍以上面的例子说明。 第一步：按最小元素法给出表的初始解，然后做下表：即在对应表的数字格处填入单位运价：,2019/1/3,29,第二步：在上
12、表中增加一行一列，在列中填入ui，在行中填入vj，得下表：,先令u1=0，然后按ui+vj=cij， i，jB相继地确定ui，vj。由上表可见，当u1=0时，由u1+v3=3可得v3=3，由u1+v4=10可得v4=10；在v4=10时，由u3+v4=5可得u3=-5，以此类推可确定所有的ui，vj的数值。,2019/1/3,30,第三步：按ij=cij-(ui+vj), i,jN计算所有空格的检验数。如 11=c11-(u1+v1)=3-(0+2)=1， 12=c12-(u1+v2)=11-(0+9)=2 这些计算可直接在上表中进行。 为了方便，特设计计算表，如下表所示,在此表中还有负检验数
13、，说明未得最优解，还可以改进。,2019/1/3,31,亦即：,成本表,u2+v1=1 u2+ v3 =2 u3+v2=4 u1+ v4 =10 u1+v3=3 u3+ v4 =5 令： u10,u10 v12 u2 1 v2 9 u3 5 v3 3v4 10,(ui+vj),2019/1/3,32,按ij=cij(ui+vj) 计算检验数，并以ij0 检验，或用(ui+vj) cij 0 检验。,cij,(ui+vj),表中还有负数，说明还未得到最优解，应继续调整。,ij,2019/1/3,33,三 解的改进 闭回路调整法,当非基变量的检验数出现负值时，则表明当前的基本可行解不是最优解。在这
14、种情况下，应该对基本可行解进行调整，即找到一个新的基本可行解使目标函数值下降，具体步骤为： （1）选负检验数中最小者 rk，那么 xrk 为主元，作为进基变量（上页图中 x24 ）; （2）以 xrk 为起点找一条闭回路，除 xrk 外其余顶点必须为基变量格（上页图中的回路）;,2019/1/3,34,解的改进 闭回路调整法,（3）为闭回路的每一个顶点标号， xrk 为 1，沿一个方向（顺时针或逆时针）依次给各顶点标号； （4）求=Minxijxij对应闭回路上的偶数标号格= xpq，那么确定xpq为出基变量，为调整量； （5）对闭回路的各奇标号顶点调整为：xij + ，对各偶标号顶点 调整为
15、：xij - ，特别 xpq - = 0， xpq变为非基变量。重复(2)、(3)步，直到所有检验数均非负，得到最优解。,接上例：,3,1,3,4,6,3,（1）,（1）,（1）,（1）,2019/1/3,36,2019/1/3,37,经检验,所有ij0 得到最优解， 最小运费为85元。,.无穷多最优解：产销平衡的运输问题必定存最优解。如果非基变量的ij0，则该问题有无穷多最优解。如上例：（1.1）中的检验数是 0，经过调整，可得到另一个最优解：可在最优表中以(1，1)为调入格，作闭回路(1，1)+-(1，4)-(2，4)+-(2，1)-(1，1)+ 确定=min(2,3)=2。经调整后得到另
16、一最优解，见下表：,四、需要说明的几个问题,2019/1/3,39,2019/1/3,40,退化：表格中一般要有(m+n-1)个数字格。但有时，在分配运量时则需要同时划去一行和一列，这时需要补一个0，以保证有(m+n-1)个数字格。一般可在划去的行和列的任意空格处加一个 0 即可。在用闭回路法调整时，在闭回路上出现两个和两个以上的具有(-1)标记的相等的最小值。这时只能选择其中一个作为调入格。而经调整后，得到退化解。这时另一个数字格必须填入一个0，表明它是基变量。当出现退化解后，并作改进调整时，可能在某闭回路上有标记为(-1)的取值为0的数字格，这时应取调整量=0。,2019/1/3,41,第
17、三节 产销不平衡的运输问题,前面所讲表上作业法，都是以产销平衡为前提条件的，但是实际问题中产销往往是不平衡的，这就需要设法把产销不平衡的问题化成产销平衡的问题。 一、总产量大于总销量 即此时，运输问题的数学模型可写为,2019/1/3,42,要求解此问题，需要先将原问题变成平衡问题，可以假设一个销地Bn+1(即实际问题中考虑产地的存量)，即：,2019/1/3,43,单位运价表中的单位运价 二、销大于产：同样假设一个产地即可，变化同上。,40,30,30,20,30,20,20,用最小元素法求初始方案,例题：,然后用表上作业法求解即可（略）,已知某运输问题的资料如下表所示,表中的发量、收量单位
18、为：吨，运价单位为：元/吨试求出最优运输方案.,练习：,2019/1/3,46,解：用最小元素法求初始方案,运费为108元/吨,2、用位势法判断：,成本表,2019/1/3,47,u1+v3=5 u2+v4 =1u1+v4 =3 u3+v2=2u2+v1=1 u3+v4 =4 令： u10,u10 v13 u22 v2 1 u3 1 v3 5v4 3,2019/1/3,48,(ui+vj),2019/1/3,49,cij,表中还有负数，说明没有得到最优解，调整运输方案。,ij,(ui+vj),2019/1/3,50,2,2,2,2,新的运送方案,新的成本表,(ui+vj)1,总的运费 105元/吨,2019/1/3,51,表中还有负数，说明没有得到最优解，继续调整运输方案。,cij,(ui+vj)1,(ij)1,2019/1/3,52,2019/1/3,53,cij,(ui+vj)2,(ij)2,表中没有负数，说明已经得到最优解。但有无穷多最优解。,0,2019/1/3,54,
