2[1].2_改进单纯形法.ppt

第五节改进单纯形算法

在每次迭代过程中不必要地计算了很多与迭代无关的数字,影响了计算效率.

单纯形法的迭代过程实质上是从一组基到另一组基的变换.

而每次迭代中真正有用的数字是基变量列数字、基的逆矩阵、非基变量检验数以及最大正检验数所对应的非基变量系数列向量。;引例：

maxZ=6x1＋8x2

2x1十x2≤12

x1十4x2≤20

xj≥0

标准化得：

maxZ=6x1＋8x2

2x1十x2十x3＝12

x1十4x2十x4＝20

xj≥0;求解：;第二步迭代中：

基变量XB=(x3、x2)

基矩阵B2＝（P3、P2）＝

B2－1→

B2－1＝正是x3、x4在第二步迭代表中的系数;根据矩阵理论，第二步迭代表中的任何数都可由B2－1左乘原始数据得到

如：p‘2＝B2－1p2＝

p‘3＝B2－1p3＝

p‘4＝B2－1p4＝

b‘＝B2－1b＝;检验数σj＝Cj－CBB－1Pj

σ1＝C1－CBB2－1P1＝6－（08）

＝6－2＝4

σ4＝C4－CBB2－1P4＝0－（08）

＝0－2＝－2

令Y＝CBB2－1

称单纯形乘子;第三步迭代中：

基变量XB=(x1、x2)

基矩阵B3＝（P1、P2）＝

B3－1→

B3－1＝

正是x3、x4

在第三步

迭代表中

的系数;根据矩阵理论，第三步迭代表中的任何数都可由B3－1左乘原始数据得到

如：p‘1＝B3－1p1＝

b‘＝B3－1b＝

检验数σj＝Cj－CBB－1Pj

σ1＝C1－CBB3－1P1

＝6－（68）

＝6－6＝0;在单纯形法的迭代计算中重复计算了很多与迭代过程无关的数字。

单纯形法的迭代计算实际上是对约束方程的系数矩阵实施行的初等变换。

显然有

上述公式是改进单纯形法计算的依据，也是下一章中要讲述的灵敏度分析等内容的基础，CB是基变量的目标函数中的系数向量，Pj和Cj是初始单纯形表中非基变量Xj的系数向量和它在目标函数中的系数值，Pj’是新单纯形表中非基变量的系数向量。;改进单纯形法的计算步骤可以归结为：

(1)在下一步迭代的基变量确定后，求新单纯形表中基矩阵B的逆矩阵B－1，并求基可行解XB＝B－1b

(2)计算和单纯形乘子Y=CBTB-1

并计算非基变量的检验数

σN＝CN－CBB－1N＝CN－YN

σj＝Cj－CBB－1Pj＝Cj－YPj，

如果有σj≤0，达到最优解，计算结束。否则转入下一步

(3)根据max{?l∣?l0，l∈IN}=?k,确定xk为???入变量(即为新基的基变量)

计算P‘k＝B－1Pk＝（a1k’、a2k’、……amk’）

若aik’≤0，线性规划问题有无界解，计算结束。否则转入下一步

（4）按最小比值原则来确定第r行基变量xr为换出变量;定理在单纯形法的相邻两次迭代中，没迭代前的可行基为

B＝(P1，P2，…，Pr－1，Pr，Pr＋1，…Pm)

经过换基运算后，得到另一个可行基

＝(P1，P2，…，Pr－1，Pr，Pr＋1，…Pm)

则迭代后所得基的逆矩阵为

＝ErkB-1;例1.已知:迭代前的基为:

1000

B=(p1,p2,p3,p4)=0100 =I

0001

;解:由定理4中初等变换矩阵Erk的构造方法可知,本例中的;第一次迭代:取初始基

B0=(p1,p4,p6)=I

相应的变量x1,x4,x6为初始基变量，x2,x3,x5为非基变量

B0-1=I

CB0=(C1,C4,C6)=(0,0,0),

初始基可行解

xB0=B0-1b=(7,12,10)T;求单纯形乘子CB0=(C1,C4,C6)=(0,0,0)

Y0=CB0.B0-1=(0,0,0),

求检验数σN=CN–YN,并确定进基变量:显然有

σ1=σ4=σ6=0;

σ2=－1,σ3=＋3,σ5=－2

由于σ3=3＞0k＝3

故应取X3为进基变量,进基列向量为:;确定出基变量:由

θ=min﹛－，12/4,10/3﹜=3;r＝2

可知应取X4为出基变量,