当前位置:网站首页>[Study Notes] Dual Theorem of Linear Programming

[Study Notes] Dual Theorem of Linear Programming

2022-08-11 10:54:00 OneInDark

注:Vectors are not necessarily marked in bold,尤其是 x , y x,y x,y,Because their bold is so ugly.

Definition. 向量比较 x < y x<y x<y 等价于 ∀ i ,    x i < y i \forall i,\;x_i<y_i i,xi<yi,That is, a numerical comparison is made on each component.

半空间( affine halfspaces \text{affine halfspaces} affine halfspaces):形如 { x : a T x ⩽ b } \{x:\mathbf a^{\sf T}x\leqslant b\} { x:aTxb},其中 x ∈ R dim ⁡ a x\in\Reals^{\dim\bf a} xRdima .

polyhedron( polyhedra/polyhedron \text{polyhedra/polyhedron} polyhedra/polyhedron):A limited half-space intersection.它等价于 P = { x : A x ⩽ b } P=\{x:\mathscr Ax\leqslant\bf b\} P={ x:Axb},其中 A \mathscr A A m × n m\times n m×n 矩阵, b ∈ R m \mathbf b\in\Reals^m bRm x ∈ R n x\in\Reals^n xRn .

多胞体( polytope \text{polytope} polytope):The convex hull generated by a finite number of vectors.Equivalent to finite size p o l y h e d r o n \rm polyhedron polyhedron .

Comment. I'm not sure what the Chinese translations of the above things should be.

Proposition. 若 x ∉ P x\notin P x/P 其中 P P P p o l y h e d r a \rm polyhedra polyhedra,There is a half space C C C 使得 P ⊆ C P\subseteq C PC x ∉ C x\notin C x/C .

The proof of the theorem is difficult to separate from the support of geometry:需要找到 y ∈ P y\in P yP 使得 ∥ x − y ∥ \Vert x-y\Vert xy 最小,而 y y y The existence of is guaranteed by its geometrically intuitive properties.But if you admit the use of geometry in it,Then the theorem is obvious.Therefore, no proof is given here.

Lemma.( Farkas’ lemma \text{Farkas' lemma} Farkas’ lemma)  A x = b \mathscr Ax=\mathbf b Ax=b x ⩾ 0 x\geqslant\mathbf 0 x0 的解,当且仅当 A T y ⩾ 0 ∧ b T y < 0 \mathscr A^{\sf T}y\geqslant\mathbf 0\land \mathbf b^{\sf T}y<0 ATy0bTy<0 y ∈ R n y\in\R^n yRn 的解.

Proof. 当 A x = b \mathscr Ax=\mathbf b Ax=b 有解时,Counter-evidence y y y 有解,则
0 > b T y = ( A x ) T ⋅ y = x T ( A T y ) ⩾ 0 0>\mathbf b^{\sf T}y=(\mathscr Ax)^{\sf T}\cdot y=x^{\sf T}(\mathscr A^{\sf T}y)\geqslant 0 0>bTy=(Ax)Ty=xT(ATy)0

A x = b \mathscr Ax=\mathbf b Ax=b 无解时,设 a 1 , a 2 , … , a n a_1,a_2,\dots,a_n a1,a2,,an A \mathscr A A 的列向量,令 C C C for their convex hull,则 x x x Unsolved is equivalent to b ∉ C \mathbf b\notin C b/C,So there is a half plane division b \mathbf b b C C C,That is, there is a vector y y y 使得 b T y < 0 \mathbf b^{\sf T}y<0 bTy<0 x T y ⩾ 0    ( x ∈ C ) x^{\sf T}y\geqslant 0\;(x\in C) xTy0(xC),The latter contains A T y ⩾ 0 \mathscr A^{\sf T}y\geqslant 0 ATy0 . ■ \blacksquare

Corollary. 设 P = { x : A x ⩽ b } ≠ ∅ P=\{x:\mathscr Ax\leqslant\mathbf b\}\ne\varnothing P={ x:Axb}=,则 ∀ x ∈ P ,    c T x ⩽ δ \forall x\in P,\;\mathbf c^{\sf T}x\leqslant\delta xP,cTxδ 等价于 ∃ y ⩾ 0 ,    A T y = c ∧ b T y ⩽ δ \exists y\geqslant 0,\;\mathscr A^{\sf T}y=c\land\mathbf b^{\sf T}y\leqslant\delta y0,ATy=cbTyδ .

Proof. 充分性:若 y y y 存在,there is immediately
A x ⩽ b    *    y T A x ⩽ y T b    *    c T x ⩽ y T b    *    c T x ⩽ δ \mathscr Ax\leqslant\mathbf b \implies y^{\sf T}\mathscr Ax\leqslant y^{\sf T}\mathbf b \implies\mathbf c^{\sf T}x\leqslant y^{\sf T}\mathbf b \implies\mathbf c^{\sf T}x\leqslant\delta Axb*yTAxyTb*cTxyTb*cTxδ

必要性:It's essentially solving equations
( y T λ ) ( A b 0 1 ) = ( c T δ ) \begin{pmatrix}y^{\sf T} & \lambda\end{pmatrix} \begin{pmatrix} A & b\\ 0 & 1\end{pmatrix} =\begin{pmatrix}\mathbf c^{\sf T} & \delta\end{pmatrix} (yTλ)(A0b1)=(cTδ)

y T , λ ⩾ 0 y^{\sf T},\lambda\geqslant 0 yT,λ0 的解.由 Farkas’ lemma \text{Farkas' lemma} Farkas’ lemma 可知存在 z , μ ∈ R z,\mu\in\Reals z,μR 使得
( A b 0 1 ) ( z μ ) ⩾ 0 and ( c T δ ) ( z μ ) < 0 \begin{pmatrix} A & b\\ 0 & 1\end{pmatrix} \begin{pmatrix}z \\ \mu\end{pmatrix}\geqslant 0 \quad\text{and}\quad \begin{pmatrix}\mathbf c^{\sf T}&\delta\end{pmatrix} \begin{pmatrix}z \\ \mu\end{pmatrix}<0 (A0b1)(zμ)0and(cTδ)(zμ)<0

Known by the former μ ⩾ 0 \mu\geqslant 0 μ0 .讨论 μ = 0 \mu=0 μ=0 μ > 0 \mu>0 μ>0 的情况,分别取 x = x 0 − γ z    ( γ → + ∞ ) x=x_0-\gamma z\;(\gamma\rightarrow+\infty) x=x0γz(γ+) x = − μ − 1 z x=-\mu^{-1}z x=μ1z Contradictions can be drawn,完成证明.Leave it as an exercise. ■ \blacksquare

Theorem.( Duality theorem of linear programming \text{Duality theorem of linear programming} Duality theorem of linear programming) 令 A \mathscr A A m × n m\times n m×n 矩阵, b ∈ R m ,    c ∈ R n \mathbf b\in\Reals^m,\;\mathbf c\in\Reals^n bRm,cRn,则
max ⁡ { c T x : A x ⩽ b } = min ⁡ { y T b : y ⩾ 0 ∧ A T y = c } \max\{\mathbf c^{\sf T}x:\mathscr Ax\leqslant\mathbf b\}=\min\{y^{\sf T}\mathbf b:y\geqslant\mathbf 0\land \mathscr A^{\sf T}y=\mathbf c\} max{ cTx:Axb}=min{ yTb:y0ATy=c}

Proof. Let the left form be δ \delta δ,根据 Corollary \text{Corollary} Corollary 立刻有 ∃ y ⩾ 0 \exists y\geqslant 0 y0 使得 A T y = c ∧ b T y ⩽ δ \mathscr A^{\sf T}y=\mathbf c\land\mathbf b^{\sf T}y\leqslant\delta ATy=cbTyδ .This shows that the right-hand form is at most δ \delta δ .

The left does not exceed the right,因为 c T x = ( y T A ) ⩽ y T b \mathbf c^{\sf T}x=(y^{\sf T}\mathscr A)\leqslant y^{\sf T}\mathbf b cTx=(yTA)yTb . ■ \blacksquare

In the usual linear programming duality x ⩾ 0 x\geqslant\mathbf 0 x0,If unified to A x ⩽ b \mathscr Ax\leqslant\mathbf b Axb 中,不难发现有
max ⁡ { c T x : x ⩾ 0 ∧ A x ⩽ b } = min ⁡ { b T y : y ⩾ 0 ∧ A T y ⩾ c } \max\{\mathbf c^{\sf T}x:x\geqslant\mathbf 0\land\mathscr Ax\leqslant\mathbf b\}=\min\{\mathbf b^{\sf T}y:y\geqslant\mathbf 0\land\mathscr A^{\sf T}y\geqslant\mathbf c\} max{ cTx:x0Axb}=min{ bTy:y0ATyc}

This is a special form of the duality theorem,It is also the usual duality. crashed \textsf{crashed} crashed 指出:The mutual relaxation theorem is obvious.Therefore omitted.

原网站

版权声明
本文为[OneInDark]所创,转载请带上原文链接,感谢
https://yzsam.com/2022/223/202208111038183927.html