运筹学复习整理(保准管用)
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
1. 简答题
(1) 运筹学的工作步骤
提出和形成问题:即要弄清问题的目标,可能的约束,问题的可控变量以及相关的参数,搜集相关资料;
建立模型:即把问题中可控变量,参数,目标与约束之间的关系用模型表示出来;
求解:用各种手段将模型求解,解可以是最优解,次优解,满意解。
复杂模型的求解需用计算机,解得精度要求可有决策者提出;
解的检验:首先检查求解步骤和程序有无错误,然后检查解是否反映现实问题;
解的控制:通过控制解的变化过程决定对解是否做一定的改变; 解的实施:是指将解用到实际中必须考虑的实际问题,如向实际部门讲清解的用法,在实施中可能产生的问题和修改。
(2)
退化产生原因及解决办法
单纯形法计算中用θ规则确定换出变量时,有时存在两个以上相同的最小比值,这样在下一次迭代中就有一个或几个基变量等于零,这就出现退化解。
勃兰特规则:
1.选取cj-zj >0中下标最小的非基变量xk 为换入变量,即k=min(j |cj-zj >0)
2. 当按θ规则计算存在两个和两个以上最小比值时,选取下标最小的基
变量为换出变量。
(3)对偶问题的经济解释
• 这说明yi 是右端项bi 每增加一个单位对目标函数Z 的贡献。
• 对偶变量 yi 在经济上表示原问题第i 种资源的边际价值。
• 对偶变量的值 yi*所表示的第i 种资源的边际价值,称为影子价值。
∑
∑
====
=
n
j m
i i i j j y b x c Z 1
1
ω
i
i
y b Z =
∂∂
若原问题的价值系数Cj 表示单位产值,则yi 称为影子价格; 若原问题的价值系数Cj 表示单位利润,则yi 称为影子利润。
影子价格不是资源的实际价格,而是资源配置结构的反映,是在其它数据相对稳定的条件下某种资源增加一个单位导致的目标函数值的增量变化。
(4)分枝定界法步骤
a) 先求出整数规划相应的LP(即不考虑整数限制)的最优解, b) 若求得的最优解符合整数要求,则是原IP 的最优解; c) 若不满足整数条件,则任选一个不满足整数条件的变量来构造新的约束,在原可行域中剔除部分非整数解。
d) 然后,再在缩小的可行域中求解新构造的线性规划的最优解,这样通过求解一系列线性规划问题,最终得到原整数规划的最优解。
(5)树的性质
一个无圈的连通图称为树。
1 树至少有两个悬挂点。
2 一个图为树的充要条件是:不含圈,边数比点数少1.
3 一个图为树的充要条件是:连通,边数比点数少1.
4 一个图为树的充要条件是:任两点之间恰有一条链。
2. 建模题
(1)线性规划建模:
)
.(x ,,x ,x b ),(x 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 z max(min)n m
n m m m n n n n n
n 310
21112122112
2222121112121112211≥≥=≤+++≥=≤+++≥=≤++++++=
约束条件
目标函数
(2)目标规划建模:
最好等于:min d - -d +
最好不大于:min d + 最好不小于:min d -
目标的重要程度不同,用优先等级因子P k 来表示第k 等级目标。
优先等级因子P k 是正的常数,P k >> P k+1 。
同一优先等级下的目标的相对重要性,赋以不同的加权系数w
(3) 整数规划模型
Max (min) Z = Σcjxj s.t. Σaijxj ≤ bi(i=1,2,…m) xj ≥ 0 且部分或全部是整数
n
j d d x K
k E d d x c m
i b x a
d w d w
P Z k k j k
n
j k k j kj i
n
j j ij
l kl L l l kl K
k k
k
,...,2,10
,,,...,2,1,...,2,1),()
(min *
1
11
1
=≥==-+==≥≤+=+
-
=+
-
=+
+=--=∑
∑∑∑非负性约束
目标约束绝对约束
3. 证明题
(1)证明可行域为凸集
为了证明满足线性规划问题的约束条件
的所有点(可行解)组成的集合是凸集,只要证明D 中任意两点连线上的点必然在D 内即可。
设
是D 内的任意两点;X(1)≠X(2)。
(2)证明无界解的判定
构造一个新的解 X (1),它的分量为
因 σm+k >0,所以对任意的λ>0都是可行解,把x(1)代入目标函数内得
∑==≥=n
j j j
j
n
j x b x
P 1
,,2,1,0,
()()()()()()()()()()
T
n
T
n x x x X x x x X 22
2
21
2
112111,,,,,, ==则有
()
()
()
()
∑∑===≥==≥=n
j j j j
n
j j j
j
n
j x b x
P n
j x b x
P 1
22111,,2,1,0,,,2,1,0,
令X=(x 1,x 2,…,x n )T 为x (1),x (2)连线上的任意一点,即
X=αX (1)+(1-α)X (2)
(0≤α≤1) X 的每一个分量是()()
21)1(j j j x x x αα-+=,将它代入约束条件, 得到
()
()()
[
]
()
()
()
b
b b b x
P x
P x
P x
x P x
P n
j n
j j
j j
j n
j j
j n
j n
j j
j
j j
j =-+=-+=--=
∑∑
∑
∑∑
=====ααα
α
αα1
1
221
11
1
211又因()()01,0,0,21>->≥ααj j x x ,所以x j ≥0,j=1,2,…,n 。
由此可见X ∈D ,D 是凸集。
证毕。
()
()
()()k
m j n m j x x a b x j k
m k
m i i i
+≠+===
>-
=++并且,,,1;0011','
1 λ
λλ
z=z0+λσm+k ;
因σm+k >0,故当λ→+∞,则z →+∞,故该问题目标函数无界。
(3)证明弱对偶性
4. 计算题
(1) 标准型,单纯行法计算
∑∑∑==++=+++++-
-
-
---------→-m
i in
i
n m
i m i i
m m
i i
i
m
mn
m m m
m
m
n m n m n m m B B i
n m m j a
c c a
c c b
c z a a b x c a a b x c a a b x c x x x x b X C c c c c c 1
1
1
,11
1
,2
21
,2222111,111111110
01
0001
θθθθ.
.
,0
;;min :,证毕于是得到
得到右乘上式将所以满足
是对偶问题的可行解,因原问题的对偶问题是
左乘上式,得到
将是对偶问题的可行解,
若即
以满足约束条件
是原问题的可行解,所
因设原问题是b
Y X A Y X C C X A Y X C
A Y Y Y C YA Yb b
Y X A Y Y Y b
X A X 0
X b;AX CX;z max ≤≤≥≥≥≥=≤≤≥≤=ω
(2) 对偶型计算
原问题(LP )
对偶问题
0,,,max
211
212
1112112211≥⎪⎪⎪⎭
⎫
⎝⎛=⎪⎪
⎪⎪
⎪⎭
⎫
⎝
⎛⎪⎪⎪⎭⎫
⎝⎛+++=n m
n mn
m m n n
n x x x b b x x x a a a a a a x c x c x c z
()()0
,,,,,,,,,min 21212
1
11211
21m 2211≥≥⎪⎪⎪
⎭
⎫
⎝⎛+++=n n mn m m n m m
y y y c c c a a a a a a y y y b y b y b y
ω
(3)运输问题
1.初始解确定
最小元素法:
–从单位运价表中逐次挑选最小元素,安排运量min{a i,b j}。
–然后,划去该元素所在行或列:
•当产大于销,划去该元素所在列;
•当产小于销,划去该元素所在行。
伏格尔法:
第一步:求出每行次小运价与最小运价之差,记为ui,i=1,2,…,m ;同时求出每列次小运价与最小运价之差,记为vj,j=1,2,…,n ;
第二步:找出所有行、列差额的最大值,即L=max{ui,vi},差额L对应行或列的最小运价处优先调运;
第三步:这时必有一列或一行调运完毕,在剩下的运价中再求最大差额,进行第二次调运,依次进行下去,直到最后全部调运完毕,就得到一个初始调运方案。
2.判定最优解
闭合回路法:
从每一空格出发找一条闭回路。
它是以某空格为起点,用水平或垂直线向前划,当碰到一数字格时可以转90°后,继续前进,直到回到起始空格为止。
当检验数还存在负数时,说明原方案不是最优解。
位势法:
变量的检验数σij=cij –ui –vj=0, 即cij =ui +vj ,且令u1 =0,计算位势量ui 和vj 计算非基变量的检验数σij = c ij –u i – v j 3. 改进方法
确定进基变量
a) 检查非基变量xij 的检验数σij ,按 min{σij| σij <0}= σlk 确定xlk 进基。
确定离基变量
b) 非基变量xlk 进基之后,能让它的运量增加多少呢?
就要求它所在行和列的运量保持产销平衡。
保持产销平衡的方法是闭回路法。
c) 闭回路法:以进基变量xlk 所在格为始点和终点,其余顶点均为基变
量的封闭回路。
d) 闭回路的画法:从进基变量xlk 所在格开始,用水平或垂直线向前划,
每碰到一个基变量格转90º,继续前进,直到返回始点。
e) 奇偶点: 始点是偶点,依次奇偶相间标注;偶点标“+” ,表示运量
增加量;奇点标“-” ,表示运量减少量。
f) 调整量:最小可减少的运量,即奇点运量的最小值。
奇点运量的最小值所在格的基变量离基。
4. 产销不平衡
产大于销:只要增加一个假想的销地j=n+1(实际上是储存),该销地总需
要量为
而在单位运价表中从各产地到假想销地的单位运价为
销大于产:可以在产销平衡表中增加一个假想的产地i=m+1,该地产量为 在单位运价表上令从该假想产地到各销地的运价为
(4) 指派问题
1. 各行各列出现零元素
a. 每行元素减去该行最小元素
b. 每列元素减去该行最小元素 2. 进行试指派,寻求最优解
a.
给只有一个0元素的行(列)的0加圈,然后划去0元素所在行的其他0元素,记作φ
b. 加圈0元素数目等于矩阵的阶数,则指派问题达到最优解
3. 做最少的直线覆盖所有0元素,已确定该系数矩阵中能找到最多的独立元素数
a. 对没有圈的行打√号
b. 对已打√的行所有含φ元素的列打√
c. 在对打有√的列中含圈的元素打√
d. 对没有打√的行及打√的列画一纵线,这就是覆盖所有0元素的最少直线
e.
∑∑==-n
j j
m i i b a 110
;
1,=+n i c ∑∑==-n j m i j j a b 11
;
,1=+j m c
(5) 最小支撑树
破圈法:任取一圈,去掉权重最大的边,重复进行直到无圈可破。
避圈法:选取权重最小的边,重复进行直到形成部分树,并保证不构成圈。
(6) 最短路径问题 Dijkstra 算法:
• S: 已确定最短路的(即具有P 标号)的节点集合。
• P :最短路径长度信息; • T: 目前路径长度信息。
•
λ: 相关长度的路径信息
(7) 最大流问题
标号过程:
(1)给vs 标号(-,+∞),vs 成为已标号未检查的点,其余都是未标号点。
(2)取一个已标号未检查的点vi ,对一切未标号点vj :若有非饱和弧(vi,vj),则vj 标号(vi,l(vj)),其中l(vj)=min[l(vi),cij-fij],vj 成为已标号未检查的点;若有非零弧(vj,vi),则vj 标号(-vi,l(vj)),其中l(vj)=min[l(vi),fji],vj 成为已标号未检查的点。
vi 成为已标号已检查的点。
(3)重复步骤(2),直到vt 成为标号点或所有标号点都检查过。
若vt 成为标号点,表明得到一条vs 到vt 的增广链,转入调整过程;若所有标号点都检查过,表明这时的可行流就是最大流,算法结束。
(4) 调整过程:在增广链上,前向弧流量增加l(vt),后向弧流量减少l(vt)。
截集:
给定容量网络D =(V,A,C),若点集V 被剖分为两个非空集合V1和V2,使 vs ∈V1 ,vt ∈V2,则把弧集(V1,V2)称为(分离vs 和vt 的)截集。
显然,若把某一截集的弧从网络中去掉,则从v s 到v t 便不存在路。
所以,直观上说,截集是从v s 到v t 的必经之路。
截集的容量(简称截量) 最小截集。
(8) 网络计划
工序最早开始时间: 工序最早可能开始时间。
T ES (i,j)=T E (i) =max{T E (i)+T(i, j)} 标号的上界
{}
ij i j j w v P v T v T +=)(),(min )(
工序最迟结束时间:工序最迟必须结束的时间。
T LF(i,j)=T L(j) =min{T L(j)-T(i, j)} 标号的下界
工序最早结束时间:工序最早可能结束的时间。
T EF(i,j)=T ES(i,j)+T(i,j) 标号上界+持续时间
工序最迟开始时间:工序最迟必须开始的时间。
TLS(i,j)=TLF(j)-T(i,j) 标号下界-持续时间
•工序总时差:
•TE(i,j)= TLF(i,j)-TEF(i,j)= TLS(i,j)-TES(i,j)
TE(i,j)=0,则为关键工序
•工序自由时差:
FF(i,j)= ES(j,k)-EF(i,j) 标号上界+持续时间-标号下界(紧后)紧后工序最早开始时间-工序最早完成时间。