灵敏度分析与对偶理论及应用举例.ppt

灵敏度分析 图解法下的灵敏度分析 定义：研究LP问题中目标函数系数C与约束条件中系数矩阵A与常数向量b的变化（局部或全部变动）对LP的最优解 与最优值 的影响程度分析的问题称为LP的灵敏度分析。 由于C，b，A通常表述产品需求，原材料价格，未来的商品售价及企业的加工能力、资源消耗，故在数学建模时的估计有可能出现误差，且未来的环境亦是不确定的，有可能出现变化，因此有必要对应用问题的解作灵敏度分析。 图解法下的灵敏度分析 灵敏度分析通常讨论如下内容： C，b，A变动（局部或全部）后LP最优解是否仍存在？若存在的话，最优解有多大程度的变化？ C，b，A在什么范围内变动时能保证LP的原最优解不变？ 若C，b，A的变动超出上述范围后，如何利用LP的原最优解 来求解变动后的 最优解。 这些C，b，A的变动，在出现情况下对决策者有利？在什么情况下对决策者不利？此称为对偶价格研究对象（C的变动对目标最优解的影响）。 各种问题求解 问题建模 求解算法（理论、算法步骤、计算复杂性） 解的存在性、唯一性、可靠性（误差分析、灵敏度分析） LPa（LP原规划） LPb（C的变动对LP最优解的影响） LPc（b1的变动对LP最优解的影响） LPd（b2的变动对LP最优解的影响） 命题Ⅰ：在例1 的优化模型中，A，b不变，讨论C变动的如下问题： 欲使LP的最优解不变，仍为 ，C1与C2应满足什么条件？ 当C2=100 不变，欲使LP 的最优解不变，仍为 ，C1应满足什么条件？ 当C1=50 不变，欲使LP 的最优解仍为 ,C2应满足什么条件？ 当C1，C2均变动，例如C1=60，C2=50时，最优解是否变动？ 欲使最优解变动为 ，求C1，C2应满足的条件？ 解： 由图(a)，原最优解 对应B点，设LE，LF，LG直线之斜率分别为kE，kF，kG，则由对应直线方程可得kE= -1，kF=0，kG= -2，kz= - G1/G2，此中kz为等值线z=50x1+100x2之斜率。由图(a)中知，欲使最优点仍为B，则应有 当C2=100时，由*式知，需要 方使最优点不变。 当C1=50时，由*式知，需要 方使最优点不变。 解： 当C1=60，C2=50时，由图(b)知最优点为C (100，200）,故最优解 且由*式知：此时显然不满足条件 欲使最优解变动为 即最优点为C时，由图(b)知应有 此时有最优值 z(100,200)=100C1+200C2。 问题Ⅱ：在例1的优化模型中，C不变，讨论b与A分别变动时的如下问题： b1=300?310，b2，b3，A不变时，求设备(台时）的对偶价格Pb1。 b2=400?410，b1，b3，A不变时，求原料(克）的对偶价格Pb2。 解： 。 。 计算机解法下的灵敏度分析(P102~P112) 注1 注2 注3 设G为R3中的凸集，它有8个顶点,每一个顶点对应一个基本可行解，当采用单纯形表迭代时，其基本原则是以一个初始顶点开始，然后沿相邻顶点行进(迭代)，直到终点(最优解)为止。 注意到在R3中的凸集D的每一个顶点的相邻顶点有三个（也可有五个如（b)等）凸集D中，每一顶点或基本可行解对应一个基矩阵，则相邻顶点对应的基矩阵只有一列不同），因此从一个起点B0(对应一个初始基本可行解或初始基)开始，沿相邻顶点行进（迭代），直到终点（对应最优解）B2的迭代或运行路径可以有多条，如B0?B1?B2，B0?B5?B4?B3?B2等，这多条路径尽管有相同的起点（初始可行解），相同的终点（最优解），但迭代步数可以不等。 例9 解1 解2 对上述LP的b1作灵敏度分析 解 结论：当 时，不仅LP与LP有相同最优解，见迭代过程亦相同。一般迭代过程可不同，但最终单纯形表结果相同，由于迭代过程为初等行变换，故变换后之方程与原方程同解，故不影响最终结果，此二数应在0点的两侧。 例9 在原例1中，该厂除生产Ⅰ，Ⅱ产品外，现又试制成一种新产品Ⅲ，并知生产单位产品Ⅲ需设备2台时，消耗原料A需0.5Kg，原料B需1.5Kg，此时可获利150元，试问该厂是否应该生产产品Ⅲ，或此时保持原方案是否合适？ 若该厂对产品Ⅲ的工艺