对偶理论与灵敏度分析.ppt
《对偶理论与灵敏度分析.ppt》由会员分享,可在线阅读,更多相关《对偶理论与灵敏度分析.ppt(90页珍藏版)》请在三一办公上搜索。
1、1,第二章 线性规划的对偶理论,窗含西岭千秋雪,门泊东吴万里船,本章主要内容:线性规划的对偶问题概念、理论及经济意义线性规划的对偶单纯形法线性规划的灵敏度分析,2,第一节 对偶问题的提出,假设工厂考虑不进行生产而把全部资源都转让,问如何定价这些资源,既能使其获利不低于安排生产所获得的收益,又能使资源租让具有竞争力。,一、引例,Max Z=2000 x1+4000 x2+3000 x3 s.t.3x1+4x2+2x3600 2x1+x2+2x3400 x1+3x2+3x3300 x1+2x2+4x3200 x1,x2,x30,Min W=600y1+400y2+300y3+200y4 s.t.3
2、y1+2y2+y3+y42000 4y1+y2+3y3+2y44000 2y1+2y2+3y3+4y43000 y1,y2,y3,y40,x1x2x3,y1 y2 y3 y4,3,比较上述模型,可以得出两者之间的一些关系:1.两个问题,一个是极大化,另一个是极小化;2.一个问题的变量数等于另一问题的方程数,反之亦然;3.一个问题的目标函数系数是另一个问题的约束方程右端常数,反之亦然;4.两个问题的约束方程系数矩阵互为转置。,称变量yi为第一个LP的第i个对偶变量,或第一个LP的第i约束相应的对偶变量,4,二、对偶问题,(1)对称LP问题的定义,(2)对称LP问题的对偶问题,第一类对称形式,第二
3、类对称形式,(1)变量为非负;(2)约束条件为不等式。对于max,约束为“”;对于min,约束为“”,5,例1 写出下列LP问题的对偶问题,Max Z=2x1+3x2 s.t.x1+2x28 4x1 16 4x2 12 x1,x2 0,Min W=8y1+16y2+12y3 s.t.y1+4y2 2 2y1+4y3 3 y1,y2,y3 0,6,(3)对偶问题的对偶是原问题,推导过程,变形,对偶,Max Z=CXs.t.AXb X0,Min W=Ybs.t.YAC Y0,Max W=-Ybs.t.-YA-C Y0,Min Z=(-C)Xs.t.(-A)X(-b)X0,Max W=Y(-b)s.
4、t.Y(-A)(-C)Y0,7,例2 写出下列LP问题的对偶问题,解:上述LP问题的 对偶问题为:,8,三、非对称LP问题的对偶问题,例3 写出下列LP问题的对偶问题,解:用x2=-x2,x3=x3-x3代入上述LP问题,并将其化为第一类对称形式,Max Z=x1+2x2+x3 x1+x2-x3 2 s.t.x1-x2+x3=1 2x1+x2+x3 2 x10,x20,x3无约束,Max Z=x12x2+x3 x3 x1 x2 x3+x3 2 x1+x2+x3 x3 1s.t.x1 x2 x3+x3 1 2x1+x2 x3+x3 2 x1,x2,x3,x3 0,x1-x2+x3 1x1-x2+
5、x3 1,x1-x2+x3 1-x1+x2-x3-1,-2x1-x2-x3-2,9,上述第一类对称形式LP问题的对偶问题为:,则上述问题变为:,Min W=2y1+y2 y32y4 y1+y2 y32y4 1 y1+y2 y3+y4-2s.t.y1+y2 y3y4 1 y1 y2+y3+y4-1 y1,y2,y3,y4 0,Min W=2u1+u2+2u3 u1+u2+2u3 1 s.t.u1 u2+u3 2 u1+u2+u3=1 u10,u30,u2无约束,令 u1=y1 u2=y2y3 u3=y4,Max Z=x12x2+x3 x3 x1 x2 x3+x3 2 x1+x2+x3 x3 1s
6、.t.x1 x2 x3+x3 1 2x1+x2 x3+x3 2 x1,x2,x3,x3 0,-y1+y2-y3-y4 1,y1-y2+y3-y4 2,10,对偶关系:一个问题第i个变量的约束情况决定另一问题第i个约束不等式的方向,反之亦然。,正常的对正常的,不正常的对不正常的,原问题约束条件是等式对应于对偶问题决策变量无约束,反之亦然,11,例3 直接写出LP问题的对偶问题,12,13,原问题,对偶问题,目标函数,Max z,Min w,变量,n个,约束条件,n个,=,0,0,无约束,约束条件,m个,=,变量,m个,0,0,无约束,约束条件右端项,目标函数变量的系数,目标函数变量的系数,约束条
7、件右端项,原问题,对偶问题,14,第二节 LP问题的对偶理论,若X(0),Y(0)分别为(L),(D)的可行解,则有CX(0)Y(0)b,定理1(弱对偶定理):极大化原问题目标函数值总是不大于其对偶问题的目标函数值。,证明:max z=CX;AXb;X0(L)min w=Yb;YAC;Y0(D)由于X(0)是(L)的可行解,有AX(0)b,X(0)0.由于Y(0)是(D)的可行解,有Y(0)0.Y(0)左乘不等式组AX(0)b的两边得:Y(0)AX(0)Y(0)b(1),15,min w=Yb;YAC;Y0(D)又Y(0)AC用X(0)(X(0)0)右乘上式,得Y(0)AX(0)CX(0)(2
8、)由(1),(2)得:CX(0)Y(0)AX(0)Y(0)b所以 CX(0)Y(0)b,16,推论1(P119),若对偶问题有无界解,则其原问题无可行解;若对偶问题无可行解,则其原问题或有无界解或无可行解。,原问题,对偶问题,若LP问题有无界解,则其对偶问题无可行解(弱对偶性);若LP问题无可行解,则对偶问题或有无界解或无可行解。,17,推论2,极大化问题(L)的任何一个可行解所对应的目标函数值都是其对偶问题(D)目标函数值的下界。,极小化问题(D)的任何一个可行解所对应的目标函数值都是其对偶问题(L)目标函数值的上界。,推论3,YbCX(0),CX Y(0)b,max z=CX;AXb;X0
9、(L)min w=Yb;YAC;Y0(D),18,例4 考虑下面一对LP问题,其对偶问题为:,由于X(0)=(1,1,1,1)T,Y(0)=(1,1)T分别为(L),(D)的可行解,,Max Z=x1+2x2+3x3+4x4 x1+2x2+2x3+3x4 20 s.t.2x1+x2+3x3+2x4 20 x1,x2,x3,x4 0,Min W=20y1+20y2 y1+2y2 1 2y1+y2 2 s.t.2y1+3y2 3 3y1+2y2 4 y1,y20,Z(0)=10,W(0)=40,W10,推论2,Z40,推论3,故Z40,W10.,19,例,试用对偶理论证明上述线性规划问题无最优解,
10、证明:首先看到该问题存在可行解,如X=(0,0,0)上述问题的对偶问题为,由第一约束条件知对偶问题无可行解,因原问题有可行解,由推论1,原问题无最优解。,若LP问题无可行解,则对偶问题或有无界解或无可行解。,20,定理2(最优性准则)当LP问题目标函数值与其对偶问题目标函数值相等时,各自的可行解即为最优解。(P57),若X(0),Y(0)分别为(L),(D)的可行解,且CX(0)Y(0)b,则X(0),Y(0)分别为(L),(D)的最优解。,证明:由定理1可知,对于(L)的任意可行解X,有CX Y(0)b 由于CX(0)=Y(0)b,CX CX(0),故X(0)为(L)的最优解;同理,由定理1
11、可知,对于(D)的任意可行解Y,有YbCX(0),由于CX(0)=Y(0)b,YbY(0)b,Y(0)为(D)的最优解。,21,例5,Max Z=x1+2x2+3x3+4x4 x1+2x2+2x3+3x4 20 s.t.2x1+x2+3x3+2x4 20 x1,x2,x3,x4 0,Min W=20y1+20y2 y1+2y2 1 2y1+y2 2 s.t.2y1+3y2 3 3y1+2y2 4 y1,y20,由于X(0)=(0,0,4,4)T,Y(0)=(6/5,1/5)T是(L),(D)的可行解且 CX(0)=bTY(0)=28,所以X(0),Y(0)分别为(L),(D)的最优解。,22,
12、线性规划的矩阵表示,Max Z=CXs.t.AX b X0,令A=(B,N),X=XB,C=(CB,CN)XN,Max Z=CX+0Xss.t.AX+IXs=b X,Xs0,用非基变量表示基变量,用非基变量表示目标函数,23,=(0,CN-CBB-1N,-CB B-1)=(C-CBB-1A,-CB B-1),(C-CBB-1A,-CB B-1)=(C-CBB-1(B,N),-CB B-1),=(CB-CBB-1B,CN-CBB-1N,-CB B-1),=(0,CN-CBB-1N,-CB B-1),若CN-CBB-1N0,-CB B-10,则Z为最优解,令XN=0,Xs=0,可以得到基B的基可行
13、解和目标函数值,XB=B-1b-B-1NXN-B-1Xs,Z=CBB-1b+(CN-CBB-1N)XN-CB B-1Xs,X=(XB,XN,Xs)=(B-1b,0,0),Z=CBB-1b,24,CB,XB,B-1b,CB,CN,0,I,CN-CBB-1N,B-1N,B-1,0,-CBB-1,-CBB-1b,25,初始单纯形表,0 x3 0 x4 0 x5,2 3 0 0 0,8 1612,1 2 1 0 0 4 0 0 1 0 0 4 0 0 1,2 3 0 0 0,0,4-3,经过一次迭代,初始基,第一轮迭代后的基,26,初始单纯形表,0 x3 0 x4 0 x5,2 3 0 0 0,8 1
14、612,1 2 1 0 0 4 0 0 1 0 0 4 0 0 1,2 3 0 0 0,0,4-3,经过一次迭代,B-1 b=,Z=CBB-1b=9,27,初始单纯形表,0 x3 0 x4 0 x5,2 3 0 0 0,8 1612,1 2 1 0 0 4 0 0 1 0 0 4 0 0 1,2 3 0 0 0,0,4-3,经过一次迭代,B-1A=,-CBB-1,=(0,0,-3/4),A,28,定理3(强对偶定理)若(L),(D)均有可行解,则(L),(D)均有最优解,且目标函数最优值相等。,证明:设X(0),Y(0)分别为(L),(D)的可行解,由弱对偶定理(推论3),对于(L)的任意可行
15、解X,有CX Y(0)b,所以CX在可行域内有上界,故(L)有最优解。同理,对于(D)的任意可行解Y,由弱对偶定理(推论2),有Y bCX(0),所以Yb在可行域内有下界,故(D)有最优解。,设(L)的最优解为X(0),且X(0)所对应的最优基为B,X(0)可以表示为X(0)=XB(0)=B-1b XN(0)0,max z=CX;AXb;X0(L),min w=Yb;YAC;Y0(D),29,则(A,S)=(C,0)CBB-1(A,I)=(C CBB-1A,CBB-1)0由于CCBB-1A0,所以CBB-1A C(1)又CBB-10,故CBB-10.(2),令Y(0)=CBB-1,由(1)、(
16、2)知Y(0)是(D)的可行解.因为CX(0)=(CB,CN)XB(0)=CBXB(0)+CNXN(0)=CBB-1b XN(0)而Y(0)b=CBB-1b 则CX(0)=Y(0)b,由最优性准则知,X(0),Y(0)分别为(L),(D)的最优解,且目标函数最优值相等。,30,推论:在用单纯形法求解LP问题(L)的最优单纯形表中,松弛变量的检验数的相反数就是其对偶问题(D)的最优解。,例5 求下列问题对偶问题的最优解,Max Z=2x1+3x2 s.t.x1+2x28 4x1 16 4x2 12 x1,x2 0,解:化为标准型,Max Z=2x1+3x2+0 x3+0 x4+0 x5 s.t.
17、x1+2x2+x3=8 4x1+x4=16 4x2+x5=12 x1,x2,x3,x4,x50,CB,XB,B-1b,CB,CN,0,I,CN-CBB-1N,B-1N,B-1,0,-CBB-1,-CBB-1b,31,此时达到最优解X*=(4,2)T,Max Z=14。对偶问题的最优解为Y*=(3/2,1/8,0)T,运用单纯形法计算得原问题的最终表如下:,32,推论:对偶问题的最优解对应于原问题的最优单纯型表中松弛变量检验数的相反数或剩余变量的检验数。,33,Min z=2x1+3x2 x1+x2 350 x1 1252x1+x2 600 x1,x2 0,Max w=350y1+125y2+6
18、00y3 s.t.y1+y2+2y32 y1+y3 3 y1,y2 0;y3 0,原问题的最优单纯型表如下所示:,0 0 4 0 1-4+M M,Y*=(4,0,-1)T,w*=800,例,34,定理4(互补松弛定理)在最优情况下,原问题的第i个决策变量与其对偶问题第i个约束中的松弛变量的乘积恒为零。,设X(0),Y(0)分别为(L),(D)的可行解,则X(0),Y(0)分别为(L),(D)的最优解的充要条件为,有,35,证明(P120),36,例7 考虑下面问题,Max Z=x1+2x2+3x3+4x4 x1+2x2+2x3+3x4 20 s.t.2x1+x2+3x3+2x4 20 x1,x
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 对偶 理论 灵敏度 分析
链接地址:https://www.31ppt.com/p-4904558.html