[高等教育]《管理运筹学》第2版6-10章教案.doc

第六章 单纯形法的灵敏度分析与对偶 §1 单纯形表的灵敏度分析 §2 线性规划的对偶问题 §3 对偶规划的基本性质 §4 对偶单纯形法 §1 单纯形表的灵敏度分析 一、目标函数中变量Ck系数灵敏度分析 1. 在最终的单纯形表里，X k是非基变量 由于约束方程系数增广矩阵在迭代中只是其本身的行的初等变换与Ck没有任何关系，所以当Ck变成Ck+ΔCk时，在最终单纯形表中其系数的增广矩阵不变，又因为Xk是非基变量，所以基变量的目标函数的系数不变，即CB不变，可知Zk也不变，只是Ck变成了Ck+ΔCk。这时δK= Ck-Zk就变成了Ck+ΔCk- Zk=δK +ΔCk。要使原来的最优解仍为最优解，只要δK+ ΔCk≤0即可，也就是Ck的增量ΔCk≤—δK。 2. 在最终的单纯形表中， X k是基变量 当Ck变成Ck+ ΔCk时，最终单纯形表中约束方程的增广矩阵不变，但是基变量的目标函数的系数CB变了，则ZJ(J=1,2,…,N)一般也变了，不妨设CB=(CB1, CB2…, Ck,…，CBm)，当CB变成=(CB1, CB2…,Ck+ ΔCk,…,CBm)，则： ZJ=(CB1, CB2…, Ck,…，CBm)(a’1j , a’2j ,…, a’Kj ,…, a’mj) Z’J=(CB1, CB2…, Ck+ΔCk,…，CBm)(a’1j , a’2j ,…, a’Kj ,…, a’mj) T = ZJ + ΔCk a’Kj 根据上式可知 检验数δJ (J=1,2,…..,M)变成了δ’J，有δ’ J=CJ-Z’J=δJ+ΔCK a’Kj 。要使最优解不变，只要当J≠K时，δ’J =0 例： 目标函数：Max Z=50X1+100X2 约束条件：X1+X2≤300 2X1+X2≤400 X2≤250 X1，X2≥0 最优单纯形表如下 表6—1 我们先对非基变量S1的目标函数的系数C3进行灵敏度分析。 这里δ3= -50，所以当c3的增量Δc3≤50，最优解不变。 再对基变量x1的目标函数的系数c1进行灵敏度分析。 在a11? ,a12? ,a13? ,a14? ,a15? 中，除了知道a11? 和 a13? 大于0， a15? 小于0，可知 ， 有 。同样有 。这样可以知道当-50≤Δc1≤ 50时，也就是x1的目标函数c1＇在0≤c1＇≤100时最优解不变。 在最终的单纯形表中，用C’1代替原来的C1=50，计算得表 表6—2 从δ3≤0，得到-c1’≤0,即c1’≥0，并且从δ5≤0，得到c1’≤100。 那么如果c1’取值超出这个范围，必然存在一个检验数大于0，我们可以通过迭代来得到新的最优解。 二、约束方程中常数项的灵敏度分析 表6—3 从上表我们可以发现各个松弛变量的值，正好等于相应变量的对偶价格。在最优解中S2 =50是基变量，即为，原料A有50千克没用完，再增加A原料是不会增加利润的，A的对偶价格为0。对于任何为基变量的松弛变量所对应的约束条件的对偶价格为0。 可以看出，上题中对于设备台时数约束来说，当其松弛变量在目标函数中从0变到Z3=50时，也就是只要当前余下一台时数设备从不能获利变成获利50元时，譬如有人愿意出50元买一个设备时，我们就不必为生产Ι、П产品而使用完所有的设备台时了，这说明了设备台时数的对偶价格就是Z3=50元。 对于含有大于等于号的约束条件，添加剩余变量化为标准型。这时这个约束条件的对偶价格就和这个剩余变量的 有关了。这将使得最优目标值特别“恶化”而不是改进，故这时约束条件的对偶价格应取 值的相反数- 。 对于含有等于号的约束条件，其约束条件的对偶价格就和该约束方程的人工变量有关了。其约束条件的对偶价格就等于此约束方程的人工变量的 值。 下表给出了一个由最终单纯形表对于不同约束类型的对偶价格的取值。 表6-4 从对偶价格的定义,可以知道当对偶价格为正时它将改进目标函数的值，当对偶价格为负时它将使得目标函数朝着与最优化相反的方向前进。 下面我们研究当右端项bj发生变化时，在什么范围内其对偶价