两阶段法求解线性规划求解

用两阶段法求解
min f=2x1-x2+x3 
s.t x1+2x2- x3=1 
    2x1+ x2+ x3=5 
    x1- x2+2x3=4 
    xi>=0,i=1,2,3 
引入人工变量x4,x5,x6,使
min g=x4+x5+x6
s.t x1+2x2-x3+x4=1
    2x1+x2+x3+x5=5
    x1-x2+2x3+x6=4
    xi>=0,i=1,...,6
其单纯形表是
基变量 x1    x2   x3   x4   x5   x6   右项
-g         0      0     0    1     1     1        0
x4         1      2    -1    1     0     0       1
x5         2      1    1     0     1     0        5
x6         1     -1    2     0     0     1      4
-f           2     -1    1     0     0    0       0
将x4,x5,x6这三行乘以-1加到-g这行，得表
基变量 x1    x2   x3   x4   x5   x6   右项
-g        -4    -2   -2     0    0     0    -10
x4       1      2    -1    1    0      0     1
x5       2      1     1     0    1     0     5
x6        1     -1     2     0    0     1     4
-f         2     -1    1     0    0    0      0
因为-4<-2<0，所以x1入基，又1/1<5/2<4/1，所以x4出基，
并把x4所在行的主元化成1，且把x1所在列非零项化成零,得表
基变量 x1    x2   x3    x4   x5   x6   右项
-g          0      6    -6    4    0     0    -6
x1         1      2     -1    1    0     0    1
x5         0     -3     3    -2    1     0    3
x6        0      -3    3    -1    0     1    3
-f         0      -5    3    -2    0     0    -2
因为-6<0，所以x3入基，又3/3=3/3选取x5出基，
并把x5所在行的主元化成1，且把x3所在列非零项化成零,得表
基变量 x1    x2   x3   x4   x5    x6   右项
-g         0      0    0    0     2     0    0
x1        1      1    0    1/3 1/3   0    2
x3        0     -1    1    -2/3 1/3   0    1
x6        0     0     0    1     -1     1    0
-f          0     -2     0    0     -1    0    -5
目标gmin=0说明原问题有最优解，把上面所得的结果作为求解
原问题的初始基本可行解进行单纯形法迭代。接上表有
基变量 x1    x2   x3      右项
-f          0     -2   0        -5
x1        1      1    0        2
x3        0     -1    1       1
因为-2<0，所以x2入基，x1出基，并把x1所在行的主元化成1，
且把x2所在列非零项化成零,得表
基变量 x1    x2   x3      右项
-f         2      0     0        -1
x2        1     1     0         2
x3        1     0      1         3
f行的左边全非负，迭代完毕。
最优解是x1=0,x2=2,x3=3,fmin=1=2x1-x2+x3=2*0-2+3