（十二）支持向量机（Support Vecor Machine）2

支持向量机（二）

  在支持向量机原理(一) 线性支持向量机中，我们对线性可分SVM的模型和损失函数优化做了总结。最后我们提到了有时候不能线性可分的原因是线性数据集里面多了少量的异常点，由于这些异常点导致了数据集不能线性可分，本篇就对线性支持向量机如何处理这些异常点的原理方法做一个总结。

1. 线性分类SVM面临的问题

  有时候本来数据的确是可分的，也就是说可以用 线性分类SVM的学习方法来求解，但是却因为混入了异常点，导致不能线性可分，比如下图，本来数据是可以按下面的实线来做超平面分离的，可以由于一个橙色和一个蓝色的异常点导致我们没法按照上一篇线性支持向量机中的方法来分类。

  另外一种情况没有这么糟糕到不可分，但是会严重影响我们模型的泛化预测效果，比如下图，本来如果我们不考虑异常点，SVM的超平面应该是下图中的红色线所示，但是由于有一个蓝色的异常点，导致我们学习到的超平面是下图中的粗虚线所示，这样会严重影响我们的分类模型预测效果。

  如何解决这些问题呢？SVM引入了软间隔最大化的方法来解决。

2. 线性分类SVM的软间隔最大化

  所谓的软间隔，是相对于硬间隔说的，我们可以认为上一篇线性分类SVM的学习方法属于硬间隔最大化。回顾下硬间隔最大化的条件：

min12||w||22 s.tyi(wTxi+b) ≥1(i=1,2,...m) $$min\frac{1}{2}||w|{|}_2^{2}s.t{y}_i(w^T{x}_i+b)\ge 1(i=1,2,...m)$$

  接着我们再看如何可以软间隔最大化呢？ SVM对训练集里面的每个样本(xi,yi)$(x_i,y_i)$引入了一个松弛变量ξi≥0${\xi }_i\ge 0$,使函数间隔加上松弛变量大于等于1，也就是说：

yi(w∙xi+b)≥1−ξi$$y_i(w\bullet x_i+b)\ge 1-{\xi }_i$$

  对比硬间隔最大化，可以看到我们对样本到超平面的函数距离的要求放松了，之前是一定要大于等于1，现在只需要加上一个大于等于 0 的松弛变量能大于等于1就可以了。当然，松弛变量不能白加，这是有成本的，每一个松弛变量ξi${\xi }_i$, 对应了一个代价ξi${\xi }_i$，这个就得到了我们的软间隔最大化的SVM学习条件如下：

min12||w||22+C∑i=1mξi$$min\frac{1}{2}||w|{|}_2^2+C\sum _{i=1}^m{\xi }_i$$

s.t. yi(wTxi+b)≥1−ξi(i=1,2,...m)$$s.t.y_i(w^T{x}_i+b)\ge 1-{\xi }_i(i=1,2,...m)$$

ξi≥0(i=1,2,...m)$${\xi }_i\ge 0(i=1,2,...m)$$

  这里 C>0$C>0$ 为惩罚参数，可以理解为我们一般回归和分类问题正则化时候的参数。C$C$越大，对误分类的惩罚越大，C$C$越小，对误分类的惩罚越小。也就是说，我们希望12||w||22$\frac{1}{2}||w|{|}_2^2$尽量小，误分类的点尽可能的少。C是协调两者关系的正则化惩罚系数。在实际应用中，需要调参来选择。
  这个目标函数的优化和上一篇的线性可分SVM的优化方式类似，我们下面就来看看怎么对线性分类SVM的软间隔最大化来进行学习优化。

3. 线性分类SVM的软间隔最大化目标函数的优化

  和线性可分SVM的优化方式类似，我们首先将软间隔最大化的约束问题用拉格朗日函数转化为无约束问题如下：

L(w,b,ξ,α,μ)=12||w||22+C∑i=1mξi−∑i=1mαi[yi(wTxi+b)−1+ξi]−∑i=1mμiξi$$L(w,b,\xi ,\alpha ,\mu )=\frac{1}{2}||w|{|}_2^2+C\sum _{i=1}^m{\xi }_i-\sum _{i=1}^m{\alpha }_i[y_i(w^T{x}_i+b)-1+{\xi }_i]-\sum _{i=1}^m{\mu }_i{\xi }_i$$

  其中 μi≥0,αi≥0${\mu }_i\ge 0,{\alpha }_i\ge 0$,均为拉格朗日系数。
  也就是说，我们现在要优化的目标函数是：

minw,b,ξmaxαi≥0,μi≥0,L(w,b,α,ξ,μ)$$\underset{w,b,\xi }{\underbrace{min}}\underset{{\alpha }_i\ge 0,{\mu }_i\ge 0,}{\underbrace{max}}L(w,b,\alpha ,\xi ,\mu )$$

  这个优化目标也满足KKT条件，也就是说，我们可以通过拉格朗日对偶将我们的优化问题转化为等价的对偶问题来求解如下：

maxαi≥0,μi≥0,minw,b,ξL(w,b,α,ξ,μ)$$\underset{{\alpha }_i\ge 0,{\mu }_i\ge 0,}{\underbrace{max}}\underset{w,b,\xi }{\underbrace{min}}L(w,b,\alpha ,\xi ,\mu )$$

  我们可以先求优化函数对于w,b,ξ$w,b,\xi$的极小值, 接着再求拉格朗日乘子α$\alpha$和 μ$\mu$的极大值。首先我们来求优化函数对于w,b,ξ$w,b,\xi$的极小值，这个可以通过求偏导数求得：

∂L∂w=0⇒w=∑i=1mαiyixi$$\frac{\mathrm{\partial }L}{\mathrm{\partial }w}=0\Rightarrow w=\sum _{i=1}^m{\alpha }_i{y}_i{x}_i$$

∂L∂b=0⇒∑i=1mαiyi=0$$\frac{\mathrm{\partial }L}{\mathrm{\partial }b}=0\Rightarrow \sum _{i=1}^m{\alpha }_i{y}_i=0$$

∂L∂ξ=0⇒C−αi−μi=0$$\frac{\mathrm{\partial }L}{\mathrm{\partial }\xi }=0\Rightarrow C-{\alpha }_i-{\mu }_i=0$$

  好了，我们可以利用上面的三个式子去消除w$w$和b$b$了。

L(w,b,ξ,α,μ)=12||w||22+C∑i=1mξi−∑i=1mαi[yi(wTxi+b)−1+ξi]−∑i=1mμiξi　=12||w||22−∑i=1mαi[yi(wTxi+b)−1+ξi]+∑i=1mαiξi=12||w||22−∑i=1mαi[yi(wTxi+b)−1]=12wTw−∑i=1mαiyiwTxi−∑i=1mαiyib+∑i=1mαi=12wT∑i=1mαiyixi−∑i=1mαiyiwTxi−∑i=1mαiyib+∑i=1mαi=12wT∑i=1mαiyixi−wT∑i=1mαiyixi−∑i=1mαiyib+∑i=1mαi=−12wT∑i=1mαiyixi−∑i=1mαiyib+∑i=1mαi=−12wT∑i=1mαiyixi−b∑i=1mαiyi+∑i=1mαi=−12(∑i=1mαiyixi)T(∑i=1mαiyixi)−b∑i=1mαiyi+∑i=1mαi=−12∑i=1mαiyixTi∑i=1mαiyixi−b∑i=1mαiyi+∑i=1mαi=−12∑i=1mαiyixTi∑i=1mαiyixi+∑i=1mαi=−12∑i=1,j=1mαiyixTiαjyjxj+∑i=1mαi=∑i=1mαi−12∑i=1,j=1mαiαjyiyjxTixj(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(13)(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(13)$$\begin{array}{}\text{(1)}& L(w,b,\xi ,\alpha ,\mu )& =\frac{1}{2}||w|{|}_2^2+C\sum _{i=1}^m{\xi }_i-\sum _{i=1}^m{\alpha }_i[y_i(w^T{x}_i+b)-1+{\xi }_i]-\sum _{i=1}^m{\mu }_i{\xi }_i & (1)\text{(2)}& & =\frac{1}{2}||w|{|}_2^2-\sum _{i=1}^m{\alpha }_i[y_i(w^T{x}_i+b)-1+{\xi }_i]+\sum _{i=1}^m{\alpha }_i{\xi }_i& (2)\text{(3)}& & =\frac{1}{2}||w|{|}_2^2-\sum _{i=1}^m{\alpha }_i[y_i(w^T{x}_i+b)-1]& (3)\text{(4)}& & =\frac{1}{2}{w}^{T}w-\sum _{i=1}^m{\alpha }_i{y}_i{w}^T{x}_i-\sum _{i=1}^m{\alpha }_i{y}_{i}b+\sum _{i=1}^m{\alpha }_i& (4)\text{(5)}& & =\frac{1}{2}{w}^T\sum _{i=1}^m{\alpha }_i{y}_i{x}_i-\sum _{i=1}^m{\alpha }_i{y}_i{w}^T{x}_i-\sum _{i=1}^m{\alpha }_i{y}_{i}b+\sum _{i=1}^m{\alpha }_i& (5)\text{(6)}& & =\frac{1}{2}{w}^T\sum _{i=1}^m{\alpha }_i{y}_i{x}_i-w^T\sum _{i=1}^m{\alpha }_i{y}_i{x}_i-\sum _{i=1}^m{\alpha }_i{y}_{i}b+\sum _{i=1}^m{\alpha }_i& (6)\text{(7)}& & =-\frac{1}{2}{w}^T\sum _{i=1}^m{\alpha }_i{y}_i{x}_i-\sum _{i=1}^m{\alpha }_i{y}_{i}b+\sum _{i=1}^m{\alpha }_i& (7)\text{(8)}& & =-\frac{1}{2}{w}^T\sum _{i=1}^m{\alpha }_i{y}_i{x}_i-b\sum _{i=1}^m{\alpha }_i{y}_i+\sum _{i=1}^m{\alpha }_i& (8)\text{(9)}& & =-\frac{1}{2}(\sum _{i=1}^m{\alpha }_i{y}_i{x}_i{)}^T(\sum _{i=1}^m{\alpha }_i{y}_i{x}_i)-b\sum _{i=1}^m{\alpha }_i{y}_i+\sum _{i=1}^m{\alpha }_i& (9)\text{(10)}& & =-\frac{1}{2}\sum _{i=1}^m{\alpha }_i{y}_i{x}_i^T\sum _{i=1}^m{\alpha }_i{y}_i{x}_i-b\sum _{i=1}^m{\alpha }_i{y}_i+\sum _{i=1}^m{\alpha }_i& (10)\text{(11)}& & =-\frac{1}{2}\sum _{i=1}^m{\alpha }_i{y}_i{x}_i^T\sum _{i=1}^m{\alpha }_i{y}_i{x}_i+\sum _{i=1}^m{\alpha }_i& (11)\text{(12)}& & =-\frac{1}{2}\sum _{i=1,j=1}^m{\alpha }_i{y}_i{x}_i^T{\alpha }_j{y}_j{x}_j+\sum _{i=1}^m{\alpha }_i& (12)\text{(13)}& & =\sum _{i=1}^m{\alpha }_i-\frac{1}{2}\sum _{i=1,j=1}^m{\alpha }_i{\alpha }_j{y}_i{y}_j{x}_i^T{x}_j& (13)\end{array}$$

  其中，(1)式到(2)式用到了C−αi−μi=0​$C-{\alpha }_i-{\mu }_i=0$, (2)式到(3)式合并了同类项，(3)式到(4)式用到了范数的定义||w||22=wTw​$||w|{|}_2^2=w^{T}w$, (4)式到(5)式用到了上面的w=∑i=1mαiyixi​$w=\sum _{i=1}^m{\alpha }_i{y}_i{x}_i$， (5)式到(6)式把和样本无关的wT​$w^T$提前，(6)式到(7)式合并了同类项，(7)式到(8)式把和样本无关的b​$b$提前，(8)式到(9)式继续用到w=∑i=1mαiyixi​$w=\sum _{i=1}^m{\alpha }_i{y}_i{x}_i$，（9）式到(10)式用到了向量的转置。由于常量的转置是其本身，所有只有向量xi​$x_i$被转置，（10）式到(11)式用到了上面的∑i=1mαiyi=0​$\sum _{i=1}^m{\alpha }_i{y}_i=0$，（11）式到(12)式使用了(a+b+c+…)(a+b+c+…)=aa+ab+ac+ba+bb+bc+…​$(a+b+c+\dots )(a+b+c+\dots )=aa+ab+ac+ba+bb+bc+\dots$的乘法运算法则，（12）式到(13)式仅仅是位置的调整。

  仔细观察可以发现，这个式子和我们上一篇线性可分SVM的一样。唯一不一样的是约束条件。现在我们看看我们的优化目标的数学形式：

maxα∑i=1mαi−12∑i=1,j=1mαiαjyiyjxTixj$$\underset{\alpha }{\underbrace{max}}\sum _{i=1}^m{\alpha }_i-\frac{1}{2}\sum _{i=1,j=1}^m{\alpha }_i{\alpha }_j{y}_i{y}_j{x}_i^T{x}_j$$

s.t.∑i=1mαiyi=0$$s.t.\sum _{i=1}^m{\alpha }_i{y}_i=0$$

C−αi−μi=0$$C-{\alpha }_i-{\mu }_i=0$$

αi≥0(i=1,2,...,m)$${\alpha }_i\ge 0(i=1,2,...,m)$$

μi≥0(i=1,2,...,m)$${\mu }_i\ge 0(i=1,2,...,m)$$

  对于C−αi−μi=0，αi≥0，μi≥0$C-{\alpha }_i-{\mu }_i=0，{\alpha }_i\ge 0，{\mu }_i\ge 0$这3个式子，我们可以消去μi${\mu }_i$，只留下αi${\alpha }_i$，也就是说0≤αi≤C$0\le {\alpha }_i\le C$。 同时将优化目标函数变号，求极小值，如下：

minα 12∑i=1,j=1mαiαjyiyjxTixj−∑i=1mαi$$\underset{\alpha }{\underbrace{min}}\frac{1}{2}\sum _{i=1,j=1}^m{\alpha }_i{\alpha }_j{y}_i{y}_j{x}_i^T{x}_j-\sum _{i=1}^m{\alpha }_i$$

s.t.∑i=1mαiyi=0$$s.t.\sum _{i=1}^m{\alpha }_i{y}_i=0$$

0≤αi≤C$$0\le {\alpha }_i\le C$$

  这就是软间隔最大化时的线性可分SVM的优化目标形式，和上一篇的硬间隔最大化的线性可分SVM相比，我们仅仅是多了一个约束条件0≤αi≤C$0\le {\alpha }_i\le C$。我们依然可以通过SMO算法来求上式极小化时对应的α$\alpha$向量就可以求出w和b$w和b$了。

4. 软间隔最大化时的支持向量

  在硬间隔最大化时，支持向量比较简单，就是满足yi(wTxi+b)−1=0$y_i(w^T{x}_i+b)-1=0$就可以了。根据KKT条件中的对偶互补条件α∗i(yi(wTxi+b)−1)=0${\alpha }_i^{\ast }(y_i(w^T{x}_i+b)-1)=0$，如果α∗i>0${\alpha }_i^{\ast }>0$则有yi(wTxi+b)=1$y_i(w^T{x}_i+b)=1$ 即点在支持向量上，否则如果α∗i=0${\alpha }_i^{\ast }=0$则有yi(wTxi+b)≥1$y_i(w^T{x}_i+b)\ge 1$，即样本在支持向量上或者已经被正确分类。
  在软间隔最大化时，则稍微复杂一些，因为我们对每个样本(xi,yi)$(x_i,y_i)$引入了松弛变量ξi${\xi }_i$。我们从下图来研究软间隔最大化时支持向量的情况，第i个点到对应类别支持向量的距离为ξi||w||2$\frac{{\xi }_i}{||w|{|}_2}$。根据软间隔最大化时KKT条件中的对偶互补条件α∗i(yi(wTxi+b)−1+ξ∗i)=0${\alpha }_i^{\ast }(y_i(w^T{x}_i+b)-1+{\xi }_i^{\ast })=0$我们有：

1. 如果α=0$\alpha =0$,那么yi(wTxi+b)−1≥0$y_i(w^T{x}_i+b)-1\ge 0$,即样本在支持向量上或者已经被正确分类。如图中所有远离支持向量的点。
2. 如果0<α<C$0<\alpha <C$,那么ξi=0,yi(wTxi+b)−1=0${\xi }_i=0,y_i(w^T{x}_i+b)-1=0$,即点在支持向量上。如图中在虚线支持向量上的点。
3. 如果α=C$\alpha =C$，说明这是一个可能比较异常的点，需要检查此时ξi${\xi }_i$
   a. 如果0≤ξi<1$0\le {\xi }_i<1$,那么点被正确分类，但是却在超平面和自己类别的支持向量之间。如图中的样本2和4.
   b. 如果ξi=1${\xi }_i=1$,那么点在分离超平面上，无法被正确分类。
   c. 如果ξi>1​${\xi }_i>1$,那么点在超平面的另一侧，也就是说，这个点不能被正常分类。如图中的样本1和3.

5. 软间隔最大化的线性可分SVM的算法过程

  这里我们对软间隔最大化时的线性可分SVM的算法过程做一个总结。输入是线性可分的m个样本(x1,y1),(x2,y2),...,(xm,ym),$(x_1,y_1),(x_2,y_2),...,(x_m,y_m),$,其中x为n维特征向量。y为二元输出，值为1，或者-1。输出是分离超平面的参数w∗和b∗$w^{\ast }和b^{\ast }$和分类决策函数。

算法过程如下：

1. 选择一个惩罚系数C>0$C>0$, 构造约束优化问题
   
   minα 12∑i=1,j=1mαiαjyiyjxTixj−∑i=1mαi$$\underset{\alpha }{\underbrace{min}}\frac{1}{2}\sum _{i=1,j=1}^m{\alpha }_i{\alpha }_j{y}_i{y}_j{x}_i^T{x}_j-\sum _{i=1}^m{\alpha }_i$$
   
   
   s.t.∑i=1mαiyi=0$$s.t.\sum _{i=1}^m{\alpha }_i{y}_i=0$$
   
   
   0≤αi≤C$$0\le {\alpha }_i\le C$$
2. 用SMO算法求出上式最小时对应的α$\alpha$向量的值α∗${\alpha }^{\ast }$向量.
3. 计算w∗=∑i=1mα∗iyixi$w^{\ast }=\sum _{i=1}^m{\alpha }_i^{\ast }{y}_i{x}_i$

4. 找出所有的S个支持向量,即满足0<αs<C$0<{\alpha }_s<C$对应的样本(xs,ys)$(x_s,y_s)$，通过ys(∑i=1SαiyixTixs+b)=1$y_s(\sum _{i=1}^S{\alpha }_i{y}_i{x}_i^T{x}_s+b)=1$，计算出每个支持向量(xx,ys)$(x_x,y_s)$对应的b∗s$b_s^{\ast }$,计算出这些b∗s=ys−∑i=1SαiyixTixs$b_s^{\ast }=y_s-\sum _{i=1}^S{\alpha }_i{y}_i{x}_i^T{x}_s$. 所有的b∗s$b_s^{\ast }$对应的平均值即为最终的b∗=1S∑i=1Sb∗s$b^{\ast }=\frac{1}{S}\sum _{i=1}^S{b}_s^{\ast }$

   • 这样最终的分类超平面为：w∗∙x+b∗=0$w^{\ast }\bullet x+b^{\ast }=0$，最终的分类决策函数为：f(x)=sign(w∗∙x+b∗)$f(x)=sign(w^{\ast }\bullet x+b^{\ast })$

6. 合页损失函数

  线性支持向量机还有另外一种解释如下：

minw,b[1−yi(w∙x+b)]++λ||w||22$$\underset{w,b}{\underbrace{min}}[1-y_i(w\bullet x+b){]}_{+}+\lambda ||w|{|}_2^2$$

  其中L(y(w∙x+b))=[1−yi(w∙x+b)]+$L(y(w\bullet x+b))=[1-y_i(w\bullet x+b){]}_{+}$称为合页损失函数(hinge loss function)，下标+表示为：

[z]+={z0z>0z≤0$$[z{]}_{+}=\{\begin{array}{ll}z& z>0\\ 0& z\le 0\end{array}$$

  也就是说，如果点被正确分类，且函数间隔大于1，损失是0，否则损失是1−y(w∙x+b)$1-y(w\bullet x+b)$,如下图中的绿线。我们在下图还可以看出其他各种模型损失和函数间隔的关系：对于0-1损失函数，如果正确分类，损失是0，误分类损失1， 如下图黑线，可见0-1损失函数是不可导的。对于感知机模型，感知机的损失函数是[−yi(w∙x+b)]+$[-y_i(w\bullet x+b){]}_{+}$，这样当样本被正确分类时，损失是0，误分类时，损失是−yi(w∙x+b)$-y_i(w\bullet x+b)$，如下图紫线。对于逻辑回归之类和最大熵模型对应的对数损失，损失函数是log[1+exp(−y(w∙x+b))]$log[1+exp(-y(w\bullet x+b))]$, 如下图红线所示。

  线性可分SVM通过软间隔最大化，可以解决线性数据集带有异常点时的分类处理，但是现实生活中的确有很多数据不是线性可分的，这些线性不可分的数据也不是去掉异常点就能处理这么简单。那么SVM怎么能处理中这样的情况呢？我们在下一篇就来讨论线性不可分SVM和核函数的原理。