运筹学05-灵敏度分析课件.ppt

第五章 线性规划灵敏度分析 5.1 目标函数系数的灵敏度分析 5.2 右端项的灵敏度分析 5.3 约束系数的灵敏度分析 5.4 参数规划 5.1 目标函数系数的灵敏度分析 解：最优单纯形表 (2) 若 ck 是基变量的系数 例 5.2 右端项的灵敏度分析 设分量 br 变化为 br + ?br ，根据前面的讨论： 最优解的基变量 xB = B-1b，那么只要保持 B-1(b + ?b) ≥ 0 ， 则最优基不变，即基变量保持不变，只有值的变化； 否则，需要利用对偶单纯形法继续计算。 例 5.3 约束系数的灵敏度分析 (1) 增加一个变量 增加一个变量，相当于系数矩阵增加一列。 增加变量 xn+1 则有相应的pn+1 ,cn+1 。 那么计算出B-1pn+1 , ?n+1=cn+1-cBB-1 pn+1 填入最优单纯形表, 若 ?n+1 ≤ 0 则 最优解不变； 否则，进一步用单纯形法求解。 例:前例增加x6 , p6=( 2, 6, 3 )T, c6=5 计算得到 (2) 增加一个约束条件 增加一个约束条件相当于系数矩阵中增加一行。 增加一个约束条件之后，应把最优解带入新的约束，若满足则最优解不变，否则填入最优单纯形表作为新的一行，引入一个新的非负变量（原约束若是小于等于形式可引入非负松弛变量，否则引入非负人工变量），并通过矩阵行变换把对应基变量的元素变为0，进一步用单纯形法或对偶单纯形法求解。 例:前例增加3x1+ 2x2≤15，原最优解不满足这个约束。于是 (3) 技术系数改变 （计划生产的产品工艺结构改变） ① 非基变量xj工艺改变 只影响单纯形表Pj 列， ? j . 关键看? j ? 0? 还是0? . 用增加新变量类似方法解决。 ② 基变量xj工艺改变，重新开始单纯型表的整个过程，复杂，在此暂不予讨论。 5.4 参数线性规划 在线性规划的实际应用中，由于某种原因，线性规划问题的目标函数的价值系数C和约束条件的右端常数b会随着某个参数而连续变动。 当数据随着某个参数连续变化时，研究其对最优解的影响，即为参数线性规划问题。 目标函数的价值系数含有参数的线性规划问题 右端常数含有参数的线性规划问题 * * 上表中6个常数 a1 , a2 , a3 , b , ?1 , ?2 取值在什么范围可使 1、现可行解最优，且唯一？何时不唯一？ 2、现基本解不可行； 3、问题无可行解； 4、无有限最优解； 5、现基本解可行，由 x1 取代 x6 目标函数可改善。 C ? ? ? ? ? ? CB XB b x1 x2 x3 x4 x5 x6 ? x3 b 4 a1 1 0 a2 0 ? x4 2 -1 -5 0 1 -1 0 ? x6 3 a3 -3 0 0 -4 1 ? ? ? σ1 σ2 0 0 -3 0 线性规划标准形式 (1)、参数A，b，C在什么范围内变动，对当前方案无影响？ (2)、参数A，b，C中的一个(几个)变动，对当前方案影响？ (3)、如果最优方案改变，如何用简便方法求新方案？ 当线性规划问题中的一个或几个参数变化时，可以用单纯形法从头计算，看最优解有无变化，但这样做既麻烦又没有必要。 灵敏度分析一词的含义是指对系统或事物因周围条件变化显示出来的敏感程度的分析。 考虑检验数 (1) 若ck是非基变量的系数： 例 试求 c3 在多大范围内变动时，原最优解保持不变。 CI -2 -3 -4 0 0 CB XB b x1 x2 x3 x4 x5 -3 X2 2/5 0 1 -1/5 -2/5 1/5 -2 X1 11/5 1 0 7/5 -1/5 -2/5 -z 0 0 -9/5 -8/5 -1/5 28/5 从表中看到σ3= -9/5+Δc3 可得到Δc3 ≤ 9/5 时，原最优解不变。 Ci -2 -3 -4+Δc3 0 0 CB XB b x1 x2 x3 x4 x5 -3 X2 2/5 0 1 -1/5 -2/5 1/5 -2 X1 11/5 1 0 7/5 -1/5 -2/5 -z 0 0 -9/5+Δc3 -8/5 -1/5 28/5 推导过程见教材P98 在单纯型表中： 用检验数行中的因素除以ck对应的行的元素 取大于零的值中的最小值为ck 增量的上界，小于零的值中的最大值为ck 增量的下界 经验方法 求c2在什么范围内变动时，原最优解保持不变。 从表中可得到 -3≤Δc2≤1 时，原最优解不变。 例: 下表为最优单纯形表,考虑基变量系数c2发生变