机器学习笔记6-支持向量机(SVM)(1/3)
支持向量机(SVM)是一种二类分类模型,包括线性可分支持向量机、线性支持向量机及非线性支持向量机。当训练数据线性可分时,通过硬间隔最大化,学习一个线性的分类器,即线性可分支持向量机;当训练数据近似线性可分时,通过软间隔最大化,也学习一个线性分类器,即线性支持向量机;当训练数据线性不可分时,通过使用核技巧及软间隔最大化,学习非线性支持向量机。
- 线性可分支持向量机
__
(1)函数间隔:对于给定的训练数据集 T T T和超平面 ( w , b ) (w,b) (w,b),定义超平面 ( w , b ) (w,b) (w,b)关于样本点 ( x i , y i ) (x_i,y_i) (xi,yi)的函数间隔为 γ ^ i = y i ( w x i + b ) {\hat \gamma _i} = {y_i}(w{x_i} + b) γ^i=yi(wxi+b),定义超平面 ( w , b ) (w,b) (w,b)关于训练数据集 T T T的函数间隔为超平面 ( w , b ) (w,b) (w,b)关于 T T T中所有样本点 ( x i , y i ) (x_i,y_i) (xi,yi)的函数间隔之最小值,即 γ ^ = min i = 1 , 2... N γ ^ i \hat \gamma = \mathop {\min }\limits_{i = 1,2...N} {\hat \gamma _i} γ^=i=1,2...Nminγ^i。
__
(2)几何间隔:对于给定的训练数据集 T T T和超平面 ( w , b ) (w,b) (w,b),定义超平面 ( w , b ) (w,b) (w,b)关于样本点 ( x i , y i ) (x_i,y_i) (xi,yi)的几何间隔为 γ i = y i ( w ∥ w ∥ x i + b ∥ w ∥ ) {{ \gamma }_i} = {y_i}{\rm{(}}\frac{w}{{\left\| w \right\|}}{x_i} + \frac{b}{{\left\| w \right\|}}{\rm{)}} γi=yi(∥w∥wxi+∥w∥b), ∣ ∣ w ∣ ∣ ||w|| ∣∣w∣∣为 w w w的 L 2 L_2 L2范数。定义超平面 ( w , b ) (w,b) (w,b)关于训练数据集 T T T的几何间隔为超平面 ( w , b ) (w,b) (w,b)关于 T T T中所有样本点 ( x i , y i ) (x_i,y_i) (xi,yi)的几何间隔之最小值,即 γ = min i = 1 , 2... N γ i \gamma = \mathop {\min }\limits_{i = 1,2...N} {\gamma _i} γ=i=1,2...Nminγi。几何间隔的目的是为了对分离超平面的法向量规范化,防止 w w w和 b b b成比例改变时,超平面不变但函数间隔发生改变的现象的出现。
__
(3)间隔最大化:支持向量机的基本思想是求解能够正确划分训练数据集并且几何间隔最大的分离超平面。这一点和感知机的不同在于感知机只求能够正确划分训练数据集,因此感知机的解有无数个。而支持向量机还要要求几何间隔最大,因此支持向量机的解只有一个。间隔最大化的直观解释是:以充分大的确信度对训练数据进行分类。不仅要正确分类,而且对最难分的实例点也有足够大的确信度将其分开。这个问题可以表述成如下约束最优化问题:
m a x w , b γ \mathop {{\rm{max}}}\limits_{w,b} {\rm{ }}\gamma w,bmaxγ y i ( w ∥ w ∥ x i + b ∥ w ∥ ) ≥ γ , i = 1 , 2 . . . N {y_i}{\rm{(}}\frac{w}{{\left\| w \right\|}}{x_i} + \frac{b}{{\left\| w \right\|}}{\rm{)}} \ge \gamma {\rm{, i = 1,2}}...{\rm{N}} yi(∥w∥wxi+∥w∥b)≥γ,i=1,2...N以上可改写成函数间隔的形式:
m a x w , b γ ^ ∥ w ∥ \mathop {{\rm{max}}}\limits_{w,b} {\rm{ }}\frac{{\hat \gamma }}{{\left\| w \right\|}} w,bmax∥w∥γ^ y i ( w x i + b ) ≥ γ ^ , i = 1 , 2 . . . N {y_i}{\rm{(}}w{x_i} + b{\rm{)}} \ge \hat \gamma {\rm{, i = 1,2}}...{\rm{N}} yi(wxi+b)≥γ^,i=1,2...N函数间隔 γ ^ {\hat \gamma } γ^的取值并不影响最优化问题的解,因为它是随 w w w和 b b b成比例的变化的。取 γ ^ = 1 {\hat \gamma}=1 γ^=1,注意到 1 ∥ w ∥ \frac{1}{{\left\| w \right\|}} ∥w∥1与 1 2 ∥ w ∥ 2 \frac{1}{2}{\left\| w \right\|^2} 21∥w∥2是等价的,可得如下最优化问题
m i n w , b 1 2 ∥ w ∥ 2 \mathop {{\rm{min}}}\limits_{w,b} {\rm{ }}\frac{{\rm{1}}}{{\rm{2}}}{\left\| w \right\|^2} w,bmin21∥w∥2 y i ( w x i + b ) − 1 ≥ 0 , i = 1 , 2 . . . N {y_i}{\rm{(}}w{x_i} + b{\rm{) - 1}} \ge {\rm{0, i = 1,2}}...{\rm{N}} yi(wxi+b)−1≥0,i=1,2...N这是一个凸二次规划问题。(凸优化问题包含的内容很多,这里就不展开了,因为我了解的也不深。)
在线性可分情况下,训练数据集的样本点中与分离超平面距离最近的样本点的实例称为支持向量。支持向量都满足 y i ( w x i + b ) − 1 = 0 {y_i}{\rm{(}}w{x_i} + b{\rm{) - 1}} = {\rm{0}} yi(wxi+b)−1=0。在决定分离超平面时只有支持向量起作用,其它实例点不起作用,所以才将这种分类模型称为支持向量机。支持向量的个数一般很少,所以支持向量机由很少的训练样本确定。
(4)学习的对偶算法:为了求解以上凸二次规划问题,将它作为原始最优化问题,应用拉格朗日对偶性,通过求解对偶问题得到原始问题的最优解。这样做的优点,一是对偶问题往往更容易求解;而是可以引入核函数,进而推广到非线性分类问题。首先构建拉格朗日函数:
L ( w , b , α ) = 1 2 ∥ w ∥ 2 − ∑ i = 1 N α i y i ( w x i + b ) + ∑ i = 1 N α i L(w,b,\alpha ) = \frac{{\rm{1}}}{{\rm{2}}}{\left\| w \right\|^2} - \sum\limits_{i = {\rm{1}}}^N {{\alpha _i}{y_i}(w{x_i} + b)} + \sum\limits_{i = {\rm{1}}}^N {{\alpha _i}} L(w,b,α)=21∥w∥2−i=1∑Nαiyi(wxi+b)+i=1∑Nαi其中, α = ( α 1 , α 2 , . . . , α N ) T \alpha = {({\alpha _1},{\alpha _2},...,{\alpha _N})^T} α=(α1,α2,...,αN)T为拉格朗日乘子向量。根据拉格朗日对偶性,原始问题的对偶问题是极大极小问题:
m a x α min w , b L ( w , b , α ) \mathop {max}\limits_\alpha \mathop {\min }\limits_{w,b} L(w,b,\alpha ) αmaxw,bminL(w,b,α)所以为了得到对偶问题的解,需要先求 L ( w , b , α ) L(w,b,\alpha ) L(w,b,α)对 w , b w,b w,b的极小,再求对 α \alpha α的极大。
A. 求 min w , b L ( w , b , α ) \mathop {\min }\limits_{w,b} L(w,b,\alpha ) w,bminL(w,b,α)。将拉格朗日函数 L ( w , b , α ) L(w,b,\alpha ) L(w,b,α)分别对 w , b w,b w,b求偏导数并令其等于零。可得 w = ∑ i = 1 N α i y i x i w = \sum\limits_{i = {\rm{1}}}^N {{\alpha _i}{y_i}{x_i}} w=i=1∑Nαiyixi, ∑ i = 1 N α i y i = 0 \sum\limits_{i = {\rm{1}}}^N {{\alpha _i}{y_i}} = 0 i=1∑Nαiyi=0,进而可得
min w , b L ( w , b , α ) = − 1 2 ∑ i = 1 N ∑ i = 1 N α i α j y i y j ( x i ⋅ x j ) + ∑ i = 1 N α i \mathop {\min }\limits_{w,b} L(w,b,\alpha ) = - \frac{{\rm{1}}}{{\rm{2}}}\sum\limits_{i = {\rm{1}}}^N {\sum\limits_{i = {\rm{1}}}^N {{\alpha _i}{\alpha _j}{y_i}{y_j}({x_i}\cdot{x_j})} } + \sum\limits_{i = {\rm{1}}}^N {{\alpha _i}} w,bminL(w,b,α)=−21i=1∑Ni=1∑Nαiαjyiyj(xi⋅xj)+i=1∑Nαi
B. 求 min w , b L ( w , b , α ) \mathop {\min }\limits_{w,b} L(w,b,\alpha ) w,bminL(w,b,α)对 α \alpha α的极大,即是对偶问题:
min α 1 2 ∑ i = 1 N ∑ i = 1 N α i α j y i y j ( x i ⋅ x j ) − ∑ i = 1 N α i \mathop {\min }\limits_\alpha \frac{{\rm{1}}}{{\rm{2}}}\sum\limits_{i = {\rm{1}}}^N {\sum\limits_{i = {\rm{1}}}^N {{\alpha _i}{\alpha _j}{y_i}{y_j}({x_i}\cdot{x_j})} } - \sum\limits_{i = {\rm{1}}}^N {{\alpha _i}} αmin21i=1∑Ni=1∑Nαiαjyiyj(xi⋅xj)−i=1∑Nαi ∑ i = 1 N α i y i = 0 \sum\limits_{i = {\rm{1}}}^N {{\alpha _i}{y_i}} = 0 i=1∑Nαiyi=0 α i ≥ 0 {{\alpha _i} \ge 0} αi≥0
可以证明对偶问题和原始问题的解是一致的。
若 α ∗ = ( α 1 ∗ , α 2 ∗ , . . . , α N ∗ ) T {\alpha ^*} = {(\alpha _{_1}^*,\alpha _{_2}^*,...,\alpha _{_N}^*)^T} α∗=(α1∗,α2∗,...,αN∗)T是对偶问题的解,则存在 α j ∗ > 0 \alpha _{_j}^*>0 αj∗>0,使得原始问题的解为: w ∗ = ∑ i = 1 N α i ∗ y i x i {w^*} = \sum\limits_{i = {\rm{1}}}^N {\alpha _i^*{y_i}{x_i}} w∗=i=1∑Nαi∗yixi, b ∗ = y j − ∑ i = 1 N α i ∗ y i ( x i ⋅ x j ) {b^*} = {y_j} - \sum\limits_{i = {\rm{1}}}^N {\alpha _i^*{y_i}{\rm{(}}{x_i}\cdot{x_j}{\rm{)}}} b∗=yj−i=1∑Nαi∗yi(xi⋅xj)。分离超平面可以写成 ∑ i = 1 N α i ∗ y i ( x ⋅ x i ) + b ∗ = 0 \sum\limits_{i = {\rm{1}}}^N {\alpha _i^*{y_i}{\rm{(}}x \cdot {x_i}{\rm{) + }}} {b^*} = 0 i=1∑Nαi∗yi(x⋅xi)+b∗=0,其只依赖于输入 x x x和训练样本输入的内积。由于 ( w ∗ , b ∗ ) {\rm{(}}{w^*},{b^*}{\rm{)}} (w∗,b∗)只依赖于训练数据中对应于 α j ∗ > 0 \alpha _{_j}^*>0 αj∗>0的样本点 ( x i , y i ) (x_i,y_i) (xi,yi),其它样本点对 ( w ∗ , b ∗ ) {\rm{(}}{w^*},{b^*}{\rm{)}} (w∗,b∗)无影响。因此将训练数据中对应于 α j ∗ > 0 \alpha _{_j}^*>0 αj∗>0的样本点称为支持向量。
补充
- 拉格朗日对偶性
假设 f ( x ) f(x) f(x), c i ( x ) c_i(x) ci(x), h j ( x ) h_j(x) hj(x)是定义在 R n R^n Rn上的连续可微函数。考虑约束最优化问题
m i n x ∈ R n f ( x ) \mathop {{\rm{min}}}\limits_{x \in {R^n}} {\rm{ }}f(x) x∈Rnminf(x) c i ( x ) ≤ 0 , i = 1 , 2 . . . k {c_i}{\rm{(x)}} \le {\rm{0, }}i{\rm{ = 1,2}}...{\rm{k}} ci(x)≤0,i=1,2...k h j ( x ) ≥ 0 , j = 1 , 2 . . . l {{\rm{h}}_j}(x) \ge 0,j = {\rm{1,2}}...l hj(x)≥0,j=1,2...l此为原始问题。引进广义拉格朗日函数
L ( x , α , β ) = f ( x ) + ∑ i = 1 k α i c i ( x ) + ∑ j = 1 l β j h j ( x ) L(x,\alpha ,\beta ) = f(x) + \sum\limits_{i = {\rm{1}}}^k {{\alpha _i}{c_i}{\rm{(x)}}} + \sum\limits_{j = {\rm{1}}}^l {{\beta _j}{h_j}{\rm{(x)}}} L(x,α,β)=f(x)+i=1∑kαici(x)+j=1∑lβjhj(x)这里 α i \alpha _i αi, β j \beta _j βj是拉格朗日乘子, α i ≥ 0 \alpha _i \ge 0 αi≥0,考虑 x x x的函数:
θ P ( x ) = m a x α β : α i ≥ 0 L ( x , α , β ) {\theta _P}(x) = \mathop {max}\limits_{\alpha \beta :{\alpha _i} \ge 0} L(x,\alpha ,\beta ) θP(x)=αβ:αi≥0maxL(x,α,β)对于违反约束条件的 x x x, θ P ( x ) {\theta _P}(x) θP(x)会变为无穷大;而对于满足约束条件的 x x x,则 θ P ( x ) = f ( x ) {\theta _P}(x)=f(x) θP(x)=f(x)。所以考虑极小化问题
min x θ P ( x ) = min x m a x α β : α i ≥ 0 L ( x , α , β ) \mathop {\min }\limits_x {\rm{ }}{\theta _P}(x) = \mathop {\min }\limits_x \mathop {max}\limits_{\alpha \beta :{\alpha _i} \ge 0} L(x,\alpha ,\beta ) xminθP(x)=xminαβ:αi≥0maxL(x,α,β)与原始问题是等价的,即它们拥有相同的解。这个称为广义拉格朗日函数的极小极大问题。相应的,有广义拉格朗日函数的极大极小问题,即
max α β : α i ≥ 0 θ D ( α , β ) = m a x α β : α i ≥ 0 min x L ( x , α , β ) \mathop {\max }\limits_{\alpha \beta :{\alpha _i} \ge 0} {\theta _D}(\alpha ,\beta ) = \mathop {max}\limits_{\alpha \beta :{\alpha _i} \ge 0} \mathop {\min }\limits_x L(x,\alpha ,\beta ) αβ:αi≥0maxθD(α,β)=αβ:αi≥0maxxminL(x,α,β) θ D ( α , β ) = min x L ( x , α , β ) {\theta _D}(\alpha ,\beta ) = \mathop {\min }\limits_x L(x,\alpha ,\beta ) θD(α,β)=xminL(x,α,β)将上式表示为约束最优化问题:
max α , β θ D ( α , β ) = m a x α , β min x L ( x , α , β ) \mathop {\max }\limits_{\alpha ,\beta } {\theta _D}(\alpha ,\beta ) = \mathop {max}\limits_{\alpha ,\beta } \mathop {\min }\limits_x L(x,\alpha ,\beta ) α,βmaxθD(α,β)=α,βmaxxminL(x,α,β) α i ≥ 0 , i = 1 , 2 , . . . , k {\alpha _i} \ge 0,i=1,2,...,k αi≥0,i=1,2,...,k称为原始问题的对偶问题。若原始问题与对偶问题都有最优值,则
d ∗ = m a x α , β : α i ≥ 0 min x L ( x , α , β ) ≤ min x m a x α , β : α i ≥ 0 L ( x , α , β ) = p ∗ {d^*} = \mathop {max}\limits_{\alpha ,\beta :{\alpha _i} \ge 0} \mathop {\min }\limits_x L(x,\alpha ,\beta ) \le \mathop {\min }\limits_x \mathop {max}\limits_{\alpha ,\beta :{\alpha _i} \ge 0} L(x,\alpha ,\beta ) = {p^*} d∗=α,β:αi≥0maxxminL(x,α,β)≤xminα,β:αi≥0maxL(x,α,β)=p∗当 x ∗ , α ∗ , β ∗ x^*,\alpha^*,\beta^* x∗,α∗,β∗满足KKT条件时,原始问题和对偶问题的解是一致的,都是 x ∗ , α ∗ , β ∗ x^*,\alpha^*,\beta^* x∗,α∗,β∗。