管理运筹学讲义第1章线性规划.pptVIP

第1章 线性规划 非标准型转化举例 minZ=x1+2x2-3x3 x1+x2+x3 ≤9 -x1-2x2+x3 ≥2 3x1+x2-3x3=5 x1≤0，x2≥0， x3无约束 注： 单纯形是指0维中的点，一维中的线段，二维中的三角形，三维中的四面体，n维空间中的有n+1个顶点的多面体。例如在三维空间中的四面体，其顶点分别为(0，0，0)，(1，0，0)，(0，1，0)，(0，0，1)。具有单位截距的单纯形的方程是∑xi≤1，并且xi≥0，i=1,2,…,m。 其步骤总结如下： 2. 单纯形法和单纯形表 假定所有b1，…bm≥0。显然问题有基可行解X1=（b1，b2，…，bm，0，…，0，），相应目标函数值 为了判定X1是否最优解，首先应当把目标函数z用非基变量表示出来。由约束方程组得到： 其中： 上式是目标函数通过非基变量表达的公式，每个变量的系数称为它的检验数。显然，基变量的检验数永远为0。而非基变量xj的检验数σj是引入xj一个单位时目标函数z的改变量；只有σj0时，方值得让xj入基。 （2）单纯形表 这个表格称为单纯形表，它的右半部第二行指明各个变量的位置，第一行是它们的价值系数，其下方m行是方程组的系数矩阵，它包括一个m阶单位阵。左半部3列中，b列是右端常数，XB列依次标明各方程的基变量，CB是相应的价值系数。表的最后一行称为检验数行，填写各个变量的检验数。单纯形表的主要部分是约束方程组的增广矩阵和检验数行。应当注意的是，系数矩阵中必须包含一个m阶单位阵。 单纯形表的便利之处首先在于，从表中可以直接读出基可行解X：xi=bi（i=1，…，m），其它xj=0；相应目标函数值是CB列与b列元素乘积之和；其次，每个变量xj（包括基变量在内）的检验数σj，等于cj减去CB列元素与表中xj的系数列向量元素乘积之和 更有益的是，单纯形法的全部分析和计算过程，可以比较方便地在单纯形表中完成。下面四个法则总结了运算各个环节的具体操作方法。 3. 单纯形法的基本法则 法则1: 最优性判定法则 若对基可行解X1，所有检验数σj≤0，则X1为最优解。 证: 由公式 可知，当所有σj ≤0时， 总有Z≤Z1。而当X= X1 =(b1，…bm，0，…，0)T时，Z=Z1，因此X1是最优解。 换个角度来说，由于某个非基变量检验数非正，引入任意非基变量都不能使Z值上升，故X1是最优解。对应于最优解的单纯形表称为最优表（或最终表）。 反之，如果存在某个σj0，则X1不是最优解，计算过程要继续下去。首先应该选择换入变量。 法则2: 换入变量确定法则 设 ，则xk为换入变量。 按照检验数的意义，任何具有正检验数的非基变量均可入基，都能使目标函数值上升；选择具有最大正检验数的变量入基，目的是为了使目标函数尽快上升。 法则3: 换出变量确定法则 设xk为换入变量，变量xm+1，xm+2，…，xn中除xk外均保持为0，约束方程组实际成为： 设xk为换入变量，并设 则最小比值所在行的原基变量xl为换出变量。 强调：这个法则的目的是，保证下一个基本解的可行性，违背这一法则，下一个基本解一定包含负分量，即不是可行解。 在单纯形表中，换入变量所在列称为主元列，换出变量所在行称为主元行，两者交叉处的元素称为主元。主元在迭代过程中起重要作用，通常用括号括起来。 法则4: 换基迭代运算法则 确定换入变量xk、换出变量xl以后，下一个环节是以xk替换xl，把约束方程组变换为对新基变量的解出形式。前面已经提及，这种变换等价于方程组的增广矩阵实施行的初等变换，结果使得变量xk的系数列向量成为单位向量，且其第l个分量为1。因此得到以下法则： 对方程组的增广矩阵实施行的初等变换，使主元alk化为1，主元列其它元素化为0。具体地说，第l行乘以1/alk，并将第l行的－aik/alk倍加到第i行上去，i=1，2，…，m，i≠l。 与此同时，对单纯形表的CB列和XB列作适当调整，再按公式计算检验数，就得到下一张单纯形表。然后返回法则1，检验新表的最优性。循环使用上述法则，直到求出最优解或判定无最优解。 单纯形计算表 c 50 30 0 0 cB xB b x1 x2 x3 x4 ? 0 x3 120