机器学习笔记之支持向量机(一)模型构建思路

引言

从本节开始，将介绍支持向量机，主要包括模型思路构建和模型最优参数的求解过程。

支持向量机介绍

支持向量机(Support Vector Machine,SVM)是一种基于监督学习的对数据进行 二元分类 的广义线性分类器(Generalized Linear Classifier)。

从线性分类的角度观察，支持向量机是一个硬分类模型，在机器学习笔记之线性分类——感知机算法中介绍过，硬分类模型使用的激活函数是非连续的。下面对支持向量机的模型 进行表示。

基于分类任务的性质，线性分类本质上就是使用直线(超平面)对样本空间进行划分，从而实现各样本在对应样本子空间的分类效果。

定义该直线(超平面) 的表达式为${\mathcal{W}}^{T}x+b=0(x\in \left\{x^{(i)}{|}_{i=1,2,\cdots ,N}\right\})$。基于硬分类模型的性质，定义支持向量机的模型表达如下：
$$f(\mathcal{W})=sign({\mathcal{W}}^{T}x+b)(x\in \left\{x^{(i)}{|}_{i=1,2,\cdots ,N}\right\})$$
其中$sign$表示符号函数：
$$sign(a)=\{\begin{array}{l}1ifa\ge 0\\ -1else\end{array}$$
因此，支持向量机是一个纯粹的判别模型。(注意：这里要和线性分类中软分类的概率判别模型区分开)：

• 判别模型 在这里的意思主要指 模型映射结果与真实标签的特征空间相同，区别于软分类的概率结果；
  判别模型的定义和硬分类的定义相类似，‘模型预测过程’和‘概率’没有关联关系。
• 概率判别模型 指的是对后验概率$P(\mathcal{Y}\mid \mathcal{X})$直接进行求解，并直接对结果进行比较。两者不是一个概念。

可以看出，支持向量机的模型描述和感知机(Perceptron)、线性判别分析(Linear Discriminant Analysis,LDA)相同，就是一个 基于硬分类的线性模型。

支持向量机构建的朴素思想

从几何角度观察，线性分类的本质即样本点在被划分后的样本子空间中被分类正确。但分类正确的直线(超平面)并不唯一。具体图像实例如下(图源：机器学习(周志华著))：

从图中可以看出，上述所有直线均可以将样本正确划分，如何从众多直线中找出一条最优直线，使得样本划分效果最优？

上述问题中隐含了一个评判标准：满足什么标准的直线是划分效果优的直线，如果存在这么一个标准(策略)，将这个标准满足到极致——就可以得到划分效果最优的直线。

从机器学习的角度出发，在机器学习笔记——极大似然估计与最大后验概率估计中介绍到，样本集合 X = { x ( 1 ) , x ( 2 ) , ⋯   , x ( N ) } \mathcal X = \{x^{(1)},x^{(2)},\cdots,x^{(N)}\} X={x(1),x(2),⋯,x(N)}可以理解为基于概率模型$P(\mathcal{X};\theta )$生成样本$x^{(i)}$组成的集合。因此，给定概率模型$P(\mathcal{X};\theta )$条件下，可以 源源不断地生成样本——相反，想要估计概率模型$P(\mathcal{X};\theta )$，却只能基于有限的数据集合进行估计。

因此，在线性分类中，基于有限的数据集合，我们不仅要关注模型是否会将样本正确划分，还要关注模型自身的鲁棒性：鲁棒性差的模型可能在训练集上可能得到一个不错的结果，但如果对一个不在训练集内部的新数据进行预测，其结果可能并不会满意。
过拟合(Overfitting)同样也是模型鲁棒性差的一种体现。

言归正传，如何通过模型的某些性质衡量模型的鲁棒性呢？观察上图中的直线1和直线2：

相比于直线2，直线1对于 数据集合样本空间边界的样本点的距离更加靠近，这种现象会导致直线1对应的模型鲁棒性更差。解释：如果在样本空间边界附近新增新的样本点(图中的红色样本点)，该样本点和原始样本之间可能只存在极小的噪声差距，但最终却被直线1错误分类。由此看出：相比于直线2，直线1对于样本噪声的容忍度更低，泛化误差更大。

支持向量机 策略构建的的朴素思想是：找出一条直线，使得位于数据集合的各样本子空间的样本点到该直线的距离均最大即可。

我们将满足上述思想的策略 称作 最大间隔分类器，其对应模型称为硬间隔SVM模型(hard-margin SVM)。

硬间隔SVM的模型构建

场景描述

在模型构建之前，照例对场景进行部署：
数据集合$Data={\left\{(x^{(i)},y^{(i)})\right\}}_{i=1}^N$，其中任意$x^{(i)}(i=1,2,\cdots ,N)$是$p$维向量，对应$y^{(i)}$是一个标量，并针对二分类任务性质进行如下表示：
x ( i ) = { x 1 ( i ) , x 2 ( i ) , ⋯   , x p ( i ) } T y ( i ) ∈ { − 1 , 1 } x^{(i)} = \left\{x_1^{(i)},x_2^{(i)},\cdots,x_p^{(i)}\right\}^{T} \\ y^{(i)} \in \{-1,1\} x(i)={x1(i)​,x2(i)​,⋯,xp(i)​}Ty(i)∈{−1,1}

模型构建过程

基于最大间隔分类器朴素思想，对其进行解析：

• 首先要满足的条件：基于当前样本点，需要分类正确；
  $$\{\begin{array}{l}{\mathcal{W}}^T{x}^{(i)}+b>0\to y^{(i)}=1\\ {\mathcal{W}}^T{x}^{(i)}+b<0\to y^{(i)}=-1\end{array}$$
  对于上述分段条件，进行进一步化简：只要满足${\mathcal{W}}^T{x}^{(i)}+b$与对应$y^{(i)}$同号(同正、同负)，即可分类正确。数学符号表达如下：
  $$y^{(i)}({\mathcal{W}}^T{x}^{(i)}+b)>0(i=1,2,\cdots ,N)$$

• 在样本点分类正确的基础上，使各样本子空间的样本点到该直线的距离均最大。

  但实际上，样本点是数据集合给定的，即样本点在$p$维样本空间的位置固定不变。因此，如果在$p$维样本空间中任意产生一条直线，在该直线确定的条件下，那么空间中的 $N$个样本点到该直线的距离均固定不变。

因此，可以将上述思想转化为：从能够将样本点分类正确的直线中找出这样一条直线，该直线与$N$个样本点对应的$N$个距离中找出长度最小的距离，而基于该直线找出最小距离比其他直线的都要大，该直线就是我们要找的直线。

基于上述思想，对符号进行定义：

• 定义样本点$x^{(i)}$到直线${\mathcal{W}}^{T}x+b=0$的距离$\mathcal{D}$表示如下：
  点到直线的距离公式~
  $$\mathcal{D}(\mathcal{W},b,x^{(i)})=\frac{1}{||\mathcal{W}||}\left|{\mathcal{W}}^T{x}^{(i)}+b\right|$$
• 长度最小距离$\mathcal{M}(\mathcal{W},b)$表示如下：
  $$\mathcal{M}(\mathcal{W},b)=\underset{\mathcal{W},b,x^{(i)}}{\min }\mathcal{D}(\mathcal{W},b,x^{(i)})$$
• 至此，长度最小距离已经得到，最终通过调整$\mathcal{W},b$使得$\mathcal{M}(\mathcal{W},b)$达到最大：
  $$\max \mathcal{M}(\mathcal{W},b)$$

对上面步骤进行整理，可以得到如下表示：
$$\underset{\mathcal{W},b}{\max }\underset{x^{(i)}{|}_{i=1}^N}{\min }\frac{1}{||\mathcal{W}||}\left|{\mathcal{W}}^T{x}^{(i)}+b\right|$$

至此，我们可以将步骤2视为目标函数，目标1视为约束条件，表达结果如下：
$$\{\begin{array}{l}\underset{\mathcal{W},b}{\max }\underset{x^{(i)}{|}_{i=1}^N}{\min }\frac{1}{||\mathcal{W}||}\left|{\mathcal{W}}^T{x}^{(i)}+b\right|\\ s.t.y^{(i)}({\mathcal{W}}^T{x}^{(i)}+b)>0\end{array}$$

观察目标函数：

• 由于 y ( i ) ∈ { − 1 , 1 } y^{(i)} \in \{-1,1\} y(i)∈{−1,1}且$y^{(i)}({\mathcal{W}}^T{x}^{(i)}+b)$结果恒正。因此，将$|{\mathcal{W}}^T{x}^{(i)}+b|$进行如下变换：
  $$|{\mathcal{W}}^T{x}^{(i)}+b|=y^{(i)}({\mathcal{W}}^T{x}^{(i)}+b)$$
• 由于$\frac{1}{||\mathcal{W}||}$和$x^{(i)}$无关，因此可以将$\frac{1}{||\mathcal{W}||}$提前，最终表达为如下形式：
  $$\underset{\mathcal{W},b}{\max }\frac{1}{||\mathcal{W}||}\underset{x^{(i)}{|}_{i=1}^N}{\min }{y}^{(i)}\left({\mathcal{W}}^T{x}^{(i)}+b\right)$$

观察约束条件：
$y^{(i)}\left({\mathcal{W}}^T{x}^{(i)}+b\right)$是由$|{\mathcal{W}}^T{x}^{(i)}+b|$转化而来，原因是：

• y ( i ) ∈ { − 1 , 1 } y^{(i)} \in \{-1,1\} y(i)∈{−1,1}
• $y^{(i)}$与${\mathcal{W}}^T{x}^{(i)}+b$同号。

但在超平面${\mathcal{W}}^T\mathcal{X}+b=0$确定的条件下，$|{\mathcal{W}}^T{x}^{(i)}+b|$可以相对地表示： 样本点$x^{(i)}$距离超平面${\mathcal{W}}^T\mathcal{X}+b=0$远近程度的一种描述。
需要注意的点，$|{\mathcal{W}}^T{x}^{(i)}+b|$它不是距离，它只是一种对‘距离’趋势的描述。如果$|{\mathcal{W}}^T{x}^{(i)}+b|$较大，那么样本点$x^{(i)}$距离超平面距离较远，反之同理。

小插曲：函数间隔

抛开 最大间隔分类器的条件，$y^{(i)}\left({\mathcal{W}}^T{x}^{(i)}+b\right)$结果正负均有可能：

• $y^{(i)}\left({\mathcal{W}}^T{x}^{(i)}+b\right)>0$，样本$(x^{(i)},y^{(i)})$分类正确；
• $y^{(i)}\left({\mathcal{W}}^T{x}^{(i)}+b\right)<0$，样本$(x^{(i)},y^{(i)})$分类错误；

从$y^{(i)}\left({\mathcal{W}}^T{x}^{(i)}+b\right)$结果连续性的角度观察，它的 数值部分(忽略它的符号) 包含实际意义：

• 当$y^{(i)}\left({\mathcal{W}}^T{x}^{(i)}+b\right)$数值部分越大，意味着样本点$(x^{(i)},y^{(i)})$在超平面${\mathcal{W}}^T\mathcal{X}+b=0$分类结果的确信度 越高；
  • $y^{(i)}\left({\mathcal{W}}^T{x}^{(i)}+b\right)$大于0时，数值越大意味着不仅分类正确，并且样本点距离分类超平面足够远，使得该分类结果越有说服力；
  • $y^{(i)}\left({\mathcal{W}}^T{x}^{(i)}+b\right)$小于0时，数值越大意味着不仅分类错误，并且样本点距离分类超平面足够远，使得该分类结果越离谱；
• 相反，当$y^{(i)}\left({\mathcal{W}}^T{x}^{(i)}+b\right)$数值结果越小，意味着样本点$(x^{(i)},y^{(i)})$在超平面${\mathcal{W}}^T\mathcal{X}+b=0$分类结果的确信度 越低；

这个表示确信度的数值，我们称之为函数间隔(Functional Margin)，定义其符号为$\mathcal{H}$。
$${\mathcal{H}}^{(i)}=y^{(i)}\left({\mathcal{W}}^T{x}^{(i)}+b\right)$$

假设对超平面进行等比例缩放：
$${\mathcal{W}}^T\mathcal{X}+b=0\to k{\mathcal{W}}^T\mathcal{X}+kb=0$$
两边对$k$约分，超平面没有任何变化。
但它的函数间隔却进行了缩放：
$$y^{(i)}\left(k\cdot {\mathcal{W}}^T{x}^{(i)}+k\cdot b\right)=k\cdot {\mathcal{H}}^{(i)}$$
$k$的取值方式有无穷多种，因此缩放方式有无穷多种，对应的函数间隔同样有无穷多种，此时，函数间隔最大化就无意义了。

继续观察最大间隔分类器的约束条件，由于${\mathcal{H}}^{(i)}=y^{(i)}\left({\mathcal{W}}^T{x}^{(i)}+b\right)>0$，则有：
$$\exists \gamma >0\to \underset{(x^{(i)},y^{(i)}{)}_{i=1}^N}{\min }{y}^{(i)}({\mathcal{W}}^T{x}^{(i)}+b)=\underset{(x^{(i)},y^{(i)}{)}_{i=1}^N}{\min }{\mathcal{H}}^{(i)}=\gamma$$

由于超平面的等比例缩放对函数间隔产生的影响，因此$\gamma$可能是大于0的任意值。
为了消除等比例缩放的影响，我们需要对$\gamma$设定一个具体值，为了方便计算，设定$\gamma =1$。
在确定$\gamma$后意味着无论缩放成什么情况，都可以映射为最小值是1的那个函数间隔。

至此，对上述结论进行整理：

• $y^{(i)}({\mathcal{W}}^T{x}^{(i)}+b)>0$，为消除等比例缩放对函数间隔的影响，转化为：
  $$\underset{(x^{(i)},y^{(i)}{)}_{i=1}^N}{\min }{y}^{(i)}({\mathcal{W}}^T{x}^{(i)}+b)=1$$
  或者：
  $$y^{(i)}({\mathcal{W}}^T{x}^{(i)}+b)\ge 1$$
• 基于步骤1，将目标函数改写如下：
  $$\underset{\mathcal{W},b}{\max }\frac{1}{||\mathcal{W}||}\underset{x^{(i)}{|}_{i=1}^N}{\min }{y}^{(i)}\left({\mathcal{W}}^T{x}^{(i)}+b\right)=\underset{\mathcal{W},b}{\max }\frac{1}{||\mathcal{W}||}\cdot 1=\underset{\mathcal{W},b}{\max }\frac{1}{||\mathcal{W}||}$$
  消除范式格式、分数形式：
  $\frac{1}{2}$是用于简化计算添加的一个系数，由于和$\mathcal{W}$无关。因此不影响$\mathcal{W}$的最终取值。
这里同样省略了一个根号，因为根号在当前定义域中式‘单调递增函数’，可能最后的‘最小值结果’是不同的，但不影响其结果是‘最小值’这个信息。
  $$\underset{\mathcal{W},b}{\max }\frac{1}{||\mathcal{W}||}\Rightarrow \underset{\mathcal{W},b}{\min }||\mathcal{W}||\Rightarrow \underset{\mathcal{W},b}{\min }\sqrt{{\mathcal{W}}^T\mathcal{W}}\Rightarrow \underset{\mathcal{W},b}{\min }\frac{1}{2}{\mathcal{W}}^T\mathcal{W}$$

最终，化简后公式表示如下：
$$\{\begin{array}{l}\underset{\mathcal{W},b}{\min }\frac{1}{2}{\mathcal{W}}^T\mathcal{W}\\ s.t.y^{(i)}({\mathcal{W}}^T{x}^{(i)}+b)\ge 1(\forall i=1,2,\cdots ,N)\end{array}$$

通过观察，该公式本质上是凸优化问题，并包含$N$个约束条件。下一节将介绍如何基于上式求解最优模型参数。

相关参考：
机器学习-支持向量机1-硬间隔SVM-模型定义(最大间隔分类器)
函数间隔与几何间隔