第6章单纯形法的灵敏度分析与对偶研究报告.ppt

管 理 运 筹 学第六章 单纯形法的灵敏度分析与对偶 1 单纯形表的灵敏度分析 2 线性规划的对偶问题 3 对偶规划的基本性质 4 对偶单纯形法1管 理 运 筹 学1 单纯形表的灵敏度分析一、目标函数中变量Ck系数灵敏度分析1.在最终的单纯形表里，X k是非基变量 由于约束方程系数增广矩阵在迭代中只是其本身的行的初等变换与Ck没有任何关系，所以当Ck变成Ck+ Ck时，在最终单纯形表中其系数的增广矩阵不变，又因为Xk是非基变量，所以基变量的目标函数的系数不变，即CB不变，可知Zk也不变，只是Ck变成了Ck+ Ck这时 K= Ck-Zk就变成了Ck+ Ck- Zk= K+ Ck要使原来的最优解仍为最优解，只要 K+ Ck0即可，也就是Ck的增量 Ck- K2.在最终的单纯形表中， X k是基变量 当Ck变成Ck+ Ck时，最终单纯形表中约束方程的增广矩阵不变，但是基变量的目标函数的系数CB变了，则ZJ(J=1,2,.,N)一般也变了，不妨设CB=(CB1, CB2 Ck,， CBm),当CB变成=(CB1, CB2Ck+ Ck,CBm),则： ZJ=(CB1, CB2 Ck,，CBm)(a1j , a2j , aKj , amj) ZJ=(CB1, CB2。

Ck+ Ck,，CBm)(a1j , a2j , aKj , amj) = ZJ + Ck aKj 2管 理 运 筹 学1 单纯形表的灵敏度分析根据上式可知 检验数 J (J=1,2,.,M)变成了 J,有 J=CJ-ZJ= J+ CK aKj 要使最优解不变，只要当J K时， J 0，就可求出 的取值范围，也就是使得第K个约束条件的对偶价格不变的bk的变化范围 ，11管 理 运 筹 学1 单纯形表的灵敏度分析下面我们仍以第二章例1在最终单纯形表上对bj 进行灵敏度分析最终单纯形表如下所示：迭代次数基变量CBX1 X2 S1 S2 S3b50 100 0 0 02X1501 0 1 0 -150 S200 0 -2 1 150 X21000 1 0 0 1250 ZJ50 100 50 0 5027500CJ -ZJ0 0 -50 0 -5012管 理 运 筹 学1 单纯形表的灵敏度分析 我们对b1进行灵敏度分析，因为在第一个约束方程中含有松弛变量S1, 实际意义可以描述为：当设备台时数的对偶价格不变，都为每设备台时数在250与325之间变化，则设备台时数的对偶价格不变，都为每台设备台时50元。

13管 理 运 筹 学1 单纯形表的灵敏度分析三、约束方程系数矩阵A灵敏度分析下面分两种情况讨论 1.在初始单纯形表上的变量Xk的系数列Pk改变为Pk经过迭代后，在最终单纯形表上Xk是非基变量由于单纯形表的迭代是约束方程的增广矩阵的行变换，Pk变成Pk仅仅影响最终单纯形表上第k列数据，包括Xk的系数列、Zk以及 k，这时最终单纯形表上的Xk的系数列就变成了B-1Pj,而Zk就变成CBB-1Pk，新的检验数 k=Ck-CBB-1Pk若 k0，则原最优解仍然为最优解若 k 0，则继续进行迭代以求出最优 例 以第二章例1为基础，设该厂除了生产，种产品外，现在试制成一个新产品，已知生产产 品,每件需要设备 2台时，并消耗A原料0.5公斤B原料1.5公斤，获利150元，问该 厂应该 生产该产 品多少？解：这是一个增加新变量的问题 我们可以把它认为 是一个改变变 量X3在初始表上的系数列的问题 ，14管 理 运 筹 学1 单纯形表的灵敏度分析接上页迭代次数基变量CBX1 X2 S1 S2 S3 X3 b50 100 0 0 0 150X1501 0 1 0 -1 0.550 S200 0 -2 1 1 -250 X21000 1 0 0 1 1.5250 ZJ50 100 50 0 50 17527500CJ -ZJ0 0 -50 0 -50 -2515管 理 运 筹 学1 单纯形表的灵敏度分析例 假设上例题中产品的工艺结 构有了改进，这时 生产1件产品需要使用1.5台设备 ，消耗原料A为2千克，原料B为1千克，每件产品的利润为 160元，问该 厂的生产计 划是否要修改。

解：首先求出X3在最终表上的系数列 迭代次数基变量CBX1 X2 S1 S2 S3 X3 b50 100 0 0 0 1502X1501 0 1 0 -1 0.55050/0.5 S200 0 -2 1 1 050 X21000 1 0 0 1 1250250/1 ZJ50 100 50 0 50 12527500CJ -ZJ0 0 -50 0 -50 3516管 理 运 筹 学1 单纯形表的灵敏度分析接下来又可以有新的迭代S3进基，迭代次数基变量CBX1 X2 S1 S2 S3 X3 b50 100 0 0 0 1503X31602 0 2 0 -2 1100- S200 0 -2 1 1 05050/1 X2100-20 1 -2 0 3 0150250/3 ZJ120 100 120 0 -20 16031000CJ -ZJ-70 0 -120 0 20 017管 理 运 筹 学1 单纯形表的灵敏度分析接上页 可知此规模的最优解X1=0, X2=0, S1=0, S2=0, S3=50, X3=200,此时，最大目标函数为32000元也就是说，该厂的新的生产计划为不生产、产品，生产产品200件, 可获得最大利润32000元。

迭代次数基变量CBX1 X2 S1 S2 S3 X3 b50 100 0 0 0 1504X31602 0 2 0 -2 1200- S300 0 -2 1 1 05050/1 X2100-2 1 4 -3 0 00250/3 ZJ120 100 80 20 0 16032000CJ -ZJ-70 0 -80 -20 0 018管 理 运 筹 学1 单纯形表的灵敏度分析 2.在初始表上的变量XK的系数PK改变为PK，经过迭代后，在最终表上XK是基变量，在这种情况下原最优解的可行性和最优解都可能被破坏，问题十分复杂，一般不去修改原表而是直接计算19管 理 运 筹 学1 单纯形表的灵敏度分析四、增加一个约束条件的灵敏度分析 在原线性规划中增加一个约束条件时,先将原问题的最优解的变量值代入新增的约束条件，如满足则说明新增的条件没有起到限制作用，故原最优解不变，否则将新增的约束添入原最终单纯形表上进一步求解 下面仍以第三章例1为例来加以说明 例：假如该工厂除了在设备台时，原材料A、B上对该厂的生产有限制外，还有电力供应上的限制最高供应电量为5000度，而生产一个产品需要用电10度，而生产一个产品需要用电30度。

试分析此时该厂获得最大利润的生产计划？20管 理 运 筹 学1 单纯形表的灵敏度分析 解：先将原问题的最优解=50，=250代入用电量的约束条件得：1050+30250=500+75005000,所以原题的最优解不是本题的最优解在用电量的约束条件中加入松驰变量S4后得：把这个约束条件加入到原最终单纯形表上，其中S4为基变量，得表如下：迭代次数基变变量b比值值501000000501010-1050000-2110501000100102500103000015000501005005002750000-500-50021管 理 运 筹 学1 单纯形表的灵敏度分析 在上表中的X1，X2不是单位向量，故进行行的线性变换，得迭代次数基变变量CBx1x2s1s2s3s4b比值值501000000 x1501010-1050s2000-211050 x2100010010250s4000-100-201-3000zj501005005002750000-500-500把上表中的S4行的约束可以写为：上式两边乘以（-1），再加上人工变量a1得：用上式替换上表中的S4行，得下表：22管 理 运 筹 学1 单纯形表的灵敏度分析迭代次数基变变量x1x2s1s2s3s4a1b比值值501000000-Mx1501010-10050s2000-21(1)0050 x21000100100250s4-M00-100-20113000zj5010050-10M050-20M0-M0010M-50020M-5000 x15010-11000100s3000-2110050 x2100012-1000200s4-M0050-200-112000zj50100150-50M20M-500M-M050M-15050-20M0-M0 x1501003/50-1/501/50140s300001/51-2/502/50130 x2100010-1/502/50-2/50120s40001-2/50-1/501/5040zj5010001003-300-100-3-M+323管 理 运 筹 学1 单纯形表的灵敏度分析 由上表可知，最优解为： 即该工厂在添加了用电限量以后的最优生产计划为产品生产140件，产品生产120件。

24管 理 运 筹 学 每一个线性规划问题，都存在每一个与它密切相关的线性规划的问题，我们称其为原问题，另一个为对偶问题例题1 某工厂在计划期内安排、两种产品，生产单位产品所需设备A、B、C台时如表所示 该工厂每生产一单位产品 可获利50元，每生产一单位产品可获利100元，问工厂应分别生产多少 产品和产品，才能使工厂获利最多？解：设 为产品 的计划产量， 为产品的计划产量，则有目标函数： Max z=50 +100约束条件： ，2 线性规划的对偶问题25管 理 运 筹 学 现在我们从另一个角度来考虑这个问题假如有另外一个工厂要求租用该厂的设备A、B、C，那么该厂的厂长应该如何来确定合理的租金呢？ 设 分别为设备A、B、C的每台时的租金为了叙述方便，这里把租金定义为扣除成本后的利润作为出租者来说，把生产单位 产品所需各设备的台时各总租金不应低于原利润50元，即 ，否则就不出租还是用于生产 产品以获利50元；同样把 生产一单位 产品所需各设备的台时的总租金也不应当低于原利润100元， 即，否则这些设备台时就不出租，还是用于生产 产品以获利100元但对于租用者来说，他要求在满足上述要求的前提下，也就是在出租者愿意出租的前提下尽量要求全部设备台时的总租金越低越好，即min ，这样我们得到了该问题的数学模型： 目标函数： 约束条件： 这样从两个不同的角度来考虑同一个工厂的最大利润（最小租金）的问题，所建立起来的两个线性模型就是一对对偶问题，其中一个叫做原问题，而另外一个叫对偶问题。

2 线性规划的对偶问题26管 理 运 筹 学 如果我们把求目标函数最大值的线性规划问题看成原问题，则求目标函数最小值的线性规划问题看成对偶问题下面来研究这两个问题在数学模型上的关系 1 求目标函数最大值的线性规划问题中有n 个变量 m个约束条件，它的约束条件都是小于等于不等式而其对偶则是求目标函数为最小值的线性规划问题，有m个变量n个约束条件，其约束条件都为大于等于不等式 2 原问题的目标函数中的变量系数为对偶问题中的约束条件的右边常数项，并且原问题的目标函数中的第i个变量的系数就等于对偶问题中的第i个约束条件的右边常数项 3 原问题的约束条件的右边常数项为对偶问题的目标函数中的变量的系数并且原问题的第i个约束条件的右边常数项就等于零对偶问题的目标函数中的第i个变量的系数 4 对偶问题的约束条件的系数矩阵A是原问题约束矩阵的转置 设 A=则 2 线性规划的对偶问题27管 理 运 筹 学如果我们用矩阵形式来表示，则有原问题： 其中A是 矩阵m*n，该问题有m个约束条件n个变量，x= ,b= , c= 对偶问题： 其中 是A的转置， 是b的转置, 是c的转置, y= 现在我们用单纯形法求对偶问题的解。

2 线性规划的对偶问题28管 理 运 筹 学 加上剩余变量 和人工变量 ，把此问题化成标准型如下：把上述数据填入单纯形表计算2 线性规划的对偶问题29管 理 运 筹 学迭代。