1.线性规划及单纯形法(第二部分)
合集下载
第二章单纯形法
s.t.
32x1x133x2x22xx33
100 120
x1, x2, x3 0
CB XB b 0 x4 100 0 x5 120
OBJ = 0 zj cj-zj
40 45 24 0 0
x1 x2 x3 x4 x5
23110 33201 00000 40 45 24 0 0
求解过程：
序
40 45 24 0 0
基本步骤：
始
确定初试基础可行解
检查是否为
是
最优解？
否 确定改善方向
求新的基础可行解
求最优解的目标函数值
1、初始基本可行解的确定
对目标函数为(MAX≤)形式的线性规划背景模型，通过标准化， 每一个约束方程引入一个松弛变量，松弛变量为基变量，其 他变量为非基变量，得到一个初始基本可行解。
n
max f (x) cj xj j 1
cB
p
' j
c ia
' ij
,
j
m
1,
n;
机会成本
i1
j c j z j, j m 1, n ;
检验数
线性规划问题的典式展开式：
max z cBB1b (cm1 cB pm' 1)xm1 (cj cB p'j )xj
(cn cB pn' )xn
x1 x2
a x ' 1,m1 m1
a1', j xj
a' 1,n
xn
b1'
a x ' 2,m1 m1a2' , xja' 2,n
xn
b2'
xm
a x ' m,m1 m1
第二章线性规划及单纯形法总结
第一章
工厂需要的原棉存放在三个仓库中，现将原棉运往工 厂以满足工厂生产的需求。已知原棉运到各个工厂的单位 运费如表所示。问使总运费最小的运输方案？
仓库\工厂
1 2 3 需求
1
2 2 3 40
2
1 2 4 15
3
3 4 2 35
库存
50 30 10
2．线性规划数学模型
解：设xij为i 仓库运到 j工厂的原棉数量(i =1,2,3
1．线性规划介绍
第一章
线性规划研究的主要问题： 有一定的人力、财力、资源条件下，如何 合理安排使用，效益最高?
某项任务确定后，如何安排人、财、物， 使之最省?
2．线性规划数学模型
第一章
例1 美佳公司计划制造I，II两种家电产品。已知各 制造一件时分别占用的设备A、B的台时、调试时间及A、 B设备和调试工序每天可用于这两种家电的能力、各售出 一件时的获利情况如表I—l所示。问该公司应制造A、B两 种家电各多少件，使获取的利润为最大？
第一章
j =1,2,3)
minZ= 2x11 + x12+3x13+2x21 +2x22 +4x23 +3x31 +4x32 +2x33 x11 +x12+x13 x21+x22+x23 x31+x32+x33 50 30 10 40
st.
x11 +x21+x31 =
x12 +x22+x32 =
x13 +x23+x33 = xij 0
15
35
2．线性规划数学模型
第一章
练习4 连续投资10万元 A：从第1年到第4年每年初投资，次年末回收本利1.15； B：第3年初投资，到第5年末回收本利1.25，最大投资4万元； C：第2年初投资，到第5年末回收本利1.40，最大投资3万元； D：每年初投资，每年末回收本利1.11。 求：使5年末总资本最大的投资方案。 分析： A 1 x1A 2 x2A x2C x1D x2D x3D x4D x5D 3 x3A 4 x4A 5
第二章 图解法与单纯形法
表1-4 XB
基变量 x1 x2
进基列 x3
bi /ai2，ai2>0 x4 b
将3化为1
(1)
θi 40 10
出 基 行
x3
x4
2
1 3
1
3 4
1
0 0
0
1 0
40
30
σj
x3
乘 以 1/3 后 得 到
5/3
0 1 0 0 1
1 0 0 3/5 -1/5
-1/3 1/3 -4/3 -1/5 2/5
x2
40
例题
2 x1 x2 40 x1 1.5x2 30
(15,10)
max Z 3x1 4x2 2 x1 x2 40
30
x1 1.5 x2 30 x1 0, x2 0
20
最优解X=(15,10) 最优值Z=85
10
O
10
20
30
40
x1
2.1 线性规划问题的图解法
θ M 20
0 λj
0 2 λj 1 2 λj
x5
x4 x2 x1 x2
1/3 1
3 1/3 1/3 1 0 0
1 2
0 1 0 0
5 1
17 5 －9 17/3
0 0
1 0 0 1/3
1 0
3 1 －2 1
20
75 20 25
25 60
1 0
28/9 －1/9 2/3 －98/9 －1/9 －7/3
1.通过图解法了解线性规划有几种解的形式 2.作图的关键有三点 (1)可行解区域要画正确 (2)目标函数增加的方向不能画错 (3)目标函数的直线怎样平行移动
线性规划单纯形法(清华2)
增加单位产品甲（x1）比乙对目标函数 的贡献大（检验数最大），把非基变量 x1换成基变量，称x1为进基变量，而把 基变量x4换成非基变量，称x4为出基变 量。 （在选择出基变量时，一定保证消去系 统为正消去系统）（最小比值原则）
增加单位产品甲（x1）比乙对目标函数 的贡献大（检验数最大），把非基变量 x1换成基变量，称x1为进基变量，而把 基变量x4换成非基变量，称x4为出基变 量。 （在选择出基变量时，一定保证消去系 统为正消去系统）（最小比值原则）
因为B为一个基, det(B)<>0
有 XB = B-1b- B-1N XN
S = CB B-1b + (CN- CB B-1N ) XN
令非基变量XN = 0 则
Xt = (XB , XN) =（ B-1b , 0）为基础解， 其目标函数值为 S = CB B-1b 只要XB = B-1b 0, Xt =（ B-1b , 0） 0
X为基础可行解, B就是可行基。
另外，若满足 CN- CB B-1N 0 则对任意的 x 0 有 S = CX CB B-1b
即对应可行基B的可行解x为最优解。
定理1-5（最优解判别准则）
对于可行基B ,若
C - CB B-1A 0
则对应于基B的基础可行解x就是基础最 优解，此时的可行基就是最优基。 C - CB B-1A为检验数。 由于基变量的检验数：CB - CB B-1B = 0
c2
Ct= …… cn X=
x2
0= …… xn
0
….. 0
并且
r(A)=m<n.
1.最优解判别定理：
不妨假设 A=（B , N）（B为一个基）
相应地有 Xt= (XB , XN)
单纯形法(第三章线性规划2)
-f 3 –6M -1+M -1+3M 0
0 -M -M x5 x6 x7 B-1b
0 0 0 11 -1 1 0 3 001 1 -M 0 0 4M
3 -1 -1 0 0 -M -M
xj x1 x2 x3 x4 x5 x6 x7 B-1b
0 x4 -M x6 -1 x3
3 -2 0 1 0100 -2 0 1 0
0 1 0 0 0.5 12
40 0 0 0 -25 -600
6/1=6 36/3=12 __
第二步迭代
40 50 0 0 0
xj
基变量
x1 x2
x3 x4 x5
b
40 x1
1 0 1 0 -1 6
0 x4 0 0 -3 1 2 18
50 x2
0 1 0 0 0.5 12
0 0 -40 0 15 -840
f 428 1.36 x4 0.52 x5
X 3 (20 24 84 0 0)T 目标函数值 f 3 = 428。
X3为最优解
即当A产品生产20kg，B产品生产24kg，工厂才能获得最大利 润428百元。x3=84代表煤的剩余量为84t，x4 = x5 = 0表示电力 和劳动日完全利用，没有剩余。
2.单纯形法的主要步骤
Step1. 标准化，找初始基可行解，建立初始的单纯形表；
对于(max , )，松弛变量对应的列构成一个单位阵 Step2.检验当前基可行解是否为最优解
所有检验数 λj 0，则得到最优解（若存在λk ＞0，且pk 0,则该问题
无最优解，停止计算） 否则进行下一步。
Step3.换基迭代（改进基可行解）
例2 用单纯形法求解下列LP问题
第2章线性规划建模及其单纯形法
目标函数 Max z =1500x1+2500x2 约束条件 s.t. 3x1+2x2≤65
2x1+x2≤40 3x2≤75 x1 ,x2 ≥ 0
7
这是一个典型的利润最大化的生产计划问题。 其中，“Max”是英文单词“Maximize”的缩写，含 义为“最大化”； “s.t.”是“subject to”的缩写，表示“满足于…”。 因此，上述模型的含义是：在给定条件限制下，求 使目标函数z达到最大的x1 ,x2的取值
a21x1 + a22x...2 + … + a2nxn≤ b2 am1x1 + am2x2 + … + amnxn≤ bm
x1 , x2 , … , xn≥0
20
•标准形式 •目标函数： Max z=c1x1 + c2x2 + … + cnxn
•约束条件： a11x1 + a12x2 + … + a1nxn=b1 a21x1 + a22x...2 + … + a2nxn=b2 am1x1 + am2x2 + … + amnxn=bm x1 , x2 , … , xn≥0
4 x1 + 2 x2 + 3 x3 - 9 x4≥ 39 6 x2 + 2 x3 + 3 x4≤-58 x1 , x3 , x4 ≥0
31
解8x3:+7首x4先；，将目标函数转换成极大化:令z=-f=3x1–5x2– 其次考虑约束，有3个不等式约束，引进松弛变量 x5 ,x6 ,x7 ≥0 ； 由于x2无非负限制，可令x2=x2’-x2”,其中x2’≥0 x2”≥0 由于第3个约束右端项系数为-58，于是把该式两端乘 以-1。 于是，我们可以得到以下标准形式的线性规划问题：
2x1+x2≤40 3x2≤75 x1 ,x2 ≥ 0
7
这是一个典型的利润最大化的生产计划问题。 其中，“Max”是英文单词“Maximize”的缩写，含 义为“最大化”； “s.t.”是“subject to”的缩写，表示“满足于…”。 因此，上述模型的含义是：在给定条件限制下，求 使目标函数z达到最大的x1 ,x2的取值
a21x1 + a22x...2 + … + a2nxn≤ b2 am1x1 + am2x2 + … + amnxn≤ bm
x1 , x2 , … , xn≥0
20
•标准形式 •目标函数： Max z=c1x1 + c2x2 + … + cnxn
•约束条件： a11x1 + a12x2 + … + a1nxn=b1 a21x1 + a22x...2 + … + a2nxn=b2 am1x1 + am2x2 + … + amnxn=bm x1 , x2 , … , xn≥0
4 x1 + 2 x2 + 3 x3 - 9 x4≥ 39 6 x2 + 2 x3 + 3 x4≤-58 x1 , x3 , x4 ≥0
31
解8x3:+7首x4先；，将目标函数转换成极大化:令z=-f=3x1–5x2– 其次考虑约束，有3个不等式约束，引进松弛变量 x5 ,x6 ,x7 ≥0 ； 由于x2无非负限制，可令x2=x2’-x2”,其中x2’≥0 x2”≥0 由于第3个约束右端项系数为-58，于是把该式两端乘 以-1。 于是，我们可以得到以下标准形式的线性规划问题：
1-3 单纯形法第2部分
（2）类型二：目标要求是“Max”，约束条 件是“=”类型——左边引入非负的人工变量， 并将引入的人工变量作为初始基变量，则初 始可行基是一个单位阵，然后用大M法或两 阶段法求解。
（3）类型三：目标要求是“Max”，约束条 件是“≥”类型——约束条件标准化，左边 减去非负的剩余变量，变成等式约束，化为 类型二。
nm
n
m
n
m
n
Z c j x j c j x j cni xni c j x j cni (bi' ai'j x j )
j1
j1
i1
j1
i1
j1
m
n
mn
cnibi' c j x j
cni ai'j x j
i1
j1
i1 j1
m
n
m
cnibi' (c j cniai'j )x j
（注意：用非基变量表示基变量的表达式）
讨论
①如果限制条件中既有“≤”类型的约束， 又有“≥”或“=”类型的约束，怎么办？
构造“不完全的人造基”！
②为什么初始可行基一定要选单位阵？
b列正好就是基变量的取值，检验数行
和b列交叉处元素也正好对应目标函数值，
因此称b列为解答列
（2）写出初始基本可行解——
用表格单纯形法求解该问题，其过程如下：
CB XB
Cj -1 -2 0 0 -M j
b xj x1 x2 x3 x4 x5
0 X4 3
10010
-M X5 2
-1 2 -1 0 1 1
-Z 2M -(M+1) 2(M-1) -M 0 0
0 X4 3 -2 X2 1
《运筹学》课堂作业及相应答案解析
第一部分绪论第二部分线性规划与单纯形法1 判断下列说法是否正确：(a)图解法同单纯形法虽然求解的形式不同,但从几何上理解,两者是一致的；(b)线性规划模型中增加一个约束条件,可行域的范围一般将缩小,减少一个约束条件,可行域的范围一般将扩大；(c)线性规划问题的每一个基解对应可行域的一个顶点；(d)如线性规划问题存在可行域,则可行域一定包含坐标的原点；(e)对取值无约束的变量x i,通常令其中,在用单纯形法求得的最优解中有可能同时出现(f)用单纯形法求解标准型的线性规划问题时,与对应的变量都可以被选作换入变量；(g)单纯形法计算中,如不按最小比值原则选取换出变量,则在下一个解中至少有一个基变量的值为负；(h)单纯形法计算中,选取最大正检验数δk对应的变量x k作为换入变量,将使目标函数值得到最快的增长；(i)一旦一个人工变量在迭代中变为非基变量后,则该变量及相应列的数字可以从单纯形表中删除,而不影响计算结果；(j)线性规划问题的任一可行解都可以用全部基可行解的线性组合表示；(k)若x1,x2分别是某一线性规划问题的最优解,则也是该线性规划问题的最优解,其中λ1,λ2可以为任意正的实数；(1)线性规划用两阶段法求解时,第一阶段的目标函数通常写为X ai为人工变量),但也可写为,只要所有k i均为大于零的常数；(m)对一个有n个变量、m个约束的标准型的线性规划问题,其可行域的顶点恰好为个；(n)单纯形法的迭代计算过程是从一个可行解转转换到目标函数值更大的另一个可行解；(o)线性规划问题的可行解如为最优解,则该可行解一定是基可行解；(p)若线性规划问题具有可行解,且其可行域有界,则该线性规划问题最多具有有限个数的最优解；(q)线性规划可行域的某一顶点若其目标函数值优于相邻的所有顶点的目标函数值,则该顶点处的目标函数值达到最优；(r)将线性规划约束条件的“≤”号及“≥”号变换成“＝”号,将使问题的最优目标函数值得到改善；(s)线性规划目标函数中系数最大的变量在最优解中总是取正的值；(t)一个企业利用3种资源生产4种产品,建立线性规划模型求解得到的最优解中,最多只含有3种产品的组合；(u)若线性规划问题的可行域可以伸展到无限,则该问题一定具有无界解；(v)一个线性规划问题求解时的迭代工作量主要取决于变量数的多少,与约束条件的数量关系相对较小。
第2章 线性规划与单纯形法(2)
max z = 20 x1 + 30 x2 3 x1 + 10 x2 ≤ 150, x1 ≤ 30, x1 + x2 ≥ 40, x1 , x2 ≥ 0.
管
max z = 20 x1 + 30 x2 − Mx6 3x1 + 10 x2 + x3 = 150, x1 + x4 = 30, x1 + x2 − x5 + x6 = 40, xi ≥ 0, i = 1, 2,..., 6
管
理
运
筹
学
4
要注意到人工变量是与松弛、剩余变量不同的。 松弛变量、剩余变量它们可以取零值，也可以取 正值，而人工变量只能取零值。一旦人工变量取 正值，那么有人工变量的约束方程和原始的约束 方程就不等价了，这样所求得的解就不是原线性 规划的解了。为了竭尽全力地要求人工变量为零， 我们规定人工变量在目标函数中的系数为－M， M 这里M为任意大的数。这样为了使目标函数实现 最大就必须把人工变量从基变量中换出。如果一 直到最后，人工变量仍不能从基变量中换出，也 就是说人工变量仍不为零，则该问题无可行解。 以下讨论如何解含有人工变量的线性规划问题
• 由于不存在单位矩阵,在第1,2个约束条件加上 一个人工变量x6,x7,并在目标函数中加上-Mx6Mx7得到的线性规划问题：
max f ' = −2 x1 − 3 x2 + 0 x3 + 0 x4 + 0 x5 − Mx6 − Mx7
x1 + x 2 − x 3 + x 6 = 3 5 0, x1 − x 4 + x 7 = 1 2 5, 2 x1 + x 2 + x 5 = 6 0 0, x1 , x 2 , x 3 , x 4 , x 5 , x 6 , x 7 ≥ 0 .
第三章-线性规划的单纯形算法2
x1
这样便结束了求全部最优基本可行解的过程，共得 四个基本最优解：
表1的：（2，3，0，1，0，0，0，0）T 表2的：（3 / 2，2，0，0，1/ 2，0，0，0）T 表3的：（0，1，0，1，0，0，2，0）T 表5的：（0，1/ 2，0，0，1/ 2，0，3 / 2，0）T
小结与复习提要：
x3
9x7 1/ 4 4 60 5 1/ 25 6 9 7 3x7 2 1/ 2 4 90 5 `1/ 50 6 3 7 x6 1 3 6
x j 0( j 17)
易见摄动问题的约束条件Ax=b(ε）中右端 的j 系数与左边
系数x j 相同，这是由b(ε）的构造决定的。
代公式知，应从下面两式中找θ，即：
min
x10 ( ) 1/ 4
x20 ( ) 1/ 2
min
1/
4
4
60 5 1/ 1/ 4
25
6
9
7
2 1/ 2 4 90 5 1/ 50 6 3 7
1/ 2
ε足够小，多项式取值主要取决于ε的较低次幂。
这里，0次项：
0 1/ 4
0 1/ 2
,1次项：1 1/ 4
X2 0 0 0 1 0 1/2 -90 -1/50 3
X3
0
1
0
0
1
0
0
1
0
σj
0 0 0 -3/4 150 -1/50 0
（注XB处只列出 0的系数，XB的取值为对应的 0系数及
j 行与该行中元素积之和。）
k min j 4 0 k 4，如何找l？在k 4这一列，
14 1/ 4 0，24 1/ 2 0 ，ε>0足够小时，由单纯形法迭
这样便结束了求全部最优基本可行解的过程，共得 四个基本最优解：
表1的：（2，3，0，1，0，0，0，0）T 表2的：（3 / 2，2，0，0，1/ 2，0，0，0）T 表3的：（0，1，0，1，0，0，2，0）T 表5的：（0，1/ 2，0，0，1/ 2，0，3 / 2，0）T
小结与复习提要：
x3
9x7 1/ 4 4 60 5 1/ 25 6 9 7 3x7 2 1/ 2 4 90 5 `1/ 50 6 3 7 x6 1 3 6
x j 0( j 17)
易见摄动问题的约束条件Ax=b(ε）中右端 的j 系数与左边
系数x j 相同，这是由b(ε）的构造决定的。
代公式知，应从下面两式中找θ，即：
min
x10 ( ) 1/ 4
x20 ( ) 1/ 2
min
1/
4
4
60 5 1/ 1/ 4
25
6
9
7
2 1/ 2 4 90 5 1/ 50 6 3 7
1/ 2
ε足够小，多项式取值主要取决于ε的较低次幂。
这里，0次项：
0 1/ 4
0 1/ 2
,1次项：1 1/ 4
X2 0 0 0 1 0 1/2 -90 -1/50 3
X3
0
1
0
0
1
0
0
1
0
σj
0 0 0 -3/4 150 -1/50 0
（注XB处只列出 0的系数，XB的取值为对应的 0系数及
j 行与该行中元素积之和。）
k min j 4 0 k 4，如何找l？在k 4这一列，
14 1/ 4 0，24 1/ 2 0 ，ε>0足够小时，由单纯形法迭
第二章 线性规划与单纯形法
max z =2x1+3x2 2x1+2x2+x3=12 x1+2x2+x4=8 4x1+x5=16 4x2+x6=12 x1-6≥0
图解法
第2节 解
例8：求下述线性规划的所有基解、基可行 解及最优解。 max z =3x1+x2+3x3 x1+x2+x3＝2 x1+2x2+4x3＝6 x1，x2，x3≥0
第1节 线性规划问题及其数学模型
要求：将下列线性规划问题转化为标准型。 例4：min z =x1+2x2+3x3 -2x1+x2+x3≤9 -3x1+x2+2x3≥4 3x1-2x2-3x3=-6 x1≤0，x2≥0，x3取值无约束
第1节 线性规划问题及其数学模型
例4： x x , x x xz , z 解：令 1 13 3 3 max z′=x1′-2x2-3x3′+3x3〞+0x4+0x5 2x1′+x2+x3′-x3〞+x4=9 3x1′+x2+2x3′-2x3〞-x5=4 3x1′+2x2+3x3′-3x3〞=6 x1′，x2，x3′，x3〞，x4，x5≥0
第3节 图解法
三、图解法解的类型 唯一最优解：仅有一点使目标函数值取得最大 （小）值 无穷多（多重）最优解：线段（射线）上任意 一点都使目标函数值取得相同的最大（小）值 无界解：可行域无界，目标函数值可以增大到 无穷大 无可行解：可行域为空集 无界解和无可行解统称为无最优解
目标函数求最大值 函数约束条件全为等式 决策变量全为非负 函数约束条件右端项全为非负
图解法
第2节 解
例8：求下述线性规划的所有基解、基可行 解及最优解。 max z =3x1+x2+3x3 x1+x2+x3＝2 x1+2x2+4x3＝6 x1，x2，x3≥0
第1节 线性规划问题及其数学模型
要求：将下列线性规划问题转化为标准型。 例4：min z =x1+2x2+3x3 -2x1+x2+x3≤9 -3x1+x2+2x3≥4 3x1-2x2-3x3=-6 x1≤0，x2≥0，x3取值无约束
第1节 线性规划问题及其数学模型
例4： x x , x x xz , z 解：令 1 13 3 3 max z′=x1′-2x2-3x3′+3x3〞+0x4+0x5 2x1′+x2+x3′-x3〞+x4=9 3x1′+x2+2x3′-2x3〞-x5=4 3x1′+2x2+3x3′-3x3〞=6 x1′，x2，x3′，x3〞，x4，x5≥0
第3节 图解法
三、图解法解的类型 唯一最优解：仅有一点使目标函数值取得最大 （小）值 无穷多（多重）最优解：线段（射线）上任意 一点都使目标函数值取得相同的最大（小）值 无界解：可行域无界，目标函数值可以增大到 无穷大 无可行解：可行域为空集 无界解和无可行解统称为无最优解
目标函数求最大值 函数约束条件全为等式 决策变量全为非负 函数约束条件右端项全为非负
第二章 线性规划与单纯形法14节
2、标准形式的特征？？？
2018/10/11 10
二、 线性规划的标准形
3、线性规划的标准化方法
（1）把最小化目标函数转化为求最大化问题。 （2）约束条件右端项为负时两边同乘以－1 （3）把约束方程中的不等式转化为等式。具体做法是：对于≤的 情况，引进松弛变量，对于≥的情况，引进剩余变量。 （4）将变量中的非正限制或无限制转化为非负限制。其中，对 于无限制变量的处理：一是同时引进两个非负变量，然后用它 们的差代替无限制变量，即令 二是从约束方程 ' " xk x k xk 中任取一个包含无限制变量的等式约束，解出该变量，并把它 代入目标函数和其他约束方程中去，以消除该无限制变量。
2018/10/11
13
小
结
1.什么是线性规划，掌握线性规划在管理中的几个应用 例子 2.线性规划数学模型的组成及其特征 3.线性规划数学模型的一般表达式和标准形式。 4.
返回首页
2018/10/11
图解法
Exit
14
第二节 线性规划的图解法
1.图解法的含义 在直角坐标系中，描绘出约束条件和变量限制的公 共区域，然后通过观察确定符合目标要求的变量的取值。 2.几个概念 （ 1 ）可行解 : 由约束条件和变量取值限制围成的公共 区域中的每一个点都称为线性规划问题的可行解。 （2)可行域:所有可行解的集合，构成线性规划问题的 可行域。 (3)等值线:使目标函数取相等值的所有点的集合，称 为目标函数的等值线。 (4)法向量: 与等值线垂直的向量。分为正法向量和负 2018/10/11 15 法向量。
基：约束系数矩阵A中，m个线性无关的列向量，称为
m维实空间中的一个基。其中，每个列向量称为基向 量，全部基向量构成基矩阵（也可简称为基），剩下 的n-m个列向量称为非基向量，所有的非基向量构成 非基矩阵。
第二章 线性规划及单纯形法
标准形式
目标函数： 目标函数： 约束条件： 约束条件： Max z = c1 x1 + c2 x2 + … + cn xn s.t. a11 x1 + a12 x2 + … + a1n xn = b1 a21 x1 + a22 x2 + … + a2n xn = b2 …… …… am1 x1 + am2 x2 + … + amn xn = bm x1 ，x2 ，… ，xn ≥ 0，bi ≥0 ，
（一）一般式
Max(min)Z=C1X1+ C2X2+…+CnXn a11X1+ a12X2+…+ a1nXn ≥(=, ≤)b1 a21X1+ a22X2+…+ a2nXn ≥(=, ≤)b2 … … … am1X1+ am2X2+…+ amnXn ≥(=, ≤)bm Xj ≥0(j=1,…,n) 0( )
三、线性规划问题的标准形式 线性规划问题的标准形式
2、约束条件不是等式的问题： 约束条件不是等式的问题： 设约束条件为
ai1 x1+ai2 x2+ … +ain xn ≤ bi
可以引进一个新的变量s ，使它等于约束右边与左 边之差
s=bi–(ai1 x1 + ai2 x2 + … + ain xn ) (
一、问题提出
Ⅰ 设备A 设备 设备B 设备 调试工序 利润 0 6 1 2
例1生产计划问题
Ⅱ 5 2 1 1 每天可用能力 15 24 5
两种家电各生产多少, 可获最大利润? 两种家电各生产多少, 可获最大利润
运筹学(第二版)课后答案
式中 x4，x5 为松弛变量，x5 可作为一个基变量，第一、三约束分别加 入人工变量 x6 x7，目标函数中加入-Mx6-Mx7 一项，得到大 M 单纯形法数 学模型
max z 4 x1 5 x 2 x3 3 x1 2 x 2 x3 x 4 x6 18 2 x x x 4 1 2 5 st x1 x 2 x3 x7 5 x j 0, ( j 1, ,7)
４０６
附录四习题参考答案
1 0 1 σj 1 0 1 σj
X6 X1 X7 X6 X2 X7
0 1 0 0 -1 2 -1 2
1/2 1/2 1/2 -1 0 1 0 0
1 0 -1 0 1 0 -1 0
-1 0 0 1 -1 0 0 1
-3/2 1/2 -1/2 2 -2 2 -1 3
1 0 0 0 1 0 0 0
0 0 1 0 0 0 1 0
12 2 3 10 4 1
表 1.4-1.2 在第一阶段的最优解中人工变量不为零，则原问题无可行解。 注：在第二阶段计算时，初始表中的检验数不能引用第一阶段最优表的检 验数，必须换成原问题的检验数。 (2) 无穷多最优解，如 X1＝（4，0，0） ；X2＝（0，0，8） (3) 无界解 (4) 唯一最优解 X*＝（5/2，5/2，5/2，0） (5) 唯一最优解 X*＝（24，33） (6) 唯一最优解 X*＝（14，0，－4） 1.5 （1） X*仍为最优解，max z＝λCX； （2） 除 C 为常数向量外，一般 X*不再是该问题的最优解； （3） 最优解变为λX*，目标函数值不变。 1.6 （1） d≥0,c1＜0, c2＜0 （2） d≥0,c1≤0, c2≤0,但 c1，c2 中至少一个为零 （3） d＝0 或 d＞0，而 c1＞0 且 d/4=3/a2 （4） c1＞0,d/4＞3/a2 （5） c2＞0,a1≤0 （6） x5 为人工变量，且 c1≤0, c2≤0 1.7 解： 设 xj 表示第 j 年生产出来分配用于作战的战斗机数； yj 为第 j 年已 培训出来的驾驶员； （aj－xj）为第 j 年用于培训驾驶员的战斗机数；zj 为第 j 年用于培训驾驶员的战斗机总数。则模型为 max z = nx1+(n-1)x2+…+2xn-1+xn s.t. zj=zj-1+(aj-xj) yj=yj-1+k(aj-xj) x1+x2+…+xj≤yj xj,yj,zj≥0 (j=1,2, …,n) 1.8
max z 4 x1 5 x 2 x3 3 x1 2 x 2 x3 x 4 x6 18 2 x x x 4 1 2 5 st x1 x 2 x3 x7 5 x j 0, ( j 1, ,7)
４０６
附录四习题参考答案
1 0 1 σj 1 0 1 σj
X6 X1 X7 X6 X2 X7
0 1 0 0 -1 2 -1 2
1/2 1/2 1/2 -1 0 1 0 0
1 0 -1 0 1 0 -1 0
-1 0 0 1 -1 0 0 1
-3/2 1/2 -1/2 2 -2 2 -1 3
1 0 0 0 1 0 0 0
0 0 1 0 0 0 1 0
12 2 3 10 4 1
表 1.4-1.2 在第一阶段的最优解中人工变量不为零，则原问题无可行解。 注：在第二阶段计算时，初始表中的检验数不能引用第一阶段最优表的检 验数，必须换成原问题的检验数。 (2) 无穷多最优解，如 X1＝（4，0，0） ；X2＝（0，0，8） (3) 无界解 (4) 唯一最优解 X*＝（5/2，5/2，5/2，0） (5) 唯一最优解 X*＝（24，33） (6) 唯一最优解 X*＝（14，0，－4） 1.5 （1） X*仍为最优解，max z＝λCX； （2） 除 C 为常数向量外，一般 X*不再是该问题的最优解； （3） 最优解变为λX*，目标函数值不变。 1.6 （1） d≥0,c1＜0, c2＜0 （2） d≥0,c1≤0, c2≤0,但 c1，c2 中至少一个为零 （3） d＝0 或 d＞0，而 c1＞0 且 d/4=3/a2 （4） c1＞0,d/4＞3/a2 （5） c2＞0,a1≤0 （6） x5 为人工变量，且 c1≤0, c2≤0 1.7 解： 设 xj 表示第 j 年生产出来分配用于作战的战斗机数； yj 为第 j 年已 培训出来的驾驶员； （aj－xj）为第 j 年用于培训驾驶员的战斗机数；zj 为第 j 年用于培训驾驶员的战斗机总数。则模型为 max z = nx1+(n-1)x2+…+2xn-1+xn s.t. zj=zj-1+(aj-xj) yj=yj-1+k(aj-xj) x1+x2+…+xj≤yj xj,yj,zj≥0 (j=1,2, …,n) 1.8
运筹学第四版第二章线性规划及单纯形法
方案的制定受到那些现实条件制约：
确定约束条件
人力资源（劳动力）的限制： 9x1 4x2 360
设备工时的限制：
4x1 5x2 200
原材料资源的限制：
3x1 10x2 300
此外，决策变量的取值不应为负值即 x1 0, x2 0
6
综上所述，我们得到了这个问题的数学模型
目标函数 约束条件
大？
项目
Ⅰ
设备A (h)
0
设备B (h)
6
调试工序(h) 1
利润（元） 2
Ⅱ 每天可用能力
5
15
2
24
表1-2
1
5
1
12
其数学模型为：
max Z 2x1 x2
5x2 15
6xx11
2x2 x2
24 5
x1, x2 0
13
例3：捷运公司在下一年度的1～4月份的4个月内拟租用仓库
堆放物资。已知各月份所需仓库面积列于下表1-3。仓库租
借费用随合同期而定，期限越长，折扣越大，具体数字见表
1-4。租借仓库的合同每月初都可办理，每份合同具体规定
租用面积和期限。因此该厂可根据需要，在任何一个月初办
理租借合同。每次办理时可签一份合同，也可签若干份租用
面积和租用期限不同的合同。试确定该公司签订租借合同的
最优决策，目的是使所租借费用最少。
14
max Z 70 x1 120 x2
9x1 s.t. 43xx11
x1,
4x2 5x2 10x2 x2 0
360 200 300
资源约束
非负约束
其中 约束条件可记 s.t (subject to), 意思为“以… 为条件“、”假定“、”满足“之意。
第2章 线性规划及单纯形法1-2节
2x1+ x2 400
A
B
最优解 (50, 250)
x2 250
x1 + x2 300
C
100 —
50 —
可 行 域
x1+ x2=300 x2=250
O0 50x1 + 100x2 =0
| | | |D | | | | 50 100 150 200 250 300 350 400
x1
二、线性规划问题解的存在情况：
例5:
Max z =1500x1+2500x2 s.t. 3x1+2x2≤ 65 (A) 2x1+x2≤ 40 (B) 3x2≤ 75 (C)
x1 ,x2 ≥0 (D、E)
B
x2
(5,25)T
A
40
C
25
目标函数 等值线
Z
0
20
Z
x1
存在唯一最优解
例6: 目标函数变为： Max z = 1500 x1 + 1000 x2
线性规划问题的规范形式和标准形式
规范形式：
Max Z =c1x1+c2x2+…+cnxn s.t. a11x1 +a12x2 +…+a1nxn ≤b1 a21x1 +a22x2 +…+a2nxn ≤b2 am1x1 +am2x2 +…+amnxn ≤bm x1 , x2 ,… , xn ≥0
矩阵型式:
§2.3 线性规划的图解法
一、线性规划的图解法
概念
线性规划的图解法(解的几何表示)对于只有 两个变量的线性规划问题，可以在二维直角 坐标平面上作图表示线性规划问题的有关概 念，并求解。 图解法有助于理解LP问题的求解原理。
A
B
最优解 (50, 250)
x2 250
x1 + x2 300
C
100 —
50 —
可 行 域
x1+ x2=300 x2=250
O0 50x1 + 100x2 =0
| | | |D | | | | 50 100 150 200 250 300 350 400
x1
二、线性规划问题解的存在情况：
例5:
Max z =1500x1+2500x2 s.t. 3x1+2x2≤ 65 (A) 2x1+x2≤ 40 (B) 3x2≤ 75 (C)
x1 ,x2 ≥0 (D、E)
B
x2
(5,25)T
A
40
C
25
目标函数 等值线
Z
0
20
Z
x1
存在唯一最优解
例6: 目标函数变为： Max z = 1500 x1 + 1000 x2
线性规划问题的规范形式和标准形式
规范形式：
Max Z =c1x1+c2x2+…+cnxn s.t. a11x1 +a12x2 +…+a1nxn ≤b1 a21x1 +a22x2 +…+a2nxn ≤b2 am1x1 +am2x2 +…+amnxn ≤bm x1 , x2 ,… , xn ≥0
矩阵型式:
§2.3 线性规划的图解法
一、线性规划的图解法
概念
线性规划的图解法(解的几何表示)对于只有 两个变量的线性规划问题，可以在二维直角 坐标平面上作图表示线性规划问题的有关概 念，并求解。 图解法有助于理解LP问题的求解原理。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
MinZ = CX AX = b s.t. X ≥ 0 C = (3,1,1,0,0), X = ( x1 , x2 , x3 , x4 , x5 )T 1 2 1 1 0 11 A = 4 1 2 0 1, b = 3 2 0 1 1 0 0
添加人工变量x6, x8—大M法
凸集
极点
凸集
不是凸集
3.2 几个基本定理
定理1 若线性规划问题存在可行解,则问题的可行域 是凸集. 引理 线性规划问题的可行解X=(x1, … , x1)T为基可行 解的充要条件是X的正分量所对应的系数列向量是线 性独立的. 定理2 线性规划总理的基可行解X对应线性规划问题 可行域(凸集)的顶点. 定理3 若线性规划问题有最优解,一定存在一个解是 最优解.
MinZ = 3 x1 + x2 + x3 + 0 × ( x4 + x5 ) + M ( x6 + x8 ) x1 2 x2 + x3 + x4 = 11 4 x + x + 2 x x + x = 3 1 2 3 5 6 s.t. 2 x1 + x3 + x8 = 1 x j ≥ 0 M 可以理解为一个量级很大的正数( + ∞)
初始表—阶段1
C→ 0 0 x2 -2 1 0 -1 0 x3 1 2 1 -3 0 x4 1 0 0 0 0 x5 0 -1 0 1 1 x6 0 1 0 0 1 x8 0 0 1 0 θj 11/1 3/2 1/1 CB XB B-1b x1 0 1 1 x4 x6 x8 σj 11 3 1 Z0= 4 1 -4 -2 6
B
可行域OABCD
C:最优解
2 x2 1 0
0 4 8
-1
O
-2 -3
D
x1
单纯形法的基本步骤
第一步:求出线性规划的初始基可行解 第二步:进行最优性检验--如果表中所有检验数 σk ≤ 0 则为最优解,否则继续. 第三步:从一个基可行解转换到另一个目标函数值更大的基可 行解,列出新的单纯形表. 确定换入基的变量 σ k = max{σ j σ j > 0} 确定换出基的变量 用换入的变量替换基变量中的掏出变量,得到一个新的基 第四步:……
第一阶段求解的两种可能结果
1.
2.
最优时,还有人工变量在最优基变量中:原问题 无可行解 最优时,所有人工变量都是非基变量:进入求解 的第二阶段
第二阶段:目标函数换回原来的目标函数, 在第一阶段得到的最优基基础上,继续求解 (检验数,基转换)
初始表—阶段2
C→ -3 1 x2 0 1 0 0 1 x3 0 0 1 0 0 x4 1 0 0 0 0 x5 -2 -1 0 1 0 x6 2 1 0 -1 0 x8 -5 -2 1 1 θj 12/3 — — CB XB B-1b x1 0 1 1 x4 x2 x3 σj 12 1 1 Z0= 2 3 0 -2 -1
迭代2—阶段1
C→ 0 0 x2 0 1 0 0 0 x3 0 0 1 0 0 x4 1 0 0 0 0 x5 -2 -1 0 0 1 x6 2 1 0 1 1 x8 -5 -2 1 1 θj CB XB B-1b x1 0 0 0 x4 x2 x3 σj 12 1 1 Z0= 0 3 0 -2 0
出基变量选择:最小比值法(保证新的基解还 是可行的)
单纯形表
C→ CB c1 c2 … cm XB x1 x2 … xm σj B-1b b1 b2 … bm Z0 c1 x1 a11 a21 … am1 σ1 c2 x2 a12 a22 … am2 σ2 … … … … … … … cn xn a1n a2n … amn σn θj θ1
3.1 凸集和顶点 凸集:如果集合C中任意两个点x1, x2,其连线上 的点也都是集合C中的点,称C为为凸集. 顶点:如果集合C中不存在任意两个不同的点 x1, x2的点,使x成为其连线上的点.
可行域的性质
● ●
线性规划的可行域是凸集 线性规划的最优解在极点上
σj Z0=2
添加人工变量x6, x8—两阶段法
第一阶段:构筑所有人工变量之和最小化的目标 函数,用单纯形法求解
Minω = x6 + x8 x1 2 x2 + x3 + x4 = 11 4 x + x + 2 x x + x = 3 1 2 3 5 6 s.t. 2 x1 + x3 + x8 = 1 x j ≥ 0
初始表
C→ -3 1 x2 -2 1 0 1 x3 1 2 1 0 x4 1 0 0 0 0 x5 0 -1 0 M M x6 0 1 0 0 M x8 0 0 1 0 θj 11/1 3/2 1/1 CB XB B-1b x1 0 M M x4 x6 x8 σj 11 3 1 Z0= 4M 1 -4 -2
θ = min{ bi b a j > 0} = l a ik a lk
基本结论
线性规划问题的可行域(如果存在)是一个凸 集 该凸集有有限个顶点(基可行解) 线性规划问题的最优解(如果存在),一定在 某个顶点(基可行解)达到
求解的可能结果
唯一最优解 多最优解,即无穷多最优解 无界解,即目标函数值趋近无穷(通常是遗漏 了约束) 无可行解(通常是存在矛盾的约束)
B 可行基 BX B + NX N = b X B = B 1b B 1 N 代入目标函数 Z = C B X B + C N X N = C B B 1b + (C N C B B 1 N ) X N
σ N = C N C B B 1 N = (σ m +1 , σ m + 2 ,..., σ n ) σ j = c j C B B 1 Pj
4.5 4 3.5 3 2.5 x2 2 1.5 1 0.5 0
0 4 8
(4,2)
x1
欲速则不达 进基变量选择讨论
例1种,经过三次顶点转换才达到最优点;从 图示看出,实际可以二次转换就达到 没有更好的指导标准
初始基可行解的确定 人工变量法
对于任意一个线性规划问题,通常难以确定初 始顶点 思路:强行添加人工变量,使得这些人工变量 与模型中原有的某些变量(如松弛变量)的系 数矩阵,能够构成一个单位矩阵 求解:对含有人工变量的模型进行求解,在求 解过程中,逐渐用非人工变量替代基变量中的 人工变量,当人工变量全部出基后,得到的就 是原问题的初始基可行解
迭代2
C→ -3 1 x2 0 1 0 0 1 x3 0 0 1 0 0 x4 1 0 0 0 0 x5 -2 -1 0 1 M x6 2 1 0 M x8 -5 -2 1 θj 12/3 — — CB XB B-1b x1 0 1 1 x4 x2 x3 σj 12 1 1 Z0= 2 3 0 -2 -1
-3 1 1 + - - 6M M 3M
迭代1
C→ -3 1 x2 -2 1 0 1 - M 1 x3 0 0 1 0 0 x4 1 0 0 0 0 x5 0 -1 0 M M x6 0 1 0 0 M x8 -1 -2 1 3M - 1 θj — 1/1 — CB XB B-1b x1 0 M 1 x4 x6 x3 σj 10 1 1 Z0= M +1 3 0 -2 -1
如何让人工变量尽快出基?
大M法
MAX:给每个人工变量一个量级很小的系数(-M) MIN:给每个人工变量一个量级很大的系数(M)
两阶段法
第一阶段:MIN Z =∑人工变量 第二阶段:原来的目标函数
示例
MinZ= 3x1 + x2 + x3 x1 2x2 + x3 ≤ 11 4x + x + 2x ≥ 3 1 2 3 s.t. 2x1 + x3 = 1 x j ≥ 0
0 4 8
(0,0)
x1
最小比值法
2 x2 + x3 = 8 x3 = 8 2 x2 s.t 0 x2 + x4 = 16 x4 = 16 0 x2 4 x + x = 12 x = 12 4 x 2 2 5 5
迭代1
C→ CB 0 0 3 XB x3 x4 x2 σj B-1b 2 16 3 Z0=9 2 x1 1 4 0 2↑ 3 x2 0 0 1 0 0 x3 1 0 0 0 0 x4 0 1 0 0 0 x5 -1/2 0 1/4 -3/4 θj 2/1 → 16/4 —
基可行解
基:Am×n的非奇异子阵Bm×m 基变量与非基变量: Bm×m各列对应的变量为 基变量XB,其余变量为非基变量XN 基解:X=(XB,XN),XB=B-1b, XN =0 基可行解:基解中,XB=B-1b≥0 目标函数值Z0=CBB-1b 可行基:基可行解对应的基
最优性判断—非基变量检验数
迭代1—阶段1
C→ 0 0 x2 -2 1 0 -1 0 x3 0 0 1 0 0 x4 1 0 0 0 0 x5 0 -1 0 1 1 x6 0 1 0 0 1 x8 -1 -2 1 3 θj — 1/1 — CB XB B-1b x1 0 1 0 x4 x6 x3 σj 10 1 1 Z0= 1 3 0 -2 0
初始表
C→ CB 0 0 0 XB x3 x4 x5 σj B-1b 8 16 12 Z0=0 2 x1 1 4 0 2 3 x2 2 0 4 3↑ 0 x3 1 0 0 0 0 x4 0 1 0 0 0 x5 0 0 1 0 θj 8/2 — 12/4 →
4.5 4 3.5 3 2.5 x2 2 1.5 1 0.5 0
M M -1 +1
迭代3
C→ -3 1 x2 0 1 0 0 1 x3 0 0 1 0 0 x4 0 x5 M x6 M x8 θj CB XB B-1b x1 -3 1 1 x1 x2 x3 4 1 9 1 0 0 0
添加人工变量x6, x8—大M法
凸集
极点
凸集
不是凸集
3.2 几个基本定理
定理1 若线性规划问题存在可行解,则问题的可行域 是凸集. 引理 线性规划问题的可行解X=(x1, … , x1)T为基可行 解的充要条件是X的正分量所对应的系数列向量是线 性独立的. 定理2 线性规划总理的基可行解X对应线性规划问题 可行域(凸集)的顶点. 定理3 若线性规划问题有最优解,一定存在一个解是 最优解.
MinZ = 3 x1 + x2 + x3 + 0 × ( x4 + x5 ) + M ( x6 + x8 ) x1 2 x2 + x3 + x4 = 11 4 x + x + 2 x x + x = 3 1 2 3 5 6 s.t. 2 x1 + x3 + x8 = 1 x j ≥ 0 M 可以理解为一个量级很大的正数( + ∞)
初始表—阶段1
C→ 0 0 x2 -2 1 0 -1 0 x3 1 2 1 -3 0 x4 1 0 0 0 0 x5 0 -1 0 1 1 x6 0 1 0 0 1 x8 0 0 1 0 θj 11/1 3/2 1/1 CB XB B-1b x1 0 1 1 x4 x6 x8 σj 11 3 1 Z0= 4 1 -4 -2 6
B
可行域OABCD
C:最优解
2 x2 1 0
0 4 8
-1
O
-2 -3
D
x1
单纯形法的基本步骤
第一步:求出线性规划的初始基可行解 第二步:进行最优性检验--如果表中所有检验数 σk ≤ 0 则为最优解,否则继续. 第三步:从一个基可行解转换到另一个目标函数值更大的基可 行解,列出新的单纯形表. 确定换入基的变量 σ k = max{σ j σ j > 0} 确定换出基的变量 用换入的变量替换基变量中的掏出变量,得到一个新的基 第四步:……
第一阶段求解的两种可能结果
1.
2.
最优时,还有人工变量在最优基变量中:原问题 无可行解 最优时,所有人工变量都是非基变量:进入求解 的第二阶段
第二阶段:目标函数换回原来的目标函数, 在第一阶段得到的最优基基础上,继续求解 (检验数,基转换)
初始表—阶段2
C→ -3 1 x2 0 1 0 0 1 x3 0 0 1 0 0 x4 1 0 0 0 0 x5 -2 -1 0 1 0 x6 2 1 0 -1 0 x8 -5 -2 1 1 θj 12/3 — — CB XB B-1b x1 0 1 1 x4 x2 x3 σj 12 1 1 Z0= 2 3 0 -2 -1
迭代2—阶段1
C→ 0 0 x2 0 1 0 0 0 x3 0 0 1 0 0 x4 1 0 0 0 0 x5 -2 -1 0 0 1 x6 2 1 0 1 1 x8 -5 -2 1 1 θj CB XB B-1b x1 0 0 0 x4 x2 x3 σj 12 1 1 Z0= 0 3 0 -2 0
出基变量选择:最小比值法(保证新的基解还 是可行的)
单纯形表
C→ CB c1 c2 … cm XB x1 x2 … xm σj B-1b b1 b2 … bm Z0 c1 x1 a11 a21 … am1 σ1 c2 x2 a12 a22 … am2 σ2 … … … … … … … cn xn a1n a2n … amn σn θj θ1
3.1 凸集和顶点 凸集:如果集合C中任意两个点x1, x2,其连线上 的点也都是集合C中的点,称C为为凸集. 顶点:如果集合C中不存在任意两个不同的点 x1, x2的点,使x成为其连线上的点.
可行域的性质
● ●
线性规划的可行域是凸集 线性规划的最优解在极点上
σj Z0=2
添加人工变量x6, x8—两阶段法
第一阶段:构筑所有人工变量之和最小化的目标 函数,用单纯形法求解
Minω = x6 + x8 x1 2 x2 + x3 + x4 = 11 4 x + x + 2 x x + x = 3 1 2 3 5 6 s.t. 2 x1 + x3 + x8 = 1 x j ≥ 0
初始表
C→ -3 1 x2 -2 1 0 1 x3 1 2 1 0 x4 1 0 0 0 0 x5 0 -1 0 M M x6 0 1 0 0 M x8 0 0 1 0 θj 11/1 3/2 1/1 CB XB B-1b x1 0 M M x4 x6 x8 σj 11 3 1 Z0= 4M 1 -4 -2
θ = min{ bi b a j > 0} = l a ik a lk
基本结论
线性规划问题的可行域(如果存在)是一个凸 集 该凸集有有限个顶点(基可行解) 线性规划问题的最优解(如果存在),一定在 某个顶点(基可行解)达到
求解的可能结果
唯一最优解 多最优解,即无穷多最优解 无界解,即目标函数值趋近无穷(通常是遗漏 了约束) 无可行解(通常是存在矛盾的约束)
B 可行基 BX B + NX N = b X B = B 1b B 1 N 代入目标函数 Z = C B X B + C N X N = C B B 1b + (C N C B B 1 N ) X N
σ N = C N C B B 1 N = (σ m +1 , σ m + 2 ,..., σ n ) σ j = c j C B B 1 Pj
4.5 4 3.5 3 2.5 x2 2 1.5 1 0.5 0
0 4 8
(4,2)
x1
欲速则不达 进基变量选择讨论
例1种,经过三次顶点转换才达到最优点;从 图示看出,实际可以二次转换就达到 没有更好的指导标准
初始基可行解的确定 人工变量法
对于任意一个线性规划问题,通常难以确定初 始顶点 思路:强行添加人工变量,使得这些人工变量 与模型中原有的某些变量(如松弛变量)的系 数矩阵,能够构成一个单位矩阵 求解:对含有人工变量的模型进行求解,在求 解过程中,逐渐用非人工变量替代基变量中的 人工变量,当人工变量全部出基后,得到的就 是原问题的初始基可行解
迭代2
C→ -3 1 x2 0 1 0 0 1 x3 0 0 1 0 0 x4 1 0 0 0 0 x5 -2 -1 0 1 M x6 2 1 0 M x8 -5 -2 1 θj 12/3 — — CB XB B-1b x1 0 1 1 x4 x2 x3 σj 12 1 1 Z0= 2 3 0 -2 -1
-3 1 1 + - - 6M M 3M
迭代1
C→ -3 1 x2 -2 1 0 1 - M 1 x3 0 0 1 0 0 x4 1 0 0 0 0 x5 0 -1 0 M M x6 0 1 0 0 M x8 -1 -2 1 3M - 1 θj — 1/1 — CB XB B-1b x1 0 M 1 x4 x6 x3 σj 10 1 1 Z0= M +1 3 0 -2 -1
如何让人工变量尽快出基?
大M法
MAX:给每个人工变量一个量级很小的系数(-M) MIN:给每个人工变量一个量级很大的系数(M)
两阶段法
第一阶段:MIN Z =∑人工变量 第二阶段:原来的目标函数
示例
MinZ= 3x1 + x2 + x3 x1 2x2 + x3 ≤ 11 4x + x + 2x ≥ 3 1 2 3 s.t. 2x1 + x3 = 1 x j ≥ 0
0 4 8
(0,0)
x1
最小比值法
2 x2 + x3 = 8 x3 = 8 2 x2 s.t 0 x2 + x4 = 16 x4 = 16 0 x2 4 x + x = 12 x = 12 4 x 2 2 5 5
迭代1
C→ CB 0 0 3 XB x3 x4 x2 σj B-1b 2 16 3 Z0=9 2 x1 1 4 0 2↑ 3 x2 0 0 1 0 0 x3 1 0 0 0 0 x4 0 1 0 0 0 x5 -1/2 0 1/4 -3/4 θj 2/1 → 16/4 —
基可行解
基:Am×n的非奇异子阵Bm×m 基变量与非基变量: Bm×m各列对应的变量为 基变量XB,其余变量为非基变量XN 基解:X=(XB,XN),XB=B-1b, XN =0 基可行解:基解中,XB=B-1b≥0 目标函数值Z0=CBB-1b 可行基:基可行解对应的基
最优性判断—非基变量检验数
迭代1—阶段1
C→ 0 0 x2 -2 1 0 -1 0 x3 0 0 1 0 0 x4 1 0 0 0 0 x5 0 -1 0 1 1 x6 0 1 0 0 1 x8 -1 -2 1 3 θj — 1/1 — CB XB B-1b x1 0 1 0 x4 x6 x3 σj 10 1 1 Z0= 1 3 0 -2 0
初始表
C→ CB 0 0 0 XB x3 x4 x5 σj B-1b 8 16 12 Z0=0 2 x1 1 4 0 2 3 x2 2 0 4 3↑ 0 x3 1 0 0 0 0 x4 0 1 0 0 0 x5 0 0 1 0 θj 8/2 — 12/4 →
4.5 4 3.5 3 2.5 x2 2 1.5 1 0.5 0
M M -1 +1
迭代3
C→ -3 1 x2 0 1 0 0 1 x3 0 0 1 0 0 x4 0 x5 M x6 M x8 θj CB XB B-1b x1 -3 1 1 x1 x2 x3 4 1 9 1 0 0 0