5第五节 线性规划.ppt

二、单纯形表 将(2.19)式与目标函数组成n+1个变量，m+1个方程的方程组。 为了便于迭代运算，可将上述方程组写成增广矩阵 若将z看作不参与基变换的基变量，它与x1,x2,…,xm的系数构成一个基，这时可采用行初等变换将c1,c2,…,cm变换为零，使其对应的系数矩阵为单位矩阵，得到 可根据上述增广矩阵设计计算表，见表2-2。 表2-2 cj → c1 …… cm cm+1 …… cn CB XB b x1 …… xm xm+1 …… xn …… …… …… …… …… …… -z 0 …… 0 …… XB列中填入基变量，这里是x1,x2,…,xm； CB列中填入基变量的价值系数，这里是c1,c2,…,cm； b列中填入约束方程组右端的常数； cj行中填入变量的价值系数c1,c2,…,cn； i列的数字是在确定变换入变量后，按 最后一行称为检验数行，对应各非基变量xj的检验数是 规则计算后填入。 表2-2称为初始单纯形表，每迭代一步构造一个新单纯形表。 三、单纯形法的计算步骤 (1)找出初始可行基，确定初始基可行解，建立初始单纯形表。 (2)检验各非基变量xj的检验数是 ，若 则己得到最优解，可停止计算。否则转入下一步。 (3)在 中，若有某个 对应xk的系数列向量 则此问题是无界，停止计算。否则，转入下一步。 (4)根据 ，确定xk为换入变量，按 规则计算 可确定xl为换出变量。转入下一步。 (5)以a1k为主元素进行迭代(即用高斯消去法或称为旋转运算)，把xk所对应的列向量 ←第l行 将XB列中的xl换为xk，得到新的单纯形表。重复(2)－(5)，直到终止。现用例l的标准型来说明上述计算步骤： (1)根据例l的标准型，取松弛变量x3,x4,x5为基变量，它对应的单位矩阵为基。这就得到初始基可行解 X(0)=(0，0，8，16，12)T 将有关数字填入表中，得到初始单纯形表，见表2-3。 表2-3 cj → 2 3 0 0 0 CB XB b x1 x2 x3 x4 x5 0 0 0 x3 x4 x5 8 16 12 1 4 0 2 0 [4] 1 0 0 0 1 0 0 0 1 4 - 3 -z 0 2 3 0 0 0 (2)因检验数都大于零，且P1，P2有正分量存在，转入下一步。 (3) ；对应的变量x2为换入变量。计算 它所在行对应的x5为换出变量。x2所在列和x5所在行的交叉处[4]为主元素。 (4)以[4]为主元素进行旋转运算，即初等行变换，使P2变换为(0，0，1)T，在 xB列中将x2替换x5，于是得到新表2-4。 表2-4 cj → 2 3 0 0 0 CB XB b x1 x2 x3 x4 x5 0 0 3 x3 x4 x2 2 16 3 [1] 4 0 0 0 1 1 0 0 0 1 0 -1/2 0 1/4 2 4 - -z -9 2 0 0 0 -3/4 b列的数字是 于是得到新的基可行解 目标函数的取值z=9。 (5)检查表2-4的所有cj-zj，这时有c1-z1=2；说明x1应为换入变量。重复 (2)-(4) 的计算步骤，得表2-5。 表2-5 cj → 2 3 0 0 0 CB XB b x1 x2 x3 x4 x5 2 0 3 x1 x4 x2 2 8 3 1 0 0 0 0 1 1 -4 0 0 1 0 -1/2 [2] 1/4 - 4 12 -z -13 0 0 -2 0 1/4 2 0 3 x1 x5 x2 4 4 2 1 0 0 0 0 1 0 -2 1/2 1/4 1/2 -1/8 0 1 0 -z -14 0 0 -3/2 -1/8 0 (6)表2-5最后一行的所有检验数都己为负或零，这表示目标函数值已不可能再增大，于是得到最优解 目标函数值 第四节 人工变量法 一、大M法 例8 现有线性规划问题 试用大M法求解。 解 在上述问题的约束条件中加入松弛变量，剩余变量，人工变量，得到 这里M是一个任意大的正数。 用单纯形法进行计算时，见表2-6。因本例是求Min，所以用所有cj-zj 0来判别目标函数是否实现了最小化。表2-6中的最终表表明得到最优解是 表2-6 cj → -3 1 1 0 0 M M CB XB b x1 x2 x3 x4 x5 x6 x7 0 M M x4 x6 x7 11 3 1 1 -4 -2 -2 1 0 1 2 [1] 1 0 0 0 -1 0 0 1 0 0 0 1 11 3/2 1 -z -4M -3+6M 1-M 1-3M 0 M 0 0 0 M