不求甚解學經濟-一些簡單的綫性規劃

(Brandon) Song Li Lv4

一些簡單的綫性規劃

綫性規劃問題的表述

Standard 綫性規劃問題

給定係數矩陣 Am×nA_{m\times n},相應的向量 bm×1b_{m\times 1}cn×1c_{n\times 1},尋找 xn×1x_{n\times 1}ym×1y_{m\times 1} 使得

最大化問題 最小化問題
maxcTx\max c^{T}x minbTy\min b^{T}y
s.t.AxbAx\leq b s.t.y0m×1y\geq 0_{m\times 1}
x0n×1x\geq 0_{n\times 1} ATycA^{T}y\geq c

為一組對偶的 Standard 綫性規劃問題。

Canonical 綫性規劃問題

給定係數矩陣 Am×nA_{m\times n},相應的向量 bm×1b_{m\times 1}cn×1c_{n\times 1},尋找 xn×1x_{n\times 1}ym×1y_{m\times 1} 使得

最大化問題 最小化問題
maxcTx\max c^{T}x minbTy\min b^{T}y
s.t.Ax=bAx = b s.t.yfreey\text{free}
x0n×1x\geq 0_{n\times 1} ATycA^{T}y\geq c

為一組對偶的 Canonical 綫性規劃問題。

General 綫性規劃問題

給定係數矩陣 A11,A12,A21,A22A_{11},A_{12},A_{21},A_{22},相應的向量 b1,b2b_{1},b_{2}c1,c2c_{1},c_{2},尋找 x1,x2x_{1},x_{2}y1,y2y_{1},y_{2} 使得

最大化問題 最小化問題
maxc1Tx2+c2Tx2\max c_{1}^{T}x_{2}+c_{2}^{T}x_{2} minb1Ty1+b2Ty2\min b_{1}^{T}y_{1}+b_{2}^{T}y_{2}
s.t.A11x1+A12x2b1A_{11}x_{1}+A_{12}x_{2} \leq b_{1} s.t.y10y_{1}\geq 0
s.t.A21x1+A22x2=b2A_{21}x_{1}+A_{22}x_{2} = b_{2} s.t.y2freey_{2}\text{free}
x10x_{1}\geq 0 A11Ty1+A21Ty2c1A_{11}^{T}y_{1}+A_{21}^{T}y_{2}\geq c_{1}
x2freex_{2}\text{free} A12Ty1+A22Ty2=c2A_{12}^{T}y_{1}+A_{22}^{T}y_{2}= c_{2}

為一組對偶的 Canonical 綫性規劃問題。

Standard 與 Canonical 問題的等價性

(暫時掠過)

綫性等式系統解的存在性

定理1

任意給定 rr 個實數 α1,...,αr\alpha_{1},...,\alpha_{r} 滿足 (a1,...,αr)1×r01×r( a_{1},...,\alpha_{r} )_{1\times r} \neq 0_{1\times r},給定 Rn\mathbb{R}^{n} 中的一組綫性無關的行向量 a1×n(1),...a1×n(r)a_{1\times n}(1),...a_{1\times n}(r) (這意味 rnr\leq n,因爲 nn 維空間中不會有 m>nm>n 個綫性無關的向量族),總可以找到一個向量 yn×1y_{n\times 1},使得 aT(i)n×1y1×n=αi,i=1,..,ra^{T}(i)_{n \times 1}y_{1\times n}=\alpha_{i},\forall i=1,..,r

證明:將 rr 個行向量 a1×n(1),...a1×n(r)a_{1\times n}(1),...a_{1\times n}(r) 堆叠成一個矩陣 Ar×nA_{r\times n},矩陣的行 rank 與列 rank 相同,由於這 rr 個行向量被假設為綫性無關,因此矩陣 A 的行 rank 與列 rank 均爲 rr,從而可以找到矩陣 Ar×nA_{r\times n}rr 個綫性無關的列向量 ar×1[1],...,ar×1[r]a_{r\times 1}[1],...,a_{r\times 1}[r] 構成 Rr\mathbb{R}^{r} 的基。

既然列向量 ar×1[1],...,ar×1[r]a_{r\times 1}[1],...,a_{r\times 1}[r] 構成 Rr\mathbb{R}^{r} 的基,那麽任何一個向量 (a1,...,αr)1×r( a_{1},...,\alpha_{r} )_{1\times r} 都可以表示成 ar×1[1],...,ar×1[r]a_{r\times 1}[1],...,a_{r\times 1}[r] 的綫性組合,即存在 η1,...,ηr\eta_{1},...,\eta_{r} 使得

(a1,...,αr)1×r=j=1rηjαr×1[j]( a_{1},...,\alpha_{r} )_{1\times r}=\sum_{j=1}^{r}\eta_{j}\alpha_{r\times 1}[j]

則構造向量 y1×n=(η1,...,ηr,0,...,0)y_{1\times n}=(\eta_{1},...,\eta_{r},0,...,0) (共有 nrn-r 個 0 ),這個向量就滿足

aT(i)n×1y1×n=αi+0=αi,i=1,..,ra^{T}(i)_{n \times 1}y_{1\times n}=\alpha_{i}+0=\alpha_{i},\forall i=1,..,r

定理2. 存在性:齊次綫性系統的解

定理3. 存在性:綫性等式系統的解

如下兩個問題,恰好有一個有解,不會都無解,也不會同時有解:

  1. Am×nxn×1=bm×1A_{m\times n}x_{n\times 1}=b_{m\times 1}

  2. An×mTym×1=0n×1,b1×mTym×1=α1×1,α01×1A^{T}_{n\times m}y_{m\times 1}=0_{n\times 1},b^{T}_{1\times m}y_{m\times 1}=\alpha_{1\times 1},\alpha\neq 0_{1\times 1}

證明:

(1) 假設這兩個問題同時有解,則對於問題 2 的某個解 ym×1y_{m\times 1},將其乘在問題 1 上可得

y1×mTAm×nxn×1=y1×mTbm×1y^{T}_{1\times m}A_{m\times n} x_{n\times 1}=y^{T}_{1\times m}b_{m\times 1}

根據 yy 是問題 2 的解,上式進一步變成

y1×mTAm×nxn×1=y1×mTbm×1=αy^{T}_{1\times m}A_{m\times n} x_{n\times 1}=y^{T}_{1\times m}b_{m\times 1}=\alpha

但由於 yy 是問題 2 的解,左側

y1×mTAm×nxn×1=01×nTxn×1=y1×mTbm×1=αy^{T}_{1\times m}A_{m\times n} x_{n\times 1}=0^{T}_{1\times n}x_{n\times 1}=y^{T}_{1\times m}b_{m\times 1}=\alpha

但這就意味著

0n×1Txn×1=01×1=α1×101×10^{T}_{n\times 1}x_{n\times 1}=0_{1\times 1}=\alpha_{1\times 1}\neq 0_{1\times 1}

這就矛盾了。因此,只剩下兩種可能:要麽兩個問題同時無解,要麽恰好有一個有解。

現在假設問題 1 無解。取 Am×nA_{m\times n} 的一個基 a1,...,ara_{1},...,a_{r},其中每個 aia_{i}m×1m\times 1。既然 Ax=bAx=b 無解,那就意味著 bb 不在基 a1,...,ara_{1},...,a_{r} 所張成的空間裏,即 a1,...,ar,ba_{1},...,a_{r},b 綫性無關。因此,根據定理 1,存在一個向量 ym×1y_{m\times 1} 使得

aiy=0,i=1,...,rb1×mTym×1=α,對於某個α0\begin{align} a_{i}y=0&,&\forall i=1,...,r\\ b^{T}_{1\times m}y_{m\times 1}=\alpha&,&\text{對於某個}\alpha\neq 0 \end{align}

既然 a1,...,ara_{1},...,a_{r}AA 的一個基,則 AA 的任何一個列向量 aka_{k} 都可以表示爲 a1,...,ara_{1},...,a_{r} 的綫性組合,即

ak=i=1rλi,kai,k=1,...,ny1×mTak=i=1rλi,k(y1×mTai)=01×1\begin{align} a_{k}=&\sum_{i=1}^{r}\lambda_{i,k}a_{i},\forall k=1,...,n\\ y^{T}_{1\times m}a_{k}=&\sum_{i=1}^{r}\lambda_{i,k}(y^{T}_{1\times m}a_{i})=0_{1\times 1} \end{align}

因此 y1×mTAm×n=01×ny^{T}_{1\times m}A_{m\times n}=0_{1\times n} 以及 b1×mTym×1=αb^{T}_{1\times m}y_{m\times 1}=\alpha

但這就意味著問題 2 是有解的。這就説明這兩個命題不可能同時不成立。

綫性不等式系統解的存在性

Tucker’s Lemma

給定矩陣 Am×nA_{m\times n},對偶系統

系統 I 系統 II
An×mTym×10n×1A^{T}_{n\times m}y_{m\times 1}\geq 0_{n\times 1} Am×nxn×1=0m×1A_{m\times n}x_{n\times 1}=0_{m\times 1}
- xn×10n×1x_{n\times 1}\geq 0_{n\times 1}

分別存在解 ym×1y_{m\times 1}xn×1x_{n\times 1} 使得

a1×mTym×1+x1>0a_{1\times m}^{T}y_{m\times 1}+x_{1}>0

其中 a1×ma_{1\times m}^{\prime} 為矩陣 AA 的第一行的 transpose,x1x_{1} 為解 xx 的第一個分量。

Farkas’ Lemma - version 1

給定矩陣 Am×nA_{m\times n},給定向量 bm×1b_{m\times 1};對於系統 An×mTym×10n×1A^{T}_{n\times m}y_{m\times 1}\geq 0_{n\times 1},如果它的每個解 ym×1y_{m\times 1} 都滿足 b1×mTym×10b^{T}_{1\times m}y_{m\times 1}\geq 0。那麽存在解 xx^{\star} 使得 Ax=bAx^{\star}=bx0x^{\star}\geq 0 同時成立。

Farkas’ Lemma - version 2

給定矩陣 Am×nA_{m\times n},給定向量 bm×1b_{m\times 1},如下兩個命題恰好只有一個成立,不會都不成立,也不會都成立

系統 I 系統 II
Ax=bAx=b 有非負解 x0x\geq 0 yTA0y^{T}A\geq 0yTb<0y^{T}b<0 有解 yy 可正可負

Farkas’ Lemma - version 3

Ax=bAx=b 有非負解 x0x\geq 0 當且僅當 yTA0y^{T}A\geq 0yTb0y^{T}b\geq 0 有解 yy 可正可負。

推論 1

如果 yTb0y^{T}b\geq 0 對於每個滿足 yTA0y^{T}A\geq 0 的非負解 y0y\geq 0 都成立,那麽存在非負向量 x0x\geq 0 使得 AxbAx\geq b

推論 2

定理 4:Tucker 第一存在定理

定理 5:Tucker 第二存在定理

推論 3

定理 6

定理 7

定理 8

定理 9

綫性規劃

綫性規劃解的性質

退化解

超平面

綫性不等式的幾何性質

凸性

分離超平面

  • 標題: 不求甚解學經濟-一些簡單的綫性規劃
  • 作者: (Brandon) Song Li
  • 撰寫于 : 2021-01-01 00:06:30
  • 更新于 : 2024-12-06 04:30:39
  • 連結: https://brandonsli.com/p/2b881b54.html
  • 版權宣告: 保留所有權利 © (Brandon) Song Li