运筹学 第二章 灵敏度分析
《运筹学》胡运权 第4版 第二章 线性规划的对偶理论及灵敏度分析

b2 bm
x1, x2 , , xn 0
对 称 形 式 的
的 定 义
m W ib 1 n y 1 b 2 y 2 b m y m 对
s.t.
a11 a12 a1n
a21 a22 a2n
am1 y1 c1
am2 y2 amn ym
c2 cn
偶 问 题
y1, y2 , , ym 0
a23 x3 a33 x3
b2 b3
x1 0, x2 0, x3无 约 束
(2.4a) (2.4b) (2.4c) (2.4d)
先转换成对称形式,如下:
的 的一个变量,其每个变量对应于对偶问题 的一个约束。
定
义
m Z a c 1 x 1 x c 2 x 2 c n x n 一
对 偶
a11x1 a12x2 a1n xn (,)b1
a2
1x1
a22x2
a2n xn
(, )b2
般 线 性
问 题 的 定 义
am1x1 am2 x2 amnxn (,)bm xj 0( 0,或符号不限) j 1 ~ n
问题。
对
对偶问题是对原问题从另一角度进
偶
行的描述,其最优解与原问题的最 优解有着密切的联系,在求得一个
原
线性规划最优解的同时也就得到对 偶线性规划的最优解,反之亦然。
理
对偶理论就是研究线性规划及其对 偶问题的理论,是线性规划理论的
重要内容之一。
问 题 的 导 出
例2-1
我们引用第一章中美佳公司的例子,如表1
的
x1, x2, , xn 0
对
m W ib 1 n y 1 b 2 y 2 b m y m
灵敏度分析(运筹学)

最优基不变,即在最终表中求得的经过变化后 的b列的所有元素要求不小于0
目标函数 m ax z 2 x1 3x2 x1 2 x2 8 4x 16 1 约束条件 : 4 x2 12 x1 , x2 0
0 x3 1 -2 1/2 -3/2 0 x4 1/4 1/2 -1/8 -1/8 0 x5 0 1 0 0 θ
(5)按照下表所列情况得出结论或继续计算的步 骤。
原问题 可行解 可行解 非可行解 非可行解 对偶问题 可行解 非可行解 可行解 非可行解 结论或继续计算的步骤 原最优基不变 用单纯形法继续迭代 用对偶单纯形法继续迭 代 引入人工变量 ,扩大原 单纯形表继续计算
资源数量变化是指资源中某系数 br 发生变化,即 br′=br+Δ br。并假设规划问题的其他系数都不变。 这样使最终表中原问题的解相应地变化为 XB′=B-1(b+Δ b) 这里 Δ b=(0,… , Δ br,0,… , 0)T 。只要 XB′≥0 , 因最终表中检验数不变,故最优基不变,但最优 解的值发生了变化,所以 XB′ 为新的最优解。新 的最优解的值可允许变化范围用以下方法确定。
(d) (e) -2
· · ·
1 0 0
0 1 0
cj - zj
XB x1 x5 cj - zj
b (f) 4
x1
x2
x3
x4
x5
(g) (h) 0
2 (i) 7
-1 1 (j)
1/2 1/2 (k)
0 1 (l)
--7--
--第2章 对偶问题--
以前讨论线性规划问题时,假定αij,bi,cj都是常数。 但实际上这些系数往往是估计值和预测值。如市场 条件一变,cj值就会变化;αij往往是因工艺条件的 改变而改变;bi是根据资源投入后的经济效果决定 的一种决策选择。显然,当线性规划问题中某一个 或几个系数发生变化后,原来已得结果一般会发生 变化。 因此,所谓灵敏度分析,是指当线性规划问题中的 参数发生变化后,引起最优解如何改变的分析。
运筹学灵敏度分析

只需由 j 0解得c j的范围。
(2) c j 是基变量x j的价格系数 这时要影响所有的检验 数
i ci (c1 ci ci cm ) B Pi , 应由所有的 i 0解得公共的c j。
1
p11-2
例1:在(煤电油例)中,其单纯形终表如下:
0 x 7 x 12 x
3 1
运筹学
2
84 20 24
0 1 0
0
0 0 1
0
1 0 0
0
- 0.32 0.4 - 0.12
- 1.36
1.16 - 0.2 0.16
- 0.52
z 428
(1)甲产品的价格在何范围内变化时,现最优解不变?
解:甲产品的价格c1是基变量的价格系数。 0.32 由 4 0 0 7 c1 12 0.4 2.8 0.4c1 1.44 0 0.12 得 c 3.4, 1.16 由 5 0 0 7 c1 12 - 0.2 1.4 0.2c1 1.92 0 0.16 得 c 2.6,
2
运筹学
例1:在(煤电油例)中,其单纯形终表如下:
0 x 7 x 12 x
3
1
2
84 20 24
0 1
0
0 0 1
1 0
0
- 3.12 1.16 0.4 - 0.2
- 0.12 0.16
z 428
0
0
0
- 1.36
- 0.52
(3)若有人愿以每度1元的价格向该厂供应25度电,是 否值得接受?
§3.4 灵敏度分析
灵敏度分析——研究系数变化对最优解的影响.
运筹(第二章对偶与灵敏度分析)(1)

5x2 3x3 30
x1 0, x2无约束,x3 0
2023/2/22
17
解:将原问题模型变形, 令x1 x1
min z 7x1 4x2 3x3
4x1 2x2 6x3 24
3x1 6x2 4x3 15 5x2 3x3 30
y1 y2 y3
x1 0, x2无约束,x3 0
则对偶问题是
max w 24 y1 15y2 30 y3
4 y1 3y2
7
x1
2 y1 6 y2 5 y3 4
x2
6 y1 4 y2 3x3 3
x3
y1, y2 0, x3无约束
2023/2/22
18
小结:对偶问题与原问题的关系:
目标函数:MAX
原 约束条件:m个约束
对
问
y1 y2
ym
2023/2/22
12
类似于前面的资源定价问题,每一个约束条件对 应一个“ 对偶变量”,它就相当于给各资源的单 位定价。于是我们有如下的对偶规划:
min W b1 y1 b2 y2 bm ym
a11 y1 a21 y2 am1 ym c1 a12y1 a22y2 am2ymc2 a1n y1 a2n y2 amn ym cn y1, y2 ,, ym 0
分别是原问题和对偶问题的可行解,则恒有
n
m
c j x j bi yi
j 1
i 1
m
n
考虑利用 c j aij yi 及
aij x j bi
i 1
j 1
代入。
2、无界性 如果原问题(对偶问题)有无界解,则
其对偶问题(原问题)无可行解。
2023/2/22
运筹学线性规划灵敏度分析教学案例

多个资源系数同时变动分析
例如,将 1 个小时的用工时间从3车间移到2车间,对总利润 产生什么影响?
总利润增加 3650 - 3600 = 50 元, 而目标系数未变,所以最优解肯定 发生变化,
2020/8/1
百分之百法则
如果约束右端值同时变动,计算出每一变动占允许变动量的 的百分比,如果所有的百分比之和不超过100%,那么,影子 价格依然有效;否则,就无法确定。
2020/8/1
灵敏度分析的概念
LP 问题的系数有 aij、bi 、 cj,这些系数往往是估计值 或预测值。
市场条件变化, cj 值就会变化;工艺条件和技术水平改 变, aij 就变化; bi 是根据资源投入后的经济效果决定的一种 选择,市场供应条件发生变化时,亦会改变。
提出问题:
• 当 LP 问题的系数有一个或几个发生变化时,已求得的最优 解会有什么变化; • 这些系数在什么范围内变化时,LP 问题的最优解不会变化。
再改变参数
最优解变了
2020/8/1
那么,保持最优解不变的价值系数允许 变化范围?
改变最优解的临界值是什么呢?
敏感性报告
在“规划求解结果”中 选定“敏感性报告”。 得到一个工作表:
2020/8/1
敏感性报告
最优解
目标函数系数
“递减成本” --- 表示目标函数的系数必须改变多少,才能使 决策变量有正数解。 “允许的增量”和“允许的减量” --- 给出最优解不变的范围。 如门的系数范围: 0≤c1≤750;窗的系数范围:c2≥200
2020/8/1
资源数量变化的分析
考虑只有一个右段值 bi 改变:2 车间可用工时由原来的 12小 时增加到 13 小时,最优解如何变化呢?再变化呢?
运筹学习题解答(chap2)(1)(1)

第二章 对偶问题与灵敏度分析一、写出下列线性规划的对偶问题1、P89,(a)321422m in x x x Z ++=⎪⎪⎩⎪⎪⎨⎧≥=++≤++≥++.,0,;534;332;243321321321321无约束x x x x x x x x x x x x解:原模型可化为321422m in x x x Z ++=⎪⎪⎩⎪⎪⎨⎧≥=++≥≥++.,0,;534;3-3--2-;243321321321321321无约束x x x y y y x x x x x x x x x 于是对偶模型为321532m ax y y y W +-=⎪⎪⎩⎪⎪⎨⎧≥≤+-≤+-≤+-.,0,;4334;243;22321321321321无约束y y y y y y y y y y y y2、P89,(b)321365m ax x x x Z ++=⎪⎪⎩⎪⎪⎨⎧≤≥≤++≥-+-=++.0,0,;8374;35;522321321321321x x x x x x x x x x x x 无约束解:令033≥-='x x 原模型可化为321365m ax x x x Z '-+=⎪⎪⎩⎪⎪⎨⎧≥'≥≤'+≤'='+.0,0,;83-74;3--5-;52-2321321321321321x x x y y y x x x x x x x x x 无约束于是对偶模型为321835m in y y y W +-=⎪⎪⎩⎪⎪⎨⎧≥-≥---≥+-=++.0,,;332;6752;54321321321321y y y y y y y y y y y y 无约束 或⎪⎪⎩⎪⎪⎨⎧≥≤++≥+-=++.0,,;332;6752;54321321321321y y y y y y y y y y y y 无约束二、灵敏度分析1、P92, 线性规划问题213m ax x x Z += ⎪⎩⎪⎨⎧≥≤+≤+0,1025;74212121x x x x x x最优单纯形表如下试用灵敏度分析的方法,分析:(1) 目标函数中的系数21,c c 分别在什么范围内变化,最优解不变(2) 约束条件右端常数项21,b b 分别在什么范围内变化,最优基保持不变解:(1) 1c 的分析:要使得最优解不变,则需⎪⎪⎩⎪⎪⎨⎧≤⨯-⨯+=≤⨯+⨯-=034131003513201413c c σσ 即 ⎪⎩⎪⎨⎧≤≥42511c c 所以:4251≤≤c 时可保持最优解不变。
运筹学第二章24灵敏度分析

(3)其他情况讨论: 某个产品工艺参数改变; 新品代替原产品等;
(2) N =?
舍弃中间计算过程
只考察初始表和最终表
B-1 = AB-1
2、价值系数C发生变化的情况: (1)当cj是非基变量的价值系数——它的变 化只影响 j 一个检验数。 ≤0 1 j c j CB B Pj ≥0 要进行基变换码?
j c j c j CB B Pj ≤ 0
' 1
c j ≤ CB B1 Pj c j
非基变量的价格系数变化,在原最优解 不变的条件下,确定的变化范围。
( 2 )当cj是基变量的价值系数 —— 它的变化 将影响所有非基变量的检验数. 1 N C N CB B N 当cj变化时,如能保持 0 ,则当前解仍 N 为最优解,否则可用单纯形法继续迭代求出 新的最优解。 1 C C B N 0 将cj看作待定参数,令 N N B 解这n-m个不等式,可算出保持最优解不变 时cj的变化范围 ! 基变量的系数变化,仍用c2代表x2的价值系 数(看成待定参数),原最优表格即为:
(2) 增加1个约束条件: 相当于系数阵A增加1行 首先将原最优解代入新增约束检查是 否满足?是,则说明新增约束不影响最 优解。否则再作下面的讨论:
将新增约束标准化,添加到原最优表 格中(相当于约束矩阵新增1行);
进行规格化处理 —— 用矩阵的行变换 将当前基变成单位阵; 用适当方法(通常是对偶单纯形法) 进行迭代求出新的最优解。
(1)增加1个新变量:相当于系数阵A增加1列 如开发出一种新产品,已知其有关工艺参数 (或消耗的资源量)和单位产品利润,设该种 产 品 的 产 量 为 xk , 则 ck 和 Pk 已 知 , 需 要 进 行 “是否投产”的决策。
《运筹学》胡运权第4版线性规划的对偶理论及灵敏度分析省名师优质课赛课获奖课件市赛课一等奖课件

13
2
y3
2 3
题
y1符号不限, y 2 0, y3 0
非 对 偶 形 式 旳 原对 偶 问 题
例2-4 写出下列问题旳对偶问题
max z c1x1 c2 x2 c3x3
a11x a12 x a13x3 b1
s.t.
a21x1 a31x1
a22 x2 a32 x2
a23 x3 a33 x3
出让自己旳资源?
问 题 旳 导 出
例2-1
条件:出让代价应不低于用同等数量资源由自己组织生 产活动时获取旳获利。
y1,y2,y3分别代表单位时间(h)设备A、设备B和调试工 序旳出让代价。 y1,y2,y3旳取值应满足:
6y 2
y 3
2
5y 1
2y 2
y 3
1
美佳企业用6h设备B和1h调试可 生产一件家电I,获利2元
y1, y2 , y3 0
LP1和LP2两个线性规划问题,一般称LP1为原问题, LP2为前者旳对偶问题。
max Z c1x1 c2 x2 cn xn
对 偶 问 题
s.t.
a11 a21
am1
a12 a22
am2
a1n x1 b1
a2n
x2
b2
amn xn bm
规 划 问
minW b1 y1 b2 y2 bm ym
a11 y1 a21 y2 am1 ym (, )c1
a12y1
a22 y2
am2
ym
(,
)c2
题 旳 对 偶 问
a1n y1 a2n y2 amn ym (, )cn
题
y j 0(符号不限,或 0)i 1 ~ m
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
目标式系数: 指目标函数中的系数 允许增量、允许减量:表示目标函数中的系数在允许的增
量与减量范围内变化时,原问题的最优解不变。
450和1E+30的含义是什么?
2.2.2 图解法
x2
8 7 6 5 4 3 2 1
0<=c1<=750
(2,6)是最优解
2.4.2 图解法——改变车间2的约束
x2 改变车间1的约束又会是如何的?
2x2=18
8 7 6 5 4 3 2 1
(2,6)是最优解
2x2=12
可行域
2x2=6
1 23 4 5 6 7 8
x1
2.5 多个约束右端值同时变化的灵敏度分析
分析1小时的工时从车间3移到车间2,对总利润所产生的 影响。 那么,根据影子价格,可知总利润变化量如下: 车间2: 12-->13,利润增加?元 车间3: 18-->17,利润减
课本P50,例2.3,回答五个问题
1. 产品甲的单位利润将会在3.8万元~5.2万元之间波动,公司该 如何应对这种情况,提前对生产格局做好调整预案?
2. 当资源A的限额(储备量)在42~46之间变化时,对线性规划 的影响? 3. 材料B在最优生产格局中出现了12.5单位的剩余,那么应如何 重新制定限额,做好节约工作? 4. 若公司停止生产,把各种原材料变卖。该如何决策?
max z 300 x1 500 x2 x1 4 2 x 12 2 s.t. 3 x1 2 x2 18 x1 , x2 0
现从另一角度提出问题。假定某A公司想把该工厂的资源收购过 来,它至少应付出多大代价,才能使该工厂愿意放弃生产活动, 出让自己的资源?显然该工厂愿意出让自己资源的条件是:出让 代价应不低于用同等数量资源由自己组织生产活动时获取的赢利。 设分别用y1、y2、y3代表单位时间车间1、车间2、车间3的出让代 价,因该工厂用1小时车间1和3小时车间3可生产1扇门,赢利300 元;分别用2小时车间2和车间3可生产1扇窗,赢利500元,由此, y1、y2、y3的取值应满足: y1 + 3y3 ≥ 300 2y2 + 2y3 ≥ 500
该问题为原问题的对偶问题
该问题的最优解的意义: 代表资源最优利用条件下对单位第i种资源的估价。这种 估价不是资源的市场价格,而是根据资源在生产中作出的贡 献而做的估价,为区别起见,称为影子价格(阴影价格)
影子价格
资源的市场价格是其价值的客观体现,相对比较
稳定,而它的影子价格则有赖于资源的利用情况, 是未知数。
第二章 线性规划的
灵敏度分析
在根据一定数据求得最优解后,当这些数据中某 一个或某几个发生变化时,对最优解会产生什么影响。 或者说,要使最优解保持不变,各个数据可以有多大 的幅度的变动。这种研究线性规划模型的原始数据变 化对最优解产生的影响就叫做线性规划的灵敏度分析。
回答两个问题:
①这些系数在什么范围内发生变化时,最优基不变(即最 优解或最优解结构不变)? ②系数变化超出上述范围时,如何用最简便的方法求出新 的最优解?
c1=0(z=0x1+500x2)
可行域 c1=300(z=300x1+500x2)
1 23 4 5 6 7 8
x1
c1=750(z=750x1+500x2)
2.3 多个目标函数系数同时变化
多个系数 c j 发生变化,即其他条件均不变,把300改成450, 把500改成400
max z 300 x1 500 x2 x1 4 2 x 12 2 s.t. 3 x1 2 x2 18 x1 , x2 0
2.4 单个约束右端值变化的灵敏度分析
只有一个约束右端值 bi 发生变化,即如果其他条
件均不变,把12改成13
max z 300 x1 500 x2 x1 4 2 x 13 2 s.t. 3 x1 2 x2 18 x1 , x2 0
max z 300 x1 500 x2 x1 4 2 x 12 2 s.t. 3 x1 2 x2 18 x1 , x2 0
A公司希望用最小的代价把该工厂的全部资源收买过来,故有
min z = 4y1 + 12 y2 + 18y3
显然,yi ≥ 0 (i=1,2,3)
综上,现在的问题为:
min z 4 y1 12 y2 18 y3 y1 3 y3 300 s.t. 2 y2 2 y3 500 y , y , y 0 1 2 3
规划求解得到
2.7 增加一个新变量
在例1.1中,如果增加一个变量x3,比如考虑增加一种新产品: 防盗门,单位利润400元,生产一扇防盗门会占用车间1、车间 2、车间3各2、1、1小时,此时,新的线性规划模型变为:
max z 300 x1 500 x2 x1 4 2 x 12 2 s.t. 3 x1 2 x2 18 x1 , x2 0
max z 300 x1 500 x2 x1 4 1.5 x 12 2 s.t. 3 x1 2 x2 18 x1 , x2 , x3 0
max z 300 x1 500 x2 x1 4 2 x 12 2 s.t. 3 x1 2 x2 18 x1 , x2 0
如果所得的变化的百分比总和不超过100%,则最优 解不会改变;如果超过了100%,则不能确定最优解 是否改变,可通过重新运行“规划求解”命令来判 断
敏感性报告——百分之百法则
作用 1. 可用于确定在保持最优解不变的条件下,目标函数系 数的变化范围; 2. 百分之百法则通过将允许的增加量或减少量在各个系 数之间进行分摊,从而可以直接显示出每个系数允许的变化 值; 3. 线性规划求解后,如果将来条件变化,致使目标函数 中一部分或所有系数都发生变化,那么百分之百法则可以直 接表明最初最优解是否保持不变
但此时敏感性报告是否有效?
2.5.1 敏感性报告——百分之百法则
含义:如果约束右端值同时变化,计算每一右端值变化量 占该约束右端值允许变化量的百分比,然后将每个约束右 端值变化的百分比相加。如果所得的变化的百分比总和不
超过100%,那么影子价格依然有效;如果超过了
100%,那就无法确定影子价格是否依然有效,可通过 重新运行“规划求解”命令来判断。
影子价格与线性规划的对偶问题
课本P54
max z 300 x1 500 x2 400 x3 x1 2 x3 4 2 x x 12 2 3 s.t. 3 x1 2 x2 x3 18 x1 , x2 , x3 0
规划求解得到
2.8 增加一个约束条件
增加一个约束条件,比如增加电力供应限制时, 最优解是否会发生变化? 假设生产一扇门和窗需要消耗电力分别为20kw和 10kw,工厂可供电量最多为90kw,此时应该在原 有的模型中加入新的约束条件:
例1.1 下列模型中,对最优值有 影响的因素有哪些?
max z = 300x1 + 500x2 x1 4 s.t. 2x2 12 3x1 + 2x2 18 x1 , x2 0
灵敏度分析的内容
目标函数系数 c j 变化的灵敏度分析
约束右端值 bi 变化的灵敏度分析
(单个变化和多个变化) (单个变化和多个变化)
max z 500 x1 500 x2 x1 4 2 x 12 2 s.t. 3 x1 2 x2 18 x1 , x2 0
用Excel求解
最优解没有改变
2.2.1 敏感性报告
几个基本概念
递减成本: 它的绝对值表示目标函数中决策变量的系数必须
变化的百分比总和为133.34%, 所以最优解变化与否不确定,需通过“规划求解”重新计算
重新“规划求解”得到
最优解改变,最大利润改变
例1.1 百分比总和=100% ——最优解不变
门的单位利润 c1=300 c1=525,占允许增加量的百分比
窗的单位利润 c2=500 c2=350 ,占允许减少量的百分比
最优解改变,利润增加
如果车间2的可用工时不断增加,会出现什么情况?
为什么?
最优解不改变,利润不变
2.4.1 敏感性报告
阴影价格,显示了约束右端值每增加(或减少) 1个单位,目标函数值(或最优值)相应的增加 量(减少量) 可从敏感性报告中获得的信息??
补充说明: 影子价格
影子价格的解释: 1. 线性规划的对偶问题 例1.1的线性规划问题为
例1.1 百分比总和>100% ——最优解变化与否不确定
门的单位利润 c1=300 c1=600,占允许增加量的百分比 窗的单位利润 c2=500 c2=300 ,占允许减少量的百分比
600 300 100% 66.67% 450
500 300 100% 66.67% 300
max z 450 x1 400 x2 x1 4 2 x 12 2 s.t. 3 x1 2 x2 18 x1 , x2 0
最优解没有改变
敏感性报告
不能反映多个目标函数系数改变时带来的影响
敏感性报告——百分之百法则
定义 如果目标函数系数同时变化,计算出每一系数变化量占 该系数允许变化量的百分比,然后将各个系数变化的百分比 相加。
525 300 100% 50% 450
500 350 100% 50% 300
变化的百分比总和为100%, 所以最优解不变
2.6 约束条件系数变化的灵敏度分析