单纯形法的灵敏度分析与对偶

上传人:夏** 文档编号:589201276 上传时间:2024-09-10 格式:PPT 页数:66 大小:934KB
返回 下载 相关 举报
单纯形法的灵敏度分析与对偶_第1页
第1页 / 共66页
单纯形法的灵敏度分析与对偶_第2页
第2页 / 共66页
单纯形法的灵敏度分析与对偶_第3页
第3页 / 共66页
单纯形法的灵敏度分析与对偶_第4页
第4页 / 共66页
单纯形法的灵敏度分析与对偶_第5页
第5页 / 共66页
点击查看更多>>
资源描述

《单纯形法的灵敏度分析与对偶》由会员分享,可在线阅读,更多相关《单纯形法的灵敏度分析与对偶(66页珍藏版)》请在金锄头文库上搜索。

1、管管 理理 运运 筹筹 学学第六章第六章 单纯形法的灵敏度分析与对偶单纯形法的灵敏度分析与对偶1 1 单纯形表的灵敏度分析单纯形表的灵敏度分析2 2 线性规划的对偶问题线性规划的对偶问题3 3 对偶规划的基本性质对偶规划的基本性质4 4 对偶单纯形法对偶单纯形法1管管 理理 运运 筹筹 学学1 1 单纯形表的灵敏度分析单纯形表的灵敏度分析一、目标函数中变量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)一般也变了,不妨设CB=(CB1, CB2。, Ck,, CBm),当CB变成=(CB1, CB2。,Ck+ Ck,CBm),则: ZJ=(CB1, CB2。, Ck,,CBm)(a1j , a2j , aKj

3、 , amj) ZJ=(CB1, CB2。, Ck+ Ck,,CBm)(a1j , a2j , aKj , amj) = ZJ + Ck aKj 2管管 理理 运运 筹筹 学学 CB XB cjCBCN xj b XBTXNTCBTXBB-1b B-1BB-1N- CB B-1b CB- CB B-1BCN- CB B-1N若B是最优基,则最优表形式如下灵敏度分析总是在最优表上进行 3管管 理理 运运 筹筹 学学灵敏度分析的步骤灵敏度分析的步骤将参数的改变计算反映到最终表上来。具体计将参数的改变计算反映到最终表上来。具体计算公式可以使用算公式可以使用检查原问题是否仍为可行解检查原问题是否仍为可

4、行解检查对偶问题是否仍为可行解检查对偶问题是否仍为可行解对检查情况按下表进行处理对检查情况按下表进行处理4管管 理理 运运 筹筹 学学原问题原问题对偶问题对偶问题结论或继续计算步骤结论或继续计算步骤可行解可行解可行解可行解问题的最优解或最优基不变问题的最优解或最优基不变可行解可行解非可行解非可行解用单纯形法继续迭代求最优用单纯形法继续迭代求最优解解非可行解非可行解可行解可行解用对偶单纯形法继续迭代求用对偶单纯形法继续迭代求最优解最优解非可行解非可行解非可行解非可行解引进人工变量,编制新的单引进人工变量,编制新的单纯形表重新计算纯形表重新计算5管管 理理 运运 筹筹 学学1 1 单纯形表的灵敏度

5、分析单纯形表的灵敏度分析根据上式可知 检验数 J (J=1,2,.,M)变成了 J,有 J=CJ-ZJ= J+ CK aKj 。要使最优解不变,只要当J K时, J 0,就可求出 的取值范围,也就是使得第K个约束条件的对偶价格不变的bk的变化范围。 ,14管管 理理 运运 筹筹 学学1 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

6、 100 50 0 5027500CJ -ZJ0 0 -50 0 -5015管管 理理 运运 筹筹 学学1 1 单纯形表的灵敏度分析单纯形表的灵敏度分析 我们对b1进行灵敏度分析,因为在第一个约束方程中含有松弛变量S1, 实际意义可以描述为:当设备台时数的对偶价格不变,都为每设备台时数在250与325之间变化,则设备台时数的对偶价格不变,都为每台设备台时50元。16管管 理理 运运 筹筹 学学1 1 单纯形表的灵敏度分析单纯形表的灵敏度分析三、约束方程系数矩阵A灵敏度分析下面分两种情况讨论 1.在初始单纯形表上的变量Xk的系数列Pk改变为Pk经过迭代后,在最终单纯形表上Xk是非基变量。由于单纯

7、形表的迭代是约束方程的增广矩阵的行变换,Pk变成Pk仅仅影响最终单纯形表上第k列数据,包括Xk的系数列、Zk以及 k,这时最终单纯形表上的Xk的系数列就变成了B-1Pj,而Zk就变成CBB-1Pk,新的检验数 k=Ck-CBB-1Pk。若 k0,则原最优解仍然为最优解。若 k 0,则继续进行迭代以求出最优。 例例 以第二章例1为基础,设该厂除了生产,种产品外,现在试制成一个新产品,已知生产产品,每件需要设备2台时,并消耗A原料公斤。B原料公斤,获利150元,问该厂应该生产该产品多少?解:这是一个增加新变量的问题。我们可以把它认为是一个改变变量X3在初始表上的系数列的问题,17管管 理理 运运

8、筹筹 学学1 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 -2518管管 理理 运运 筹筹 学学1 1 单纯形表的灵敏度分析单纯形表的灵敏度分析例 假设上例题中产品的工艺结构有了改进,这时生产1件产品需要使用台设备 ,消耗原料A为2千克,原料B为1千克,每件产品的利润为160元,问该厂的生产

9、计划是否要修改。 解:首先求出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 3519管管 理理 运运 筹筹 学学1 1 单纯形表的灵敏度分析单纯形表的灵敏度分析接下来又可以有新的迭代S3进基,迭代次数基变量CBX1 X2 S1 S2 S3 X3 b50 100 0 0 0 1503X31602 0 2

10、 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 020管管 理理 运运 筹筹 学学1 1 单纯形表的灵敏度分析单纯形表的灵敏度分析接上页 可知此规模的最优解X1=0, X2=0, S1=0, S2=0, S3=50, X3=200,此时,最大目标函数为32000元。也就是说,该厂的新的生产计划为不生产、产品,生产产品200件, 可获得最大利润32000元。迭代次数基变量CBX1 X2 S1 S2 S3 X3 b50 10

11、0 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 021管管 理理 运运 筹筹 学学1 1 单纯形表的灵敏度分析单纯形表的灵敏度分析 2.在初始表上的变量XK的系数PK改变为PK,经过迭代后,在最终表上XK是基变量,在这种情况下原最优解的可行性和最优解都可能被破坏,问题十分复杂,一般不去修改原表而是直接计算。22管管 理理 运运 筹筹 学学1 1 单纯形表的灵敏度分析单纯形表的灵敏度分析四

12、、增加一个约束条件的灵敏度分析 在原线性规划中增加一个约束条件时,先将原问题的最优解的变量值代入新增的约束条件,如满足则说明新增的条件没有起到限制作用,故原最优解不变,否则将新增的约束添入原最终单纯形表上进一步求解。 下面仍以第三章例1为例来加以说明。 例:假如该工厂除了在设备台时,原材料A、B上对该厂的生产有限制外,还有电力供应上的限制。最高供应电量为5000度,而生产一个产品需要用电10度,而生产一个产品需要用电30度。试分析此时该厂获得最大利润的生产计划?23管管 理理 运运 筹筹 学学1 1 单纯形表的灵敏度分析单纯形表的灵敏度分析 解:先将原问题的最优解=50=50,=250代入用电

13、量的约束条件得:1050+30250=500+75005000,所以原题的最优解不是本题的最优解。在用电量的约束条件中加入松驰变量S4后得:把这个约束条件加入到原最终单纯形表上,其中S4为基变量,得表如下:迭代迭代次数次数基变量基变量b b比值比值50501001000 00 00 00 050501 10 01 10 0-1-10 050500 00 00 0-2-21 11 10 050501001000 01 10 00 01 10 02502500 0101030300 00 00 01 150005000505010010050500 050500 0275027500 00 00

14、0-50-500 0-50-500 024管管 理理 运运 筹筹 学学1 1 单纯形表的灵敏度分析单纯形表的灵敏度分析 在上表中的X1,X2不是单位向量,故进行行的线性变换,得迭代迭代次数次数基变量基变量C CB Bx x1 1x x2 2s s1 1s s2 2s s3 3s s4 4b b比比值值50501001000 00 00 00 0x x1 150501 10 01 10 0-1-10 05050s s2 20 00 00 0-2-21 11 10 05050x x2 21001000 01 10 00 01 10 0250250s s4 40 00 00 0-10-100 0-2

15、0-201 1-3000-3000z zj j505010010050500 050500 027500275000 00 0-50-500 0-50-500 0把上表中的S4行的约束可以写为:上式两边乘以(-1),再加上人工变量a1得:用上式替换上表中的S4行,得下表:25管管 理理 运运 筹筹 学学1 1 单纯形表的灵敏度分析单纯形表的灵敏度分析迭代迭代次数次数基基变量量x x1 1x x2 2s s1 1s s2 2s s3 3s s4 4a a1 1b b比比值50501001000 00 00 00 0-M-Mx x1 150501 10 01 10 0-1-10 00 05050s

16、 s2 20 00 00 0-2-21 1(1)(1)0 00 05050x x2 21001000 01 10 00 01 10 00 0250250s s4 4-M-M0 00 0-10-100 0-20-201 11 130003000z zj j505010010050-10M50-10M0 050-20M50-20M0 0-M-M0 00 010M-5010M-500 020M-5020M-500 00 0x x1 150501 10 0-1-11 10 00 00 0100100s s3 30 00 00 0-2-21 11 10 00 05050x x2 21001000 01

17、12 2-1-10 00 00 0200200s s4 4-M-M0 00 05050-20-200 0-1-11 120002000z zj j5050100100150-50M150-50M20M-5020M-500 0M M-M-M0 050M-15050M-15050-20M50-20M0 0-M-M0 0x x1 150501 10 00 03/53/50 0-1/50-1/501/501/50140140s s3 30 00 00 00 01/51/51 1-2/50-2/502/502/50130130x x2 21001000 01 10 0-1/5-1/50 02/502/5

18、0-2/50-2/50120120s s4 40 00 00 01 1-2/5-2/50 0-1/50-1/501/501/504040z zj j50501001000 010100 03 3-3-30 00 0-10-100 0-3-3-M+3-M+326管管 理理 运运 筹筹 学学1 1 单纯形表的灵敏度分析单纯形表的灵敏度分析 由上表可知,最优解为: 即该工厂在添加了用电限量以后的最优生产计划为产品生产140件,产品生产120件。27管管 理理 运运 筹筹 学学 每一个线性规划问题,都存在每一个与它密切相关的线性规划的问题,我们称其为原问题,另一个为对偶问题。例题例题1 1 某工厂在计

19、划期内安排、两种产品,生产单位产品所需设备A、B、C台时如表所示 该工厂每生产一单位产品 可获利50元,每生产一单位产品可获利100元,问工厂应分别生产多少 产品和产品,才能使工厂获利最多?解:设 为产品 的计划产量, 为产品的计划产量,则有目标函数: Max z z=50 +100约束条件: ,2 2 线性规划的对偶问题线性规划的对偶问题28管管 理理 运运 筹筹 学学 现在我们从另一个角度来考虑这个问题。假如有另外一个工厂要求租用该厂的设备A、B、C,那么该厂的厂长应该如何来确定合理的租金呢? 设 分别为设备A、B、C的每台时的租金。为了叙述方便,这里把租金定义为扣除成本后的利润。作为出租

20、者来说,把生产单位 产品所需各设备的台时各总租金不应低于原利润50元,即 ,否则就不出租还是用于生产 产品以获利50元;同样把 生产一单位 产品所需各设备的台时的总租金也不应当低于原利润100元, 即,否则这些设备台时就不出租,还是用于生产 产品以获利100元。但对于租用者来说,他要求在满足上述要求的前提下,也就是在出租者愿意出租的前提下尽量要求全部设备台时的总租金越低越好,即min ,这样我们得到了该问题的数学模型: 目标函数: 约束条件: 这样从两个不同的角度来考虑同一个工厂的最大利润(最小租金)的问题,所建立起来的两个线性模型就是一对对偶问题,其中一个叫做原问题原问题,而另外一个叫对偶问

21、题对偶问题。2 2 线性规划的对偶问题线性规划的对偶问题29管管 理理 运运 筹筹 学学 如果我们把求目标函数最大值的线性规划问题看成原问题,则求目标函数最小值的线性规划问题看成对偶问题。下面来研究这两个问题在数学模型上的关系。 1 求目标函数最大值的线性规划问题中有n 个变量 m个约束条件,它的约束条件都是小于等于不等式。而其对偶则是求目标函数为最小值的线性规划问题,有m个变量n个约束条件,其约束条件都为大于等于不等式。 2 原问题的目标函数中的变量系数为对偶问题中的约束条件的右边常数项,并且原问题的目标函数中的第i个变量的系数就等于对偶问题中的第i个约束条件的右边常数项。 3 原问题的约束

22、条件的右边常数项为对偶问题的目标函数中的变量的系数。并且原问题的第i个约束条件的右边常数项就等于零对偶问题的目标函数中的第i个变量的系数。 4 对偶问题的约束条件的系数矩阵A是原问题约束矩阵的转置。 设 A=则 2 2 线性规划的对偶问题线性规划的对偶问题30管管 理理 运运 筹筹 学学如果我们用矩阵形式来表示,则有原问题: 其中A是 矩阵m*n,该问题有m个约束条件n个变量,x= ,b= , c= 对偶问题: 其中 是A的转置, 是b的转置, 是c的转置, y= 2 2 线性规划的对偶问题线性规划的对偶问题31管管 理理 运运 筹筹 学学规范形式的线性规划问题规范形式的线性规划问题原问题(L

23、P)对偶问题(DLP) 32管管 理理 运运 筹筹 学学非规范形式线性规划的对偶问题1 变量取值范围不符合非负要求的情况33管管 理理 运运 筹筹 学学非规范形式线性规划的对偶问题2 约束方程不是“”的情况 34管管 理理 运运 筹筹 学学求线性规划的对偶规划规律总结方程对变量,变量对方程;正常对正常,不正常对不正常;变量正常是非负,方程正常看目标(max ,min )。 35管管 理理 运运 筹筹 学学例 求解下面线性规划的对偶规划LPDLP36管管 理理 运运 筹筹 学学 加上剩余变量 和人工变量 ,把此问题化成标准型如下:把上述数据填入单纯形表计算。2 2 线性规划的对偶问题线性规划的对

24、偶问题现在我们用单纯形法求对偶问题的解。37管管 理理 运运 筹筹 学学迭代变量基变量 b-300-400-25000-M 1-M1 0 -1 0 15050/2-250 1 1 1 0 -1 0 100100/1-M-250-2M-250-250M250-M-50M-25000M-2502M-1500-M-25002-4001/210-1/201/225-2501/2011/2-1-1/275-325-400-25075250-75-287502500-75-250-M+753-300120-10150-2500-111-1-150-300-350-25050250-50-275000-50

25、0-50-250-M+502 2 线性规划的对偶问题线性规划的对偶问题38管管 理理 运运 筹筹 学学 由上表,最优解: =50, -f 的最大值为-27500,即目标函数f的最大值为f=27500元。 从上面可知租金:A设备为50元,B设备为0元,C设备为50元。这样把工厂的所有设备出租可共得租金27500元。对出租者来说这租金是出租者愿意出租设备的最小费用,因为这是目 标函数的最小值。 通过比较,我们发现:对偶问题的最优解即最佳租金恰好等于原问题各种设备的对偶价格,这在道理上也能讲得通。 对于两个有对偶关系的线性规划的问题,我们只要求得了其中一个最优解,就可以从这个问题的对偶价格而求得其对

26、偶问题的最优解,知道其中一个最优值也就找到了其对偶问题的最优值,因为这两个最优值相等。 2 2 线性规划的对偶问题线性规划的对偶问题39管管 理理 运运 筹筹 学学 下面来阐述如何写出一个线性规划问题的对偶问题。为了便于阐述,我们不妨以下面的线性规划为例,写出它的对偶问题。 S.T. 2 2 线性规划的对偶问题线性规划的对偶问题40管管 理理 运运 筹筹 学学 这是一个求最大值的线性规划问题,为了写出它的对偶问题,我们不妨把它的约束条件都变换成取小于号的不等式。显然第一个约束条件已符合要求,不要做任何变动,而第二个约束条件,我们只要两边都乘以(-1),使不等号方向改变即可,得 这样第二个约束条

27、件也就符合要求。对于第三个约束条件,我们可以用小于等于和大于等于两个约束条件来替代它。即有 显然,这两个约束条件与原来第三个约束条件是等价的,我们再把其中的两边都乘以(-1),得 2 2 线性规划的对偶问题线性规划的对偶问题41管管 理理 运运 筹筹 学学 通过上面的一些变换,我们得到了一个和原线性规划等价的线性规划问题: s.t. 2 2 线性规划的对偶问题线性规划的对偶问题42管管 理理 运运 筹筹 学学 这个求最大值的线性规划问题的约束条件都取小于等于号,我们马上可以写出其对偶问题: .2 2 线性规划的对偶问题线性规划的对偶问题43管管 理理 运运 筹筹 学学 这里 和 一样都是不同的

28、决策变量,为了表示这两个决策变量都来源于原问题的第三个约束条件,记为 。 因为在该对偶问题中 和 的系数只相差一个符号,我们可以把上面的对偶问题化为: s.t.2 2 线性规划的对偶问题线性规划的对偶问题44管管 理理 运运 筹筹 学学 进一步,我们可以令 ,这时当 时, ,当 时, 。这也就是说,尽管 但 的取值可以为正,可以为0,可以为负,即 没有非负限制。 这样我们把原规划的对偶问题化为 s.t. 没有限制。 对照原线性规划问题,我们可以知道: 当原线性规划问题的第i个约束条件取等号时,则其对偶问题的 i个决策变量没有非负限制。 如果当原线性规划问题中的第 i个决策变量 没有非负限制时,

29、我们也可以用 进行替换,这里 , ,用类似的方法知道其对偶问题中第 i个约束条件取等号。2 2 线性规划的对偶问题线性规划的对偶问题45管管 理理 运运 筹筹 学学 另外,用大于等于0的两个决策变量之差来代替无非负限制的决策变量也是求解含有无非负限制的决策变量的线性规划问题的一种方法。 原线性规划问题为: s.t. 无非负限制。2 2 线性规划的对偶问题线性规划的对偶问题46管管 理理 运运 筹筹 学学3 3 对偶规划的基本性质对偶规划的基本性质对偶规划的基本性质对偶规划的基本性质1对称性对称性。即对偶问题的对偶是原问题。2弱对偶性弱对偶性。即对于原问题()和对偶问题()的可行解 都有C bT

30、 。 由弱对偶性,可得出以下推论:(1)原问题任一可行解的目标函数值是其对偶问题目标函数值的下界;反之对偶问题任一可行解的目标函数值是其原问题目标函数值的上界。(2)如原问题有可行解且目标函数值无界(或具有无界解),则其对偶问题无可行解;反之对偶问题有可行解且目标函数值无界,则其原问题无可行解(注意:本点性质的逆不成立,当对偶问题无可行解时,其原问题或具有无界解或无可行解,反之亦然)。(3)若原问题有可行解而其对偶问题无可行解,则原问题目标函数值无界;反之对偶问题有可行解而其原问题无可行解,则对偶问题的目标函数值无界。47管管 理理 运运 筹筹 学学对称性定理:对称性定理:线性规划的对偶问题的

31、对偶问题是原问题。线性规划的对偶问题的对偶问题是原问题。证明:证明: 对偶定义令w=-w;约束方程左右同乘“-1”对偶定义令z=-z;约束方程左右同乘“-1” 48管管 理理 运运 筹筹 学学弱对偶性定理:弱对偶性定理: 如果X、Y分别是原问题和对偶问题的一个可行解,则其对应的原问题的目标函数值不大于对偶问题的目标函数值,也即证明:证明:因为X、Y分别是原问题()与对偶问题()的可行解,故: 所以49管管 理理 运运 筹筹 学学3 3 对偶规划的基本性质对偶规划的基本性质3最优性最优性。如果 是原问题()的可行解, 是对偶问题()的可行解,并且 C = bT ,则 和 分别为原问题()和对偶问

32、题()的最优解。4强对偶性强对偶性。即若原问题()及其对偶问题()都有可行解,则两者都有最优解;且它们的最优解的目标函数都相等。5互补松弛性互补松弛性。在线性规划问题的最优解中,如果对应某一约束条件的对偶变量值为非零,则该约束条件取严格等式;反之,如果约束条件取严格不等式,则其对应的对偶变量一定为零也即 若yi*0,则有 若 ,则有yi*=050管管 理理 运运 筹筹 学学互补松弛定理51管管 理理 运运 筹筹 学学性质性质 设原问题是设原问题是 Max Z=C X AX+Xs=b X ,Xs 0 其对偶问题为其对偶问题为 Min w=Y b YA-Ys=C Y,Ys 0 则原问题单纯形表的则

33、原问题单纯形表的松弛变量的检验数对应对偶问题的一个解松弛变量的检验数对应对偶问题的一个解证明:证明: MaxZ=CX 非基变量的检验数非基变量的检验数: AX+Xs=b B XB + NXN +XS = b X ,Xs 0 对偶:对偶: Min w=Y b YA-Ys=C Y,Ys 0 XB = B-1 b - B-1 NXN - B-1 XS 代入代入Z=C X = CB XB + CN XN = CB B-1 b - CB B-1 NXN - CB B-1 XS + CN XN = CB B-1 b + ( CN - CB B-1 N) XN - CB B-1 XS 注注意意非非基基变变量

34、量的的检检验验数数和和松松弛弛变变量量的的检检验验数数,松松弛弛变变量量XS的的检检验验数数 - CB B-1对对应应对对偶问题的一个解偶问题的一个解 CB B-1 52管管 理理 运运 筹筹 学学第四节第四节 对偶单纯形法对偶单纯形法按对偶问题与原问题之间的关系,对最大化问按对偶问题与原问题之间的关系,对最大化问题,在用单纯形法求解原问题时,最终表不但题,在用单纯形法求解原问题时,最终表不但给出了原问题的最优解,而且其检验数的相反给出了原问题的最优解,而且其检验数的相反数就是对偶问题的最优解。数就是对偶问题的最优解。53管管 理理 运运 筹筹 学学单纯形法求解的基本思路单纯形法求解的基本思路

35、基可行解基可行解检验数非正检验数非正保持解的可行性保持解的可行性对偶单纯形法的基本思路对偶单纯形法的基本思路对偶问题基可行解对偶问题基可行解(检验数非正)(检验数非正)原问题基可行原问题基可行解解保持对偶问题解的保持对偶问题解的可行性(检验数非可行性(检验数非正正(对偶问题可行解)(对偶问题可行解)54管管 理理 运运 筹筹 学学 对偶单纯形法的基本思想对偶单纯形法的基本思想 对对偶偶单单纯纯形形法法的的基基本本思思想想是是:从从原原规规划划的的一一个个基基本本解解出出发发,此此基基本本解解不不一一定定可可行行,但但它它对对应应着着一一个个对对偶偶可可行行解解(检检验验数数非非正正),所所以以

36、也也可可以以说说是是从从一一个个对对偶偶可可行行解解出出发发;然然后后检检验验原原规规划划的的基基本本解解是是否否可可行行,即即是是否否有有负负的的分分量量,如如果果有有小小于于零零的的分分量量,则则进进行行迭迭代代,求求另另一一个个基基本本解解,此此基基本本解解对对应应着着另另一一个个对对偶偶可可行行解解(检验数非正)。(检验数非正)。55管管 理理 运运 筹筹 学学 如果得到的基本解的分量皆非负则该基如果得到的基本解的分量皆非负则该基本解为最优解。也就是说,对偶单纯形法在本解为最优解。也就是说,对偶单纯形法在迭代过程中始终保持对偶解的可行性(即检迭代过程中始终保持对偶解的可行性(即检验数非

37、正),使原规划的基本解由不可行逐验数非正),使原规划的基本解由不可行逐步变为可行,当同时得到对偶规划与原规划步变为可行,当同时得到对偶规划与原规划的可行解时,便得到原规划的最优解。的可行解时,便得到原规划的最优解。56管管 理理 运运 筹筹 学学对偶单纯形法在什么情况下使用:对偶单纯形法在什么情况下使用: 应用前提:有一个基,其对应的基满足应用前提:有一个基,其对应的基满足: : 单纯形表的检验数行全部非正(对偶单纯形表的检验数行全部非正(对偶可行);可行); 变量取值可有负数(非可行解)。变量取值可有负数(非可行解)。 注:注:通过矩阵行变换运算,使所有相应变通过矩阵行变换运算,使所有相应变

38、量取值均为非负数即得到最优单纯性表。量取值均为非负数即得到最优单纯性表。57管管 理理 运运 筹筹 学学对偶单纯形法的适用范围对偶单纯形法的适用范围 对对偶偶单单纯纯形形法法适适合合于于解解如如下下形形式式的的线线性性规规划划问题问题58管管 理理 运运 筹筹 学学对偶单纯形法计算步骤对偶单纯形法计算步骤适应于求解这样的适应于求解这样的LP问题:问题:标准化后不含初标准化后不含初始基变量,但将某些约束条件两端乘以始基变量,但将某些约束条件两端乘以“-1”后,即可找出初始基变量。后,即可找出初始基变量。要求:要求:初始单纯形表中的检验数满足最优性初始单纯形表中的检验数满足最优性条件条件59管管

39、理理 运运 筹筹 学学4 4 对偶单纯形法对偶单纯形法 对偶单纯形法也是解决线性规划问题的一种方法。对偶单纯形法是在保持原有问题的所有检验数都小于0的情况下,通过迭代使得所有的约束都大于等于0,最后求得最优解。 简化计算是对偶单纯形法的优点,但是它在使用上有很大的局限,这主要是大多数线性规划问题很难找到初始解使得其所有检验数都小于0。但是在灵敏度分析中,有时需要对偶单纯形法,这样可以简化处理。下面以第二节例一为例。 上节分析中已知当250b1325时第一个约束条件的对偶价格不变,现在b1=300变成b1=350,请问这时第一个约束方程的对偶价格应为多少呢? 解:求出在第二次迭代表上的常数列 6

40、0管管 理理 运运 筹筹 学学4 4 对偶单纯形法对偶单纯形法迭代次数基变量CBX1 X2 S1 S2 S3b50 100 0 0 02X1501 0 1 0 -1100 S200 0 -2 1 1-50 X21000 1 0 0 1250 ZJ50 100 0 50 50CJ -ZJ0 0 -50 0 -501.确定出基变量,在常数列中找一个最小的负常量,把这个常量所在行作为出基变量61管管 理理 运运 筹筹 学学4 4 对偶单纯形法对偶单纯形法 4.检查常数列值,若已经都非负结束迭代,即为最优,如果还有负数重复1-4步。迭代次数基变量CBX1 X2 S1 S2 S3b50 100 0 0

41、02X1501 0 1 1/2 -1/275 S200 0 -2 -1/2 -1/225 X21000 1 0 0 1250 ZJ50 100 0 25 7528750CJ -ZJ0 0 0 -25 -7562管管 理理 运运 筹筹 学学例例 用对偶单纯形法求解如下的用对偶单纯形法求解如下的LP问题问题化成标准形式63管管 理理 运运 筹筹 学学将各约束条件两端同乘“-1”得用对偶单纯形法求解得64管管 理理 运运 筹筹 学学最优解:最优解:最优解:最优解:x x1 1=0, x=0, x2 2=1/4, x=1/4, x3 3=1/2, x=1/2, x4 4=0, x=0, x5 5=0=0最优目标函数值:()最优目标函数值:()最优目标函数值:()最优目标函数值:()注:通常很少直接使用对偶单纯形法求解线性规划问题。65管管 理理 运运 筹筹 学学单纯形法和对偶单纯形法步骤是是是是否否否否所有所有得到最优解计算计算典式对应原规划的基本解是可行的典式对应原规划的基本解的检验数所有所有计算计算以为中心元素进行迭代以为中心元素进行迭代停没有最优解没有最优解单纯形法对偶单纯形法66

展开阅读全文
相关资源
正为您匹配相似的精品文档
相关搜索

最新文档


当前位置:首页 > 高等教育 > 研究生课件

电脑版 |金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号