【公务员考试】灵敏度分析.pptVIP

灵敏度分析 灵敏度分析是在建立数学模型和求得最优解之后，研究线性规划的一些系数Cj、aij、bi变化时，对最优解产生什么影响？ 当线性规划问题中的一个或几个参数变化时．可以用单纯形法从头计算，看最优解有无变化．但这样做既麻烦又没有必要。因为单纯形法的迭代计算是从一组基向量变换为另一组基向量，表中每步迭代得到的数字只随基向量的不同选择而改变，因此可把个别参数的变化直接在计算得到最优解的单纯形表上反映出来。这样就不需要从头计算，而直接对计算得到最优解的单纯形表进行审查，看一些数字变化后，是否仍满足最优解的条件，如果不满足的话，再从这个表开始进行迭代计算，求得最优解。 1. 灵敏度分析的基本原理 maxZ CX; AX b; LP s.t. X≥0. 对于选定的基B，不妨设 B P 1, P 2,…, P m . 将问题（LP）化为关于基B的典式 maxZ CBB-1b + CN — CBB-1N XN ; XB + B-1NXN B-1b ; s.t. XB, XN ≥0. 列出单纯形表 灵敏度分析是： 1 为了保持现有的最优解或最优基不变，找出这些数据变化的范围．即所谓数据的稳定性区间． 2 当这些数据的变化超出了 1 的范围时，如何在原有最优解或最优基的基础上，作微小的调整，尽快求出新的最优解或最优基． 由表可以看出，某些数据只和表中的某些块有关，因而当这些数据发生变化时，只需对上表中的某些块进行修改，便可得到新问题的单纯形表，从而能够进行判别和迭代，而不必从头开始计算线性规划问题，这正是单纯形法的优点之一． 下面分别讨论这些变化对最优解或最优解的影响． 一、目标函数系数的灵敏度 1. 非基变量xj的价值系数cj的变化 cj改变为c’j cj +Δcj 变化后的检验数为 σ‘j cj +Δcj — CBB-1Pj. 要保持原有最优解不变，则必须 σ’j cj + Δcj — CBB-1Pj ≤0 即 cj — CBB-1Pj cj + Δcj σj + Δcj ≤0 所以即： Δcj ≤ —σj　　　 2.基变量xj的价值系数cj的变化 某基变量xr的价值系数cr 改变为 c’r cr + Δcj，因cr ∈CB，则 CB +ΔCB B-1A CB B-1A+ 0,…, Δcr,…,0 B-1A CBB-1A + Δcr a‘r1, a‘r2,…, a‘r n 其中 a‘r1, a‘r2,…, a‘r n 是矩阵B-1A的第r行 变化后的检验数为： σ’j cj — CBB-1Pj —Δcra‘rj σj —Δcra‘rj （j 1，2，…，n）。 若要最优解不变，则必须满足，即 σ’j σj —Δcra‘rj ≤0 （ j 1，2，…，n）。 由此导出 当a‘rj 0时，有Δcr≤σj / a‘rj； 当a‘rj 0时，有Δcr≥σj / a‘rj； 当a‘rj 0时，有Δcr≤σj / a‘rj； 当a‘rj 0时，有Δcr≥σj / a‘rj； 因此，Δc允许变化范围是 尽管我们可以推导确定cj可变范围上下限的一些公式，但不需要去记忆此类公式。 例1 　　　 求其中目标函数系数C2的可变范围。 解：原问题加入松弛变量X3 ， X4 ， X5 ，化为标准型后，可将该问题的初始单纯形表列出如表示。 经过求解的最后单纯形表如下。 在最优单纯形表中，x2的目标函数系数c2的当前值为120。现把这个系数的值用符号c2代替。此时，相应的单纯形表如表所示。 为使问题的最优解不变，即c2要满足不等式组： 　　　 -28+0.12c2≤0 　　　 14-0.16c2≤0 由此可得到c2的可变范围： 87.5≤c2≤233.33 或 87.5≤120＋Δc2≤233.33 原c2＝120 即 87.5－120≤Δc2≤233.33－120 -32.5≤Δc2≤113.33 其他