SVM 白板推导| 由最大间隔化目标演化的损失函数推导过程

SVM 英文名为 Support Vector Machine,叫支持向量机,SVM 是一个二分类模型。

感知机学习算法会因为初始值不同而得到不同的超平面,而 SVM 试图寻找一个最佳超平面,来划分数据,怎么才算最佳呢?我们自然能想到,距离样本点距离尽可能的远,模型的泛化性能和鲁棒性更强为最好,而且可以找到一条唯一的超平面。

如下图可以看出,感知机的超平面可以有无数条,SVM只要找到一条距离样本点最大间隔的超平面即可。

SVM 目标

我们知道,感知机是一个现在网上已经有很多最初被提出是为了解决二分类问题,假设如下有正负样本。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-fITE6hAO-1585487278443)(https://imgkr.cn-bj.ufileos.com/14067463-82d2-4afb-8219-ed367fdddef9.png)]

由此,我们知道,SVM 的目标就是要找到一条间隔最大化的超平面,使得样本点距离超平面尽可能的远,那么,我们如何找到间隔最大的超平面呢?

首先,我们要先想到,间隔最大,即是求样本点到超平面的距离

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-9MfW4IKZ-1585487278447)(https://imgkr.cn-bj.ufileos.com/b8e5891b-462c-427e-885d-420dbd3be1a5.png)]

而现在,我们有 N 个样本点,则我们可能有 N 个距离,那我们最终想要的就是,离超平面距离最小的样本点的最大值。

既然是最大间隔分类器,我们拆开来解释:

最大:找到距离样本点最大的间隔

间隔:N 个样本点距离超平面的最小值

分类器:既然是二分类问题,那么我们的超平面就是要将样本分开,对于超平面 y = w T x + b y=w^Tx+b 而言,我们可以假设有两条超平面,一条是 w T x + b > 0 w^Tx+b>0 w T x + b < 0 w^Tx+b<0 ,我们可以引入函数 y i = { 1 + 1 yi=\{_{-1}^{+1} ,将两个超平面用一个函数表示,如图。

现在就需要我们想办法计算样本点到超平面的距离,如下图,表示点 ( x i , y i ) (x_i,y_i) 到超平面 y = w T x + b y=w^Tx+b 的距离的计算,最终将目标转化为找到 N 个样本点距离超平面的最小距离。

m i n d i s t a n c e ( w , b , x ) min-distance(w,b,x) 代入 m a x m a r g i n ( w , b ) max-margin(w,b) ,条件是 y i = w T x i + b > 0 y_i=w^Tx_i+b>0 ,并且是恒大于0,由此,我们可以将下图的式子转换一下,将模去掉。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-zJU21QEm-1585487278450)(https://imgkr.cn-bj.ufileos.com/e5e1d1b2-b0e3-462d-80c9-4bbe6faa9fe4.png)]

并且我们还能看出, 1 w \frac{1}{||w||} x i x_i 无关,我们可以将它提到前面去。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-Xx5ogZPC-1585487278451)(https://imgkr.cn-bj.ufileos.com/66eff5b1-c794-45a0-897e-a0b520ee209a.png)]

这样一来,我们的问题就转化成求解一个带约束的优化问题。对于条件 y i ( w T x i + b ) > 0 y_i(w^Tx_i+b)>0 而言,一定存在一个 γ,使得 m i n ( y i ( w T x i + b ) > 0 ) = γ min(y_i(w^Tx_i+b)>0)=γ ,为了简化计算,我们可以令 y i ( w T x i + b ) = γ y_i(w^Tx_i+b)=γ

假设令 γ=1,则我们可以得到最终的式子只与 W 有关,我们要求损失函数,则将其转化成最小化问题。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-hbAatu1C-1585487278454)(https://imgkr.cn-bj.ufileos.com/2700d10c-74b4-43bf-b227-dbdaeacd8f92.png)]

可能有的人会疑惑,为什么 γ=1可行,我们知道, y = w T x + b y=w^Tx+b 和$ y=2w^Tx+2b$ 表示的是同一个超平面,我们如果将 2-范数 固定下来,这样去指定 y = w T x + b y=w^Tx+b 的值是可以确定下来,否则有无穷多个,因为可以随机缩放,因为 γ > 0,所以无论 γ 等于多少,都可以任意比例缩放成 1,对整个等式无影响。

我们可以看到,目标函数是二次的,约束条件是线性的,所以它是一个带约束的凸二次规划问题,这个问题可以用现成的 QP 优化包进行求解,一言以蔽之:在一定约束下,目标最优,损失最小。

引入拉格朗日乘子

通过拉格朗日乘子,我们已经将带约束的原问题转化成了无约束的问题。

那么什么是拉格朗日对偶性呢?简单来讲,通过给每个约束条件加上拉格朗日乘子 λ,定义拉格朗日函数 (通过拉格朗日函数将约束条件融合到目标函数中,从而只用一个函数表达式便能清楚的表达出我们的问题)。

我们可以看到,原式子是一个关于 w 的凸二次函数,约束条件是 N 个关于 x 的不等式。通过引入拉格朗日乘子之后,我们可以将有约束问题转化为一个无约束的,只与 λ 有关的式子。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-NCp6U9NB-1585487278467)(https://imgkr.cn-bj.ufileos.com/07ad9910-b401-4e9f-bd99-bba4ffa54071.png)]

我们可以简单的理解一下,为什么可以这样定义,假设如下是关于 w,b 的集合,我们可以将其表示为两个部分,一部分是 △>0,一部分是 △<=0。

当 △>0 时,即 1 y i ( w T x i + b ) > 0 1-y_i(w^Tx_i+b)>0 ,显然 m a x L ( w , b , λ ) = 1 2 w T w + = maxL(w,b,λ)=\frac{1}{2}w^Tw+∞=∞ ,此时,对于我们求解无意义,我们也得不到最小化的值。

当 △<=0 时,即 1 y i ( w T x i + b ) < = 0 1-y_i(w^Tx_i+b)<=0 ,显然 m a x L ( w , b , λ ) = 1 2 w T w + 0 = 1 2 w T w maxL(w,b,λ)=\frac{1}{2}w^Tw+0=\frac{1}{2}w^Tw ,我们可以看到,与原问题相符。

所以,当我们谈到满足最优化的条件时,通常情况都是 △<=0 的情况。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-3celtfL0-1585487278470)(https://imgkr.cn-bj.ufileos.com/e946fbbe-ed53-4aa3-a057-d9d00595cf44.png)]

此外,由于这个问题的特殊结构,还可以通过拉格朗日对偶性变换到对偶变量的优化问题,即通过求解与原问题等价的对偶问题得到原始问题的最优解,这就是线性可分条件下支持向量机的对偶算法,这样做的优点在于:一是对偶问题往往更容易求解,二是可以自然的引入核函数进而推广到非线性分类问题。

对偶问题与原问题

这里用 p p^* 表示原问题的最优解,且和最初的问题是等价的。如果直接求解,那么一上来便得面对 w , b w,b 两个参数,而 λ 又是不等式约束,这个求解过程不好做。不放把最小和最大的位置交换下。

交换以后的新问题是原问题的对偶问题,对偶问题的最优解用 d d^* 来表示,而且有 d < = P d^* <= P^* ,在满足某些条件的情况下,这两者相等,这个时候,口可以通过求解对偶问题来间接求解原问题。

换言之,之所以从 m i n min m a x max 原始问题 p p^* 转化为 m a x max m i n min 的对偶问题 d d^* ,一者是因为 d d^* p p^* 是近似解,二者,转化成对偶问题更容易求解。

下面可以先求 L L w , b w,b 的极小,再求 L L 对 λ 的极大。

KKT

上面提到 d < = P d^* <= P^* 在满足某些条件的情况下,两者相等,某些条件就是 KKT 条件。

对偶问题求解的 3 个步骤

1,固定 λ,让 L L 关于 w , b w,b 最小化,对 w , b w,b 求偏导,令 δ L δ w = 0 \frac{δL}{δw}=0 δ L δ b = 0 \frac{δL}{δb}=0

将以上结果代入 L,最终结果我们可以看到,最终的结果只与 x 有关。

如何理解这一结果呢?我们可以看看 KKT 条件中的其中一个。

当 λ=0 时,此时的 L 只与 w 有关,对超平面无限制。

当 λ>0 时,此时 1 y i ( w T x i + b ) = 0 1-y_i(w^Tx_i+b)=0 ,因为 y i = { 1 + 1 yi=\{_{-1}^{+1} ,此时,L 与在边界上的点有关,即只与支持向量有关。

2, 对 λ 求极大,即是关于对偶问题的最优化问题,经过上面的第一个步骤的求 w 和 b,得到拉格朗日函数式子已经没有了变量 w,b,只有 λ。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-kizmGXpW-1585487278493)(https://imgkr.cn-bj.ufileos.com/389ea7e3-c7ad-4713-9407-75f423f0fab2.png)]

并且我们也在之前求导过程中求出了 w。

将上述式子代入 L 中,可以求出 w 和 b

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-yD11gXZv-1585487278504)(https://imgkr.cn-bj.ufileos.com/d9fec992-30cd-4c99-8bdc-9b3c9d5ce914.png)]

最终,我们可以得到新的分类决策函数。

3, 在求得 L ( w , b , λ ) L(w,b,λ) 关于 w w b b 最小化,以及对 λ λ 极大之后,最后一步便是利用 SMO 算法求解对偶问题中的拉格朗日乘子 λ。

我们的拉格朗日函数中,在 λ i = λ 1 , λ 2 , . . . , λ n λ_i={λ_1,λ_2,...,λ_n} 上求上述目标函数的最小值,为了求解这些乘子,每次从中任意抽取两个乘子 λ 1 , λ 2 λ_1,λ_2 ,然后固定 λ 1 , λ 2 λ_1,λ_2 以外的乘子 λ i = λ 3 , λ 4 , . . . , λ n λ_i={λ_3,λ_4,...,λ_n} ,使得目标函数只关于 λ 1 , λ 2 λ_1,λ_2 的函数,这样,不断的从一堆乘子中任意抽取两个求解,不断的迭代求乘子问题,最终达到求解原问题的目的。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-ZtLgum3G-1585487278505)(https://imgkr.cn-bj.ufileos.com/0cb2d5eb-7618-4921-bf6f-b8d25d60d588.png)]

如果您觉得文章对你有帮助,欢迎点赞转发关注。

发布了30 篇原创文章 · 获赞 12 · 访问量 1万+

猜你喜欢

转载自blog.csdn.net/WBST5/article/details/105186329