强对偶性、弱对偶性以及KKT条件的证明(对偶问题的几何证明)[通俗易懂]

强对偶性、弱对偶性以及KKT条件的证明(对偶问题的几何证明)[通俗易懂]本文主要探讨优化问题中强、弱对偶性以及KKT条件的证明。

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE稳定放心使用

1.原问题

首先给出问题的一般形式:
在这里插入图片描述
  上式表明我们一共有M+N个约束条件,对于不是求最小值或者约束条件大于等于0的情况,我们添加一个负号就可以变成上面这种形式。
  上述问题我们一般称之为带约束的原问题。

  利用拉格朗日乘子法,我们构造一个新的函数以及约束条件如下:
在这里插入图片描述
其中:
在这里插入图片描述

  我们称上面的问题为无约束的原问题(对x不再有约束)。上述L是拉格朗日乘子法的基本形式,这个就不再证明。

2.对偶问题

  对于无约束的原问题,我们先直接给出它的对偶问题形式(其实就是简单交换min和max):
在这里插入图片描述
  上述问题我们称之为原问题的对偶问题。

2.1弱对偶性的一般证明

  所谓弱对偶性,指的是:
在这里插入图片描述
  在再谈SVM(hard-margin和soft-margin详细推导、KKT条件、核技巧)中,我们大致口头证明了弱对偶性的成立,即“凤尾”>=“鸡头”。何谓“凤尾”?我先选出最强的一批人( max ⁡ f \max f maxf),然后组成实验班,实验班倒数第一就是 min ⁡   max ⁡ f \min \ \max f min maxf;何谓“鸡头”?我先选出最弱的一批人( min ⁡ f \min f minf),然后在这批比较弱的人当中选出最强的那个人,也即是 max ⁡   min ⁡ f \max \ \min f max minf,那么“鸡头”与“凤尾”孰强孰弱,是显而易见的。
  现在我们利用数学推导来大致证明一下弱对偶性。
  对于 L ( x , λ , η ) L(x,\lambda,\eta) L(x,λ,η)这个函数,我们知道下面这个不等式一定成立:
在这里插入图片描述
上面这个不等式很好理解,中间 L ( x , λ , η ) L(x,\lambda,\eta) L(x,λ,η)我们可以理解为L的值域,值域里面的任何一个数,必然是大于等于它的最小值,小于等于它的最大值。其实这一步已经证明出弱对偶性了,不过为了更容易理解,我们可以进一步说明。
  上述不等式最左边的表达式最后是关于 λ , η \lambda,\eta λ,η的一个函数,而最右边是一个关于 x x x的函数,因此我们又令:
在这里插入图片描述
因此我们有:
在这里插入图片描述
证毕。

2.2弱对偶性的几何证明

  为了使问题简化,同时方便证明,我们去掉原问题中等式的约束条件,同时不等式约束条件只保留一个,即原问题变成:
在这里插入图片描述
那么拉格朗日函数就变成:
在这里插入图片描述
我们又令:
在这里插入图片描述
p ∗ p^* p是原问题的最优解, d ∗ d^* d是对偶问题的最优解。证明弱对偶性实际上就是证明 d ∗ ≤ p ∗ d^*\leq p^* dp
  我们令区域G的表达形式为:
在这里插入图片描述

D是原问题的定义域,G表示一个个点的集合,点的横坐标是约束条件 u = m 1 ( x ) u=m_{1}(x) u=m1(x),纵坐标是原函数 t = f ( x ) t=f(x) t=f(x)
  有了上述集合G的定义之后,我们就可以对 p ∗ , d ∗ p^*,d^* p,d进行变式。
  首先对 p ∗ p^* p进行变式:
在这里插入图片描述
因为 t = f ( x ) t=f(x) t=f(x),所以 p ∗ p^* p实际上就是t的最小值,反映到集合G中去就是指一个点的纵坐标,这个点要满足两个条件:一是肯定要在G中,二是 m 1 ( x ) ≤ 0 m_{1}(x)\leq0 m1(x)0也就是该点的横坐标小于等于0。
如下图所示:
在这里插入图片描述
  我们对 u ≤ 0 u\leq0 u0那部分,也就是图中阴影部分上的每个点,找到一个最低的点,它的纵坐标就是 p ∗ p^* p
  接着对 d ∗ d^* d进行变形:
在这里插入图片描述
上述 t + λ u t+\lambda u t+λu的来源为:
在这里插入图片描述
对变形后 d ∗ d^* d我们令:
在这里插入图片描述
我们先找到 g ( λ ) g(\lambda) g(λ)在图中的位置:我们知道 t + λ u t+\lambda u t+λu实际上也是一个值,我们不妨令 t + λ u = k t+\lambda u=k t+λu=k,该式表示一条斜率为 − λ -\lambda λ并过(0,k)的直线,我们要找的是 t + λ u t+\lambda u t+λu的最小值,实际上就是k的最小值,实际上就是该直线与纵轴交点的最小值。而在求 min ⁡ x   t + λ u \min \limits_{x}\ t+\lambda u xmin t+λu的最小值时, λ \lambda λ固定的,因此斜率 − λ -\lambda λ也是固定的:
在这里插入图片描述
  我们保持斜率不变移动直线,不断往上移动,则该直线与纵轴交点的纵坐标k也不断增大,因为限制条件还有一个就是 ( u , t ) ∈ G (u,t)\in G (u,t)G,因此该直线必须经过区域G,我们一直往上移动,直到直线第一次与G相交,记下相应的k值为 k 1 k_{1} k1,再继续往上也都满足条件,知道该直线与G不再相交,但是现在我们求得是最小值,那么最小值其实就是 k 1 k_{1} k1,即 g ( λ ) g(\lambda) g(λ)就等于 k 1 k_{1} k1
  进一步,我们要求:
在这里插入图片描述
这里要重点注意:上一步我们求得了 g ( λ ) g(\lambda) g(λ)就等于 k 1 k_{1} k1,但是这种情况只是一种情况,在上一步求 g ( λ ) g(\lambda) g(λ)时,我们假若改变斜率 − λ -\lambda λ,那么 k 1 k_{1} k1的值是会变的,如下所示:
在这里插入图片描述
  我们换一个 λ \lambda λ固定时, g ( λ ) g(\lambda) g(λ)也就是 k 1 k_{1} k1自然也就在变。
  第二步要干的其实就是让我们求这个 g ( λ ) g(\lambda) g(λ)的最大值。那什么时候是最大的?实际上就是以G的最低点为轴,我们旋转直线,直到与左上方最低点相交时, g ( λ ) g(\lambda) g(λ)是最大的。如下所示:
在这里插入图片描述

  可能很多人就有疑问了:我为什么不可以让斜率继续增大?让直线穿过G?这里有疑问的同学不妨回忆一下第一个步骤:
在这里插入图片描述
我们在确定 g ( λ ) g(\lambda) g(λ)的时候,第一次相切我们就停止了,而后改变斜率继续相切:如上图所示,假如你继续增大斜率,那么该直线就不跟G的最低点相切了。我们是先推第一步再推第二步,而推第二步的时候肯定必须满足第一步的条件,所以 d ∗ d^* d只能是在那个位置。
在这里插入图片描述
  从上图也可以看出来: d ∗ ≤ p ∗ d^*\leq p^* dp,也就是对偶问题的解是小于等于原问题的解的,也即是说满足弱对偶性。因此这里我们进一步证明了弱对偶性。

2.3强对偶性的几何表示以及条件

  什么是强对偶性?就是指原问题的解与对偶问题的解是相同的,也即是: d ∗ = p ∗ d^*=p^* d=p
  画个图:
在这里插入图片描述
  假设G是一个凸集,那么根据上面找 d ∗ d^* d p ∗ p^* p的思路,我们很容易知道这个时候二者是相等的,也就是满足强对偶关系。
  那上面这句话的意思就是说:只要是凸集就一定满足强对偶关系。这句话不是正确的,不是所有的凸集都满足强对偶关系,但是加上slater条件就一定满足。

2.4 slater condition

  先直接给出slater条件的定义:对于x的定义域D,如果它存在一个内点(不是边界上的点) x ∗ x^* x满足对任意的 m i ( x ) < 0 , i = 1 , 2 , . . . , M m_{i}(x)\lt0,i=1,2,…,M mi(x)<0,i=1,2,...,M,则说明该问题满足slater条件。
  对slater条件做两点说明:

  1. 对于大多数的凸优化问题来说,slater condition都是成立的
  2. 放松的slater条件:如果约束函数 m i ( x ) , i = 1 , 2 , . . . , M m_{i}(x),i=1,2,…,M mi(x),i=1,2,...,M中有K个是仿射函数,则我们只需要那M-K个约束函数满足第一个条件,我们也说该问题满足slater条件。
  3. 什么是仿射函数?仿射函数,即最高次数为1的多项式函数。常数项为零的仿射函数称为线性函数。简单来说,就是比较简单的函数。

  对第一个条件进一步说明:为什么我们要满足这个条件?第一个条件放在G中的意思就是:G在u<0部分必须有点存在
在这里插入图片描述
假设纵坐标左边没有点,那么在我们寻找 g ( λ ) g(\lambda) g(λ)时,那条直线实际上就会是纵坐标轴。

3.KKT条件的证明

  通过上面的推导我们知道了:
在这里插入图片描述
满足强对偶关系之后我们就得到一个结论: d ∗ = p ∗ d^*=p^* d=p,但是也到此为止了,我们肯定得解出那些未知最优参数(带 * 的变量),KKT条件就是干这件事。
  KKT条件有三部分:可行条件、互补松弛条件以及偏导为0条件,我们一个一个推导。

3.1可行条件

  所谓可行条件,指的是一开始就满足的一些条件:
在这里插入图片描述

这三个条件肯定得满足,这个没啥可说的,天然满足

3.2互补松弛条件

  我们知道:
在这里插入图片描述
带星号的都是最优解的意思。根据可行条件我们知道: λ i ≥ 0 , m i ≤ 0 \lambda_{i}\geq0,m_{i}\leq0 λi0,mi0,所以 λ i m i ≤ 0 \lambda_{i}m_{i}\leq0 λimi0,所以上面继续变换:
在这里插入图片描述
因为最后一步等于第一步,所以中间推导步骤中的 ≤ \leq 都应该变成=,倒数第三步等于倒数第二步,所以我们有:
在这里插入图片描述
而前面我们又知道 λ i m i ≤ 0 \lambda_{i}m_{i}\leq0 λimi0,所以互补松弛条件如下:
在这里插入图片描述

3.3偏导为0条件

在这里插入图片描述
继续看这个推导,第二、三步之间应该用等号连接,即:
在这里插入图片描述
意思就是说L在 x = x ∗ x=x^* x=x处有最小值,于是偏导为0条件就出来了:
在这里插入图片描述
于是KKT条件为:
在这里插入图片描述
证毕。

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/182318.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)
blank

相关推荐

  • 一个简单完整的网页密码_简单的个人网页

    一个简单完整的网页密码_简单的个人网页获得源码链接,点击这里网页头部+banner和信息部分+新闻部分+底部一头部效果:先对css进行初始化分析:头部有一张图片和一个input输入框还有一个按钮+下面的通栏因为用到左浮,右浮的地方不同我们可以写一个通类这里的logo图片如果不定义宽高会影响下面的通栏的设置,影响其中的第一个为首的顺序无法对齐二、通栏(宽度为适应屏幕所以是10…

    2022年10月13日
  • linux下杀死全部进程,linux下批量杀死进程[通俗易懂]

    linux下杀死全部进程,linux下批量杀死进程[通俗易懂]psaux|greppython|grep-vgrep|cut-c9-15|xargskill-15管道符“|”用来隔开两个命令,管道符左边命令的输出会作为管道符右边命令的输入。下面说说用管道符联接起来的几个命令:“psaux”是linux里查看所有进程的命令。这时检索出的进程将作为下一条命令“greppython”的输入。“greppython”的输出结果是,所有含有关键…

  • idea激活码2021.10_最新在线免费激活「建议收藏」

    (idea激活码2021.10)好多小伙伴总是说激活码老是失效,太麻烦,关注/收藏全栈君太难教程,2021永久激活的方法等着你。IntelliJ2021最新激活注册码,破解教程可免费永久激活,亲测有效,下面是详细链接哦~https://javaforall.cn/100143.html1STL5S9V8F-eyJsaWNlbnNlSWQi…

  • dirsearch使用方法_查看es某个索引下的所有数据

    dirsearch使用方法_查看es某个索引下的所有数据下载地址:https://github.com/maurosoria/dirsearch要求安装并配置好:python3环境,python3环境,python3环境!!!使用方法,在cmd切换到安装目录pydirsearch.py-u网址-e语言(我一般用*)运行结果截图-u指定url-e指定网站语言例如pydirsearch.py-uwww.xxx.com…

  • https httppost ssl单向验证,json格式传输 java案例

    https httppost ssl单向验证,json格式传输 java案例https httppost ssl单向验证,json格式传输 java案例

  • mybatis 创建SqlSessionFactory

    mybatis 创建SqlSessionFactory创建SqlSessionFactory有两种方式:一种是通过XML,另一种是通过Java代码。Stringresource="org/mybatis/example/mybatis-config.xml";InputStreaminputStream=Resources.getResourceAsStream(resource);SqlSessionFactorysqlSess…

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号