学习笔记【机器学习重点与实战】——6 支持向量机原理

1 支持向量机

支持向量机(support vector machines，SVM)是一种二类分类模型。它的基本模型是定义在特征空间上的间隔最大的线性分类器，间隔最大使它有别于感知机；支持向量机还包括核技巧，这使它成为实质上的非线性分类器。支持向量机的学习策略就是间隔最大化，可形式化为一个求解凸二次规划(convex quadratic programming)的问题，也等价于正则化的合页损失函数的最小化问题。支持向量机的学习算法是求解凸二次规划的最优化算法。

支持向量机学习方法包含构建由简至繁的模型：线性可分支持向量机(linear support vector machine in linearly separable case)、线性支持向量机(linear support vector machine)及非线性支持向量机(non-linear vector machine)。简单模型是复杂模型的基础，也是复杂模型的特殊情况。当训练数据线性可分时，通过硬间隔最大化(hard margin maximization)，学习一个线性的分类器，即线性可分支持向量机，又称为硬间隔支持向量机；当训练数据近似线性可分时，通过软间隔最大化(soft margn maximization)，也学习一个线性的分类器，即线性支持向量机，又称为软间隔支持向量机；当训练数据线性不可分时，通过使用核技巧(kernel trick)及软间隔最大化，学习非线性支持向量机。

——《统计学习方法》P95$P_{95}$

优点：泛化错误率低，计算开销不大，结果易解释。

缺点：对参数调节和核函数的选择敏感，原始分类器不加修改仅适用于处理二分类问题。

适用数据类型：数值型和标称型数据。

——《机器学习实战》P89$P_{89}$

2 线性可分支持向量机

给定线型可分训练数据集，通过间隔最大化或等价地求解相应的凸二次规划问题学习得到的分离超平面为：

ωT⋅Φ(x)+b=0$${\omega }^T·\Phi (x)+b=0$$

分类决策函数为：

f(x)=sign(ωT⋅Φ(x)+b)$$f(x)=sign({\omega }^T·\Phi (x)+b)$$

该决策函数称为线型可分支持向量机。

Φ(x)$\Phi (x)$是某个确定的特征空间转换函数，作用是将x映射到（更高的）维度。最简单的为Φ(x)=x$\Phi (x)=x$。

2.1 原理及定义

对于上图所示的二分类问题，线性可分支持向量机学习的目标是在特征空间中找到一个分离超平面，然而存在无穷个分离超平面可将两类数据正确分开；线性可分支持向量机则是利用间隔最大化求最优分离超平面，解是唯一的，图中的红色直线即为求解得到的最优分离超平面，将实例分到不同的类。

线性可分支持向量机中的间隔最大化又称为硬间隔最大化（与下面小节的软间隔最大化相对应）。

间隔最大化的直观解释是：对训练数据集找到几何间隔最大的超平面意味着以充分大的确信度对训练数据进行分类。也就是说，不仅将正负实例点分开，而且对最难分的实例点（离超平面最近的点）也有足够大的确信度将它们分开。这样的超平面应该对未知的新实例有很好的分类预测能力。

支持向量（support vector）：训练数据集的样本点中与分离超平面距离最近的样本点的实例。上图中两条虚线上的三个点即为支持向量。约束条件为：

yi(ωT⋅Φ(x)i+b)−1=0$$y_i({\omega }^T·\Phi (x{)}_i+b)-1=0$$

间隔边界：支持向量所在的两个超平面。

间隔（margin）：间隔边界之间的距离。其值为2||ω||$\frac{2}{||\omega ||}$。

在决定分离超平面时只有支持向量起作用，而其他实例点并不起作用、如果移动支持向量将改变所求的解；但是如果在间隔边界以外移动其他实例点，甚至去掉这些点，则解是不会改变的。由于支持向量在确定分离超平面中起着决定性作用，所以将这种分类模型称为支持向量机。支持向量的个数一般很少．所以支持向量机由很少的“重要的”训练样本确定。

2.2 算法步骤

输入：线型可分训练集T={(x1,y1),(x2,y2),...,(xN,yN)}$T=\{(x_1,y_1),(x_2,y_2),...,(x_N,y_N)\}$，
其中xi∈X=Rn,yi∈Y={−1,+1},i=1,2,...,N;$x_i\in X=R^n,y_i\in Y=\{-1,+1\},i=1,2,...,N;$

输出：分离超平面和分类决策树

（1）构造并求解约束最优化问题

minαs.t.αi≥0,12∑i=1N∑j=1Nαiαjyiyj(xi⋅xj)−∑i=1Nαi∑i=1Nαiyi=0i=1,2,...,N(3)(4)(5)$$\begin{array}{}\text{(3)}& \underset{\alpha }{min}& \frac{1}{2}\sum _{i=1}^N\sum _{j=1}^N{\alpha }_i{\alpha }_j{y}_i{y}_j(x_i·x_j)-\sum _{i=1}^N{\alpha }_i\text{(4)}& s.t.& \sum _{i=1}^N{\alpha }_i{y}_i=0\text{(5)}& {\alpha }_i\ge 0,& i=1,2,...,N\end{array}$$

求得最优解α∗=(α∗1,α∗2,...,α∗N)T${\alpha }^{\ast }=({\alpha }_1^{\ast },{\alpha }_2^{\ast },...,{\alpha }_N^{\ast }{)}^T$

（2）计算

ω∗=∑i=1Nα∗iyixi$${\omega }^{\ast }=\sum _{i=1}^N{\alpha }_i^{\ast }{y}_i{x}_i$$

并选择α∗${\alpha }^{\ast }$的一个正分量α∗>0${\alpha }^{\ast }>0$，计算

b∗=yi−∑i=1Nα∗iyi(xi⋅xj)$$b^{\ast }=y_i-\sum _{i=1}^N{\alpha }_i^{\ast }{y}_i(x_i·x_j)$$

（3）求得分离超平面

ω∗⋅x+b∗=0$${\omega }^{\ast }·x+b^{\ast }=0$$

分类决策函数：

f(x)=sign(ω∗⋅x+b∗)$$f(x)=sign({\omega }^{\ast }·x+b^{\ast })$$

在线型可分支持向量机中，ω∗${\omega }^{\ast }$和b∗$b^{\ast }$只依赖于训练数据中对应于α∗>0${\alpha }^{\ast }>0$的样本点(xi,yi)$(x_i,y_i)$，而其它样本点对ω∗${\omega }^{\ast }$和b∗$b^{\ast }$没有影响。则训练数据中对应于α∗>0${\alpha }^{\ast }>0$的实例点xi∈Rn$x_i\in R^n$称为支持向量。

3 线性支持向量机

3.1 原理及定义

对于线性可分问题，上述线性可分支持向量机的学习（硬间隔最大化）算法是完美的。但是，训练数据集线性可分是理想的情形。在现实问题中，训练数据集往往是线性不可分的，即在样本中出现噪声或特异点。这就需要将硬间隔最大化，改为软间隔最大化，允许支持向量机在一些样本上出错。

对每个样本点(xi,yi)$(x_i,y_i)$引进松弛变量ζi≥0${\zeta }_i\ge 0$，使得函数间隔加上松弛变量≥1。约束条件变为：

yi(ω⋅xi+b)≥1−ζi$$y_i(\omega ·x_i+b)\ge 1-{\zeta }_i$$

同时，对每个松弛变量ζi${\zeta }_i$，支付一个代价ζi${\zeta }_i$。目标函数由原来的12||ω||2$\frac{1}{2}||\omega |{|}^2$变成：

12||ω||2+C∑i=1Nζi$$\frac{1}{2}||\omega |{|}^2+C\sum _{i=1}^N{\zeta }_i$$

这里，C>0称为惩罚参数，一般由应用问题决定。

C越大，训练精度会变高，有可能过拟合，过渡带宽度越来越小。

软间隔的支持向量或者在间隔边界上，或者在间隔边界与分离超平面之间，或者在分离超平面误分一侧。若α∗i<C${\alpha }_i^{\ast }<C$，则ζi=0${\zeta }_i=0$，支持向量xi$x_i$恰好落在间隔边界上；若α∗i=C${\alpha }_i^{\ast }=C$，0<ζi<1$0<{\zeta }_i<1$，则分类正确，xi$x_i$在间隔边界与分离超平面之间；若α∗i=C${\alpha }_i^{\ast }=C$，ζi=1${\zeta }_i=1$，则在分离超平面上；若α∗i=C${\alpha }_i^{\ast }=C$，ζi>1${\zeta }_i>1$，则xi$x_i$位于分离超平面误分一侧。

3.2 损失函数

然而，0/1损失函数非凸、非连续，数学性质不太好，使得ζi${\zeta }_i$不易直接求解。于是，人们通常用其他一些函数来代替 0/1损失函数， 称为”替代损失” (surrogate loss)。替代损失函数一般具有较好的数学性质，如它们通常是凸的连续函数且是 0/1损失函数的上界。下图给出了三种常用的替代损失函数：

hinge损失：lhinge(z)=max(0,1−z)指数损失（exponentialloss）：lexp(z)=exp(−z)对率损失（logisticloss）：llog(z)=log(1+exp(−z))(6)(7)(8)$$\begin{array}{}\text{(6)}& & hinge损失：l_{hinge}(z)=max(0,1-z)\text{(7)}& & 指数损失（exponentialloss）：l_{exp}(z)=exp(-z)\text{(8)}& & 对率损失（logisticloss）：l_{log}(z)=log(1+exp(-z))\end{array}$$

由图可知，”替代损失”不仅要分类正确，而且确信度足够高时损失才是0，也就是说，对学习有更高的要求。

3.3 算法步骤

输入：线型可分训练集T={(x1,y1),(x2,y2),...,(xN,yN)}$T=\{(x_1,y_1),(x_2,y_2),...,(x_N,y_N)\}$，
其中xi∈X=Rn,yi∈Y={−1,+1},i=1,2,...,N;$x_i\in X=R^n,y_i\in Y=\{-1,+1\},i=1,2,...,N;$

输出：分离超平面和分类决策树

（1）选择惩罚参数C>0，构造并求解凸二次规划问题

minαs.t.12∑i=1N∑j=1Nαiαjyiyj(xi⋅xj)−∑i=1Nαi∑i=1Nαiyi=00≤αi≤C,i=1,2,...,N(9)(10)(11)$$\begin{array}{}\text{(9)}& \underset{\alpha }{min}& \frac{1}{2}\sum _{i=1}^N\sum _{j=1}^N{\alpha }_i{\alpha }_j{y}_i{y}_j(x_i·x_j)-\sum _{i=1}^N{\alpha }_i\text{(10)}& s.t.& \sum _{i=1}^N{\alpha }_i{y}_i=0\text{(11)}& & 0\le {\alpha }_i\le C,i=1,2,...,N\end{array}$$

求得最优解α∗=(α∗1,α∗2,...,α∗N)T${\alpha }^{\ast }=({\alpha }_1^{\ast },{\alpha }_2^{\ast },...,{\alpha }_N^{\ast }{)}^T$

（2）计算

ω∗=∑i=1Nα∗iyixi$${\omega }^{\ast }=\sum _{i=1}^N{\alpha }_i^{\ast }{y}_i{x}_i$$

并选择α∗${\alpha }^{\ast }$的一个正分量C>α∗>0$C>{\alpha }^{\ast }>0$，计算

b∗=yi−∑i=1Nα∗iyi(xi⋅xj)$$b^{\ast }=y_i-\sum _{i=1}^N{\alpha }_i^{\ast }{y}_i(x_i·x_j)$$

（3）求得分离超平面

ω∗⋅x+b∗=0$${\omega }^{\ast }·x+b^{\ast }=0$$

分类决策函数：

f(x)=sign(ω∗⋅x+b∗)$$f(x)=sign({\omega }^{\ast }·x+b^{\ast })$$

4 非线性支持向量机

4.1 原理及定义

前面两个算法都是假设训练样本是线性可分的，即存在一个分离超平面能将训练样本正确分类。然而在现实任务中，原始样本空间内也许并不存在一个能正确划分两类样本的超平面。例如下图中的” 异或 ” 问题就不是线性可分的。

对这样的问题，可将样本从原始空间映射到一个更高维的特征空间，使得样本在这个特征空间内线性可分。例如在上图中，若将原始的二维空间映射到一个合适的三维空间 ，就能找到一个合适的分离超平面。幸运的是，如果原始空间是有限维，即属性数有限，那么一定存在一个高维特征空间使样本可分。

用线性分类方法求解非线性分类问题分为两步：首先使用一个变换将原空间的数据映射到新空间；然后在新空间里用线性分类学习方法从训练数据中学习分类模型。

4.2 核函数

核技巧的做法是，在学习与预测中只定义核函数K(x,z)$K(x,z)$，而不显式地定义映射函数。

即将上面优化目标函数中的(xi⋅xj)$(x_i·x_j)$替换为K(xi⋅xj)$K(x_i·x_j)$。

核技巧巧妙的利用线性分类学习方法和核函数解决非线性问题。在实际应用中，往往以各领域知识直接选择核函数，核函数的选择的有效应需要通过实验验证。

常用核函数如下表：

名称	表达式	参数
线性核	K(xi,xj)=xTixj$K(x_i,x_j)=x_i^T{x}_j$
多项式核	K(xi,xj)=(xTixj)d$K(x_i,x_j)=(x_i^T{x}_j{)}^d$	d≥1$d\ge 1$ 为多项式的次数
高斯核	K(xi,xj)=exp(−||xi−xj||22σ2)$K(x_i,x_j)=exp(-\frac{||x_i-x_j|{|}^2}{2{\sigma }^2})$	σ>0$\sigma >0$ 为高斯核的带宽 (width)
拉普拉斯核	K(xi,xj)=exp(−||xi−xj||σ)$K(x_i,x_j)=exp(-\frac{||x_i-x_j||}{\sigma })$	σ>0$\sigma >0$
Sigmoid 核	K(xi,xj)=tanh(βxTixj+θ)$K(x_i,x_j)=tanh(\beta x_i^T{x}_j+\theta )$	tanh 为双曲正切函数， β>0,θ<0$\beta >0,\theta <0$

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

• 若K1$K_1$和K2$K_2$为核函数，则对于任意正数γ1${\gamma }_1$、γ2${\gamma }_2$，其线性组合

γ1K1+γ2K2$$\gamma 1K_1+\gamma 2K_2$$

也是核函数；

• 若K1$K_1$和K2$K_2$为核函数，则而核函数的直积

K1⊙K2(x,z)=K1(x,z)K2(x,z)$$K_1\odot K_2(x,z)=K_1(x,z)K_2(x,z)$$

也是核函数；

• 若K1$K_1$为核函数，则对于任意函数g(x)$g(x)$

K(x,z)=g(x)K1(x,z)g(z)$$K(x,z)=g(x)K_1(x,z)g(z)$$

也是核函数。

4.3 算法步骤

输入：训练集T={(x1,y1),(x2,y2),...,(xN,yN)}$T=\{(x_1,y_1),(x_2,y_2),...,(x_N,y_N)\}$，其中
xi∈X=Rn,yi∈Y={−1,+1},i=1,2,...,N;$x_i\in X=R^n,y_i\in Y=\{-1,+1\},i=1,2,...,N;$

输出：分类决策树

（1）选择适当的核函数K(x,z)$K(x,z)$和适当的惩罚参数C>0，构造并求解最优化问题

minαs.t.12∑i=1N∑j=1NαiαjyiyjK(xi⋅xj)−∑i=1Nαi∑i=1Nαiyi=00≤αi≤C,i=1,2,...,N(12)(13)(14)$$\begin{array}{}\text{(12)}& \underset{\alpha }{min}& \frac{1}{2}\sum _{i=1}^N\sum _{j=1}^N{\alpha }_i{\alpha }_j{y}_i{y}_{j}K(x_i·x_j)-\sum _{i=1}^N{\alpha }_i\text{(13)}& s.t.& \sum _{i=1}^N{\alpha }_i{y}_i=0\text{(14)}& & 0\le {\alpha }_i\le C,i=1,2,...,N\end{array}$$

求得最优解α∗=(α∗1,α∗2,...,α∗N)T${\alpha }^{\ast }=({\alpha }_1^{\ast },{\alpha }_2^{\ast },...,{\alpha }_N^{\ast }{)}^T$

（2）选择α∗${\alpha }^{\ast }$的一个正分量C>α∗>0$C>{\alpha }^{\ast }>0$，计算

b∗=yi−∑i=1Nα∗iyiK(xi⋅xj)$$b^{\ast }=y_i-\sum _{i=1}^N{\alpha }_i^{\ast }{y}_{i}K(x_i·x_j)$$

（3）构造决策函数：

f(x)=sign(∑i=1Nα∗iyiK(x⋅xi)+b∗)$$f(x)=sign(\sum _{i=1}^N{\alpha }_i^{\ast }{y}_{i}K(x·x_i)+b^{\ast })$$

5 序列最小优化算法-SMO

5.1 原理及定义

支持向量机的学习问题可以形式化为求解凸二次规划问题。这样的凸二次规划问题具有全局最优解，并且有许多最优化算法可以用于这一问题的求解，但是当训练样本容量很大时，这些算法往往变得非常低效，以致无法使用。所以，如何高效地实现支持向量机学习就成为一个重要的问题。

1998年Platt提出序列最小最优化(sequential minimal optimization，SMO)算法。该算法是支持向量机学习的一种快速算法，其特点是不断地将原二次规划问题分解为只有两个变量的二次规划子问题，并对子问题进行解析求解，直到所有变量满足KKT条件为止。这样通过启发式的方法得到原二次规划问题的最优解。因为子问题有解析解，所以每次计算子问题都很快，虽然计算子问题次数很多，但在总体上还是高效的。

5.2 算法步骤

输入：训练集T={(x1,y1),(x2,y2),...,(xN,yN)}$T=\{(x_1,y_1),(x_2,y_2),...,(x_N,y_N)\}$，
其中xi∈X=Rn,yi∈Y={−1,+1},i=1,2,...,N$x_i\in X=R^n,y_i\in Y=\{-1,+1\},i=1,2,...,N$，精度 ε ；

输出：近似解α^$\hat{\alpha }$

（1）取初值α(0)=0${\alpha }^{(0)}=0$，令k=0$k=0$；

（2）选取优化变量α(k)1,α(k)2${\alpha }_1^{(k)},{\alpha }_2^{(k)}$，解析求解两个变量的最优化问题

minα1,α2W(α1,α2)=s.t.12K11α21+12K22α22+y1y2K12α1α2−(α1+α2)+y1α1∑i=3NyiαiKi1+y2α2∑i=3NyiαiKi2α1y1+α2y2=∑i=3Nαiyi=ζ0≤αi≤C,i=1,2(15)(16)(17)(18)$$\begin{array}{}\text{(15)}& \underset{{\alpha }_1,{\alpha }_2}{min}W({\alpha }_1,{\alpha }_2)=& \frac{1}{2}{K}_{11}{\alpha }_1^2+\frac{1}{2}{K}_{22}{\alpha }_2^2+y_1{y}_2{K}_{12}{\alpha }_1{\alpha }_2\text{(16)}& & -({\alpha }_1+{\alpha }_2)+y_1{\alpha }_1\sum _{i=3}^N{y}_i{\alpha }_i{K}_{i1}+y_2{\alpha }_2\sum _{i=3}^N{y}_i{\alpha }_i{K}_{i2}\text{(17)}& s.t.& {\alpha }_1{y}_1+{\alpha }_2{y}_2=\sum _{i=3}^N{\alpha }_i{y}_i=\zeta \text{(18)}& & 0\le {\alpha }_i\le C,i=1,2\end{array}$$

求得最优解α(k+1)1,α(k+1)2${\alpha }_1^{(k+1)},{\alpha }_2^{(k+1)}$，更新α$\alpha$为α(k+1)${\alpha }^{(k+1)}$；

（3）若在精度 ε 范围内满足停机条件

∑i=1Nαiyi=00≤αi≤C,i=1,2,...,Nyi⋅g(xi)=⎧⎩⎨≥1,{xi|αi=0}=1,{xi|0≤αi≤C}≤1,{xi|αi=C}(19)(20)(21)$$\begin{array}{}\text{(19)}& \sum _{i=1}^N{\alpha }_i{y}_i=0\text{(20)}& 0\le {\alpha }_i\le C,i=1,2,...,N\text{(21)}& y_i·g(x_i)=\{\begin{array}{l}\ge 1,\{x_i|{\alpha }_i=0\}\\ =1,\{x_i|0\le {\alpha }_i\le C\}\\ \le 1,\{x_i|{\alpha }_i=C\}\end{array}\end{array}$$

其中，

g(xi)=∑j=1NαjyjK(xj,xi)+b$$g(x_i)=\sum _{j=1}^N{\alpha }_j{y}_{j}K(x_j,x_i)+b$$

则转（4）；否则令k=k+1$k=k+1$，转（2）；

（4）取α^=α(k+1)$\hat{\alpha }={\alpha }^{(k+1)}$

6 支持向量回归

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

引入松弛变量ζ_i和\hatζ_i，SVR可形式化为最优化问题：

mins.t.12||ω||2+C∑i=1N(ζi+ζ^i)f(xi)−yi≤ε+ζi,yi−f(xi)≤ε+ζ^i,ζi≥0,ζ^i≥0,i=1,2,...,N(22)(23)(24)(25)$$\begin{array}{}\text{(22)}& min& \frac{1}{2}||\omega |{|}^2+C\sum _{i=1}^N({\zeta }_i+{\hat{\zeta }}_i)\text{(23)}& s.t.& f(x_i)-y_i\le \epsilon +{\zeta }_i,\text{(24)}& & y_i-f(x_i)\le \epsilon +{\hat{\zeta }}_i,\text{(25)}& & {\zeta }_i\ge 0,{\hat{\zeta }}_i\ge 0,i=1,2,...,N\end{array}$$

7 参考

1. 机器学习升级版视频 - 邹博
2. 《统计学习方法》第7章 支持向量机
3. 《机器学习实战》第6章 支持向量机
4. 《机器学习 - 周志华》第6章 支持向量机

---

===========文档信息============
学习笔记由博主整理编辑，供非商用学习交流用
如本文涉及侵权，请随时留言博主，必妥善处置
版权声明：非商用自由转载-保持署名-注明出处
署名(BY) ：dkjkls（dkj卡洛斯）
文章出处：http://blog.csdn.net/dkjkls