对偶与灵敏度分析
- 格式:doc
- 大小:578.54 KB
- 文档页数:20
对偶问题例题1:某养鸡场所用的混合饲料由n 种天然饲料配合而成。
要求在这批配合饲料中必须含有m 种不同的营养成分,且第i 种营养成分的含量不低于bi 。
已知第i 种营养成分在每单位第j 种天然饲料中的含量为a ij ,每单位第j 天然饲料的价格为c j 。
试问,应如何对这n 种饲料配方,使这批饲料的费用最小? 解 设x j 为第j 种天然饲料的用量。
显然,a ij x j 即为所用第j 种天然饲料中第i 种营养成分的含量,1nij j j a x =∑为这批混合饲料中第i 种营养成分的总含量;它不应低于bi 。
于是,我们得下列线性规划模型(1—1):1min nj jj f c x ==∑11,,..01,,nij j i j j a x b i m s t x j n=⎧≥=⎪⎨⎪≥=⎩∑现设想有一个饲料加工厂欲把这m 种营养成分分别制成m 种营养丸。
设第i 种营养丸的价格为ui(i =1,…,m)。
则养鸡场采购一个单位的第j 种天然饲料,就相当于对这m 种营养丸分别采购数量a 1j ,…a mj ,所化费用为1mij ii a u =∑养鸡场自然希望在用营养丸代替天然饲料时,在价格上能相对地比较便宜,故而饲料加工厂为了能与天然饲料供应者竞争,在制订价格时必然满足下述条件:11,,mij ij i a uc j n =≤=∑另一方面,养鸡场如果全部采购营养丸来代替天然饲料进行配料,则第i 种营养丸就需采购bi 个单位,所化费用为b i u i ,总费用为z=∑b i u i饲料加工厂面临的问题是:应把这m 种营养丸的单价ui(f=1,…,m)定为多少,才能使养鸡场乐意全部采用该厂生产的营养丸来取代这批天然饲料,且使本厂在竞争中得到最大收益。
为该问题建立数学模型,即得如下线性规划(1—2):1max mi i i z b u ==∑11,,..01,,mij i j i ia u c j n s t u i m =⎧≤=⎪⎨⎪≥=⎩∑我们称问题(1—2)为原有问题 (1—1)的对偶问题(记为(D))。
《运筹学》第三章线性规划对偶理论与灵敏度分析习题及答案一、填空题1. 在线性规划问题中,若原问题存在最优解,则其对偶问题也一定存在最优解,这是线性规划的基本性质之一,称为______。
答案:对偶性2. 在线性规划问题中,若原问题与对偶问题均存在可行解,则它们均有______。
答案:最优解3. 对于线性规划问题,若原问题约束条件系数矩阵为A,目标函数系数向量为c,则其对偶问题的目标函数系数向量是______。
答案:c的转置(c^T)二、选择题1. 线性规划的原问题与对偶问题之间的关系是:A. 原问题的最优解和对偶问题的最优解相同B. 原问题的最优解是对偶问题的最优解的负数C. 原问题的最优解与对偶问题的最优解互为对偶D. 原问题的最优解和对偶问题的最优解没有关系答案:C2. 在线性规划中,若原问题不可行,则其对应的对偶问题:A. 可行B. 不可行C. 无界D. 无法确定答案:B三、判断题1. 线性规划的原问题和对偶问题具有相同的可行解。
()答案:错误2. 若线性规划的原问题存在唯一最优解,则其对偶问题也一定存在唯一最优解。
()答案:正确四、计算题1. 已知线性规划问题:max z = 3x1 + 2x2s.t.x1 + 2x2 ≤ 42x1 + x2 ≤ 5x1, x2 ≥ 0求该问题的对偶问题,并求解原问题和对偶问题的最优解。
答案:对偶问题为:min w = 4y1 + 5y2s.t.y1 + 2y2 ≥ 32y1 + y2 ≥ 2y1, y2 ≥ 0原问题和对偶问题的最优解如下:原问题最优解:x1 = 2, x2 = 1,最大利润z = 8对偶问题最优解:y1 = 2, y2 = 1,最小成本w = 82. 某工厂生产甲、乙两种产品,生产一件甲产品需要2小时的机器时间和3小时的工人劳动时间,生产一件乙产品需要1小时的机器时间和1小时的工人劳动时间。
工厂每周最多能使用12小时的机器时间和9小时的工人劳动时间。
线性规划中的对偶问题与灵敏度分析线性规划是一种优化方法,广泛应用于各个领域的决策问题。
在线性规划中,对偶问题与灵敏度分析是两个重要的概念和工具,可以帮助我们更好地理解和解决实际问题。
1. 对偶问题在线性规划中,对偶问题是指与原始问题相对应的一个问题。
它通过转换原始问题并构造一个新的问题,以便从不同的角度来解释和解决原始问题。
对偶问题能够提供原始问题的一些有用信息,并且在某些情况下,对偶问题的解与原始问题的解是相等的。
对偶问题的构造可以通过拉格朗日对偶性理论来完成。
该理论通过构造一个拉格朗日函数,将原始问题中的约束条件转化为拉格朗日乘子,从而得到对偶问题。
对偶问题的目标函数是原始问题的约束条件的线性组合。
解决对偶问题可以通过求解拉格朗日函数的最优化问题来实现。
对于线性规划问题,对偶问题的解可以通过求解一组线性方程或线性不等式来获得。
对偶问题的解不仅可以提供原始问题的一些信息,还可以用于检验原始问题的解的可行性和最优性。
2. 灵敏度分析灵敏度分析是在线性规划中评估解决方案对问题参数变化的响应程度的方法。
它可以帮助我们了解如果问题的参数发生变化,对解决方案的影响有多大,并做出相应的调整和决策。
灵敏度分析可以通过改变单个参数或多个参数来进行。
其中,常见的灵敏度分析包括目标函数系数的变化、约束条件右侧常量的变化和新增或取消约束条件。
这些变化可以用来模拟实际情况中可能发生的条件变化,以及评估解决方案的稳定性和可行性。
在进行灵敏度分析时,我们可以通过计算变动参数对解决方案的影响程度来得到一些关键指标。
例如,参数的变化导致目标函数值的变化量称为“影子价格”,而约束条件右侧常量的变化导致解决方案中相应决策变量的变化量,则称为“机会成本”。
灵敏度分析的结果可以帮助我们确定参数的重要性,判断解决方案的可行性和稳定性,以及找到最佳的决策方案。
在实际应用中,灵敏度分析可以帮助我们应对不确定性和风险,做出更加准确和可靠的决策。
§2 对偶与灵敏度分析§2.1 LP 的对偶问题无论从理论和实践角度,对偶理论是LP 中的一个最重要和有趣的概念,支持对偶理论的基本思想是:每一个LP 问题都存在一个与其对偶的问题,在求解一个问题解的时候,也同时给出了另一问题的解。
一、问题的提出例2.1:设某工厂生产两种产品甲乙,生产过程需要4种设备ABCD 进行加工,每件产品加工所需机时数,每件产品的利润值及每种设备的可利用机时如下表:1.问:充分利用设备时,应怎样安排甲乙产品的生产数量,利润才能最大?2.问:如有另外一家公司想租用该厂设备加工生产,那么,这家公司应至少对每台设备的机时价格为多少时,才能使该厂愿意出租设备? 解:1.设甲乙产品各生产1x 2x 件LP1:⎪⎪⎩⎪⎪⎨⎧≥≤≤+≤++=0,16482122232211212121x x x x x x x x x MaxZ 2.设每台设备的机时最低价分别为:1y ,2y ,3y ,4yLP2:⎪⎩⎪⎨⎧=≥≥++≥+++++=4,3,2,1,0342224212168124213214321i y y y y y y y y y y y MinZ i二、原问题和对偶问题之间的关系: 1.对称形式下的原问题与对偶问题对称形式下原问题的一般式: 矩阵形式:⎪⎪⎪⎩⎪⎪⎪⎨⎧=≥≤+++≤+++≤++++++=n j x bx a x a x a b x a x a x a b x a x a x a x c x c x c MaxZ j mn mn m m n n n n nn .......21,0 (221)122222121112121112211⎩⎨⎧≥≤=0X b AX CX Max 若用i y 代表第i 种资源的估价,则其对偶问题的一般式为:⎪⎪⎪⎩⎪⎪⎪⎨⎧=≥≥+++≥+++≥++++++=m j y cy a y a y a c y a y a y a c y a y a y a y b y b y b MinZ j nm mn n n m mn m m mm .......21,0 (221)122222112112211112211⎩⎨⎧≥≥=0Y C Y A Yb Min T T ω 2.非对称形式下原问题与对偶问题:方法一:将非对称形式转化为对称形式,求出对偶问题,然后再还原。
例2.2写出下列LP 问题的对偶问题:⎪⎪⎩⎪⎪⎨⎧≤≥≤+--≥-++=--+-+++=无约束43214321432143214321,0,0,)3.........(..........2099912)2..(....................85376)1.....(....................53432x x x x x x x x x x x x x x x x x x x x MaxZ 为写出基对偶问题,先将其转化为对称形式,再进行变化:因目标函数为Max ,故约束条件全部转化为“0≤”,所有变量均应为“0≥”。
将(1)式变为:535343214321≤--+-≥--+-x x x x x x x x 和。
再将535343214321-≤-+-≥--+-x x x x x x x x 转化为将(2)式两端同乘以“-1”,并令:0,,;0,''4'4''4'44'33'3≥-=≥∴-=x x x x x x x x 其中得:855376''4'4'321-≤-++--x x x x x所以,例2.2可以变为:⎪⎪⎪⎩⎪⎪⎪⎨⎧≥≤-++--≤-++---≤-+--≤+-++--+-+=0,,,,20999912855376533533)(432''4'4'321''4'4'321''4'4'321''4'4'321''4'4'321''4'4'321x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x MaxZ 令对应上述四个约束条件的对偶变量为3'2''1'1,,,y y y y ,则其对偶问题为: ⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧≥-≥---≥+++--≥++-≥---≥+-+-+--=0,,,4953349533393297112'6208553'2''1'13'2''1'13'2''1'13'2''1'13'2''1'132''1'13'2''1'1y y y y y y y y y y y y y y y y y y y y y y y y y y y y Min ω令'22''1'11,y y y y y -=-=,将第4与第5个不等式合并,将第三个不等式两端同乘以“-1”得:⎪⎪⎪⎩⎪⎪⎪⎨⎧≥≤=+--≤-+-≥-+≥++-++=0,0,495339329711262085321321321321321321y y y y y y y y y y y y y y y y y y Min 无约束ω 由以上的推导可以发现有以下规律:方法二:根据原问题与对偶问题之间的关系,可将其归纳如下表:§2.2 对偶问题的基本性质一、单纯形法计算的矩阵描述设线性规划问题的矩阵表达式为:⎩⎨⎧≥≤=0X bAX CXMaxZ ,加上松弛变量s X 后为:⎩⎨⎧≥≥=++=0,00s s s X X b IX AX X CX MaxZ 式中:单位矩阵为m m I x x x X m n n n s ⨯=+++,),......,(21。
单纯形法计算时,总选取为初始基I ,对应基变量为s X 。
设迭代若干步后,基变量为B X ,B X 在初始单纯形表中的系数矩阵为B 。
将B 在初始单纯形表中单独列出,而A 中去掉B 的若干列后组成矩阵N ,同时将C 也分为两块),(N B C C ,B C 是目标函数中基变量B X 的系数行向量;N C 是目标函数中非基变量N X 的系数行向量。
这样LP 的初始单纯形表可列成如下形式:于是原LP 问题可以改写为:⎩⎨⎧≥=++++=0,,0s N B s N B sN N B B X X X bIX NX BX X X C X C MaxZ 将约束条件移项并左乘1-B 后得到B X 的表达式:s N B IX B NX B b B X 111-----=代入目标函数得:s B N B N B X B C X N B C C b B C Z 111)(-----+= 所以,经过迭代之后,得出新的单纯形表为:当为最优基时,在上述单纯形表中有:01≤--N B C C B N ,01≤--B C B 而B X 的检验数可以写为:0=-I C C B B 因此有:01≤--A B C C B ,01≤--B C B称1-B C B 为单纯形乘子,若令1-=B C Y B T则有:⎩⎨⎧≥≥0Y C Y A TT可以发现这时检验数行,若取其相反数恰好是其对偶问题的一个可行解。
将这个解代入对偶问题的目标函数值,有:Z b B C b Y B T ===-1ω即:当原问题为最优解时,这时对偶问题为可行解,且两者具有相同的目标函数值。
(其实,这时对偶问题的解也为最优解。
)二、本节讨论先假定原问题和对偶问题为对称形式的LP ,即原问题为:⎩⎨⎧≥≤=0X b AX CX MaxZ 其对偶问题为:⎩⎨⎧≥≥=0Y C Y A YbMin T T ω然后说明对偶问题的基本性质在非对称形式也适用。
1.对称性:对偶问题的对偶是原问题。
证明:原问题与其对偶问题分别为:⎩⎨⎧≥≤=0X b AX CX Max 与⎩⎨⎧≥≥=0Y C Y A Yb Min T T ω2.弱对偶性:若X 是原问题可行解,Y 是对偶问题的可行解,则存在b Y X C ≤。
证明:设原问题为:0≥≤=X bAX CX Max原问题的对偶问题为:0≥≥=Y C Y A Yb Min TT ω3.无界性:若原问题为无界解,则其对偶问题为无可行解。
证明:由弱对偶性可得,本性质成立。
4.可行解是最优解时的性质:若X 是原问题可行解,Y 是对偶问题的可行解,当b Y X C =时,X 与Y 是原问题与对偶问题的最优解。
证明:5.对偶定理:若原问题有最优解,那么其对偶问题也有最优解,且目标函数值相等。
6.互补松弛性:若X 是原问题可行解,Y 是对偶问题的可行解,那么00==X Y X Y s s 和,当且仅当为X ,Y 最优解。
(在线性规划问题的最优解中,如果对应某一约束条件的对偶变量值为非零,则该约束条件取严格等式;反之如果约束条件取严格不等式,则其对应的对偶变量值一定为零。
)0,01==∑=si i nj j ij i x b x a y 则有 若 ;00,1=∑=i si i nj j ij y x b x a 则有 若 。
将互补松弛性质应用于其对偶问题时可以这样叙述:j mi i ij j c y a x =∑=1,0则有 若 ;0,1=∑=j j mi i ij x c y a 则有 若 。
由互补松弛定理可知:若原问题最优解中的第j个变量为正,则其对偶问题中与之对应的第j个约束条件在最优情况下必呈严格等式(即其松弛变量为0);而如果对偶问题中第j个约束条件在最优情况下呈严格不等式,则原问题最优解中第j个变量必为0。
类似地,若对偶问题最优解中第i个对偶变量为正,则其原问题中与之对应的第i个约束条件在最优情况下必呈严格等式(即其剩余变量为0);而如果原问题中第i个约束条件在最优情况下呈严格不等式,则对偶问题最优解中第i个对偶变量必为0。
互补松弛定量阐明了原问题及其对偶问题最优解各分量之间的关系,当已知一对对偶问题之一的最优解时,可利用该定量求出另一个问题的最优解。
7.设原问题是:0,≥=+=s s X X bX AX CX Max对偶问题为:0,≥=-=s s Y Y CY YA YbMin ω则原问题单纯形表的检验数行对应其对偶问题的一个基解,其对应关系如下表所示:这里1S Y 是对应原问题中B X 基变量的剩余变量,2S Y 是对应原问题中非基变量N X 的剩余变量。