线性规划二阶段法

1、关于线性规划的二阶段法第一张，PPT共三十二页，创作于2022年6月 线性规划的人工变量法目前有两种方法：大M法和二阶段法。人工变量法在讨论单纯形法时，我们总是假定AXb的系数矩阵A的秩r(A)=mn,或者已有一个可行基。但是，在许多问题中，初始可行基是不容易找到的，或者A不满秩。这样单纯形法就很难进行。所以，我们要探讨如何寻找第一个可行基？第二张，PPT共三十二页，创作于2022年6月大M法（1）把原问题化为下列形式：其中M是任意大的正数第三张，PPT共三十二页，创作于2022年6月大M法（2）关于大M法的几点注释：(1)在引入人工变量之前，约束条件已是等式，为了这些等式得到满足，因此在最优

2、解中人工变量取值必须为零；为此在目标函数中人工变量的取值为充分小的负数，用“-M”代表；(2)若在单纯形表中有j0，且存在非零人工变量，则原问题无可行解；若基变量中不再含有非零的人工变量，这表示原问题有解；(3)引入的人工变量个数越少越好，只要出现单位矩阵作为基阵即可。第四张，PPT共三十二页，创作于2022年6月大M法举例（1）例解：将原问题化为标准形为：第五张，PPT共三十二页，创作于2022年6月大M法举例（2）引入的人工变量个数越少越好引入人工变量y2,y30,由大M法得辅助问题为：其中大M为任意大的正数得上述辅助问题的单纯形初表为：第六张，PPT共三十二页，创作于2022年6月大M法

3、举例（3） T(B1)XB b x1 x2 x3 x4 x5 y2 y3x4 y2 y3419 -z10M 1 1 1 1 0 0 0 -2 1 -1 0 -1 1 0 0 3 1 0 0 0 1-2M-3 4M 1 0 -M 0 0 T(B2)x4 x2 y3 -z3 3 0 2 1 1 -1 0 1 -2 1 -1 0 -1 1 0 6 6 0 4 0 3 -3 1 6M6M-30 4M+1 03M-4M0人工变量优先出基第七张，PPT共三十二页，创作于2022年6月大M法举例（4） T(B3)XB b x1 x2 x3 x4 x5 y2 y3x4 x2 x1031 -z 30 0 0 1

4、 -1/2 1/2 -1/2 0 1 1/3 0 0 0 1/31 0 2/3 0 1/2 -1/2 1/60 0 3 0 3/2 -M-3/2 -M+1/2 T(B4)x4 x2 x3 -z00 0 0 1 -1/2 1/2 -1/2 5/2-1/2 1 0 0 -1/4 1/4 1/4 3/23/2 0 1 0 3/4 -3/4 1/4 -3/2-9/20 0 0-3/4-M+3/4 -M-1/4第八张，PPT共三十二页，创作于2022年6月大M法举例（5）原线性规划问题的最优解为因为在单纯形表T(B4)中，非基变量检验数均小于等于零，且人工变量均为非基变量，取值为零，故原线性规划问题达到

5、最优。第九张，PPT共三十二页，创作于2022年6月线性规划的二阶段法(1)原线性规划问题为第一阶段:y1， y2， ym称为人工变量构造原(LP)的辅助问题第十张，PPT共三十二页，创作于2022年6月线性规划的二阶段法(2)原(LP)的辅助问题的标准形式为：辅助问题必有最优解第十一张，PPT共三十二页，创作于2022年6月线性规划的二阶段法(初始单纯形表1)辅助问题的标准形式的系数矩阵为：第十二张，PPT共三十二页，创作于2022年6月线性规划的二阶段法(初始单纯形表2)用单纯形法求解,最终得到辅助问题的最优单纯形表T(B*)第十三张，PPT共三十二页，创作于2022年6月二阶段法的计算步

6、骤:第二步 若 w*0,则原线性规划无可行解,停止求解,否则转第三步.第一步 用单纯形法求辅助问题的最优单纯形表T(B*)和最优值w*. 第三步 T(B*)中基变量中不含人工变量y,则把T(B*)中人工变量所在列划去,把检验数行用原规划的目标函数的系数替代再把基变量的检验数化为0,即得原规划的一个可行基的单纯形表.再用单纯形法迭代,直到终止.否则转第四步. 第四步 w*=0,T(B*)中基变量中含有人工变量yr,若yr所在行的对应的X系数全为0 ,则划去T(B*)中yr所在行和所在的列,转第三步。否则以某变量xS的系数brs0为轴心项进行换基迭代后转第三步。 第十四张，PPT共三十二页，创作于

7、2022年6月线性规划的二阶段法举例（1）例1. 用二阶段法求解下列(LP)问题 解: 化原问题为标准形式 则第一阶段的辅助问题为 第十五张，PPT共三十二页，创作于2022年6月线性规划的二阶段法举例（1-2） 辅助问题的标准形式为 第十六张，PPT共三十二页，创作于2022年6月线性规划的二阶段法举例（例1-3）则辅助问题的单纯表 T(B1)XB b x1 x2 x3 x4 x5 x6 y2 y3x4 y2 y3643 w 7 1 1 1 1 0 0 0 0 1 0 -1 0 -1 0 1 0 0 1 -1 0 0 -1 0 1 1 1 -2 0 -1 -1 0 0 T(B2)x4 x1

8、y3 w2 0 1 2 1 1 0 -1 04 1 0 -1 0 -1 0 1 03 0 1 -1 0 0 -1 0 1 301 -1 0 0-1-10第十七张，PPT共三十二页，创作于2022年6月二阶段法举例（例1-4） T(B2)x4 XB b x1 x2 x3 x4 x5 x6 y2 y3x1 y3 w2 0 1 2 1 1 0 -1 04 1 0 -1 0 -1 0 1 03 0 1 -1 0 0 -1 0 1 30 1 -1 00 -1-10 x2 T(B3)x1y3w2 0 1 2 1 1 0 -1 04 1 0 -1 0 -1 0 1 01 0 0 -3 -1 -1 -1 1

9、1 1 0 0 -3 -1 -1 -1 0 0 因为辅助问题的最优值W*=10，则原问题无可行解。第十八张，PPT共三十二页，创作于2022年6月线性规划的二阶段法 例2-1例2 用二阶段法解线性规划 解:引进人工变量y20，建立第一阶段的辅助问题为 第十九张，PPT共三十二页，创作于2022年6月线性规划的二阶段法举例（例2-2） 辅助问题的标准形式为 则第一阶段的初始单纯形表 为T（B1）第二十张，PPT共三十二页，创作于2022年6月 XB b x1 x2 x3 x4 y2 x2 y2 2 31/2 1 1/2 -2/3 0 3/2 0 3/4 0 1 w 33/2 0 3/4 0 0例

10、2-3 T(B1) x2 x1w 1 0 1 1/4 -2/3 -1/321 0 1/2 0 2/3 00000 T(B2)-1第二十一张，PPT共三十二页，创作于2022年6月例2-4得w*=0，且基变量中不含有人工变量，则划去T(B2)中y2所在列,把检验数行用原问题的目标函数的系数替代再把基变量的检验数化为0,即得第二阶段的初始单纯形表. T(B2)c-4-300 XB b x1 x2 x3 x4 y2 x2 x1 1 20 1 1/4 -2/3 -1/3 1 0 1/2 0 2/3 w 00 0 0 0 -1-z 11 0 0 11/4 -2 第二十二张，PPT共三十二页，创作于202

11、2年6月例2-5则原问题的一个初始单纯形表如下x2 x1 12 0 1 1/4 -2/3 1 0 1/2 0 -z 0 0 11/4 -2 XBb x1 x2 x3 x4 x2 x3 -z 110-1/2 1 0 -2/3 4 2 0 1 0 0-11/2 0 0 -2 原问题最优解X*=(0,0,4,0)T原问题最优值为z*=0第二十三张，PPT共三十二页，创作于2022年6月例3-1例3 用二阶段法解下列线性规划 解: 该问题的标准形式为：第二十四张，PPT共三十二页，创作于2022年6月线性规划的二阶段法举例（例3-2） 辅助问题的标准形式为 则第一阶段的单纯形初表为T（B1）引进人工变

12、量y1,y2,y3，建立第一阶段的辅助问题为：第二十五张，PPT共三十二页，创作于2022年6月例3-3 XB b x1 x2 x3 x4 y1 y2 y3 y1 y2 y3 1 2 1-1 1 1 0 1 0 0 1 1 0 1 0 1 0 2 0 -1 1 0 0 1 w 4 2 2 0 2 0 0 0 T(B1)y1 y2w3/2 0 1 1/2 1/2 1 0 1/23/20 1 1/2 1/2 0 1 -1/2 3021 10 T(B2)0 x11/2 1 0 -1/2 1/2 0 0 1/2 -1第二十六张，PPT共三十二页，创作于2022年6月例3-4 XB bx1 x2 x3

13、x4 y1 y2 y3 y1 y2 x1 3/2 3/2 1/2 0 1 1/2 1/2 1 0 1/2 0 1 1/2 1/2 0 1 -1/2 1 0 -1/2 1/2 0 0 1/2 w 3 0 2 1 1 0 0 -1 T(B2)x2 y2w3/2 0 1 1/2 1/2 1 0 1/200 0 0 0 -1 1 -1 0000 00 T(B3)0 x11/21 0 -1/2 1/2 0 0 1/2 -2第二十七张，PPT共三十二页，创作于2022年6月例3-5x2 y2w3/2 0 1 1/2 1/2 1 0 1/200 0 0 0 -1 1 -1 000 0 00 T(B3)0 x

14、11/21 0 -1/2 1/2 0 0 1/2 -2 c 2 100得w*=0，T(B3)中基变量中含有人工变量y2,且y2所在行的对应的X系数全为0 ,则划去T(B3)中y2所在行，此时，基变量中不再含有人工变量，则划去T(B3)中人工变量所在列,把检验数行用原问题的标准形式的目标函数的系数替代再把基变量的检验数化为0,即得第二阶段的初始单纯形表.XBb x1 x2 x3 x4 y1 y2 y3第二十八张，PPT共三十二页，创作于2022年6月例3-6 XB bx1 x2 x3 x4 x2 x1 3/2 1/20 1 1/2 1/2 1 0 -1/2 1/2 c 0 2 1 0 0 T(B3) T(B4)z00-5/21/2-3/2x3x1320 2 1 11 1 0 1 z-40 -1 0 -2原最优解X*=(2,0,3,0)T最优值为z*=-4第二十九张，PPT共三十二页，创作于2022年6月Yes Yes NO NO Yes NO 二阶段法的框图表示为如下:开始w*