运筹学线性规划与目标函数讲义资料.ppt

* 3.5 迭代(旋转运算) 设x1,x2,…,xm为基变量，对应的系数矩阵是m×m单位阵I，它是可行基。令非基变量xm+1,xm+2,…,xn为零，即可得到一个基可行解。若它不是最优解，则要另找一个使目标函数值增大的基可行解。这时从非基变量中确定xk为换入变量。显然这时θ为在迭代过程中θ可表示为 其中 是经过迭代后对应于 的元素值。 * 3.5 迭代(旋转运算) 按θ规则确定xl为换出变量，xk, xl的系数列向量分别为 * 3.5 迭代(旋转运算) 为了使xk与xl进行对换，须把Pk变为单位向量，这可以通过(1-33)式系数矩阵的增广矩阵进行初等变换来实现。 * 3.5 迭代(旋转运算) 变换的步骤是：(1) 将增广矩阵(1-34)式中的第l行除以al k，得到 (2) 将(1-34)式中xk列的各元素，除al k变换为1以外，其他都应变换为零。其他行的变换是将(1-35)式乘以aik(i≠l)后，从(1-34)式的第i行减去，得到新的第i行。 * 3.5 迭代(旋转运算) 由此可得到变换后系数矩阵各元素的变换关系式 是变换后的新元素。 * 3.5 迭代(旋转运算) (3) 经过初等变换后的新增广矩阵是 * 3.5 迭代(旋转运算) (4) 由(1-36)式中可以看到x1,x2,…,xk,…，xm的系数列向量构成m×m单位矩阵;它是可行基。当非基变量xm+1,…，xl,…,xn为零时，就得到一个基可行解X(1) 在上述系数矩阵的变换中，元素al k称为主元素，它所在列称为主元列，它所在行称为主元行，元素al k位置变换后为1。 * 3.5 迭代(旋转运算) 例7 试用上述方法计算例6的两个基变换。 解：例6的约束方程组的系数矩阵写成增广矩阵 当以x3,x4,x5为基变量，x1,x2为非基变量，令x1,x2=0，可得到一个基可行解 X(0)=(0,0,8,16,12)T * 3.5 迭代(旋转运算) 现用x2去替换x5，于是将x3,,x4, x2的系数矩阵变换为单位矩阵，经变换后为 令非基变量x1,x5=0，得到新的基可行解 X(1)=(0,3,2,16,0)T * 第4节 单纯型法的计算步骤 4.1 单纯型表 4.2 计算步骤 * 4.1 单纯型表 为了便于理解计算关系，现设计一种计算表，称为单纯形表，其功能与增广矩阵相似。 将(1-22)式与目标函数组成n+1个变量，m+1个方程的方程组。 * 线性规划的方程组 * 线性规划的方程组 为了便于迭代运算，可将上述方程组写成增广矩阵形式 * 线性规划的方程组 若将z看作不参与基变换的基变量，它与x1,x2,…,xm的系数构成一个基，这时可采用行初等变换将c1,c2,…,cm变换为零，使其对应的系数矩阵为单位矩阵。得到 表1-2 * 线性规划的方程组 表1-2的说明 XB列中填入基变量，这里是x1，x2,…，xm； CB列中填入基变量的价值系数，这里是c1,c2,…,cm；它们是与基变量相对应的； b列中填入约束方程组右端的常数； cj行中填入基变量的价值系数c1,c2,…,cn； θi列的数字是在确定换入变量后，按θ规则计算后填入； 最后一行称为检验数行，对应各非基变量xj的检验数是 * 4.2 计算步骤 表1-2称为初始单纯形表，每迭代一步构造一个新单纯形表。 计算步骤： (1) 按数学模型确定初始可行基和初始基可行解，建立初始单纯形表。 (2) 计算各非基变量xj的检验数， 检查检验数，若所有检验数 则已得到最优解，可停止计算。否则转入下一步。 * 4.2 计算步骤 (3)在σj＞0,j=m+1,…,n中，若有某个σk对应xk的系数列向量Pk≤0，则此问题是无界，停止计算。 否则，转入下一步。 (4) 根据max(σj＞0)=σk，确定xk为换入变量，按θ规则计算 * 4.2 计算步骤 (5) 以alk为主元素进行迭代(即用高斯消去法或称为旋转运算)，把xk所对应的列向量 将XB列中的xl换为xk，得到新的单纯形表。重复(2)～(5)，直到终止。 * 4.2 计算步骤 现用例1的标准型来说明上述计算步骤 (1) 取松弛变量x3,x4,x5为基变量，它对应的单位矩阵为基。这就得到初始基可行解X(0)=(0,0,8,16,12)T 将有关数字填入表中，得到初始单纯形表，见表1-3。表中左上角的cj是表示目标函数中各变量的价值系数。在CB列填入初始基变量的价值系数，它们都为零。 * 4.2 计算步骤 表 1-3 * 4.2 计算步骤 计算非基变量的检验数 各非基变量的检验数为 σ1=c1?z