管理运筹学单纯形法的灵敏度分析与对偶精选文档.ppt
《管理运筹学单纯形法的灵敏度分析与对偶精选文档.ppt》由会员分享,可在线阅读,更多相关《管理运筹学单纯形法的灵敏度分析与对偶精选文档.ppt(50页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、管理运筹学单纯形法的灵敏度分析与对偶1本讲稿第一页,共五十页11单纯形表的灵敏度分析单纯形表的灵敏度分析一、目标函数中变量一、目标函数中变量Ck系数灵敏度分析系数灵敏度分析1.在最终的单纯形表里,在最终的单纯形表里,X k是非基变量是非基变量 由于约束方程系数增广矩阵在迭代中只是其本身的行的初等变换与由于约束方程系数增广矩阵在迭代中只是其本身的行的初等变换与Ck没有任何关系,没有任何关系,所以当所以当Ck变成变成Ck+Ck时,在最终单纯形表中其系数的增广矩阵不变,又因为时,在最终单纯形表中其系数的增广矩阵不变,又因为Xk是非是非基变量,所以基变量的目标函数的系数不变,即基变量,所以基变量的目标
2、函数的系数不变,即CB不变,可知不变,可知Zk也不变,只是也不变,只是Ck变变成了成了Ck+Ck。这时。这时 K=Ck-Zk就变成了就变成了Ck+Ck-Zk=K+Ck。要使原来的最优解。要使原来的最优解仍为最优解,只要仍为最优解,只要 K+Ck0即可,也就是即可,也就是Ck的增量的增量 Ck-K。2.在最终的单纯形表中,在最终的单纯形表中,X k是基变量是基变量 当当Ck变成变成Ck+Ck时,最终单纯形表中约束方程的增广矩阵不变,但是基变量的目时,最终单纯形表中约束方程的增广矩阵不变,但是基变量的目标函数的系数标函数的系数CB变了,则变了,则ZJ(J=1,2,.,N)一般也变了,不妨设一般也变
3、了,不妨设CB=(CB1,CB2。,Ck,,CBm),当当CB变成变成=(CB1,CB2。,Ck+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本讲稿第二页,共五十页11单纯形表的灵敏度分析单纯形表的灵敏度分析根据上式可知根据上式可知 检验数检验数 J(J=1,2,.,M)变成了变成了 J,有有 J=CJ-ZJ=J+CK aKj。要使最优解不变,只要当。要使最优解不变,只要当J K时,时,J 5000,250=500+7500500
4、0,所以原题的最优解不是本题的最优解。所以原题的最优解不是本题的最优解。在用电量的约束条件中加入松驰变量在用电量的约束条件中加入松驰变量S S4 4后得:后得:把这个约束条件加入到原最终单纯形表上,其中把这个约束条件加入到原最终单纯形表上,其中S S4 4为基变量,得表如下为基变量,得表如下:迭代次数基变量b比值501000000501010-1050000-2110501000100102500103000015000501005005002750000-500-5008本讲稿第八页,共五十页11单纯形表的灵敏度分析单纯形表的灵敏度分析 在上表中的在上表中的X X1 1,X X2 2不是单位
5、向量,故进行行的线性变换,得不是单位向量,故进行行的线性变换,得迭代次数基变量CBx1x2s1s2s3s4b比值501000000 x1501010-1050s2000-211050 x2100010010250s4000-100-201-3000zj501005005002750000-500-500把上表中的把上表中的S S4 4行的约束可以写为:行的约束可以写为:上式两边乘以(上式两边乘以(-1-1),再加上人工变量),再加上人工变量a a1 1得:得:用上式替换上表中的用上式替换上表中的S S4 4行,得下表:行,得下表:9本讲稿第九页,共五十页11单纯形表的灵敏度分析单纯形表的灵敏度
6、分析迭代次数基变量x1x2s1s2s3s4a1b比值501000000-Mx1501010-10050s2000-21(1)0050 x21000100100250a4-M00-100-20113000zj5010050-10M050-20M0-M0010M-50020M-5000 x15010-11000100s3000-2110050 x2100012-1000200a4-M0050-200-112000zj50100150-50M20M-500M-M050M-15050-20M0-M0 x1501003/50-1/501/50140s300001/51-2/502/50130 x2100
7、010-1/502/50-2/50120s40001-2/50-1/501/5040zj5010001003-300-100-3-M+310本讲稿第十页,共五十页11单纯形表的灵敏度分析单纯形表的灵敏度分析 由上表可知,最优解为:由上表可知,最优解为:即该工厂在添加了用电限量以后的最优生产计划为即该工厂在添加了用电限量以后的最优生产计划为产产品生产品生产140140件,件,产品生产产品生产120120件。件。11本讲稿第十一页,共五十页11单纯形表的灵敏度分析单纯形表的灵敏度分析三、约束方程中常数项的灵敏度分析三、约束方程中常数项的灵敏度分析 从上表我们可以发现各个松弛变量的值,正好等于相应变
8、量的对偶价格。在最从上表我们可以发现各个松弛变量的值,正好等于相应变量的对偶价格。在最优解中优解中S2=50是基变量,即为,原料是基变量,即为,原料A有有50千克没用完,再增加千克没用完,再增加A原料是不会增原料是不会增加利润的,加利润的,A的对偶价格为的对偶价格为0。对于任何为基变量的松弛变量所对应的约束条件的。对于任何为基变量的松弛变量所对应的约束条件的对偶价格为对偶价格为0。迭代次数基变量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 502750
9、0CJ-ZJ0 0 -50 0 -5012本讲稿第十二页,共五十页11单纯形表的灵敏度分析单纯形表的灵敏度分析 可以看出,上题中对于设备台时数约束来说,当其松弛变量在目标函数可以看出,上题中对于设备台时数约束来说,当其松弛变量在目标函数中从中从0变到变到Z3=50时,也就是只要当前余下一台时数设备从不能获利变成获利时,也就是只要当前余下一台时数设备从不能获利变成获利50元时,譬如有人愿意出元时,譬如有人愿意出50元买一个设备时,我们就不必为生产元买一个设备时,我们就不必为生产、产产品品而使用完所有的而使用完所有的设备设备台台时时了,了,这说这说明了明了设备设备台台时时数的数的对对偶价格就是偶价
10、格就是Z3=50元。元。对对于含有大于等于号的于含有大于等于号的约约束条件,添加剩余束条件,添加剩余变变量化量化为标为标准型。准型。这时这时这这个个约约束条件的束条件的对对偶价格就和偶价格就和这这个剩余个剩余变变量的量的 有关了。有关了。这这将使得最将使得最优优目目标值标值特特别别“恶恶化化”而不是改而不是改进进,故,故这时约这时约束条件的束条件的对对偶价格偶价格应应取取 值值的相反的相反数数-。对对于含有等于号的于含有等于号的约约束条件,其束条件,其约约束条件的束条件的对对偶价格就和偶价格就和该约该约束方束方程的人工程的人工变变量有关了。其量有关了。其约约束条件的束条件的对对偶价格就等于此偶
11、价格就等于此约约束方程的人工束方程的人工变变量的量的 值值。13本讲稿第十三页,共五十页下表下表给给出了一个由最出了一个由最终单纯终单纯形表形表对对于不同于不同约约束束类类型的型的对对偶价格的取偶价格的取值值。从从对对偶价格的定偶价格的定义义,可以知道当可以知道当对对偶价格偶价格为为正正时时它将改它将改进进目目标标函数的函数的值值,当,当对对偶价格偶价格为负时为负时它将使得目它将使得目标标函数朝着与最函数朝着与最优优化相反的方向前化相反的方向前进进。下面我下面我们们研究当右端研究当右端项项bj发发生生变变化化时时,在什么范,在什么范围围内其内其对对偶价格不偶价格不变变。由于由于bj的的变变化并
12、不影响系数矩化并不影响系数矩阵阵的迭代,故其最的迭代,故其最终单纯终单纯形表中的系数矩形表中的系数矩阵阵没有没有变变化。由化。由此可此可见见当当bj变变化化时时,要使原来的基不,要使原来的基不变变得到的基本可行解仍然是可行解,也就是所得到的基本可行解仍然是可行解,也就是所求的基求的基变变量的量的值值一定要大于一定要大于0。所。所谓谓使其使其对对偶价格不偶价格不变变的的bj的的变变化范化范围围,也就是使,也就是使其最其最优优解的所有基解的所有基变变量不量不变变,且所得的最,且所得的最优优解仍然是可行的解仍然是可行的bj的的变变化范化范围围。约束条件影子价格的取值 等于这个约束条件对应的松弛变量的
13、 值,即为 的相反数 等于这个约束条件对应的剩余变量的 值,即为 的相反数 =等于这个约束条件对应的人工变量的 值,即为 的相反数 11单纯形表的灵敏度分析单纯形表的灵敏度分析14本讲稿第十四页,共五十页11单纯形表的灵敏度分析单纯形表的灵敏度分析 当bj中的第k项bK 变成 时,也就是原来的初始单纯形表中的b向量变成了b向量15本讲稿第十五页,共五十页11单纯形表的灵敏度分析单纯形表的灵敏度分析 这样在最终单纯形表中基变量这样在最终单纯形表中基变量XB的解就变成了的解就变成了 如要使如要使XB成为可行解,只要使上述等式的右边成为可行解,只要使上述等式的右边0,就可求出,就可求出 的取值范围,
14、也就是使得第的取值范围,也就是使得第K个约束条件的对偶价格不变的个约束条件的对偶价格不变的bk的变化范围。的变化范围。,16本讲稿第十六页,共五十页11单纯形表的灵敏度分析单纯形表的灵敏度分析下面我们仍以第二章例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 -5017本讲稿第十七页,共五十页11单纯形表的灵敏度分析单纯形表的灵敏度分
15、析 我们对b1进行灵敏度分析,因为在第一个约束方程中含有松弛变量S1,实际意义可以描述为:当设备台时数的对偶价格不变,都为每设备台时数在250与325之间变化,则设备台时数的对偶价格不变,都为每台设备台时50元。18本讲稿第十八页,共五十页11单纯形表的灵敏度分析单纯形表的灵敏度分析三、约束方程系数矩阵A灵敏度分析下面分两种情况讨论 1.在初始单纯形表上的变量Xk的系数列Pk改变为Pk经过迭代后,在最终单纯形表上Xk是非基变量。由于单纯形表的迭代是约束方程的增广矩阵的行变换,Pk变成Pk仅仅影响最终单纯形表上第k列数据,包括Xk的系数列、Zk以及 k,这时最终单纯形表上的Xk的系数列就变成了B
16、-1Pj,而Zk就变成CBB-1Pk,新的检验数 k=Ck-CBB-1Pk。若 k0,则原最优解仍然为最优解。若 k 0,则继续进行迭代以求出最优。例例 以第二章例1为基础,设该厂除了生产,种产品外,现在试制成一个新产品,已知生产产品,每件需要设备2台时,并消耗A原料0.5公斤。B原料1.5公斤,获利150元,问该厂应该生产该产品多少?解:这是一个增加新变量的问题。我们可以把它认为是一个改变变量X3在初始表上的系数列的问题,19本讲稿第十九页,共五十页11单纯形表的灵敏度分析单纯形表的灵敏度分析接上页迭代次数基变量CBX1 X2 S1 S2 S3 X3 b50 100 0 0 0 150X15
17、01 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 -2520本讲稿第二十页,共五十页11单纯形表的灵敏度分析单纯形表的灵敏度分析例 假设上例题中产品的工艺结构有了改进,这时生产1件产品需要使用1.5台设备,消耗原料A为2千克,原料B为1千克,每件产品的利润为160元,问该厂的生产计划是否要修改。解:首先求出X3在最终表上的系数列 迭代次数基变量CBX1 X2 S1 S2 S3 X3 b50 100 0 0 0 1502X1501 0 1 0
18、 -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 3521本讲稿第二十一页,共五十页11单纯形表的灵敏度分析单纯形表的灵敏度分析接下来又可以有新的迭代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 160310
19、00CJ-ZJ-70 0 -120 0 20 022本讲稿第二十二页,共五十页11单纯形表的灵敏度分析单纯形表的灵敏度分析接上页接上页 可知此规模的最优解可知此规模的最优解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
20、 -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 023本讲稿第二十三页,共五十页11单纯形表的灵敏度分析单纯形表的灵敏度分析 2.在初始表上的变量在初始表上的变量XK的系数的系数PK改变为改变为PK,经,经过迭代后,在最终表上过迭代后,在最终表上XK是基变量,在这种情况下是基变量,在这种情况下原最优解的可行性和最优解都可能被破坏,问题十分原最优解的可行性和最优解都可能被破坏,问题十分复杂,一般不去修改原表而是直接计算。复杂,一般不去修改原表而是直接计算。24本讲稿第二
21、十四页,共五十页 每一个线性规划问题,都存在每一个与它密切相关的线性规划的问题,我们称其为原问题,另一个为对偶问题。每一个线性规划问题,都存在每一个与它密切相关的线性规划的问题,我们称其为原问题,另一个为对偶问题。例题例题1 1 某工厂在计划期内安排某工厂在计划期内安排、两种产品,生产单位产品所需设备两种产品,生产单位产品所需设备A A、B B、C C台时如表所示台时如表所示 该工厂每生产一单位产品该工厂每生产一单位产品 可获利可获利5050元,每生产一单位产品元,每生产一单位产品可获利可获利100100元,问工厂应分别生产多少元,问工厂应分别生产多少 产品和产品和产品,才能使工厂获利最多?产
22、品,才能使工厂获利最多?解:设解:设 为产品为产品 的计划产量,的计划产量,为产品为产品的计划产量,则有的计划产量,则有目标函数:目标函数:Max z=50+100 Max z=50+100约束条件:约束条件:,2 2 线性规划的对偶问题线性规划的对偶问题25本讲稿第二十五页,共五十页 现在我们从另一个角度来考虑这个问题。假如有另外一个工厂要求租用该厂的设备现在我们从另一个角度来考虑这个问题。假如有另外一个工厂要求租用该厂的设备A、B、C,那么该厂的厂长应该如何来确定合理的租金呢?,那么该厂的厂长应该如何来确定合理的租金呢?设设 分别为设备分别为设备A、B、C的每台时的租金。为了叙述方便,这里
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 管理 运筹学 单纯 灵敏度 分析 对偶 精选 文档
限制150内