单纯形法解法的矩阵描述及灵敏度分析.pptVIP

迭代后，得： 如果某个系数aij发生变化，如何 进行分析？ 最终单纯形表中，xj对应的系数列 向量发生了变化，即Pj发生了变化。 增加一个约束条件的分析 增加一个约束条件相当于增添一道工序或一种资源。 分析方法： 将最优解代入新的约束中 （1）若满足要求，则原最优解不变； （2）若不满足要求，则原最优解改变， 将新增的约束条件添入最终的单纯形表 中继续分析。 灵敏度分析的步骤归纳如下： （1）将参数的改变计算反映到最终 单纯形表上； （2）检查原问题是否仍为可行解； （3）检查对偶问题是否仍为可行解； （4）按下表所列情况得出结论和决 定继续计算的步骤。 原问题 对偶问题 结论或继续计算的步骤 可行解 可行解 问题的最优解或最优基不变 可行解 非可行解 用单纯形法继续迭代 非可行解 可行解 用对偶单纯形法继续迭代 非可行解 非可行解 编制新的单纯形表重新计算 单纯形法解法的矩阵描述及灵敏度分析 张林刚 经济与管理学院 单纯形解法的矩阵描述 线性规划问题 引入松弛变量Xs,化为标准型： 用单纯形法求解以上模型，初始单纯形表为： 单纯形解法的矩阵描述 非基变量 基变量 当基变量为 时，新的单纯形表 基变量 非基变量 当 时，得到 最优解 比较两个单纯形表 非基变量 基变量 基变量 非基变量 【例】 加入松弛变量x4、x5，对上述模型进行标准化处理 初始单纯形表 6 -2 3 0 0 CB XB b x1 x2 x3 x4 x5 0 x4 2 2 -1 2 1 0 0 x5 4 1 0 4 0 1 6 -2 3 0 0 6 -2 3 0 0 CB XB b x1 x2 x3 x4 x5 6 x1 4 1 0 4 0 1 -2 x2 6 0 1 6 -1 2 0 0 -9 -2 -2 B矩阵 N矩阵 I单位阵 B的逆阵B-1 最终单纯形表 I单位阵 B-1 N 如何而来？ 知识点1 目标函数为max时，判断最优的准则为б≤0； 目标函数为min时，迭代过程与max一样，判断最优的准则为σ≥0。 性质6：线性规划的原问题与其对偶问题之间存在一对互补的基解；其中原问题的松弛变量对应对偶问题的变量，对偶问题的剩余变量对应原问题的变量；这些互相对应的变量如果在一个问题的解中是基变量，则在另一个问题的解中是非基变量；将这对互补的基解分别代入原问题和对偶问题的目标函数有z=ω。 知识点2 【例】线性规划问题为 2 3 0 0 0 CB XB b x1 x2 x3 x4 x5 2 x1 3 1 0 1/2 0 -1/5 0 x4 4 0 0 -2 1 4/5 3 x2 3 0 1 0 0 1/5 cj-zj 0 0 -1 0 -1/5 原问题变量 原问题松弛变量 对偶问题变量 对偶问题剩余变量 其中XB=B-1b为原问题的基解；而YS=CBB-1为对偶问题的基解。 基变量 非基变量 知识点3：单纯形法与对偶单纯形法的对比 单纯形法求解的基本思想：保持原问题为可行解，通过迭代，增大目标函数，当其对偶问题的解也为可行解时，就达到了目标函数的最优值。 对偶单纯形法的基本思想：保持对偶问题为可行解，通过迭代，减小目标函数，当原问题也达到可行解时，即得到了目标函数的最优值。 灵敏度分析 所谓灵敏度分析，就是考察当线性规划问题中的参数aij、bi、cj等发生变化时，线性规划问题最优解如何变化的问题。 某一个参数发生变化时，通过计算，原最终单纯形表会出现以下几种情况： 原问题 对偶问题 结论或继续计算的步骤 可行解 可行解 问题的最优解不变 可行解 非可行解 用单纯形法继续迭代 非可行解 可行解 用对偶单纯形法继续迭代 非可行解 非可行解 编制新的单纯形表重新计算 【例】某家电厂家利用现有资源生产两种产品，有关数据如下表： 原问题 最优解 对偶问题最优解 （相差负号） 最终单纯形表 分析 的变化 的变化仅影响 的变化。 设备A 设备B 调试工序 利润（元） 0 6 1 2 5 2 1 1 15时 24时 5时 Ⅰ Ⅱ D 1.5 2 问题1：当 该公司最优生 产计划有何变化？ 最终单纯形表变为： 对变化后的单纯形表继续迭代 新的最终单纯形表为 最优解 问题2：设产品II利润为 ， 求原最优解不变时 的范围。 的变化