运筹学-对偶单纯形法概要.pptVIP

本节内容安排 对偶单纯形法的求解思路 对偶单纯形法的求解步骤 对偶单纯形法是根据对偶原理和单纯形法的 原理而设计出来求解 原LP的一种方法。 采用的技术是在原问题的单纯形表格上进行对偶处理。 注意: 对偶单纯形法不是求解对偶问题的单纯 形法。 1 单纯形法中原问题(max)的最优解满足的条件: （ 1）是基本解 （2）可行解（XB =B-1 b≥0); （3 ) 检验数C－CBB-1A?0 ， -CB B-1≤0 即YA ?C， Y?0 即对偶解可行 2 普通单纯形法的求解思路: 从满足（1），（2）的一个初始基本可行解出发 (此时原LP问题中，b列保持≥0 ， 对偶的解 一般为非可行基解）， 通过逐步迭代，增大原目标函数值， 每一步迭代，都得到一个基本可行解， 并且逐步迭代实现检验数行≤0 (对偶解可行)。 1 对偶单纯形法求解思路： 换个角度考虑LP求解过程 从满足 （1）（3）的一个非可行基解（检验数行保 持≤0）出发，（此时对偶问题的解一般为可行解）， 通过逐步迭代直至（2）得到满足， 即直到实现到b列所有的值≥0， 原问题的解在迭代过程中 从非可行解变成可行解， 最终达到最优解， 此时，对偶问题也达到最优解。 原问题基可行解 最优解判断 3 对偶单纯形法的使用条件： ①原问题的初始基解的检验数全部≤0； ②b列至少一个元素 0； 4 实施对偶单纯形法的基本原则： 在保持对偶可行的前提下进行基变换—— 每一次迭代过程中取出基变量中的一个负分量 作为换出变量 去替换某个非基变量(作为换入变量), 使原始问题的非可行解向可行解靠近。 第一步： 构造初始单位阵，确定原问题（max ) 的初始基B，使所有检验数 C j - Zj = σj = Cj - CB B -1 Pj ≤ 0， 即 Y = CB B -1 （b 列的值）是对偶可行解，建立初始单纯形表。 第二步： 可行性检验。 检验 b 列 和σj 行（即检查基变量的取值） 若 bi ≥0 （XB = B -1 b ≥0）, σj ≤ 0 , 则原问题得到最优解 ，计算停; 若bi 0 , σj ≤ 0 , 则用对偶单纯形法进行换基迭代. 第三步 先确定换出变量 解答列（b 列）中的负元素对应的基变量出基， 相应的行为主元行。 一般选最小的负元素出基， 即若min { ( B -1 b )i| (B -1b )I 0 } = ( B–1 b )l 则选取 x l 为换出变量. 检验第l 行中非基变量 xj 的系数 αlj , 若所有的αlj ≥ 0，则LP 问题 无可行解， （下面进行说明），此时计算结束。 否则转下步 例6 应用对偶单纯型法求解下面的线性规 划问题 对偶单纯形法与原始单纯形法内在的对应关系 单纯形法和对偶单纯形法步骤 例 用对偶单纯形法求解 是 是 是 是 否 否 否 否 所有 所有 得到 最优解 计算 计算 典式对应原规划的基本解是可行的 典式对应原规划的基本解的检验数≤0 所有 所有 计算 计算 以aek为中心元素进行迭代 以alk为中心元素进行迭代 停 没有最优解 没有最优解 普通单纯形法 对偶单纯形法 min z = 2x1+4x2+6x3 s.t. 2x1 - x2 + x3 ≥10 x1+2x2+2x3 ≤12 2x2 - x3 ≥4 x1，x2，x3 ≥0 解： 将问题化为： max z′= -2x1- 4x2 - 6x3 s.t. -2x1 + x2 - x3+ x4=-10 x1+2x2+2x3 +x5=12 -2x2 + x3 +x6 =-4 xj ≥0 ( j=1,2,…,6) 检验数 c cB xB -2 -4 -6 0 0 0 x1 x2 x3 x4 x5 x6 0 0 0 x4 x5 x6 -2 1 -1 1 0