logo

SVM

wangzf / 2022-09-13


目录

SVM 模型概览

SVM 介绍

支持向量机是 90 年代中期发展起来的基于统计学习理论的一种有监督机器学习方法, 通过寻求结构化风险最小来提高学习器的泛化能力,实现经验风险和置信范围的最小化, 从而达到在统计样本量较少的情况下,也能获得良好的统计规律性

img

对于一个分类问题,给定样本集​,的目的是在样本空间中找到一个划分超平面,可以将两种类别的样本分开。 事实上可能存在好几个超平面可以将正负样本分开,对于上图中不同的划分超平面, 从直观上来看应该用“中间”那个作为划分超平面,因为它距离两个类别最近的点距离差不多相等, 对于样本的局部扰动“容忍性”最高

所谓支持向量机, 顾名思义, 分为两个部分了解:

SVM 优缺点

SVM 理论

硬间隔

如何来确定上文提到的距离两个类别最近的点的距离?

超平面

首先先来定义一下划分超平面:

wTx+b=0

其中:

点到超平面的距离

点到平面的距离解释: img

给定平面 Ax+By+Cz+D=0,求平面外一点到平面的距离 d: 任意取平面上一点 P(x,y,z),连接两点 P,Q, 过 P 做平面的法向量 n=(A,B,C), 可以看到点 Q 到平面的距离正好是向量 PQ 在法向量 n 上的投影长度。 所以有: d=||PQ||cosθ=||n||||n||||PQ||cosθ=PQn||n|| 由上面的分析可知,点到平面的距离公式为向量与法向量的点积除以法向量的模

将划分超平面记为 (w,b),则样本空间任意点 x 到超平面的距离为:

γ=|wTx+b|||w||

超平面分类

假设超平面可以将样本正确分类,我们可以得到:

{wTxi+b+1,yi=+1wTxi+b+1,yi=1

注意上式中的 +1,若是超平面 (w,b) 可以将样本正确分类且距离超平面最近的样本距离为 d, 则有 w=1dwb=1db,使得:

{dwTxi+db+d,yi=+1dwTxi+db+d,yi=1

两边约去 d 总能得到上面的式子

支持向量

img

距离超平面最近的点使得上式中的等号成立,它们被称为“支持向量”(Support vectors), 不同种类的支持向量到超平面的距离之和称为“间隔”(Margin),有:

Margin=2||w||

分类的目标就是要找到这样的 w 使得划分超平面具有最大间隔, 可以等价于最小化 12||w||2, 于是得到需要求得的目标函数和约束条件

minω,b12||w||2

s.t.yi(wiTxi+b)1,i=1,2,,m

拉格朗日函数

原始问题是一个凸二次规划问题,将这个问题一般化:

minxRnf(x)

s.t.{ci(x)0,i=1,2,,khj(x)=0,j=1,2,,l

对每条约束添加拉格朗日乘子,上式写成拉格朗日函数为:

L(x,α,β)=f(x)+i=1kαici+j=1lβjhj(x)

上式将对 x 有限制条件的 f(x) 的最优化问题转转化为对 xαβ 没有限制条件的极值问题, 其中 αiβi 是拉格朗日乘子,且有 αi0

对于拉格朗日函数,问题求解变为:

θp(x)=maxα,β,αi0L(x,α,β)=maxα,β,αi0(f(x)+i=1kαici(x)+j=1lβjhj(x))

我们来对这个函数深入解析:

  1. 如果有 x 不满足原始约束条件,既有 ci(x)0hj0。 由于有​ αi0,可取 αi=+βjhj(x)=+,此时 θp(x)+
  2. 如果全部​都满足原始条件约束,即有​ ci(x)0hj=0αi0, 此时有 j=1lβjhj(x)=0i=1kαici0, 要使 θp(x) 最大,那么就有 i=1kαici(x)=0, 此时 θp(x)=f(x)。注意:max 函数的参数是 αβ,此时把 x 看做常数, 所以 f(x) 也是常数,maxf(x) 自然就是 f(x) 本身了
  3. 所以原最小化问题​ minω,bf(x) 可以变为 minω,bθp(x), 然后可以推出原问题的目标函数为

minxmaxα,β,α0L(x,α,β)

至此,我们将 SVM 要求得最大间隔的凸二次规划问题转变成了拉格朗日函数求解问题

对偶问题

得到了上述的拉格朗日函数,就要考虑如何求解了

通常情况下我们求最小化问题比较方便,即求导为 0 即可,于是我们将上述拉格朗日函数转为求其对偶问题, 即将 minmax 求解顺序对调,变为 maxmin​。有一个显而易见的结论, 即对一个序列先找极小值,然后在所有极小值中找到一个最大值 a​;对一个序列先找极大值, 然后在所有极大值中找到一个最小值​ bab 是明显成立的。因此有:

maxα,β,αi0minxL(x,α,β)minxmaxα,β,αi0L(x,α,β)

可以看到对偶问题外层的优化目标参数是拉格朗日参数,然后通过求得的拉格朗日参数, 间接得到我们最终要求的超平面的参数​ x。而原始问题外层的优化目标参数直接是​ x, 无法直接去优化得到这个参数,因为还存在着其他未知的拉格朗日参数。 (注意:此处的参数​ x 而非​ w, 是因为最开始将 12||w||2 ​抽象成了​ f(x)

我们开始回到求解 SVM 优化目标的问题上,先最小化​ L(w,b,α)

L(w,b,α)=12||w||2+i=1mαi(1yi(wTxi+b))

此式中 wx ​与线性回归中一样,为 nub_features 维的向量, m 为样本个数,对每个样本有一个​ αi, 即​ α=(α1,α2,,αm)

L(w,b,α)wb 的偏导为 0,有:

L(w,b,α)w=w[12||w||2+i=1m(αiαiyiwTxi+αib)]=w+i=1m(0αiyixi+0)=wi=1mαiyixi=0

所以有:

w=i=1mαiyixi

同理对 b 求偏导,可得:

i=1mαiyi=0

w 代入拉格朗日函数 L(w,b,α)=12||w||2+i=1mαi(1yi(wTxi+b)),同时利用上式,可以得到:

L(w,b,α)=12i=1mj=1mαiαjyiyjxiTxj+i=1mαii=1mαiyi((j=1mαjyjxj)xi+b)=12i=1mj=1mαiαjyiyjxiTxj+i=1mj=1mαiαjyiyjxiTxjbi=1mαiyi+i=1mαi=12i=1mj=1mαiαjyiyjxiTxj+i=1mαi

然后再来求此时关于 α 的极大值,即 maxα,α0minω,bL(w,b,α),故有:

maxα12imj=1mαiαjyiyjxiTxj+i=1mαi

s.t.αiyi=0α0,i=1,2,,m

此式为关于 α 的极大值求解,当求出解 α 之后,求出 wb,有:

f(x)=wTx+b=i=1mαiyixiTx+b

此处开始思考求解拉格朗日乘子 α,有约束条件,实际为 KKT 条件:

{αi01yif(xi)0αi(1yif(xi))=0

此时,对于任意样本 (xi,yi)​, 总有​ αi=0yif(xi)=1。若​ αi=0, 则该样本将不会在 f(x)=i=1mαiyixiTx+b 中出现, 也就不会对 f(x) ​有任何影响;若 αi0​,则必有 yif(xi)=1​, 所对应的样本点位于最大间隔边界上,是一个支持向量。这显示出支持向量的一个重要性质: 训练完成后,大部分的训练样本都不需保留,最终模型仅与支持向量有关

SMO 算法

那上面的 ​α 是如何求解得到的呢? 这里我们采用 Platt 在 1998 年专门针对 SVM 的对偶问题求解设计的算法 SMO(Sequential Minimal Optimization)。 它的主要思想是每次挑选出两个变量 αi ​和​ αj,固定其他的 α ​参数, 由此问题就变成了类似于二次方程求最大值的问题 。其大致为重复一下步骤:

  1. 选出​ αi,α2,,αm 中“最不好”的两个参数 αiαj
  2. 固定这两个参数以外的其他​ α, 求解​ maxαminw,bL(w,b,α) 来更新 αiαj

选取的 αiαj ​只要有一个不满足 KKT 条件, 目标函数就会在迭代会减小。而 KKT 条件的违背程度越大,则更新变量后目标函数值减小的幅度越大

那我们应该如何选取这两个参数 αi ​和​ αj 呢?

  1. 先选出违反 KKT 条件最严重的样本,以其参数作为第一个参数
  2. 采用启发式的思考方式,选择与第一个样本间隔最大的样本,其参数作为第二个参数。 直观的解释是,与对两个相似的参数进行更新相比,两个参数有很大的差别, 更新会使目标函数变化更大

软间隔

以上的推导形式都是建立在样本数据线性可分的基础上,如果样本数据中线性不可分(你中有我,我中有你), 应该如何处理呢?这里我们引入软间隔(Soft Margin),意味着,允许支持向量机在一定程度上出错

img

引入常量 C>0 和 “0/1 损失函数” l0/1(z),其特性为当 z<0 时函数为 1,否则函数为 0:

l0/1(z)={1,z00,z0

允许样本在一定程度上出错,即有样本不满足约束:

yi(wTxi+b)1

我们的目标是在最大化间隔的同时,尽可能减小不满足约束的样本数量,此时的优化目标写为:

minω,b12||w||2+Ci=1ml0/1(yi(wTxi+b)1)

img

但是上式的 l0/1 损失函数非凸、非连续,数学性质不好,用其他一些函数来代替它:

lhinge(z)=max(0,1z)

lexp(z)=ez

llog(z)=log(1+ez)

以上四种损失函数图像如下所示:

img

若采用 Hinge 损失,原优化目标变为:

minω,b12||w||2+Ci=1mmax(0,1yi(wTxi)+b)

引入“松弛变量”(slack variables),可将上式重写为:

minw,b,ξi12||w||w+Ci=1mξi

s.t.yi(wTxi+b)1ξiξi0,i=1,2,,m

仔细思考松弛变量 xi 的作用,可以看到,原来的硬间隔要求:yi(wTxi+b)1, 所有的样本距离划分超平面距离都要大于 1||w||;而引入松弛变量 ξ 之后, 我们允许那些样本距划分超平面 1ξ||w|| 的样本点(即离群点)存在,其中 ξ0。 允许离群点存在即忽略它们,不因为个别离群点改变原有的可以将绝大多数样本分离的划分超平面。可以看出,当松弛变量​ ξ 为 0 时即为硬间隔,不允许离群点存在;ξ 越大,离群点偏离分类群越远。当 ξ 为1时, 离群点在划分超平面上,当 ​ξ 大于1时,离群点已经跑到另一个分类样本群去了,即被错误分类了

再来思考惩罚因子​ C 的作用,它是随松弛变量 ξ 一起引入的,代表了有多重视离群点。​ C 越大表示越重视,越不想丢掉离群点,C ​为无穷大时迫使后半式​ i=1mξi 为 0, 而​ ξi0 所以有​ ξi=0,即不允许有离群点。 在实际应用中,C ​是需要事先设定且需要调参的数

上面得到的软间隔优化目标也是个二次规划问题,我们同样转换为拉格朗日函数:

L(w,b,α,ξ,μ)=12||w||2+i=1mαi(1ξiyi(wTxi+b))i=1mμiξi

其中:αi0μi0 是拉格朗日乘子。同样地我们转换为对偶问题,将现在的:

minω,b,ξimaxα,μL(w,b,α,ξ,μ)

转换为:

maxα,μminω,b,ξiL(w,b,α,ξ,μ)

L(w,b,α,ξ,μ)wbξi 的偏导为 0,可得:

w=i=1mαiyixi

0=i=1mαiyi

C=αi+μi

将此结果代入 L(w,b,α,ξ,μ),就得到其最小值。然后其对偶问题转换为:

maxαi=1mαi12i=1mj=1mαiαjyiyjxiTxj

s.t.i=1mαiyi=0

0αiC,i=1,2,,m

对比硬间隔的对偶问题,软间隔的差别只在于对对偶变量的约束不同: 前者是​ αi0,后者是​ 0αiC

核函数

软间隔解决的实际上是允许模型忽略某些异常点来进行划分的问题,即在绝大部分样本都能被正确分类的情况下, 某几个样本无法被划分超平面分开,我们不认为应该考虑它来重新修改模型,而是认为给出的样本是被错误分类的

那如果大部分样本都线性不可分怎么办?完全的你中有我,我中有你。

此时引入核函数,简单来说就是将低维空间的样本映射到高维空间。 以“异或”问题为例,二维情况下无论如何都不能将正负样本线性分开, 但映射到三维后就可以用一个平面将样本分开

img

由 Cover 定理可证,同一份样本数据在越高维的空间中越有可能线性可分。 若设 d 维空间中 N 个点线性可分的概率为​ p(d,N),那么就有:

p(d,N)=2i=0mCN1i2N={i=1dCN1i2N1,N>d+11,Nd+1

即样本空间维度大于样本个数时,一定线性可分(证明略)

$\phi$(\mathbf{x}) 表示将 x 映射后的特征向量,于是在特征空间中的划分超平面为:

f(x)=wTϕ(x)+b

现在的问题就转化为了如何找到合适的 ϕ()x 映射到高维空间,使得样本线性可分。

与硬间隔 SVM 对偶问题求解相似,我们要优化的目标变为:

minω,b12||w||2

s.t.yi(wTϕ(xi)+b)1,i=1,2,,m

转化为对偶问题:

maxαi=1mαi12i=1mj=1mαiαjyiyjϕ(xiT)ϕ(xj)

s.t.i=1mαiyi=0

αj0,i=1,2,,m

上式的计算中涉及到了 ϕ(xi)Tϕ(xj), 这是样本 xixj 映射到特征空间之后的内积。 由于特征空间维度可能很高甚至是无穷维,计算困难,于是换了一条路来解决这个问题,设想一个函数:

κ(xi,xj)=ϕ(xi),ϕ(xj)=ϕ(xi)Tϕ(xj)

其中:

我们把 ​κ(xi,xj) 称为核函数,它的作用显而易见, 使我们不必直接去计算​ ϕ(xi),ϕ(xj), 甚至不需要知道 ϕ() ​长什么样,直接求出​ ϕ(xi)Tϕ(xj) 的结果。

将核函数代入后求解,得

f(x)=wTϕ(x)+b=i=1mαiyiϕ(xi)Tϕ(x)+b=i=1mαiyiκ(x,xi)+b

上式表明,模型的最优解可以通过训练样本的核函数展开,这一展开式称为“支持向量展式”。

国外的英文文献把Kernel(核函数)叫做Kernel Function或者是Kernel Trick,从Trick这种叫法可知核函数其实是一种运算技巧而已。

可以证明,只要一个对称函数所对应的核矩阵半正定,它就能作为核函数使用。而对于一个半正定核矩阵,总能找到一个与之对应的 ϕ() ​。事实上,核函数的选取依然是一个未决问题,核函数选择的不合适,意味着将样本映射到了一个不合适的特征空间,可能会导致样本依然无法线性可分。

下表列出了几种常用的核函数:

img

可以发现,线性核函数对应着样本空间本身,即不对样本做函数映射

其他内容

线性可分性

给定一个数据集

T={(x1,y1),(x2,y2),...,(xn,yn)}

其中:

如果存在某个超平面 S:

ωTx+b=0

能够将数据集的正实例(y=1)和负实例(y=0)完全正确地分到超平面的两侧, 即对所有的 yi=1 的实例 i,有 ωTx+b>0; 对所有的 yi=0 的实例 i,有 ωTx+b<0, 则称数据集 T 为线性可分数据集, 否则称为线性不可分

一个二分类线性分类器就是要在 Rn 特征空间中找到一个超平面 S,其方程可以表示为:

ωTx+b=0

这个超平面将特征空间划分为两个部分, 位于两部分的点分别被分为两类

函数间隔与几何间隔

一般而言, 一个数据点距离超平面(ωTx+b=0)的远近可以表示为分类预测的确信或准确程度

在超平面 ωTx+b 确定的情况下,|ωTx+b| 能够相对表示点 x 距离超平面的远近; 而 ωTx+b 的符号与类别标记 y 的符号是否一致表示分类是否正确, 可以用指标量 y(ωTx+b) 的正负性来判定或表示分类的正确性和确信度

函数间隔

Function Margin

函数间隔:

γ^=y(ωTx+b)=yf(x)

超平面 ωTx+b 关于训练数据 T 的函数间隔为超平面 ωTx+b 关于 T 中所有样本点 (xi,yi) 的函数间隔的最小值:

γ^=minγ^i,i=1,2,...,n

上面定义的函数间隔虽然可以表示分类预测的正确性和确信度,但在选择分类超平面时,只有函数间隔是不够的。 如果成比例的改变 ωb,比如将他们改变为 2ω2b, 虽然此时超平面没有改变, 但是函数间隔的值 yf(x) 却变成了原来的 4 倍

解决防范: 可以对法向量 ω 加一些约束条件,使其表面上看起来规范化

几何间隔

Geometric Margin

参考