运筹学课件单纯形法的迭代原理.pptVIP

等式约束左端引入人工变量的目的 使约束方程的系数矩阵中出现一个单位阵，用单位阵的每一个列向量对应的决策变量作为“基变量”，这样，出现在单纯形表格中的B(i)列（即约束方程的右边常数）值正好就是基变量的取值。 ３、最优性检验和解的判别 将　　　　　代入目标函数计算： 　　第四节 单纯形法计算步骤 * 运筹学教程 * 四、单纯形法的迭代原理 ? 1、确定初始基可行解 (1)初始可行基的确定 观察法——观察系数矩阵中是否含有现成的单位阵？ LP限制条件中全部是“≤”类型的约束—— 将新增的松弛变量作为初始基变量，对应的系数列向量构成单位阵； 先将约束条件标准化，再引入非负的人工变量， 以人工变量作为初始基变量，其对应的系数列向量构成单位阵，称为“人造基”； 然后用大M法或两阶段法求解； 线性规划限制条件都是“≥”或“=”类型的约束—— ①如果限制条件中既有“≤”类型的约束，又有“≥”或“=”类型的约束，怎么办？ 构造单位阵 问题 ②初始可行基一定要选单位阵？ b列正好就是基变量的取值，因此称b列为解答列 （2）写出初始基可行解—— 令非基变量取0，基变量对应b(i)，一起构成初始基可行解 此时LP的标准型为 在约束条件中的变量系数矩阵中总会有一个单位矩阵 初始可行基 ： 当线性规划的约束条件均为≤，其松弛变量的系数矩阵为单位矩阵；当线性规划的约束条件均为≥或=，为便于找到初始基可行解，构造人工变量，人为产生一个单位矩阵。 初始基本可行解： 式中p1,…,pm 为基变量,同其所对应的x1,x2,…..,xm为基变量；其它变量xm+1,xm+2，……,xn为非基变量。令所有的非基变量等于零。 定义：两个基可行解称为相邻的，如果它们之间变换 且仅变换一个基变量。 初始基可行解的前m个为基变量， 2、基变换 代入约束条件有 系数矩阵的增广矩阵 因为p1,…,pm,是一个基，其他向量pj可以这个基的线性组合表示： 相减，然后乘上一个正数θ,加上 经过整理得到： 找到满足约束方程组 的另一点： 第j个大于0 只变换1个变量； 前m个变量必须换出1个 其中θ是X（1）的第j个坐标的值，要使X（1）是一个基可行解，对所有的i=1,…,m,存在 令这m个不等式至少有一个等号成立，当 是一个可行解。因为变量 x11, x21, xl-11, xl+11,…………xm1, xj1所对应的向量， 经过重新排列后加行b列形成的增广矩阵为： L alj×(1/alj)=1 rL×(-al-1j) +rL-1 0 -(bL/aLj)+bL-1 所以，P1,P2,…,Pl-1,Pj,Pl+1,…,Pm,是一个基。 进行初等行变换，将第L行乘上1/alj,再分别乘以 -aij,(i=1,…,l-1,l+1,…,m)加到各行，增广矩阵 的左边变成一个单位矩阵， 最优性判别 １、如果所有的检验数　　　，表明现有的顶点对应的基可行解为最优解。 ２、当所有的检验数　　　　，又对某个非基变量xj的检验数等于０，并且可以找到　0,这表明可以找到一个顶点目标函数达到最大，说明LP有无穷多个最优解。 ３、如果存在某个检验数　　0,又　　≤0, 表明目标函数达到无限，说明LP有无界解。 第一步：求初始基可行解，列出初始单纯形表。 　　将LP化为标准型，并加以整理。引入适当的松驰变量、剩余变量和人工变量，使约束条件化为等式，并且约束方程组的系数阵中有一个单位阵。 … 0 … 0 cj-zj a1n a2n . amn a1j a2j . amj 0 0 . 1 1 0 . 0 b1 b2 . bm x1 x2 . xm c1 c2 . cm xn … xj … xm … x1 b 基 CB cn … cj … cm … c1 cj 第二步：最优性检验 计算检验数，检查： ?所有检验数是否≤ 0？ 是——结束，写出最优解和目标函数最优值； ?还有正检验数——检查相应系数列≤ 0？ 是——结束，该LP无界解。 ?不属于上述两种情况，转入下一步—基变换。 确定是停止迭代还是转入基变换？ 选择（最大）正检验数 对应的系数列为主元列，主元列对应的非 基变量xk 为换入变量； 最小比值 对应的行为主元行，主元行对应的基变量 xl为换出变量。 主元行和主元列的交叉元素alk称为主元素。 第三步：基变换 利用矩阵的初等行变换