计算机系线性规划运筹学16.pptVIP

* 不妨假设A＝（P1，P2，…，Pn）中前m个列向量构成一个可行基，记为B=(P1，P2，…，Pm)。矩阵A中后n－m列构成的矩阵记为N＝（Pm+1,…Pn）,则A可以写成分块矩阵A=（B，N）。对于基B，基变量为XB=（x1,x2，…，xm ）T, 非基变量为XN=(xm+1,xm+2,…xn)T。 则X可表示成 同理将C写成分块矩阵C=（CB，CN）， CB=(C1,C2,…,Cm), CN=(Cm+1Cm+2，…，cn) 则AX=b可写成 1.5 单纯形法 Simplex Method * 因为r（B）=m（或|B|≠0）所以B —1存在，因此可有 令非基变量XN=0，XB=B—1b，由 B是 可行基的假设，则得到基本可行解 X=(B－1b，0)T 将目标函数写成 1.5 单纯形法 Simplex Method * 得到下列五个计算公式： (令XN=0) 1.5 单纯形法 Simplex Method * 上述公式可用下面较简单的矩阵表格运算得到，设初始矩阵单纯形表1-15 将B化为I（I为m阶单位矩阵）,CB化为零,即求基本可行解和检验数。用B－1左乘表中第二行,得到表1-16 XB XN b XB I B-1N B-1b Cj-Zj CB CN 0 XB XN b XB B N b Cj-Zj CB CN 0 表1－15 表1－16 1.5 单纯形法 Simplex Method * 再将第二行左乘－CB后加到第三行,得到 λΝ XB －Z0 XB XN b XB I B-1N B-1b λ＝Cj-Zj 0 CN-CBB－1N －CBB－1b 表1－17 1.5 单纯形法 Simplex Method * 五个公式的应用 【例1.24】线性规划 已知可行基 求(1)单纯形乘子π; (2)基可行解及目标值; (3)求λ3; (4)B1是否是最优基,为什么; (5)当可行基为 时求λ1及λ3。 1.5 单纯形法 Simplex Method * 【解】(1)因为B1由A中第一列、第二列组成,故x1、x2为基变量，x3、x4、x5为非基变量,有关矩阵为 CB=(c1,c2)=(1,2) CN=(c3,c4,c5)=(1,0,0) 故单纯形乘子 1.5 单纯形法 Simplex Method * (2)基变量的解为 故基本可行解为 目标函数值为 1.5 单纯形法 Simplex Method * (3) 求λ3 1.5 单纯形法 Simplex Method * (4) 要判断B1是不是最优基,亦是要求出所有检验数则否满足λj≤0,j=1…,5。x1,x2是基变量, 故λ1=0,λ2=0,而 剩下来求λ4,λ5，由λN计算公式得 因λj≤0, j=1,…,5,故B1是最优基。 1.5 单纯形法 Simplex Method * (5) 因B2是A中第四列与第二列组成的,x4、x2是基变量x1、x3、x5是非基变量,这时有 即 1.5 单纯形法 Simplex Method * 【例1.26】求解线性规划 解 用大M单纯形法，加入人工变量x4、x5，构造数学模型 1.5.4 退化与循环 1.5 单纯形法 Simplex Method * Cj 1 2 1 M M b θ CB XB x1 x2 x3 x4 x5 (1) M M x4 x5 1 4 －2 －9 [4] 14 1 0 0 1 4→ 16 1 8/7 λj 1－5M 2+11M 1－18M↑ 0 0 (2) 1 M x3 x5 [1/4] 1/2 －1/2 －2 1 0 1/4 －7/2 0 1 1→ 2 4 4 λj 3/4－1/2M↑ 5/2+2M 0 －1/4+9/2M 0 (3) 1 M x1 x5 1 0 －2 [－1] 4 －2 1 －4 0 1 4 0→ λj 0 4+M↑ －3+2M －1+5M 0 (4) 1 2 x1 x2 1 0 0 1 8 [2] 9 4 －2 －1 4 0→ λj 0 0 －11↑ M－17 M－4 (5) 1 1 x1 x3 1 0 －4 1/2 0 1 1 2 2 －1/2 4 0 λj 0 15/2 0 M－17 M－3/2 表1－18 1.5 单纯形法 Simplex Method * 单纯形法迭代对于大多数退化解时是有效的，很少出现不收敛的情形．1955年Beale提出