支持向量机SVM原理及推导

本周内容总结

1.间隔与支持向量

分类学习的基本思想:基于训练集在样本空间找到一个划分超平面,将不同类别的样本划分开。

特点:寻找位于两类训练样本“正中间”的划分超平面

原因:该超平面对训练样本局部扰动的“容忍”性最好,所产生的分类结果最鲁棒,对未见示例的泛化能力最强。

在样本空间,划分超平面可通过方程来描述,其中为法向量,决定了超平面的方向。b为位移项,决定了超平面与原点的距离。

样本空间中任意点x到超平面的距离可写为:

假设超平面(w,b)(w,b)能将训练样本正确分类,即对于任意样本(x_{i},y_{i})\in D为训练样本集,

y_{i}=+1,则有

y_{i}=-1,则有

如下图所示,距离超平面最近的这几个训练样本点使上式等号成立,它们被称为“支持向量”,两个异类支持向量到超平面的距离之和称为“间隔”。

           

欲找到具有“最大间隔”的划分超平面,也就是要找到能满足中的约束参数wb,使得\gamma最大,即

显然,为了最大化间隔,仅需要最大化||w||^{-1},这等价于最小化||w||^{2}。于是,有

这就是支持向量机的基本型。

2.对偶问题

对于一般优化问题

为了得到最大间隔划分超平面,我们采用了拉格朗日乘子法得到其“对偶问题”,即对每条约束添加拉格朗日乘子\alpha_{i}\geq 0,则该问题的拉格朗日函数可写为

其中

的偏导为零可得

将求得的导数带入中,即可消去参数wb,并得到对偶问题

对偶问题的表达式

解出\alpha后,求出参数wb即可得到模型

从对偶问题中解出的\alpha _{i} 是拉格朗日函数的拉格朗日乘子,它恰对应着训练样本(x_{i},y_{i})

注意到支持向量机的基本型中有不等式约束,因此上述过程需满足KKT(Karush-Kuhn-Tucker)条件,即要求

于是,对于任意样本(x_{i},y_{i}),总有

\alpha_{i}=0,则该样本将不会出现在的求和中,也就不会对f(x)有任何影响;

\alpha_{i} >0,则必有,所对应的样本点位于最大间隔边界上,是一个支持向量。

这显示出支持向量机的一个重要性质:训练完成后,大部分的训练样本都不需要保留,最终模型仅与支持向量有关。

如何求解对偶问题?

这里采用的是SMO(Sequential Minimal Optimization)算法

SMO的基本思路是:先固定\alpha _{i}之外的所有参数,然后求\alpha _{i}的上的极值。

由于存在约束,若固定\alpha _{i}之外的其他变量,则\alpha _{i}可由其他变量导出。因此,SMO每次选择两个变量\alpha _{i}\alpha_{j},并固定其他参数。这样,在参数初始化后,不断执行如下两个步骤直至收敛:

1)选取一对需要更新对的变量\alpha _{i}\alpha_{j}

2)固定\alpha _{i}\alpha_{j}以外的参数,求解对偶问题获得更新后的\alpha _{i}\alpha_{j}

3.核函数

之前的讨论中,我们假设训练样本是线性可分的,即存在一个划分超平面能将训练样本正确分类。然而在现实任务中,原始样本空间也许并不存在一个能正确划分两个样本的超平面。

当样本在原始空间线性不可分时,可将样本从原始空间映射到一个更高维的特征空间中,使得样本在这个特征空间内线性可分。如上图所示,将二维空间映射为一个合适的三维空间,就能找到一个合适的划分超平面。同时,如果原始空间是有限维,即属性数有限,那么一定存在一个高维特征空间使样本可分。

\phi(x)表示映射后的特征向量,于是,在特征空间中划分超平面所对应的模型可表示为

同时对支持向量机的最基本形式有

其对偶问题是

而引入这样的映射后,所要求解的对偶问题的求解中,无需求解真正的映射函数,而只需要知道其核函数。

x_{i}x_{j}在特征空间的内积等于它们在原始空间中通过函数计算的结果。有了这样的函数,我们就不必直接去计算高维甚至无穷维特征空间中的内积,于是,对偶问题可重写为

求解后即可得到

这里的函数就是“核函数”(kernel function)。

若已知映射\phi(\cdot ),则可写出核函数。但在现实任务中我们通常不知道\phi(\cdot )是什么形式,那么,合适的核函数是否一点存在?

什么样的函数能做核函数?

只要一个对称函数所对应的核矩阵半正定,它就能作为核函数。

事实上,对于一个半正定核矩阵,总能找到一个与之对应的映射\phi(\cdot )

我们希望样本在特征空间内线性可分,因此特征空间的好坏对支持向量机的性能至关重要。然而在不知道特征映射的形式时,我们并不知道什么样的核函数是合适的,而核函数也就是隐式地定义了这个空间。于是,核函数的选择 成了支持向量机的最大变数。在实际处理过程的,选取不同的核函数进行尝试对比来选取核函数的。

此外,还可通过函数组合得到,例如

为核函数,则核函数的直积也是核函数;

为核函数,则对于任意正数,其线性组合也是核函数;

为核函数,则对于任意函数g(x)也是核函数。

4.软间隔与正则化

前面的讨论中,我们一直假定训练样本在样本空间或特征空间中是线性可分的,即存在一个超平面能将不同类的样本完全划分开。然而,在现实任务中往往很难确定合适的核函数使得训练样本在特征空间中线性可分;同时,即使找到哦啊了某个核函数使训练集在特征空间中线性可分,也很难断定这个貌似线性可分的结果不是由于过拟合所造成的。

缓解该问题的一个办法是允许支持向量机在一些样本上出错。为此,引入“软间隔”(soft margin)的概念。

因此若所有样本都必须划分正确,这称为“硬间隔”(hard margin),而软间隔是允许某些样本不满足约束

当然,在最大化间隔的同时,不满足约束的样本应尽可能少。于是,优化目标可写为

其中C>0是一个常数,是“0/1损失函数”

C取有限值时,允许一些样本不满足约束;当C为无穷大时,则所有样本均满足约束。

非凸、非连续,数学性质不太好,使得目标函数不易求解。于是,通常用其他一些函数来代替,称为“替代损失”(surrogate loss)。替代损失函数一般具有较好的数学性质,如他们通常是凸的连续函数且是的上界。

以下是三种常用的替代损失函数:

若采用higne损失,则目标函数变为

引入“松弛变量”(slack variable)\xi_{i}\geq 0,则目标函数重写为

这就是常用的“软间隔支持向量机”。

软间隔支持向量机中每个样本都有一个对应的松弛变量,用以表征该样本不满足约束的程度。

同理,我们可通过拉格朗日乘子法得到软间隔支持向量机的拉格朗日函数

其中\alpha_{i}\geq 0\mu_{i}\geq 0是拉格朗日乘子。

wb\xi_{i} 的偏导为零可得

其中

将上面求得的结果带入软间隔支持向量机的拉格朗日函数

所以

故其对偶问题为

该软间隔与硬间隔的队友问题对比可看出,两个唯一的差别就在对偶变量的约束不同,前者是,后者是\alpha _{i}\geq 0。同样可采用SMO算法求解该目标函数;在引入核函数也能得到与硬间隔同样的支持向量展式。

对软间隔支持向量机,KKT条件要求

对于,对于任意训练样本(x_{i},y_{i}),总有\alpha_{i}=0.

\alpha_{i}=0,则该样本不会对f(x)有任何影响

\alpha_{i}>0,则必有,即该样本是支持向量

\alpha_{i}<C,则有\mu_{i}>0,进而有\xi_{i}=0,即该样本恰在最大间隔边界上

\alpha_{i}=C,则有\mu_{i}=0,此时若\xi_{i}\leq 1,则该样本落在最大间隔内部,若\xi_{i}> 1,则该样本被错误分类。

由此可以看出,软间隔支持向量机的最终模型仅与支持向量有关,即通过采用hinge损失函数仍保持了稀疏性。

5.支持向量回归

对样本(x,y),传统回归模型通常直接基于模型输出f(x)与真实输出y之间的差别来计算损失,当且仅当f(x)y完全相同时,损失才为零。与此不同,支持向量回归(Support Vector Regression,简称SVR)假设我们能容忍f(x)y之间最多有\epsilon的偏差,即当f(x)y之前的差别绝对值大于\epsilon时才计算损失。如图6.6所示,这相当于以f(x)为中心,构建了一个宽度为2\epsilon的间隔带,若训练样本落入此间隔带,则认为是被预测正确的。

于是,SVR问题可形式化为

其中C为正则化常数,是图6.7所示的\epsilon-不敏感损失(\epsilon-insensitive loss)函数

引入松弛变量\xi_{i} 和\hat{\xi}_{i} ,则可将SVR重写为

类似地,通过引入拉格朗日乘子,由拉格朗日乘子罚可得到SVR的拉格朗日函数

带入上式,再令的偏导数为零可得

将求导结果带入拉格朗日函数中,即可得到SVR的对偶问题

上述过程中需满足KKT条件,即要求

的约束条件全部恒等变形为小于等于0的形式可得

又由

从KKT条件可以看出,当且仅当时,\alpha _{i}能取到非零值,当且仅当时,\hat{\alpha }_{i}能取到非零值。

换言之,仅当样本(x_{i},y_{i})不落入\epsilon-间隔带中,相应的\alpha _{i} 和\hat{\alpha }_{i} 才能取非零值。

此外,约束不能同时成立,因此\alpha _{i} 和\hat{\alpha }_{i}中至少有一个为零。

w带入中,则SVR的解形如

能使上式中的的样本即为SVR的支持向量,它们必落在\epsilon-间隔带之外。

参考南瓜书中推导https://datawhalechina.github.io/pumpkin-book/#/chapter6/chapter6

  • 4
    点赞
  • 42
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值