运筹学 第三章 运输问题

4、重复第二、第三步，直至得到最优解。
2020/4/3
10
一、确定初始基本可行解：
对于有m个产地n个销地的产销平衡问题，有m个关于产量 的约束方程和n个关于销量的约束方程。表面上，共有m+n个 约束方程。
但由于产销平衡，其模型最多只有m+n-1个独立的约束方 程，所以运输问题实际上有m+n-1个基变量。在m×n的产销 平衡表上给出m+n-1个数字格，其相对应的调运量的值即为 基变量的值。
3
B4
ui
3 10
0
-1 8
-1
35
-5
10
B1
3
31
7
2
B2
11 9
64
9
B3
4(+1) 3 1 (-1) 2
10
3
B4
ui
3(-1) 10
0
+1 8
-1
35
-5
10
2020/4/3
27
调整运量后的新方案：
销地
产地
B1
A1
A2
3
A3
B2
B3
5
6
销量
3
6
5
B4
产量
2
7
1
4
3
9
6
对上表用位势法进行检验如下表，可知已达最优解。
2020/4/3
8
§2 运输问题的表上作业法
从第一节的运输问题的数学模型可知，运输问题实际上 也属于线性规划，但由于运输问题的特殊性（变量个数较多， 系数矩阵的特点），如果用单纯形表格方法迭代，计算量很 大。今天介绍的 “表上作业法”，是针对运输问题的特殊求解 方法，实质还是单纯形法，但减少了计算量。
B4
产量
3
7
10
30
4 10
8
3
9
5
30
6
20
3
0
20
表中填有数字的格对应于基变量(取值即为格中数字），而空格对应 的是非基变量（取值为零）.
在求初始基本可行解时要注意的一个问题： 当我们取定xij的值之后，会出现Ai的产量与Bj的销量都改为零的情 况，这时只能划去Ai行或Bj列，但不能同时划去Ai行与Bj列。
销地 运费单价
B1
B2
B3
B4
产量 （吨）
产地
A1
3
11
3
10
7
A2
1
9
2
8
4
A3
7
4
10
5
9
销量（吨）
3
6
5
6
2020/4/3
16
销地 产地
A1
A2
A3
2 2 2
B1
3
31 76
5
B2
B3
B4
11
9
4
1 1 1 1
53
21
10
3 3 2 2
2 10 0 0 0 7
8 1 116
3 5 12
20
20
否则产销不平衡。
2020/4/3
7
说明：从上述模型可以看出：
（1）这是一个线性规划的模型； （2）变量有m×n个； （3）约束条件有 m+n 个； （4）系数矩阵非常稀疏；系数矩阵的秩一般为（m+n-1),
而非m+n 。
若直接用单纯形法求解，显然单纯形表比较庞大，于是在 单纯形法的基础上创建了表上作业法求解运输问题这一特 殊的线性规划问题
销量
3
B2
11 9 4
6
B3
B4
4 （-）3 3
10
1 （+）2
8
10 3
5
5
6
产量 7 4 9
销地
产地
B1
A1
（+）1 3
A2
（-） 3 1
A3
76
销量
3
B2
11 9 4
6
B3
B4
4 （-）3 3
10
1 （+）2
8
10 3
5
5
6
产量 7 4 9
从非基变量 x11出发，找到一个闭回路如上表所示。回路有四个 顶点，除 外x1，1 其余都为基变量。 调整调运量：x11 ，1运费增加了3元； x，13 运1费减少3元
x21 9x22 2x23 8x24
约束条件：
7x31 4x32 10x33 5x34 x11 x12 x13 x14 7
产量约束
x21 x22 x23 x24 4
x31 x32 x33 x34 9
x11 x21 x31 3
销量约束
x12 x22 x32 6 x13 x23 x33 5
则该运输问题的模型如下：
2020/4/3
6
mn
Min f
cij xij
i 1 j 1
m
s.t
xij d j
j 1,...,n
i 1
n
xij si j 1
i 1,...m
xij 0, i 1,...m,
说明：当
m
n
si d j
i1
j 1
j 1,...,n
时，称其为产销平衡的运输问题，
，运x23费增1 加2元； ，运费x减21 少11元
调整后，总运费增加：3-3+2-1=1元。
说明如果让 x1为1 基变量，运费就会增加，其增加值1作为
检验数，
的x11
2020/4/3
20
闭回路法计算检验数：就是对于代表非基变量的空格
（其调运量为零），把它的调运量调整为1，由于产销平衡的 要求,必须对这个空格的闭回路中的各顶点的调运量加上或减 少1。最后计算出由这些变化给整个运输方案的总运输费带来 的变化。以这个变化的数值，作为各空格（非基变量）的检验 数。
x14 x24 x34 6
xij 0,i 1,2,3; j 1,2,3,4
2020/4/3
20 20
5
二、一般运输问题数学模型
设有m个产地，分别为 A1 , A2 ,.... Am ;
n 个销地，分别是 B1, B2 ,.... Bn ；
从产地 Ai运往销地 Bj 的单位运价是 cij ，运量 xij si 是产地Ai 的产量；d j 是销地Bj 的销量。
销地 运费单价
B1
产地
A1
3 x11
A2
1 x21
A3
7 x31
销量（吨）
3
B2
11 x12 9 x22 4 x32
6
B3
3 x13 2 x23 10 x33
5
B4
10 x14 8 x24 5 x34
6
产量 （吨）
7 4 9
于是可建立如下的数学模型：
2020/4/3
4
目标函数：MinZ 3x11 11x12 3x13 10x14
3
B4
ui
3 10
0
-1 8
-1
35
-5
10
销地 产地
A1 A2 A3
vj
B1
13 31 10 7
2
B2
2 11 19 64
9
B3
43 12 12 10
3
B4
ui
3 10
0
-1 8
-1
35
-5
10
我们先给u1赋个任意数值，不妨设u1=0，则从基 变量x11的检验数求得 v3=c13-u1=3-0=3 。
其中 Pij (0,....1,0....1, ,0...0)T
第i个分量
第m+j个分 量
又因为基变量的检验数为0，于是由（m+n-1)个基
变 量的检验数 cij ui v j 0
可解出 (u1,...,um , v1，,...进vn而)计算其他非基变量的检
验数。
2020/4/3
25
三、改进运输方案的办法——闭回路调整法
2020/4/3
17
二、最优解的判别
判别解的最优性需要：计算检验数。方法有两种
1.闭回路法
闭回路：是在已给出的调运方案的运输表上从一个代表 非基变量的空格出发，沿水平或垂直方向前进，遇到代表基 变量的填入数字的格可转90度（当然也可以不改变方向） 继续前进，这样继续下去，直至回到出发的那个空格，由此 形成的封闭折线叫做闭回路。一个空格存在唯一的闭回路。
同理可以求得 v4=10，u2= -1，等等见上表。
检验数的求法，即用公式 ij cij u，i v j
如 11 c11 u1 v1 3 0 2 1
。
2020/4/3
24
位势法计算检验数：
检验数： ij cij CB B1Pij
cij YPij cij (u1,...,um , v1,...vn )Pij
（或者在同时划去Ai行与Bj列时，在该行或该列的任意空格处填加一 个0。）
这样可以保证填过数或零的格为m+n-1个，即保证基变量的个数为 m+n-1个。
2020/4/3
15
2.Vogel法
Vogel法的思想是:一地的产品如果不能按照最小运
费就近供应，就考虑次小运费，这就有差额，差额越大， 说明不能按最小运费调运时，运费增加得越多。因而差 额越大处，就应当采用最小运费调运。
销地 产地
A1 A2 A3
B1
03 31 97
B2
2 11 29 64
B3
53 12 12 10
vj