system: [SYSTEM]

user: 2. 主要问题讨论 2.1. 采用位势法计算检验数 产销平衡的运输问题可以通过一种特殊的单纯形法——表上作业法进行求解。当非基变量的检验数 σ i j ≥ 0 时，可确定当前基可行解为最优解。非基变量的检验数 σ i j 可以通过闭回路法或者位势法进行计算。在课程的讲解中，学生经常提出的一个问题是：用位势法计算检验数时，位势是不唯一的，为什么可以用此来计算检验数呢？本文以一个例子来说明：尽管位势不唯一，但是位势法计算得到的检验数和闭回路法计算得到的检验数一定是相同的。 2.1.1. 闭回路与非基变量的检验数 从每一个空格(非基变量)处 x i j = 0 出发一定存在且能找到唯一的闭回路。则图1中非基变量 x i j 的检验数为 图1. 闭回路示意 σ i j = c i j − ( c i k − c l k + c l s − c t s + c t j ) . (2) 2.1.2. 位势法计算检验数 满足方程 u i + v j = c i j ，   ( i = 1 ， 2 ， ⋯ ， m ; j = 1 ， 2 ， ⋯ ， n ) (3) 的解 u i ，   v j 称为基本可行解对应的位势。其中 c i j 为基变量 x i j 对应的单位运价。考虑到产销平衡的运输问题基变量为 m + n − 1 个，故位势方程中包含 m + n − 1 个方程， m + n 个变量，因此有1个变量为自由未知量，故位势不是唯一的。但是利用 σ i j = c i j − ( u i + v j ) ， (4) 计算得到的位势与前面闭回路法得到的位势却是一致的。因为 σ i j = c i j − ( u i + v j ) = c i j − ( u i + v k − v k − u l + u l + v s − v s − u u + u u + v j ) = c i j − [ ( u i + v k ) − ( u l + v k ) + ( u l + v s ) − ( u t + v s ) + ( u t + v j ) ] = c i j − ( c i k − c l k + c l s − c t s + c t j ) . (5) 2.2. 极大化问题与表上作业法 在实际的运输问题中，目标函数可以为求利润最大或营业额最大等问题。其数学模型为 max f = ∑ i = 1 m ∑ j = 1 n c i j x i j ， (6) s .t . { ∑ j = 1 n x i j = a i ， i = 1 ， 2 ， ⋯ ， m . ∑ i = 1 m x i j = b j ， j = 1 ， 2 ， ⋯ ， n . x i j ≥ 0 ， i = 1 ， 2 ， ⋯ ， m ， j = 1 ， 2 ， ⋯ ， n . 实际求解时，将极大化问题转化为极小化问题。设常数M满足 M ≥ max c i j 将目标函数更改为 min z = ∑ i = 1 m ∑ j = 1 n ( M − c i j ) x i j = ∑ i = 1 m ∑ j = 1 n ( − c i j ) x i j + M ∑ i = 1 m a i . (7) 将 M − c i j ≥ 0 作为极小化问题的单位运价，用原有的表上作业法求出最优解。 此外，亦可以直接采用表上作业法对极大化问题进行求解，具体的求解步骤为： 1) 初始化：采用最大元素法或其他方法给出初始可行解； 2) 判断当前可行解是否是最优解；采用闭回路法或是位势法对表上空格处的检验数，当所有空格处的检验数 σ i j = c i j − C B B − 1 P i j ≤ 0 时，当前基可行解为最优解。 3) 更新基变量与基解。 a) 选出 σ r q > 0 中的最大值 σ i j = max r ， q { σ r q | σ r q > 0 } . 则 x i j 为入基变量。 b) 在运输表上找一条包含 x i j 为顶点的闭回路，这条闭回路上的其他顶点皆为基变量(格)。 c) 计算 θ = min { x r q | x r q 对 应 闭 回 路 上 的 偶 点 格 } = x l s 称 为 θ 的 调 整 量 . 确定 x l s 为出基变量，所在格子调整为空格。 在闭回路上各个奇点变量值均加上 θ ，各个偶点格均减去 θ ，而不在闭回路上的各个变量值不变，这样便完成了一次迭代，并保证目标函数上升。 以图1中的闭回路为例，迭代前目标函数 z old = ∑ ∑ x r q 为 基 变 量 且 不 在 闭 回 路 上 c r q x r q + [ c i k x i k + c l k x l k + c l s x l s + c u s x u s + c u j x u j ] . 迭代后的目标函数 z new = ∑ ∑ x r q 为 基 变 量 且 不 在 闭 回 路 上 c r q x r q + [ c i k ( x i k − θ ) + c l k ( x l k − θ ) + c l s ( x l s − θ ) + c u s ( x u s + θ ) + c u j ( x u j − θ ) + c i j θ ] = ∑ ∑ x r q 为 基 变 量 且 不 在 闭 回 路 上 c r q x r q + [ c i k x i k + c l k x l k + c l s x l s + c u s x u s + c u j x u j ] + ( c i j − c i k + c l k − c l s + c t s − c t j ) θ = z old + σ i j θ ≥ z old . 2.3. 表上作业法与Karush-Kuhn-Tucker条件 本文中的运输问题本质上为一个带约束的优化问题，而针对优化问题，学生在“线性规划”，“最优化问题” [ 3 ] [ 4 ] 等课程中均已接触过相关问题的求解方法，也明白：Karush-Kuhn-Tucker条件是带不等式约束的非线性优化问题取得最优解的必要条件 [ 5 ] 。如果考虑用此方法对产销平衡的运输问题进行求解，会得到什么样的结果？我们用以下定理说明Karush-Kuhn-Tucker乘子与运输问题中位势之间的关系。 定理2.1设 L ( X ; λ ， μ ， d ) = ∑ i = 1 m ∑ j = 1 n c i j x i j + ∑ i = 1 m λ i ( ∑ j = 1 n x i j − a i ) + ∑ j = 1 n μ j ( ∑ i = 1 m x i j − b j ) + ∑ i = 1 m ∑ j = 1 n d i j ( − x i j ) . (8) 为对模型1.1中的最优化问题对应的Lagrange乘子函数，则运输问题(1)中 1) 等式约束对应的Lagrange乘子为位势的相反数： λ i = − u i ， μ j = − v j 。 2) 非负性约束对应的Lagrange乘子为非基变量的检验数： d i j = c i j + λ i + μ j = c i j − ( u i + v j ) . 证明：若 x i j * ( i = 1 ， 2 ， ⋯ ， m ， j = 1 ， 2 ， ⋯ ， n ) 为运输问题1.1的最优解，则 x i j * 满足对应的Karush-Kuhn-Tucker 条件。 平稳性条件： ∂ L ∂ x i j = c i j + λ i + μ j − d i j = 0 ， 原始可行性： ∂ L ∂ μ j = ∑ i = 1 m x i j * − b j = 0 ， ∂ L ∂ λ i = ∑ j = 1 n x i j * − a i = 0 ， x i j ≥ 0 ， 对偶可行性： d i j ≥ 0 ， 互补松弛性： d i j ( − x i j * ) = 0. 为后续描述方便，我们首先考虑运输问题1.1的对偶问题 max w = ∑ i = 1 m a i u i + ∑ j = 1 n b j v j ， (9) s .t .   { u i + v j ≤ c i j ， u i ， v j 无 符 号 限 制 . 由平稳性条件和对偶可行性条件 ( − λ i ) + ( − μ j ) = c i j − d i j ≤ c i j . (10) 则 − λ i ， − μ j 为对偶问题1.1的可行解。 另一方面，由平稳性条件和互补松弛性条件知 ( − λ i − μ j ) x i j * = ( c i j − d i j ) x i j * = c i j x i j * . (11) 从而有 ∑ i = 1 m ∑ j = 1 n ( − λ i − μ j ) x i j * = ∑ i = 1 m ∑ j = 1 n c i j x i j * ， (12) 即 ∑ i = 1 m ( − λ i ) ∑ j = 1 n x i j * + ∑ j = 1 n ( − μ j ) ∑ i = 1 m x i j * = ∑ i = 1 m ∑ j = 1 n c i j x i j * . (13) 结合原始可行性条件，我们有 ∑ i = 1 m a i ( − λ i ) + ∑ j = 1 n b j ( − μ j ) = ∑ i = 1 m ∑ j = 1 n c i j x i j * . (14) 由此，满足Karush-Kuhn-Tucker条件的 x i j * ， − λ i ， − μ j 分别是原问题(1)和对偶问题(9)的可行解，且满足目标函数相等，见等式(14)。则它们分别为原问题和对偶问题的最优解。进一步，结合(10)，我们有： 当 x i j 为基变量时， d i j = 0 ，则 − λ i − μ j = c i j ，根据位势的定义(3)， − λ i ， − μ j 为位势； 当 x i j 为非基变量时， d i j = c i j + λ i + μ j 恰为检验数。