一些簡單的綫性規劃
綫性規劃問題的表述
Standard 綫性規劃問題
給定係數矩陣 Am×n,相應的向量 bm×1 和 cn×1,尋找 xn×1 和 ym×1 使得
最大化問題 |
最小化問題 |
maxcTx |
minbTy |
s.t.Ax≤b |
s.t.y≥0m×1 |
x≥0n×1 |
ATy≥c |
為一組對偶的 Standard 綫性規劃問題。
Canonical 綫性規劃問題
給定係數矩陣 Am×n,相應的向量 bm×1 和 cn×1,尋找 xn×1 和 ym×1 使得
最大化問題 |
最小化問題 |
maxcTx |
minbTy |
s.t.Ax=b |
s.t.yfree |
x≥0n×1 |
ATy≥c |
為一組對偶的 Canonical 綫性規劃問題。
General 綫性規劃問題
給定係數矩陣 A11,A12,A21,A22,相應的向量 b1,b2 和 c1,c2,尋找 x1,x2 和 y1,y2 使得
最大化問題 |
最小化問題 |
maxc1Tx2+c2Tx2 |
minb1Ty1+b2Ty2 |
s.t.A11x1+A12x2≤b1 |
s.t.y1≥0 |
s.t.A21x1+A22x2=b2 |
s.t.y2free |
x1≥0 |
A11Ty1+A21Ty2≥c1 |
x2free |
A12Ty1+A22Ty2=c2 |
為一組對偶的 Canonical 綫性規劃問題。
Standard 與 Canonical 問題的等價性
(暫時掠過)
綫性等式系統解的存在性
定理1
任意給定 r 個實數 α1,...,αr 滿足 (a1,...,αr)1×r=01×r,給定 Rn 中的一組綫性無關的行向量 a1×n(1),...a1×n(r) (這意味 r≤n,因爲 n 維空間中不會有 m>n 個綫性無關的向量族),總可以找到一個向量 yn×1,使得 aT(i)n×1y1×n=αi,∀i=1,..,r。
證明:將 r 個行向量 a1×n(1),...a1×n(r) 堆叠成一個矩陣 Ar×n,矩陣的行 rank 與列 rank 相同,由於這 r 個行向量被假設為綫性無關,因此矩陣 A 的行 rank 與列 rank 均爲 r,從而可以找到矩陣 Ar×n 的 r 個綫性無關的列向量 ar×1[1],...,ar×1[r] 構成 Rr 的基。
既然列向量 ar×1[1],...,ar×1[r] 構成 Rr 的基,那麽任何一個向量 (a1,...,αr)1×r 都可以表示成 ar×1[1],...,ar×1[r] 的綫性組合,即存在 η1,...,ηr 使得
(a1,...,αr)1×r=j=1∑rηjαr×1[j]
則構造向量 y1×n=(η1,...,ηr,0,...,0) (共有 n−r 個 0 ),這個向量就滿足
aT(i)n×1y1×n=αi+0=αi,∀i=1,..,r
定理2. 存在性:齊次綫性系統的解
定理3. 存在性:綫性等式系統的解
如下兩個問題,恰好有一個有解,不會都無解,也不會同時有解:
-
Am×nxn×1=bm×1
-
An×mTym×1=0n×1,b1×mTym×1=α1×1,α=01×1
證明:
(1) 假設這兩個問題同時有解,則對於問題 2 的某個解 ym×1,將其乘在問題 1 上可得
y1×mTAm×nxn×1=y1×mTbm×1
根據 y 是問題 2 的解,上式進一步變成
y1×mTAm×nxn×1=y1×mTbm×1=α
但由於 y 是問題 2 的解,左側
y1×mTAm×nxn×1=01×nTxn×1=y1×mTbm×1=α
但這就意味著
0n×1Txn×1=01×1=α1×1=01×1
這就矛盾了。因此,只剩下兩種可能:要麽兩個問題同時無解,要麽恰好有一個有解。
現在假設問題 1 無解。取 Am×n 的一個基 a1,...,ar,其中每個 ai 為 m×1。既然 Ax=b 無解,那就意味著 b 不在基 a1,...,ar 所張成的空間裏,即 a1,...,ar,b 綫性無關。因此,根據定理 1,存在一個向量 ym×1 使得
aiy=0b1×mTym×1=α,,∀i=1,...,r對於某個α=0
既然 a1,...,ar 是 A 的一個基,則 A 的任何一個列向量 ak 都可以表示爲 a1,...,ar 的綫性組合,即
ak=y1×mTak=i=1∑rλi,kai,∀k=1,...,ni=1∑rλi,k(y1×mTai)=01×1
因此 y1×mTAm×n=01×n 以及 b1×mTym×1=α
但這就意味著問題 2 是有解的。這就説明這兩個命題不可能同時不成立。
綫性不等式系統解的存在性
Tucker’s Lemma
給定矩陣 Am×n,對偶系統
系統 I |
系統 II |
An×mTym×1≥0n×1 |
Am×nxn×1=0m×1 |
- |
xn×1≥0n×1 |
分別存在解 ym×1 和 xn×1 使得
a1×mTym×1+x1>0
其中 a1×m′ 為矩陣 A 的第一行的 transpose,x1 為解 x 的第一個分量。
Farkas’ Lemma - version 1
給定矩陣 Am×n,給定向量 bm×1;對於系統 An×mTym×1≥0n×1,如果它的每個解 ym×1 都滿足 b1×mTym×1≥0。那麽存在解 x⋆ 使得 Ax⋆=b 和 x⋆≥0 同時成立。
Farkas’ Lemma - version 2
給定矩陣 Am×n,給定向量 bm×1,如下兩個命題恰好只有一個成立,不會都不成立,也不會都成立
系統 I |
系統 II |
Ax=b 有非負解 x≥0 |
yTA≥0 且 yTb<0 有解 y 可正可負 |
Farkas’ Lemma - version 3
Ax=b 有非負解 x≥0 當且僅當 yTA≥0 且 yTb≥0 有解 y 可正可負。
推論 1
如果 yTb≥0 對於每個滿足 yTA≥0 的非負解 y≥0 都成立,那麽存在非負向量 x≥0 使得 Ax≥b。
推論 2
定理 4:Tucker 第一存在定理
定理 5:Tucker 第二存在定理
推論 3
定理 6
定理 7
定理 8
定理 9
綫性規劃
綫性規劃解的性質
退化解
超平面
綫性不等式的幾何性質
凸性
分離超平面