深入了解 机器学习:支持向量机(Support Vector Machine, SVM)
支持向量机(Support Vector Machine, SVM)
1. 概要
1.1 简介
自从大半年前接触到SVM以来,感觉一直没怎么把SVM整明白。直到最近上的《模式识别》课程才仿佛打通了我的任督二脉,使我终于搞清楚了SVM的来龙去脉,所以写个博客作个总结。
SVM是什么? 先来看看*上对SVM的定义:
支持向量机(英语:support vector machine,常简称为SVM,又名支持向量网络)是在分类与回归分析中分析数据的监督式学习模型与相关的学习算法。给定一组训练实例,每个训练实例被标记为属于两个类别中的一个或另一个,SVM训练算法创建一个将新的实例分配给两个类别之一的模型,使其成为非概率二元线性分类器。SVM模型是将实例表示为空间中的点,这样映射就使得单独类别的实例被尽可能宽的明显的间隔分开。然后,将新的实例映射到同一空间,并基于它们落在间隔的哪一侧来预测所属类别。
如果从未接触SVM的话,维基的这一大段解释肯定会让你一头雾水。简单点讲,SVM就是一种二类分类模型,他的基本模型是的定义在特征空间上的间隔最大的线性分类器,SVM的学习策略就是间隔最大化。
1.2 直观理解
我们先来看看下面这个图:
图1.1
图中有分别属于两类的一些二维数据点和三条直线。如果三条直线分别代表三个分类器的话,请问哪一个分类器比较好?
我们凭直观感受应该觉得答案是H3。首先H1不能把类别分开,这个分类器肯定是不行的;H2可以,但分割线与最近的数据点只有很小的间隔,如果测试数据有一些噪声的话可能就会被H2错误分类(即对噪声敏感、泛化能力弱)。H3以较大间隔将它们分开,这样就能容忍测试数据的一些噪声而正确分类,是一个泛化能力不错的分类器。
对于支持向量机来说,数据点若是pp维向量,我们用p−1p−1维的超平面来分开这些点。但是可能有许多超平面可以把数据分类。最佳超平面的一个合理选择就是以最大间隔把两个类分开的超平面。因此,SVM选择能够使离超平面最近的数据点的到超平面距离最大的超平面。
以上介绍的SVM只能解决线性可分的问题,为了解决更加复杂的问题,支持向量机学习方法有一些由简至繁的模型:
- 线性可分SVM
当训练数据线性可分时,通过硬间隔(hard margin,什么是硬、软间隔下面会讲)最大化可以学习得到一个线性分类器,即硬间隔SVM,如上图的的H3。
- 线性SVM
当训练数据不能线性可分但是可以近似线性可分时,通过软间隔(soft margin)最大化也可以学习到一个线性分类器,即软间隔SVM。
- 非线性SVM
当训练数据线性不可分时,通过使用核技巧(kernel trick)和软间隔最大化,可以学习到一个非线性SVM。
2. 线性可分SVM——硬间隔
考虑如下形式的线性可分的训练数据集:
(X1,y1),(X2,y2),…,(Xn,yn)(X1,y1),(X2,y2),…,(Xn,yn)
其中 XiXi是一个含有dd个元素的列向量, 即Xi∈RdXi∈Rd; yiyi是标量, y∈+1,−1y∈+1,−1, yi=+1yi=+1时表示XiXi属于正类别, yi=−1yi=−1时表示XiXi属于负类别。
注: 本文中, XX、XiXi、WW等都是(列)向量,有的文章一般用 xixi 表示一个向量而用 XX 表示所有 xixi 组成的一个矩阵,注意区分。
回忆一下感知机的目标: 找到一个超平面使其能正确地将每个样本正确分类。感知机使用误分类最小的方法求得超平面,不过此时解有无穷多个(例如图1.1的H2和H3以及它俩的任意线性组合)。而线性可分支持向量机利用间隔最大化求最优分离超平面,这时解是唯一的。
2.1 超平面与间隔
一个超平面由法向量WW和截距bb决定,其方程为XTW+b=0XTW+b=0, 可以规定法向量指向的一侧为正类,另一侧为负类。下图画出了三个平行的超平面,法方向取左上方向。
注意: 如果XX和WW都是列向量,即XTWXTW会得到XX和WW的点积(dot product, 是一个标量),等价于X⋅WX⋅W和W⋅XW⋅X。
图2.1
为了找到最大间隔超平面,我们可以先选择分离两类数据的两个平行超平面,使得它们之间的距离尽可能大。在这两个超平面范围内的区域称为“间隔(margin)”,最大间隔超平面是位于它们正中间的超平面。这个过程如上图所示。
2.2 间隔最大化
将高数里面求两条平行直线的距离公式推广到高维可求得图2.1中margin的ρρ:
margin=ρ=2||W||(2.2.1)(2.2.1)margin=ρ=2||W||
我们的目标是使ρρ最大, 等价于使ρ2ρ2最大:
maxW,bρ⟺maxW,bρ2⟺minW,b12||W||2(2.2.2)(2.2.2)maxW,bρ⟺maxW,bρ2⟺minW,b12||W||2
上式的1212是为了后续求导后刚好能消去,没有其他特殊意义。
同时也不要忘了有一些约束条件:
XTiW+b≥+1,yi=+1XTiW+b≤−1,yi=−1(2.2.3)(2.2.3)XiTW+b≥+1,yi=+1XiTW+b≤−1,yi=−1
总结一下,间隔最大化问题的数学表达就是
minW,bJ(W)=minW,b12||W||2s.t.yi(XTiW+b)≥1,i=1,2,…n.(2.2.4)(2.2.4)minW,bJ(W)=minW,b12||W||2s.t.yi(XiTW+b)≥1,i=1,2,…n.
通过求解上式即可得到最优超平面 W^W^ 和 b^b^ 。具体如何求解见2.4和2.5节。
2.3 支持向量
在线性可分的情况下,训练数据集的样本点中与分离超平面距离最近的数据点称为支持向量(support vector),支持向量是使(2.2.4)(2.2.4)中的约束条件取等的点,即满足
yi(XTiW+b)=1(2.3.1)(2.3.1)yi(XiTW+b)=1
的点。也即所有在直线XTW+b=1XTW+b=1或直线XTW+b=−1XTW+b=−1的点。如下图所示:
图2.2
在决定最佳超平面时只有支持向量起作用,而其他数据点并不起作用(具体推导见2.4节最后)。如果移动非支持向量,甚至删除非支持向量都不会对最优超平面产生任何影响。也即支持向量对模型起着决定性的作用,这也是“支持向量机”名称的由来。
2.4 对偶问题
如何求解式 (2.2.4)(2.2.4) 呢?
我们称式 (2.2.4)(2.2.4) 所述问题为原始问题(primal problem), 可以应用拉格朗日乘子法构造拉格朗日函数(Lagrange function)再通过求解其对偶问题(dual problem)得到原始问题的最优解。转换为对偶问题来求解的原因是:
- 对偶问题更易求解,由下文知对偶问题只需优化一个变量αα且约束条件更简单;
- 能更加自然地引入核函数,进而推广到非线性问题。
首先构建拉格朗日函数。为此需要引进拉格朗日乘子(Lagrange multiplier) αi≥0,i=1,2,…nαi≥0,i=1,2,…n。则拉格朗日函数为:
L(W,b,α)=12||W||2−∑i=1nαi[yi(XTiW+b)−1](2.4.1)(2.4.1)L(W,b,α)=12||W||2−∑i=1nαi[yi(XiTW+b)−1]
因此,给定一个WW和bb, 若不满足式(2.2.4)(2.2.4)的约束条件,那么有
maxαL(W,b,α)=+∞(2.4.2)(2.4.2)maxαL(W,b,α)=+∞
否则,若满足式(2.2.4)(2.2.4)的约束条件,有
maxαL(W,b,α)=J(W)=12||W||2(2.4.3)(2.4.3)maxαL(W,b,α)=J(W)=12||W||2
结合式(2.4.2)(2.4.2)和(2.4.3)(2.4.3)知,优化问题
minW,bmaxαL(W,b,α)(2.4.4)(2.4.4)minW,bmaxαL(W,b,α)
与式(2.2.4)(2.2.4)所述问题是完全等价的。
根据拉格朗日对偶性,式(2.4.4)(2.4.4)所述问题即原始问题的对偶问题是:
maxαminW,bL(W,b,α)(2.4.5)(2.4.5)maxαminW,bL(W,b,α)
以上具体推导细节可参见书籍《统计学习方法》或者知乎文章拉格朗日对偶性
为了求得对偶问题的解,需要先求得L(W,b,α)L(W,b,α)对WW和bb的极小再求对αα的极大。
(1) 求minW,bL(W,b,α)minW,bL(W,b,α):
对拉格朗日函数求导并令导数为0,有:
∇WL(W,b,α)=W−∑i=1nαiyiXi=0⟹W=∑i=1nαiyiXi(2.4.6)(2.4.6)∇WL(W,b,α)=W−∑i=1nαiyiXi=0⟹W=∑i=1nαiyiXi
∇bL(W,b,α)=−∑i=1nαiyi=0⟹∑i=1nαiyi=0(2.4.7)(2.4.7)∇bL(W,b,α)=−∑i=1nαiyi=0⟹∑i=1nαiyi=0
将上面两式代入L(W,b,α)L(W,b,α):
所以,
minW,bL(W,b,α)=−12∑i=1n∑j=1nαiαjyiyjXTiXj + ∑i=1nαi(2.4.8)(2.4.8)minW,bL(W,b,α)=−12∑i=1n∑j=1nαiαjyiyjXiTXj + ∑i=1nαi
(2) 求minW,bL(W,b,α)minW,bL(W,b,α) 对αα的极大:
等价于式(2.4.8)(2.4.8)对αα求极大,也等价于式(2.4.8)(2.4.8)取负数后对αα求极小,即
minα12∑i=1n∑j=1nαiαjyiyjXTiXj − ∑i=1nαi(2.4.9)(2.4.9)minα12∑i=1n∑j=1nαiαjyiyjXiTXj − ∑i=1nαi
同时满足约束条件:
∑i=1nαiyi=0αi≥0,i=1,2,…,n.(2.4.10)(2.4.10)∑i=1nαiyi=0αi≥0,i=1,2,…,n.
至此,我们得到了原始最优化问题(2.2.4)(2.2.4)和对偶最优化问题(2.4.9)(2.4.9)、(2.4.10)(2.4.10)。
由slater条件知,因为原始优化问题的目标函数和不等式约束条件都是凸函数,并且该不等式约束是严格可行的(因为数据是线性可分的), 所以存在 W^W^, b^b^ ,α^α^,使得W^W^, b^b^是原始问题的解,α^α^是对偶问题的解。这意味着求解原始最优化问题(2.2.4)(2.2.4)可以转换为求解对偶最优化问题(2.4.9)(2.4.9)、(2.4.10)(2.4.10)。
slater 条件:
原始问题一般性表达为minxf(x)s.t. ci(x)≤0,i=1,2,…khj(x)=0,j=1,2,…,lminxf(x)s.t. ci(x)≤0,i=1,2,…khj(x)=0,j=1,2,…,l
则其拉格朗日函数为L(x,α,β)=f(x)+∑i=1kαici(x)+∑j=1lβjhj(x),αi≥0L(x,α,β)=f(x)+∑i=1kαici(x)+∑j=1lβjhj(x),αi≥0
假设原始问题目标函数 f(x)f(x) 和不等式约束条件 ci(x)ci(x)都是凸函数,原始问题等式约束hj(x)hj(x)都是仿射函数,且不等式约束 ci(x)ci(x)是严格可行的,即存在 xx ,对所有 ii 都有 ci(x)<0ci(x)<0 ,则存在 x^x^, α^α^, β^β^,使 x^x^ 是原始问题的解, α^α^, β^β^是对偶问题的解。
那么如何求解优化问题(2.4.9)(2.4.9)、(2.4.10)(2.4.10)的最优解 α^α^ 呢?
不难发现这是一个二次规划问题,有现成的通用的算法来求解。
事实上通用的求解二次规划问题的算法的复杂度正比于训练数据样本数,所以在实际应用中需要寻求更加高效的算法,例如序列最小优化(Sequential Minimal Optimiation, SMO)算法。
假设我们现在求得了(2.4.9)(2.4.9)、(2.4.10)(2.4.10)的最优解 α^α^,则根据式(2.4.6)(2.4.6)可求得最优W^W^:
W^=∑i=1nα^iyiXi(2.4.11)(2.4.11)W^=∑i=1nα^iyiXi
因为至少存在一个 α^j>0α^j>0(若不存在,即 α^α^ 全为0,则 W^=0W^=0, 即 margin=2||W^||=∞margin=2||W^||=∞,显然不行), 再根据KKT条件,即
⎧⎩⎨⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪乘子非负:αi≥0(i=1,2,…n.下同)约束条件:yi(XTiW+b)−1≥0互补条件:αi(yi(XTiW+b)−1)=0{乘子非负:αi≥0(i=1,2,…n.下同)约束条件:yi(XiTW+b)−1≥0互补条件:αi(yi(XiTW+b)−1)=0
所以至少存在一个 jj , 使 yj(XTjW^+b^)−1=0yj(XjTW^+b^)−1=0, 即可求得最优 b^b^:
b^=1yj−XTjW^=yj−XTjW^=yj−∑i=1nα^iyiXTjXi(2.4.12)(2.4.12)b^=1yj−XjTW^=yj−XjTW^=yj−∑i=1nα^iyiXjTXi
至此,所以我们就求得了整个线性可分SVM的解。求得的分离超平面为:
∑i=1nα^iyiXTXi+b^=0(2.4.13)(2.4.13)∑i=1nα^iyiXTXi+b^=0
则分类的决策函数为
f(X)=sign(∑i=1nα^iyiXTXi+b^)(2.4.14)(2.4.14)f(X)=sign(∑i=1nα^iyiXTXi+b^)
再来分析KKT条件里的互补条件,对于任意样本 (Xi,yi)(Xi,yi) ,总会有 αi=0αi=0 或者 yif(Xi)=yi(XTiW^+b)=1yif(Xi)=yi(XiTW^+b)=1。则有
- 若αi=0αi=0,此样本点不是支持向量,对模型没有任何作用;
- 若αi>0αi>0,此样本点位于最大间隔边界上,是一个支持向量,如下图所示。
图2.3
此外,当样本点是非支持向量时,因为αi=0αi=0, 所以SVM的解中的求和项中第 ii 项就为0,所以SVM的解(2.4.11)(2.4.11)、(2.4.12)(2.4.12)可简化为如下形式:
W^=∑i∈SVα^iyiXi(2.4.15)(2.4.15)W^=∑i∈SVα^iyiXi
b^=yj−∑i∈SVα^iyiXTjXi(2.4.16)(2.4.16)b^=yj−∑i∈SVα^iyiXjTXi
类似的,判别函数也可转换成如下形式:
f(X)=sign(∑i∈SVα^iyiXTXi+b^)(2.4.17)(2.4.17)f(X)=sign(∑i∈SVα^iyiXTXi+b^)
所以,整个SVM的解只与支持向量SV有关,与非支持向量无关。这也就解释了2.3节的结论,即在决定最佳超平面时只有支持向量起作用,而其他数据点并不起作用。
3. 线性SVM——软间隔
在前面的讨论中,我们一直假定训练数据是严格线性可分的,即存在一个超平面能完全将两类数据分开。但是现实任务这个假设往往不成立,例如下图所示的数据。
图3.1
3.1 软间隔最大化
解决该问题的一个办法是允许SVM在少量样本上出错,即将之前的硬间隔最大化条件放宽一点,为此引入“软间隔(soft margin)”的概念。即允许少量样本不满足约束
yi(XTiW+b)≥1(3.1.1)(3.1.1)yi(XiTW+b)≥1
为了使不满足上述条件的样本点尽可能少,我们需要在优化的目标函数(2.2.2)(2.2.2)里面新增一个对这些点的惩罚项。最常用的是hinge损失:
lhinge(z)=max(0,1−z)(3.1.2)(3.1.2)lhinge(z)=max(0,1−z)
即若样本点满足约束条件损失就是0, 否则损失就是 1−z1−z ,则优化目标 (2.2.2)(2.2.2) 变成
minW,b12||W||2+C∑i=1nmax(0,1−yi(XTiW+b))(3.1.3)(3.1.3)minW,b12||W||2+C∑i=1nmax(0,1−yi(XiTW+b))
其中 C>0C>0 称为惩罚参数,CC 越小时对误分类惩罚越小,越大时对误分类惩罚越大,当 CC 取正无穷时就变成了硬间隔优化。实际应用时我们要合理选取 CC,CC 越小越容易欠拟合,CC 越大越容易过拟合。
如果我们引入“松弛变量” ξi≥0ξi≥0, 那么式 (3.1.3)(3.1.3) 可重写成
minW,b,ξ12||W||2+C∑i=1nξis.t. yi(XTiW+b)≥1−ξiξi≥0,i=1,2,…n.(3.1.4)(3.1.4)minW,b,ξ12||W||2+C∑i=1nξis.t. yi(XiTW+b)≥1−ξiξi≥0,i=1,2,…n.
上式所述问题即软间隔支持向量机。
3.2 对偶问题
式 (3.1.4)(3.1.4) 表示的软间隔支持向量机依然是一个凸二次规划问题,和硬间隔支持向量机类似,我们可以通过拉格朗日乘子法将其转换为对偶问题进行求解。
式 (3.1.4)(3.1.4) 对应的拉格朗日函数为
L(W,b,ξ,α,β)=12||W||2+C∑i=1nξi−∑i=1nαi[yi(XTiW+b)−1+ξi]−∑i=1nβiξi(3.2.1)(3.2.1)L(W,b,ξ,α,β)=12||W||2+C∑i=1nξi−∑i=1nαi[yi(XiTW+b)−1+ξi]−∑i=1nβiξi
类似2.4节,为了求得对偶问题的解,我们需要先求得L(W,b,ξ,α,β)L(W,b,ξ,α,β)对 WW、bb 和 ξξ 的极小再求对 αα 和 ββ 的极大。
以下两步和2.4节几乎完全一样,除了最后对 αα 的约束条件略有不同。
(1) 求minW,b,ξL(W,b,ξ,α,β)minW,b,ξL(W,b,ξ,α,β):
将 L(W,b,ξ,α,β)L(W,b,ξ,α,β) 分别对 WW、bb 和 ξξ 求偏导并令为0可得
W=∑i=1nαiyiXi(3.2.2)(3.2.2)W=∑i=1nαiyiXi
∑i=1nαiyi=0(3.2.3)(3.2.3)∑i=1nαiyi=0
C=αi+βi(3.2.4)(3.2.4)C=αi+βi
将上面三个式子代入式 (3.2.1)(3.2.1) 并进行类似式 (2.4.8)(2.4.8) 的推导即得
minW,b,ξL(W,b,ξ,α,β)=−12∑i=1n∑j=1nαiαjyiyjXTiXj + ∑i=1nαi(3.2.5)(3.2.5)minW,b,ξL(W,b,ξ,α,β)=−12∑i=1n∑j=1nαiαjyiyjXiTXj + ∑i=1nαi
注意其中的 ββ 被消去了。
(2) 求minW,b,ξL(W,b,ξ,α,β)minW,b,ξL(W,b,ξ,α,β)对 αα 的极大:
式(3.2.5)(3.2.5)对αα求极大,也等价于式(3.2.5)(3.2.5)取负数后对αα求极小,即
minα12∑i=1n∑j=1nαiαjyiyjXTiXj − ∑i=1nαi(3.2.6)(3.2.6)minα12∑i=1n∑j=1nαiαjyiyjXiTXj − ∑i=1nαi
同时满足约束条件:
∑i=1nαiyi=00≤αi≤C,i=1,2,…,n.(3.2.7)(3.2.7)∑i=1nαiyi=00≤αi≤C,i=1,2,…,n.
至此,我们得到了原始最优化问题(3.1.4)(3.1.4)和对偶最优化问题(3.2.6)(3.2.6)、(3.2.7)(3.2.7)。
类似2.4节地,假设我们现在通过通用的二次规划求解方法或者SMO算法求得了(3.2.6)(3.2.6)、(3.2.7)(3.2.7)的最优解 α^α^,则根据式(3.2.2)(3.2.2)可求得最优W^W^:
W^=∑i=1nα^iyiXi(3.2.8)(3.2.8)W^=∑i=1nα^iyiXi
再根据KKT条件,即
⎧⎩⎨⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪乘子非负:αi≥0,βi≥0(i=1,2,…n.下同)约束条件:yi(XTiW+b)−1≥ξi互补条件:αi[yi(XTiW+b)−1+ξi]=0,βiξi=0{乘子非负:αi≥0,βi≥0(i=1,2,…n.下同)约束条件:yi(XiTW+b)−1≥ξi互补条件:αi[yi(XiTW+b)−1+ξi]=0,βiξi=0
可求得整个软间隔SVM的解,即:
W^=∑i∈SVα^iyiXi(3.2.9)(3.2.9)W^=∑i∈SVα^iyiXi
b^=yj−∑i∈SVα^iyiXTjXi(3.2.10)(3.2.10)b^=yj−∑i∈SVα^iyiXjTXi
其中 jj 需满足 0<α^j<C0<α^j<C 。
对于任意样本 (Xi,yi)(Xi,yi) ,
- 若 αi=0αi=0,此样本点不是支持向量,该样本对模型没有任何的作用;
- 若 αi>0αi>0,此样本是一个支持向量。
若满足 αi>0αi>0 ,进一步地,
- 若 0<αi<C0<αi<C, 由式 (3.2.4)(3.2.4) 得 βi=0βi=0,即刚好 yi(XTiW+b)=1yi(XiTW+b)=1,样本恰好在最大间隔边界上;
- 若 αi=Cαi=C,有βi>0βi>0,此时若 βi<1βi<1则该样本落在最大间隔内部,若 βi>1βi>1 则该样本落在最大间隔内部即被错误分类。
如下图所示。
图3.2
因此,我们有与2.4节相同的结论,最优超平面只与支持向量有关而与非支持向量无关。
3.3 惩罚参数 CC
对于不同惩罚参数 CC,SVM结果如下图所示。
图3.3
再来看看我们的原始目标函数:
minW,b,ξ12||W||2+C∑i=1nξiminW,b,ξ12||W||2+C∑i=1nξi
对于更加一般化的问题,可将上述式子抽象成:
minfΩ(f)+C∑i=1nl(f(xi),yi)(3.3.1)(3.3.1)minfΩ(f)+C∑i=1nl(f(xi),yi)
前一项可以理解为“结构风险(structural risk)”,用来描述所求模型的某些性质(SVM就是要求间隔最大);第二项称为“经验风险(empirical risk)”,用来描述模型与训练数据的契合程度(即误差)。而参数 CC 就是用于对二者的折中,即我们一方面要求模型要满足某种性质另一方面又想使模型与训练数据很契合。
从正则化角度来讲, Ω(f)Ω(f) 称为正则化项,CC 称为惩罚参数,CC 越大即对误分类的惩罚越大(要求模型对训练模型更契合),这可能会存在过拟合;CC 越小即相对更加看重正则化项,此时可能存在欠拟合。
4. 非线性SVM——核技巧
前面介绍的都是线性问题,但是我们经常会遇到非线性的问题(例如异或问题),此时就需要用到核技巧(kernel trick)将线性支持向量机推广到非线性支持向量机。需要注意的是,不仅仅是SVM,很多线性模型都可以用核技巧推广到非线性模型,例如核线性判别分析(KLDA)。
4.1 核函数
如下图所示,核技巧的基本思路分为两步:
- 使用一个变换将原空间的数据映射到新空间(例如更高维甚至无穷维的空间);
- 然后在新空间里用线性方法从训练数据中学习得到模型。
图4.1
怎样映射到特征空间?
先来看看核函数的定义,
设 XX 是输入空间(欧式空间RnRn的子集或离散集合),又设 HH 是特征空间(希尔伯特空间),如果存在一个 XX 到 HH 的映射
ϕ(x):X→Hϕ(x):X→H
使得对所有 x,z∈Xx,z∈X,函数 K(x,z)K(x,z) 满足条件K(x,z)=ϕ(x)⋅ϕ(z)K(x,z)=ϕ(x)⋅ϕ(z)
则称 K(x,z)K(x,z) 为核函数, ϕ(x)ϕ(x) 为映射函数,式中 ϕ(x)⋅ϕ(z)ϕ(x)⋅ϕ(z) 为 ϕ(x)ϕ(x) 和 ϕ(z)ϕ(z) 的內积。
通常,直接计算 K(x,z)K(x,z) 比较容易而通过 ϕ(x)ϕ(x) 和 ϕ(z)ϕ(z) 计算 K(x,z)K(x,z) 并不容易。而幸运的是,在线性支持向量机的对偶问题中,无论是目标函数还是决策函数都只涉及到输入样本与样本之间的內积,因此我们不需要显式地定义映射 ϕ(x)ϕ(x) 是什么而只需事先定义核函数 K(x,z)K(x,z) 即可。也就是说,在核函数 K(x,z)K(x,z) 给定的情况下,可以利用解线性问题的方法求解非线性问题的支持向量机,此过程是隐式地在特征空间中进行的。
4.2 正定核
由上面的介绍可知,我们只需要定义核函数就可以了。但是如何不通过映射 ϕ(x)ϕ(x) 判断给定的一个函数 K(x,z)K(x,z)是不是核函数呢?或者说,K(x,z)K(x,z) 需要满足什么条件才是一个核函数。
通常所说的核函数就是正定核函数,下面不加证明的给出正定核的充要条件,具体证明略显复杂,有兴趣的可以参考《统计学习方法》。
设 X⊂RnX⊂Rn ,K(x,z)K(x,z) 是定义在 X×XX×X 上的对称函数,如果对任意的 xi∈X,i=1,2,…,mxi∈X,i=1,2,…,m,K(x,z)K(x,z) 对应的Gram矩阵
K=[K(xi,xj)]m×mK=[K(xi,xj)]m×m
是半正定矩阵,则 K(x,z)K(x,z) 是正定核。
虽然有了上述定义,但是实际应用时验证 K(x,z)K(x,z) 是否是正定核依然不容易,因此在实际问题中一般使用已有的核函数,下面给出一些常用的核函数。
-
多项式核函数(polynomial kernel function)
K(x,z)=(x⋅z+1)p(4.2.1)(4.2.1)K(x,z)=(x⋅z+1)p
-
高斯核函数(Guassian kernel function)
K(x,z)=exp(−||x−z||22σ2)(4.2.2)(4.2.2)K(x,z)=exp(−||x−z||22σ2)
4.3 非线性支持向量机
如前4.1、4.2所述,利用核技巧可以很简单地把线性支持向量机扩展到非线性支持向量机,只需将线性支持向量机中的內积换成核函数即可。下面简述非线性支持向量机学习算法。
-
首先选取适当的核函数 K(x,z)K(x,z) 和适当的参数 CC,构造最优化问题
minα12∑i=1n∑j=1nαiαjyiyjK(Xi,Xj) − ∑i=1nαis.t.∑i=1nαiyi=00≤αi≤C,i=1,2,…,n.(4.3.1)(4.3.1)minα12∑i=1n∑j=1nαiαjyiyjK(Xi,Xj) − ∑i=1nαis.t.∑i=1nαiyi=00≤αi≤C,i=1,2,…,n.
再利用现成的二次规划问题求解算法或者SMO算法求得最优解 α^α^ 。 -
选择 α^α^ 的一个满足 0<α^j<C0<α^j<C 的分量 α^jα^j ,计算
b^=yj−∑i∈SVα^iyiK(Xj,Xi)(4.3.2)(4.3.2)b^=yj−∑i∈SVα^iyiK(Xj,Xi)
-
构造决策函数:
f(x)=sign(∑i∈SVα^iyiK(Xj,Xi)+b^)(4.3.3)(4.3.3)f(x)=sign(∑i∈SVα^iyiK(Xj,Xi)+b^)
5. 总结
5.1 SVM优缺点
任何算法都有其优缺点,支持向量机也不例外。
支持向量机的优点是:
- 由于SVM是一个凸优化问题,所以求得的解一定是全局最优而不是局部最优。
- 不仅适用于线性线性问题还适用于非线性问题(用核技巧)。
- 拥有高维样本空间的数据也能用SVM,这是因为数据集的复杂度只取决于支持向量而不是数据集的维度,这在某种意义上避免了“维数灾难”。
- 理论基础比较完善(例如神经网络就更像一个黑盒子)。
支持向量机的缺点是:
- 二次规划问题求解将涉及m阶矩阵的计算(m为样本的个数), 因此SVM不适用于超大数据集。(SMO算法可以缓解这个问题)
- 只适用于二分类问题。(SVM的推广SVR也适用于回归问题;可以通过多个SVM的组合来解决多分类问题)
5.2 SVM调参经验
SVM很适合做分类任务,但是如果刚开始接触SVM而不知道如何进行合理的参数选择的话可能得不到满意的结果。下面简介运用SVM的基本步骤,或者说调参经验。
主要参考:
Hsu C W, Chang C C, Lin C J. A practical guide to support vector classification[J]. 2003.
- 将原始数据转换为SVM算法期待的格式;
- 将数据进行scaling(很重要);
- 一般考虑用高斯核RBF(如果特征维度太高,建议直接用线性SVM);
- 交叉验证寻找最优的RBF的参数以及参数 CC ;
- 用上面找到的最优参数在整个训练集上训练;