线性规划的单纯形法常见问题及详细的解决方法

1047521767

                                          线性规划简介
3 |1 o5 @1 y5 z) o线性规划(Linear Programming, LP)，是运筹学中研究较早、发展较快、应用广泛、方法较成熟的一个重要分支，它是辅助人们进行科学管理的一种数学方法。研究线性约束条件下线性目标函数的极值问题的数学理论和方法。英文缩写LP。它是运筹学的一个重要分支，广泛应用于军事作战、经济分析、经营管理和工程技术等方面。为合理地利用有限的人力、物力、财力等资源作出的最优决策，提供科学的依据。(来自百度百科)

2 M0 C6 X. _1 r4 ^0 K/ b4 z
线性规划问题:
- X/ J1 ?& |3 e. x2 I* h- ~1 `在一般的线性规划问题中，我们希望优化一个满足一组线性不等式约束的线性函数。一直一组实数 ，对应的一组变量。定义关于这些变量的一个线性函数:
$ d. N1 u8 u* |9 g

那么等式和不等式 以及  都称为”线性约束”。并且不久之后我们即将看到，我们可以把两种不等式的约束形式统一为 这种形式，然而之后在运行解决线性规划的单纯形法时又会将它变成等价的 形式。

$ Y: |3 s) y1 `+ D, O  c
1 ]; }3 s9 a: s% a- f, I描述线性规划问题:
  \; u4 x! ^; V4 n7 p4 [9 E线性规划问题的标准型:
' Z( i/ }- A( F- m如果我们已知n个实数与m个实数，以及n*m个实数，我们希望找到n个实数

. d1 s" h# w) B( S# A
最大化


满足约束


- \1 t2 ?8 z  l1 X! M1 u' k
# k4 e. @- @' Z9 ]- C. a; s
7 J+ `  @  e7 z' f
或者将分别将和和表示成向量b, c, x，构造m行n列矩阵。那么标准型也可以写作:
# @: \$ r0 Q0 {9 b6 {1 G% D  d" P
4 E9 t% g5 L) {6 d4 {* i& T
最大化
1 U0 b  E$ S9 x; J* c& {
6 p* K$ B4 S+ V* r* P6 c4 b3 G' q
0 z8 y! h% V* o" _1 B9 A& a) h满足约束

            

0 ]- _, g) ^0 C
) \- M0 {% x' F. |( _注意到，标准型里面目标函数的目标是最大化，所有变量都具有非负约束，除了非负约束外的约束都是具有 <= 符号的非严格不等式。对于不满足标准型形式的线性规划，我们可以想办法将它转化为标准型:

2 ]8 z) [5 @, O6 G
问题1:目标函数的目标并非最大化，而是最小化。
* ^* c4 g1 e* C9 Q* H. j* u4 O取负即可。也就是最小化 等价于最大化

7 N; K6 n7 y. D+ i% L- _
6 F; s/ L# ~4 R" f问题2:一些变量不具有非负约束
0 }5 E$ o% n2 u1 r对于不具有非负约束的变量，把每次出现的地方都替换为适当的

+ p8 V6 P" A, t4 p" p+ A
: C3 S+ u$ s8 D3 c/ W: I问题3:大于等于约束和等式约束
0 B* ~: s" e6 ?' U1 c大于等于约束两边同时取负即可。即 等价于
$ T7 O, F* e1 e" d2 Z

等式约束拆成两个不等式约束。对于 等价于 且


描述线性规划问题的松弛型
$ F" o0 F8 r/ X3 |- W% y为了利用单纯形算法高效地求解线性规划，我们更喜欢把其中的除非负约束外的约束转换成等式约束。
. v$ S+ B% I% N, ^' Q
2 a$ X5 B) Z+ U
对于一个约束 引入一个松弛变量s，约束即可写作:
7 l8 B5 ]. G- C/ D

5 H3 F  y" Y: h' y/ D! f; s一般地，为了方便常会用 作为第i个约束条件的松弛变量。因此第i约束条件就可以写作:
用一个例子来说明上述描述线性规划问题的方法:
; h8 x/ D) ~1 m& q6 {7 A最小化
2 s4 c; S. J5 K" {5 i满足约束条件
                     
: V, l' U9 c9 e2 j& `- V5 w3 \                   　　　　　　 　
最小化目标转变成最大化:
最大化
; w- b: i+ g# ^( f: J  F满足约束条件
, P- a! `4 D) X, U- F, H                  
                   　　　　　　 　
等式约束
9 T# j( N5 {1 I( s3 E拆开等式约束:
; V5 x9 c& m3 d% l7 z8 B最大化
8 U- ~: ^$ {+ r$ Q+ R2 D7 T满足约束条件
% C# Z9 y! Z/ v9 g7 O                  
                  
                     　　　　　　 　
0 A; x) y- q/ j4 C3 ?) u& `- U增加非负约束　
替换为，添加对的非负约束，得到
5 `- n! f- _5 w! n) T最大化
满足约束条件  

% D8 w  z+ r4 S, i3 `                     
/ P! W, l$ ~7 B% g- m1 F8 V- u! R
                    
4 F* B( v$ i8 c9 t2 ^+ x
                       　　　　　　 　
: s, p1 r* P' t  i) @7 J8 D- t标准型到松弛型
' Q) L8 z0 r& y: l除了非负约束外，有3个约束条件，于是增加松弛变量
- Q* t+ [+ b5 D& m6 `

/ c( n1 D( J3 g. ?1 e最大化  
& `1 N% I( h0 H) F$ J

满足约束条件  
6 m% s0 N$ ]7 r2 `) g( m) h& S3 F
4 z+ a: ~9 A5 R0 t8 A7 Q# Q: @: v2 [6 E                     
8 n& j# Z( V1 F9 o6 y) v7 j
                     

; R5 P1 |" F6 ~6 ~' t/ a9 W5 q                         　　　　　　 　
思考:如果原先有一个等式约束，直接把它写到松弛型里面，与将它转变成标准型，再变成松弛型(这样就比前者多了一个约束条件)，问题是否等价？（提示:考虑后者引入的松弛变量s1 = s2 = 0)


+ W8 P5 K$ Z- |+ m* a) U解决线性规划问题的单纯形法:
单纯形算法是求解线性规划的经典算法，它再最坏情况瞎执行时间并不是多项式级别的，但是在实际的应用中，它总是相当快速。
2 f' ], [$ r$ w7 r
& D4 R. ?/ {" r8 O, }* ^+ k) ]+ S
个人理解:


单纯形法接受一个松弛型的输入，每轮迭代尝试将松弛型L转变为一个等价的松弛型L’，直到转变成一种”易于求解的松弛型”，得到结果。这就类似化三角法求行列式值。

2 p" G2 n. y. G: u% F9 p( Y
什么是”易于求解的松弛型”呢？如果目标函数具有   
6 s; ]4 `8 G2 t) D1 u
9 g' p# o) ~# p# ^/ q1 |% D, O7 K
这样的形式，由于变量具有非负约束，因此z的最大值即为常数C。如果我们能把一般的松弛型转变为这样”易于求解的松弛型”的松弛型即可得解。
5 o8 B& |. W% L( S$ \
$ p$ X& m2 i/ f$ f  K0 \
- t  S: ]8 D. Y* u+ e具体的方法 换个便于演示的例子:
  [; h. Q  l: Z& X' {7 {. t5 e考虑下面这个松弛型:
1 l! w& `  G0 g. p* u" E  l2 j! }

最大化  
, P# F  J+ O, g6 F4 a- t$ c
. |5 _. G' \- T9 c* T  R9 o
满足约束条件
( B- t$ i) c1 k5 K7 T+ m* t
                    

: L& m. Z2 V. i1 \9 k* Y. d                    

                       　　　　　　 　

+ U* g$ m( b6 @. M
术语：我们把松弛型约束条件等号左边的变量称为基本变量，右边的称为非基本变量。


我们会选取目标函数里面某个系数为正的非基本变量，并尽可能增大的取值且不违反任何约束条件。由于的增大可能会导致其它变量值缩小，但是其它变量也存在非负约束，这样就约束了的取值的上界，我们选出对约束最紧的变量，令 与 互换(代入法)得到一个新的等价的松弛型。具体对于我们这个例子就是:


首先看变量，对它约束最紧的是第三个约束，约束了。是约束了。于是互换 与，把放到 的位置，得到

3 H6 a$ l+ B5 B1 w
* h# Z/ v9 Q5 ]* f! ]& w将约束条件不等式右边的所有用这个式子替换掉，即可得到  
( D9 Z7 u# c5 A

最大化


# \! O- u' _& R满足约束条件  

3 V6 n) n) e# O2 K$ _! q                     

                     
0 Q4 Y" ^( a6 @5 Q5 q. y
                         　　　　　　 　


1 }- J, x4 F" V9 U. s( e这样的一次操作叫做转动(pivot)


% [7 a( l5 z" u9 ~0 Q0 O(结合例子请思考:选择约束最紧的变量做互换的原因。提示:如果换的不是最紧的约束看看会违反什么性质)

$ o0 }. |! O/ b- B; y/ v+ l
然后选择变量或者，不选是因为增大会导致目标z减小。

# A. v( y; G' n- `6 W4 {2 a: r
如果选择 ，同样的方法，互换与可以得到新的但等价的约束系统:  
9 p8 a  H& r+ K- O
: r7 j' o) q& \9 V
最大化   
& ?: `9 T1 R5 _, ?% V1 L/ z
3 j7 E' x  l+ b8 @4 n" E- j满足约束条件  

  m1 f- W) f0 R  P+ O5 t& E                        

! c6 P5 ?: N% S! j3 f                        
) D1 a3 d/ O/ X* C
                        
7 b6 k+ s# y% D3 l" R
3 ^( L+ p( I8 j, |2 F0 V

* G6 t" M$ s1 a0 S) A. t- `. T接下来也只剩下这个变量可选了，与互换，得到  

5 _5 |: g0 K5 m8 i
最大化   


满足约束条件  
$ S- d' r" l/ O8 ~6 Z" ]
                     

: g+ L2 ~7 P  C3 A& A. I  c                     
) h. H* T8 G2 z, t& q  T2 N) O
                        　　　　　　 　
! g. P, T- V+ \# o% ?
7 U; q  X1 H, K4 c( h, A, D. L此时目标函数右边变量系数全部为负数，且变量具有非负约束，显然这时候能得到的z的最大值为28。当且仅当。带入解得。除去附加的松弛变量，最终这个线性规划最优解为
0 L& q9 T& B3 p& q
2 h! c' ~: g3 c* ^' n# b
0 {- j* Q5 Z4 r# G
; g% v7 v( S' g  C7 Z" u; H思考:每次转动操作交换的两个变量对单纯形法的运行时间有着怎样的影响？如何选择可以使单纯形法尽快结束？（提示:每次pivot尽可能让目标函数增大得多）

& S# l$ G8 n6 T9 h( e
单纯形算法代码实现

4 {2 U3 ?. x# L* ^, I' U
伪代码:


4 g' y5 p9 J5 V8 x0 z; {9 w& ?simplex:
8 C; O) H5 {, w* v3 I    检查是否无解
' X7 h8 A: ~* q    loop:
        找出目标函数中前面系数为正数的一个变量x
        如果找不到x
( ?7 n* a! b  C            返回目标函数的常数项(已找到目标函数最大值)
        找到对x的增长约束最紧的变量y   
        如果找不到y
" a1 Q' Z5 k* o            返回Inf
9 d- N8 g4 S4 k# _  F  q        互换变量x,y(pivot(x,y)，操作后问题与原来等价)


. q6 A8 L, h- {: P, [
' [( \/ Y* U3 L
非标准型转换为标准型魔法码(github上给出的参考代码要求输入的是标准型)
8 @* g3 T. x8 _& T% F
2 `3 ~  b# M. D# b  R
2 z( ]) k0 [7 ^* z0 p* G输出目标函数达到最大值时自变量的取值


, V: L5 {; n. y. Y3 d, r) K思考:会不会存在一种情况:虽然一个线性规划问题有解，但是单纯形算法仍会无休止迭代下去。

! J. N. T9 T8 y& H4 X5 D2 w
思考:如果判断出这样一种无解的情况:  
2 ^; t+ H+ {* d
4 ]. {6 J4 d# A1 J

8 N. `; E$ I4 _3 q

" p2 y; i* f% o3 Z　
1 h7 O+ A( U  l+ }3 D
, k5 B* L; I! m5 s请关注数学中国微信公众号和数学中国网微博号 联系QQ数学中国浅夏3243710560

" w0 e- H$ y) s6 e" L0 i1 E