第三章线性规划的对偶理论及灵敏度分析1总结
- 格式:docx
- 大小:75.83 KB
- 文档页数:20
精品文档第三章 线性规划的对偶理论及灵敏度分析主要内容:1、对偶问题及其性质; 2、对偶单纯形法;3、灵敏度分析。
重点与难点:对偶问题与原问题的对应关系,对偶问题的基本性质,对偶单纯形法的求解步骤,灵敏度分析的方法。
要 求:理解线性规划对偶问题的性质,熟练掌握对偶单纯形法的求解步骤和灵敏度分析的方法和技巧,能够用这些数学方法解决实际问题。
§1 对偶问题的对称形式一、对偶问题引例,某工厂在计划期内要安排生产甲、乙两种产品,已知生产单位产品所需要的设备台时及A 、B 两种原材料的消耗,该工厂每生产一件产品甲可获利2元,每生产一件产品乙可获利3元,问应如何安排计划才能使该工厂获利最多?解:设1x 、2x 分别为甲、乙两种产品的产量则目标函数2132m ax x x z +=约束条件 ⎪⎪⎩⎪⎪⎨⎧≥≤≤≤+0,12416482212121x x x x x x(1)假设该工厂决定不再生产甲、乙产品,而将其出租或出售。
这时要考虑每种资源的定价问题,设321,,y y y 分别为出租单位设备台时的租金和出让单位原材料A 、B 的附加额。
作一比较:若用一个单位台时和4个单位原材料A 生产一件产品甲,可获利2元,那么生产每件产品甲的设备台时和原材料出租和出让的收入应不低于生产一件甲产品的利润。
即:2421≥+y y同理,将生产每件乙产品的设备台时和原材料出租和出让的收入应不低于生产一件乙产品的利润。
即:精品文档34231≥+y y将工厂所有设备台时和资源都出租和出让,其收入为32112168y y y ++=ω对工厂来说,ω越大越好;但对接受者来说,支付的愈少愈好,所以工厂只能在满足≥所有产品的利润前提下,使其总收入尽可能小,才能实现其愿望。
为此,得到如下模型:32112168m in y y y ++=ω⎪⎩⎪⎨⎧=≥≥+≥+3,2,1,0342243121j y y y y y j(2)我们就称(2)为模型(1)的对偶问题。
第三章线性规划的对偶理论及灵敏度分析主要内容:1、对偶问题及其性质;2、 对偶单纯形法;3、 灵敏度分析。
重点与难点:对偶问题与原问题的对应关系,对偶问题的基本性质,对偶单纯形法的求解步骤,灵敏度分析的方 法。
要求:理解线性规划对偶问题的性质,熟练掌握对偶单纯形法的求解步骤和灵敏度分析的方法和技巧,能够用这些数学方法解决实际问题。
§ 1对偶问题的对称形式一、对偶问题弓侧,某工厂在计划期内要安排生产甲、乙两种产品,已知生产单位产品所需要的设备台时及 A 、B 两种原材料的消耗,该工厂每生产一件产品甲可获利 2元,每生产一件产品乙可获利 3元,问应如何安排计划才能使该工厂获利最多?解:设X i 、X 2分别为甲、乙两种产品的产量作一比较:若用一个单位台时和 4个单位原材料 A 生产一件产品甲,可获利 2元,那么生产每件产品甲的设备台 y^ 4y^ 2同理,将生产每件乙产品的设备台时和原材料出租和出让的收入应不低于生产一件乙产品的利润。
即:2力 4y 33将工厂所有设备台时和资源都出租和出让,其收入为则目标函数maxz 二2x 「3x 2x 「2x 2 岂8i4x 1 - 16 i4x 2 兰12约束条件-x 1,x^ 0(1)不再生产甲、乙产品,而将其出租或出售 3分别为出租单位设备台时的租金和出让单位原材料这时要考虑每种资源的定价问题,设A 、B 的附加额。
时和原材料出租和出让的收入应不低于生产一件甲产品的利润。
即:。
=8y 〔+ 16y 2 + 12y 3对工厂来说,••越大越好;但对接受者来说,支付的愈少愈好,所以工厂只能在满足》所有产品的利润前提下, 使其总收入尽可能小,才能实现其愿望。
为此,得到如下模型:min =8y 1 16y 212y 3"+4丫2工 2< 2y i +4y ^ 3 J j > 0 , j =1,2,3我们就称(2)为模型(1)的对偶问题。
一般地,设原问题为max z = c/ c 2x 2 ……c x n'a ii X i +a i2X 2 + … +amX n 兰 b a 2l X l +a 22X 2 +■八 +a 2n X n 兰 b 2a a a a ■■■■■s■■■■a mi Xi +a m2X2 +*a mn X n 兰 *X j _0 , j =i,2, ,n则其对偶问题为:min 二 by b ?y 2^^nNiy i +a 2〃2 + …+a mi y m A"a i2yi +a 22y2 * +a m2ym ®C 2m- a - < ■■■■a inyi +a 2ny2 + +a mnym ®C ny i 一0 , i =i,2, ,m矩阵形式:原问题 对偶问题max z = cX mi n = Yb'AX E b , 、Y A 启C (实际为A T y T ^C T )X >07 >0、原问题与对偶问题的关系例1求下列问题的对偶问题min z =2x1 3x2- 5x3x4x1 +x2-3x3 +x4 >52x1+2x3—x4兰41x2+x3+ x4=6捲_ 0, x2, x3- 0, x4无约束解:max = 5y! 4y26y3» +2y2 3 2y i *3 兰3«—3% +2y? +y3 兰一5y i -丫2*3=1y i -0』2空0小无约束§ 2对偶问题的基本性质、对称性:对偶问题的对偶是原问题。
证:设原问题为max z 二cXi X 0mi n = Yb则其对偶问题为:YA_ Ci y 0-min 二-Yb对上式两边取负号,得YA C1y0-max(-代)= min wmax(-⑷)=_Yb-YA-JCY -0上式的对偶问题为min(v)=CX-AXJ--bX-0(两边同取负号)-min(- v)二maxv maxv 二CX 二maxzAX bX 0(0)(0) cX (0)::Y(°)b二、弱对偶性:若X 是原问题的可行解,Y 是对偶问题的可行解,则存在CX 一丫b。
(0)证:X 是原问题的可行解同理Y(0)A— C,用X(0)右乘之得丫(0)AX(0)一CX(0)已知Y (0)是对偶问题的可行解,用Y(0)左乘上式得Y (0)AX(0)二Y(0)bCX (0) ’ Y (0)AX (0)’ 丫⑼b ,故CX (0X Y (0)b三、 无界性:若原问题(对偶问题)为无界解,则其对偶问题(原问题)无可行解。
注意:此性质不可逆。
(0) (0)四、 可行解是最优解时的性质最优性:设X是原问题的可行解,Y是对偶问题的可行解,当(0) . (0) (0) (0)CX - 丫 b 时,X 、丫是最优解。
五、 对偶定理(强对偶性):若原问题有最优解,那么对偶问题也有最优解,且目标函数值相等。
反之,若其一 无最优解,则另一也无最优解。
(0)(0) (0)、 ^六、 互补松弛性:若X 、Y分别是原问题和对偶问题的可行解,那么Y X s 二0和二 YAX YX s若Y s X (0)= 0, 丫(0)X s 二 0 ;则 Y (0)b 二 Y (0)AX(0)二CX (0)(0) (0)由性质4知,X 、 Y为最优解。
又如果X(0)、丫(0)为原问题和对偶问题的最优解,由性质4有CX (0)= Y (0)AX (0)= Y (0)b 即Y (0)AX (0)-Y s X (0)= 丫⑼AX (0)= Y (0)AX (0)Y (0)X s必有Y s X (0) = 0,丫(0)X s =例2已知线性规划问题(0) (0)Y S X 0,当且仅当X证:设原问题和对偶问题的标准型是max z CXAX X s = b X,X s - 0将C 二 YA - Y s , b 二 AX(0)Y 为最优解。
mi n =YbYA- Y s 二 C Y,Y s - 0X s 分别代入原问题和对偶问题目标函数得maxz= x「x2捲 + x 2 + X 3 兰2 r2x< x 2 x 3 ' 1 X i , X 2, X 3试用对偶理论证明上述线性规划问题无最优解。
上述问题的对偶问题为min 二 2y 「y 2-y< 2y^ 1 y 「y 2 - 1 yi - y^ 0 y i , y^ 0由第一个约束条件知,对偶问题无可行解,所以,由对偶定理知,原问题无最优解。
七、对偶问题的经济解释----影子价格由对偶定理可知,当达到最优解时,原问题和对偶问题的目标函数值相等,即有求z 对b 的偏导数得:其经济学意义是:在其它条件不变的情况下,单位资源变化所引起的目标函数的最优值的变化。
i 的值代表对第i 种资源的估价,这种估价是针对具体工厂的具体产品而存在的一种特殊价格,称它为“影 子价格”。
影子价格随具体情况而异,在完全市场经济条件下,当某种资源的市场价格低于影子价格时,企业应买进该 资源用于扩大生产;而当某种资源的市场价格高于企业影子价格时,则企业应把已有的资源卖掉。
可见,影子价格对证:原问题存在可行解,如(0,0,0)Tz cX ⑼二 丫⑼b 二y 10)b iy 20)b 2(0)y:z (0)石,y 2(0),y mb m*即y i*■:z b i§ 3对偶单纯形法—、基本思路对偶单纯形法是运用对偶原理求解原问题的一种方法,而不是求解对偶问题的单纯形法。
首先讨论这样一个问题:min 二 Yb; YA 丫眇 c; Y,Y s 0设B 是原问题的一个可行基,于是 A=(B|N),原问题可改写为:maxz 二 C B X B C N XBX B NX N X s = b &B ,X N ,X s = 0相应地对偶问题可以表示为minmin 二 YbY B 飞=C B ................. (1) Y N - Ys^ = C N .. (2)丫,丫0,丫5 - 0这里Y s = (Yd%)Y$ ----对应原问题中基变量 X B的剩余变量市场有调节作用。
Ys 2 ----对应原问题中非变量X N 的剩余变量设原问题: maxz=cX; AX 空 b, X - 0则其对偶问题:min=Yb; YA - c; 丫- 0化为标准型:max zcX; AX X sb;X,X s 0C B B 1。
现当求得原问题的一个基解X B二B’b,其相应的检验数为C N一C B B 1N与一分析这些检验数与对偶问题的解之间的关系:,1代入⑴、(2)得令Y = C B BY S^ 0飞=C N C B B1N由此可得出:原问题单纯形表的检验数行对应其对偶问题的一个基解,其对应关系如下:说明:在单纯形表中若在b列中得到的是原问题的基可行解,而在检验数行得到的是对偶问题的基解,当在检验数行得到对偶问题的解也是基可行解时,根据性质知,已知到最优解,即原问题与对偶问题是最优解。
根据对偶问题的对称性,可这样考虑:若保持对偶问题的解是基可行解,即C j -C B B」P j乞0,而原问题在非可行解的基础上,逐步迭代达到基可行解,这样也得到了最优解。
方法是:设原问题maxz=CX;AX =bX艺0设B是一个基,令B =(Ph,P2,…,P m),它对应的变量为X B=(为公2,…,X m)T当非基变量都为零时,可以得到X B=B,b,若在B4b中至少有一个负分量,设2北)「::0,并且在单纯形表的检验数行中的检验数都为负值,即对偶问题保持可行解,它的各分量是:1. 对应基变量x1, x2 / ,x m的检验数是-^c^C B B4P i =0, i =1,2, ,m2. 对应非基变量x m 1,x m -2/' ,x n的检验数是:■j =C j -C B B P j g j 二m 1,m 2, ,n每次迭代是将基变量中的负分量X l取出,去替换非基变量中的X k,经基变换,所有检验数仍保持负值,原问题逐步由非可行解向可行解靠近,当原问题得到可行解时,便得到了最优解。
二•计算步骤(1) 列出初始单纯形表,若所0,‘ j ' 0,则停止计算,已得到最优解。
若b中含有负元素,有b ■则需继续计算。
⑵确定换出变量叫门{( B b)i (B b)i < 0} = (B b)i ,基变量X |为换出变量。
(3)确定换入变量检查X |行的系数aj ,若所有a ij > 0,则无可行解,停止计算。
若存在 3|j <0,则继续计算。
(4)以 a lk 为主元素进行取主变换。
例3、用对偶单纯形法求解。
min z = 2x 「3x 2 4x 3x 「2X 2 X 33 r2x< x 2 3X 3 4x i , X 2, X 3解:化为标准型max z 二 2x< 3x 24x 3maxz 二 -2捲 - 3x 2 - 4x 3 0x 4 0x 5-- 2x 2 - x 3兰-3l2x 「x 23x 3 '4 心X 2, X 3— Y — O Y — Y + w 兰 —O12342x 「x 23x 3 x 5'4 X j - 0 ,厂 1,2, ,5I 貯,日9 =min t 一 按"规则, ja,,ijIf Io<_ka ,非基变量x k 为换入变量。