A Fundamental Question: 给定一个(不)等式,如何寻找其等价形式?

∙  \bullet\;∙ 对于线性或者一般的凸(不)等式 ,有 Farkas引理

∙  \bullet\;∙ 对于二次(不)等式 ,则有 S-引理 ------------------- S-Lemma本质上是Hilbert零点定理的特例?

凸锥的对偶的计算

Farkas引理

These results essentially state that a concave inequality is a (logical) consequence of some convex inequalities if and only if it is a nonnegative linear combination of those convex inequalities and an identically true inequality. This is important, since it is relatively easy to check if an inequality is a linear combination of some other inequalities.

定理. (Farkas’ lemma) 假设 A∈Rm×nA\in\mathbb{R}^{m\times n}A∈Rm×n,b∈Rmb\in\mathbb{R}^mb∈Rm。则 Ax=b  (x≥0)Ax=b \;(x\ge 0)Ax=b(x≥0) 有可行解当且仅当 ATy≤0⇒bTy≤0A^Ty\leq 0\Rightarrow b^Ty\leq 0ATy≤0⇒bTy≤0。
Proof. 定义一个凸锥 K={y ∣ ATy≤0}\mathcal{K}=\{y\,|\,A^Ty\leq 0\}K={y∣ATy≤0} 它的对偶锥是 K∗={Ax ∣ x≥0}\mathcal{K}^*=\{Ax\,|\,x\ge 0\}K∗={Ax∣x≥0} 而根据定义有 K∗={b ∣ bTy≤0,  ∀y∈K}\mathcal{K}^*=\{b\,|\,b^Ty\leq 0, \; \forall y\in\mathcal{K}\}K∗={b∣bTy≤0,∀y∈K} 注意到 Ax=b (x≥0)Ax=b\,(x\ge 0)Ax=b(x≥0) 有可行解等价于 b∈K∗b\in\mathcal{K}^*b∈K∗,引理即可得证。

定理. (Farkas’ lemma variant) 假设 A∈Rm×nA\in\mathbb{R}^{m\times n}A∈Rm×n,b∈Rmb\in\mathbb{R}^mb∈Rm。则 Ax≤bAx\leq bAx≤b 有可行解当且仅当 Ax=0  (x≥0)⇒bTx≥0Ax=0\;(x\ge 0)\Rightarrow b^Tx\ge 0Ax=0(x≥0)⇒bTx≥0。

定理. 假设 KKK 是一个闭凸锥,K∗K^*K∗ 是其对偶,则 x∈Kx\in Kx∈K 当且仅当对任意 y∈K∗y\in K^*y∈K∗ 有 yTx≥0y^T x\ge 0yTx≥0。

定理. (Farkas’ Theorem) 假设 f,g1,⋯&ThinSpace;,gm:Rn→Rf,g_1,\cdots,g_m:\mathbb{R}^n\rightarrow\mathbb{R}f,g1​,⋯,gm​:Rn→R 是一组凸函数,C⊂Rn\mathcal{C}\subset\mathbb{R}^nC⊂Rn 是一个凸集,如果 g1,⋯&ThinSpace;,gmg_1,\cdots,g_mg1​,⋯,gm​ 满足 Slater 条件,也即,存在 C\mathcal{C}C 的一个内点 xˉ\bar{x}xˉ 使得 gi(xˉ)&lt;0,i=1,⋯&ThinSpace;,mg_i(\bar{x})&lt;0, i=1,\cdots,mgi​(xˉ)<0,i=1,⋯,m。则下列两组论述等价:
(1) 不等式组
f(x)&lt;0f(x)&lt;0f(x)<0 gi(x)≤0,i=1,⋯&ThinSpace;,mg_i(x)\leq0, i=1,\cdots,mgi​(x)≤0,i=1,⋯,m x∈Cx\in\mathcal{C}x∈C 无解。
(2) 存在 y1,⋯&ThinSpace;,ym≥0y_1,\cdots, y_m\ge 0y1​,⋯,ym​≥0,使得对所有的 x∈Cx\in\mathcal{C}x∈C 有
f(x)+∑i=1mgi(x)≥0f(x)+\mathop{\sum}\limits_{i=1}^m g_i(x) \ge 0f(x)+i=1∑m​gi​(x)≥0

Proof. 定理的证明主要是基于分离技巧 (separation argument) ,一个关键的事实是:
集合 S={(u,v1,⋯&ThinSpace;,vm)∈Rm+1:∃&ThickSpace;x∈Rn,f(x)&lt;u,gi(x)≤vi,i=1,⋯&ThinSpace;,m}S=\{(u,v_1,\cdots,v_m)\in\mathbb{R}^{m+1}:\exist\;x\in\mathbb{R}^n,f(x)&lt;u,g_i(x)\leq v_i, i=1,\cdots,m\}S={(u,v1​,⋯,vm​)∈Rm+1:∃x∈Rn,f(x)<u,gi​(x)≤vi​,i=1,⋯,m} 是一个凸集,并且 (1) 中的不等式组无解当且仅当 0∉S0\notin S0∈/​S。

关于SDP的Farkas引理,可以参考 Ye 的 Linear Conic Programming 第28页。

S-引理

定理. (S-lemma, Yakubovich) 假设 f,g:Rn→Rf,g:\mathbb{R}^n\rightarrow\mathbb{R}f,g:Rn→R 是二次函数,并且存在 xˉ∈Rn\bar{x}\in\mathbb{R}^nxˉ∈Rn 使得 g(xˉ)&lt;0g(\bar{x})&lt;0g(xˉ)<0。则下列两个论述等价:
(1) 不等式组
f(x)&lt;0f(x)&lt;0f(x)<0 g(x)≤0g(x)\leq 0g(x)≤0 无解
(2) 存在 y≥0y\ge 0y≥0 使得对任意 x∈Rnx\in\mathbb{R}^nx∈Rn,有
f(x)+yg(x)≥0f(x)+yg(x)\ge 0f(x)+yg(x)≥0

Farkas'Lemma 和 S-Lemma相关推荐

  1. Lindström–Gessel–Viennot lemma

    (摘自知乎) Lindström–Gessel–Viennot lemma(Lindström-Gessel-Viennot lemma这里有详细介绍跟证明) 在一个有向无环图里,想要计算从n个起点 ...

  2. 定理(Theorem)、引理(Lemma)、推论(Corollary)的定义及LaTeX用法

    1.定理(Theorem) 1.1 定义: Theorem(定理)----a mathematical statement that is proved using rigorous mathemat ...

  3. 论文中的定理(Theorem)、引理(Lemma)、推论(Corollary)

    Theorem:就是定理,比较重要的,简称是 Thm. Lemma:小小的定理,通常是为了证明后面的定理,如果证明的篇幅很长時,可能會把证明拆成几个部分來论述,虽然篇幅可能变多,但派络却很清楚. Co ...

  4. IEEE latex 编写lemma

    \newtheorem{lemma}{Lemma} \begin{lemma}[Einstein-Podolsky-Rosenberg] \end{lemma}

  5. Differences between Theorem, Lemma and Corollary, etc.

    Definition - a precise and unambiguous description of the meaning of a mathematical term.  It charac ...

  6. 区分定理(Theorem)、引理(Lemma)、推论(Corollary)等概念

    ZZ: http://blog.sina.com.cn/s/blog_a0e53bf70101jwv1.html Theorem:就是定理,比較重要的,簡寫是 Thm. Lemma:小小的定理,通常是 ...

  7. What is the difference between a theorem, a lemma, and a corollary?

    Definition - a precise and unambiguous description of the meaning of a mathematical term.  It charac ...

  8. ADPRL - 近似动态规划和强化学习 - Note 8 - 近似策略迭代 (Approximate Policy Iteration)

    Note 8 近似策略迭代 Approximate Policy Iteration 近似策略迭代 Note 8 近似策略迭代 Approximate Policy Iteration 8.1 通用框 ...

  9. 什么样的NLP库,可以支持53种语言?

    https://www.toutiao.com/a6663746829693223427/ 在学习自然语言处理(NLP)时,常常会遇到一个挑战--我们可以为非英语语言建立模型吗?很长一段时间以来,这个 ...

最新文章

  1. Android Studio 的安装笔记
  2. mysql原理~undo
  3. python解析器_Python Python解析器
  4. shell实例第14讲:字符串截取的8种方法
  5. One Order CLOSING date修改后的执行原理
  6. 这就是搜索引擎:核心技术详解
  7. Codeforces Round #727 (Div. 2) 题解
  8. 抓取新浪的每日星座运势
  9. 【Selenium 爬爬 猪 八 戒】
  10. Python多线程编程方式1(转)
  11. html 拦截境外ip,网站怎样屏蔽所有国外IP访问? 一个小技巧可以搞定!
  12. kaldi_源码笔记 thchs-30_data_prep
  13. [leetcode] 69 Sqrt(x)
  14. 十二星座谁最不会顾及别人感受?
  15. 九九乘法表 - Java (矩形、正三角形、倒三角形)
  16. hdu 6357 Hills And Valleys——dp
  17. 路由的导航守卫过渡动效transtion导航守卫 路由懒加载 路由元信息 @stage3---wee2--day7
  18. rs232串口通讯模块 - Delphi编程
  19. ElasticSearch之HTTP索引操作和文档操作
  20. window.open 全屏展示

热门文章

  1. Codeforces Round #294 (Div. 2) -- A. A and B and Chess
  2. 华为机试题 2014
  3. php微信公众号支付实例教程,php微信支付之公众号支付功能
  4. 苹果计算机开机进入安全模式,iphone进入安全模式及退出安全模式的方法
  5. linux安装包PGP加密验证
  6. 关于热流体仿真的博客和联系方式
  7. 利用Flowchart.fun快速可视化递归过程
  8. Real-Time Rendering读书笔记——01
  9. 树莓派用c语言pwm控制电机,树莓派学习笔记之PWM控制直流电机转速
  10. 嵌入式Linux开发环境搭建-4-嵌入式编程基础知识