运筹学-对偶单纯形法.ppt

;本节内容安排; 对偶单纯形法是根据对偶原理和单纯形法的 原理而设计出来求解 原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 (对偶解可行)。; 3 普通单纯型法的求解过程： 对原问题的一个基可行解，判别是否所有检验数 非正 cj-zj?0(j=1, …,n)； 若是，又基变量中无非零人工变量，即找到问题最优解，基变量中含有非零人工变量，则无最优解； 若否，再迭代，找出相邻的目标函数值更大的基可行解，并继续判别，只要最优解存在，就一直迭代下去，直到找出最优解为止.; 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 问题 无可行解， （下面进行说明），此时计算结束。 否则转下步;当bl0，而对所有j=1,…,n，有alj?0， 则原问题无可行解。 证明：xl+al,m+1xm+1+…+al,nxn=bl 因alj?0(j=m+1,…,n)，又bl0， 故有xl0， 即不可能存在xj?0(j=1,…,n)的解， 故原问题无可行解， 此时对偶问题的目标函数值无界。;若有：Min { cj – zj / αlj|αlj 0 , x j 为非基变量} = ck – zk /αlk 则确定 xk 为换入变量，相应的列为主元列， 标出主元素αlk , 应用矩阵的初等行变换得到新的单纯形表。;alk为主元素 xk为进基变量;(1)对alj?0，因cj－zj?0，故 ，又因主元素alk0， 故有 ，所以(cj－zj)? ?0； ; 第五步：用换入变量替换换出变量，得到一个新的基，对新的基再检查是否所有 如果是，得原问题的最优解； 如果否，回到第一步再重复