当前位置:网站首页>[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:aTx⩽b},其中 x ∈ R dim a x\in\Reals^{\dim\bf a} x∈Rdima .
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:Ax⩽b},其中 A \mathscr A A 是 m × n m\times n m×n 矩阵, b ∈ R m \mathbf b\in\Reals^m b∈Rm 而 x ∈ R n x\in\Reals^n x∈Rn .
多胞体( 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 P⊆C 且 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 y∈P 使得 ∥ x − y ∥ \Vert x-y\Vert ∥x−y∥ 最小,而 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 x⩾0 的解,当且仅当 A T y ⩾ 0 ∧ b T y < 0 \mathscr A^{\sf T}y\geqslant\mathbf 0\land \mathbf b^{\sf T}y<0 ATy⩾0∧bTy<0 无 y ∈ R n y\in\R^n y∈Rn 的解.
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)T⋅y=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) xTy⩾0(x∈C),The latter contains A T y ⩾ 0 \mathscr A^{\sf T}y\geqslant 0 ATy⩾0 . ■ \blacksquare ■
Corollary. 设 P = { x : A x ⩽ b } ≠ ∅ P=\{x:\mathscr Ax\leqslant\mathbf b\}\ne\varnothing P={ x:Ax⩽b}=∅,则 ∀ x ∈ P , c T x ⩽ δ \forall x\in P,\;\mathbf c^{\sf T}x\leqslant\delta ∀x∈P,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 ∃y⩾0,ATy=c∧bTy⩽δ .
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 Ax⩽b*yTAx⩽yTb*cTx⩽yTb*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 b∈Rm,c∈Rn,则
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:Ax⩽b}=min{ yTb:y⩾0∧ATy=c}
Proof. Let the left form be δ \delta δ,根据 Corollary \text{Corollary} Corollary 立刻有 ∃ y ⩾ 0 \exists y\geqslant 0 ∃y⩾0 使得 A T y = c ∧ b T y ⩽ δ \mathscr A^{\sf T}y=\mathbf c\land\mathbf b^{\sf T}y\leqslant\delta ATy=c∧bTy⩽δ .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 x⩾0,If unified to A x ⩽ b \mathscr Ax\leqslant\mathbf b Ax⩽b 中,不难发现有
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:x⩾0∧Ax⩽b}=min{ bTy:y⩾0∧ATy⩾c}
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.
边栏推荐
- I got the P8 "top" distributed architecture manual that went viral on Ali's intranet
- 【中央任务调度系统—通信开发】
- 宝塔一键部署WordPress(含宝塔添加站点、阿里云安全组配置、阿里云子域名解析)
- LeetCode·每日一题·1417.重新格式化字符串·模拟
- 解决 Pocess finished with exit code 1 Class not found 和 Command line is too long. Shorten the command
- LeetCode · Question of the Day · 1417. Reformatting String · Simulation
- fetch请求设置请求头错误导致无法跨域
- 7 天找个 Go 工作,Gopher 要学的条件语句,循环语句 ,第3篇
- Jetpack Compose学习(9)——Compose中的列表控件(LazyRow和LazyColumn)
- 人是怎么废掉的?人是怎么变强的?
猜你喜欢
chrome无痕浏览模式中使用插件
【Mysql系列】03_系统设计
LeetCode·每日一题·1417.重新格式化字符串·模拟
LeetCode · Question of the Day · 1417. Reformatting String · Simulation
Convolutional Neural Network Gradient Vanishing, The Concept of Gradient in Neural Networks
数字钱包红海角逐,小程序生态快速引入可助力占领智慧设备入口
Database indexes and their underlying data structures
数据库内核面试中我不会的问题(4)
How to explain to my girlfriend what is cache penetration, cache breakdown, cache avalanche?
计算数组某个元素的和
随机推荐
算法---跳跃游戏(Kotlin)
SDS Observatory
php获取微信小程序码并存储到oss
【阿克曼运动控制】
数据库 SQL 优化大总结之:百万级数据库优化方案
How to build programming ideas and improve programming ideas
【Prometheus】Alertmanager告警全方位讲解
【Prometheus】 Grafana数据与可视化
打印时间的各种格式
VMWare中安装的win10,新增其他盘符,如:E盘,F盘等
国产数据库有没有在国外的应用案例?
不可思议,全靠这份Android面试题,斩获多家互联网大厂offer
运动健康服务场景事件订阅,让应用推送“更懂用户”
chrome is set to dark mode (including the entire webpage)
Cholesterol-PEG-FITC,Fluorescein-PEG-CLS,胆固醇-聚乙二醇-荧光素水溶性
困扰所有SAP顾问多年的问题终于解决了
logstash/filebeat only receives data from the most recent period
天花板级微服务大佬总结出这份451页笔记告诉你微服务就该这么学
【分享】PPT还能做成这样?你一定没见过
[Ext JS]11.14 SimXhr.js?_dc=1659315492151:65 Uncaught TypeError problem analysis and solution