机器学习复习：支持向量机（SVM）之二：线性不可分向量机

机器学习复习：支持向量机（SVM）之一：线性可分向量机
机器学习复习：支持向量机（SVM）之二：线性不可分向量机

一、硬间隔SVM和软间隔SVM
  前面简单介绍了以下线性可分SVM，本质上前面介绍的SVM属于硬间隔SVM，因为我们已经知道存在一个超平面 $S:\omega \cdot x+b=0$ 可以将所有的样本数据完美分开。但是实际上，并不是所有的样本数据都这么完美。有时候存在这种情况，给定一个样本数据集，绝大多数的样本都是线性可分的，但是存在某些特别的样本是线性不可分的。例如，如果我们将 $y=0$ 这条直线看作一个分类的超平面， $y=0$ 以上的点标记为1， $y=0$ 以下的点标记为-1，那么绝大多数的样本点符合这个规则，但是有一些点位于 $y=0$ 以上，却被标记为-1，有些点位于 $y=0$ 以下，却被标记为了1，这样就会导致整个数据集线性不可分。

  线性不可分的数据集不能使用硬间隔SVM，根本原因在于硬间隔SVM对于函数距离的使用较为严格，必须满足 $y_i(\omega \cdot x+b)\ge 1$ 这个条件。所以为了可以解决上面提到的问题，我们可以在硬间隔SVM的表达式上进行一些修改。既然条件比较严格，那么我们就对每一个条件进行一定的松弛操作，所以我们需要引入一个新的松弛变量 $\xi$ ，于是对每一个条件，我们有：
$$y_i(\omega \cdot x_i+b)\ge 1-{\xi }_i(1\le i\le m)\text{\hspace{1em}(1)}$$
  其中， ${\xi }_i\ge 0$。因为引入了多余的参数，所有我们需要对原目标函数和松弛变量之间的关系进行平衡，因此，目标函数更改为：
$$\frac{1}{2}||\omega |{|}^2+C\sum _{i=1}^m{\xi }_i(1\le i\le m)\text{\hspace{1em}(2)}$$
  这里的 $C$ 是一个人为设定的超参数，通常 $C>0$，或者可以叫做平衡系数，或者惩罚参数，主要是用来平衡两者之间的关系。当 $C$ 比较大的时候对那些对那些误分类的样本数据的惩罚作用就会比较大，反之则会比较小，所以最小化目标函数就包括了既要使得超平面 $S$ 到最近样本的 距离之和最小，又要使得误分类的数据样本的数目最小。

  所以，我们需要求解的问题就变成了如下形式（ $m$ 为样本数目）：
$$\begin{array}{rl}\underset{\omega ,b,\xi }{\min }& \frac{1}{2}||\omega |{|}^2+C\sum _{i=1}^m{\xi }_i\\ s.t.& y_i(\omega \cdot x_i+b)\ge 1-{\xi }_i(1\le i\le m)\\ & {\xi }_i\ge 0(1\le i\le m)\end{array}\text{\hspace{1em}(3)}$$

二、学习的对偶问题
  首先我们对式(3)进行一个小小的变化，如下：
$$\begin{array}{rl}\underset{\omega ,b,\xi }{\min }& \frac{1}{2}||\omega |{|}^2+C\sum _{i=1}^m{\xi }_i\\ s.t.& y_i(\omega \cdot x_i+b)-1+{\xi }_i\ge 0(1\le i\le m)\\ & {\xi }_i\ge 0(1\le i\le m)\end{array}\text{\hspace{1em}(4)}$$
  和前面一样，我们对上面的式子使用拉格朗日数乘法，因为上面由两个不同类型的约束条件，所以我们需要使用两个不同的变量进行控制。在这里我们使用 $\alpha =({\alpha }_1,{\alpha }_2,...,{\alpha }_i)$ 来对第一个条件进行约束，使用 $\mu =({\mu }_1,{\mu }_2,...,{\mu }_i)$ 来对第二个条件进行约束，于是我们的问题就转换成了求解下面的式子的极大极小值的问题（ $L(w,b,\xi ,\alpha ,\mu )$ 表示我们使用的拉格朗日函数）：
$$\underset{\alpha ,\mu }{\max }\underset{\omega ,b,\xi }{\min }L(w,b,\xi ,\alpha ,\mu )\text{\hspace{1em}(5)}$$
  我们的拉格朗日函数定义如下：
$$\begin{array}{rl}L(w,b,\xi ,\alpha ,\mu )& =\frac{1}{2}||\omega |{|}^2+C\sum _{i=1}^m{\xi }_i-\sum _{i=1}^m{\alpha }_i(y_i(\omega \cdot x_i+b)-1+{\xi }_i)-\sum _{i=1}^m{\mu }_i{\xi }_i\\ & =\frac{1}{2}||\omega |{|}^2+C\sum _{i=1}^m{\xi }_i-\sum _{i=1}^m{\alpha }_i{y}_i(\omega \cdot x_i+b)+\sum _{i=1}^m{\alpha }_i-\sum _{i=1}^m{\alpha }_i{\xi }_i-\sum _{i=1}^m{\mu }_i{\xi }_i\\ & =\frac{1}{2}||\omega |{|}^2+C\sum _{i=1}^m{\xi }_i-\omega \cdot \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-\sum _{i=1}^m{\alpha }_i{\xi }_i-\sum _{i=1}^m{\mu }_i{\xi }_i\end{array}\text{\hspace{1em}(6)}$$
  接着，我们需要先求出目标拉格朗日函数对于 $\omega$ ，$b$ ，$\xi$ 的极小值，所以我们需要对着三个变量分别求偏导，（注意，由于 $\xi$ 的每一个分量都是相同地位的，所以我们只需要对其中的一个分量求偏导即可。）如下：
$$\frac{\partial }{\partial \omega }L(w,b,\xi ,\alpha ,\mu )=\omega -\sum _{i=1}^m{\alpha }_i{y}_i{x}_i\text{\hspace{1em}(7)}$$
$$\frac{\partial }{\partial b}L(w,b,\xi ,\alpha ,\mu )=\sum _{i=1}^m{\alpha }_i{y}_i\text{\hspace{1em}(8)}$$
$$\frac{\partial }{\partial {\xi }_i}L(w,b,\xi ,\alpha ,\mu )=C-{\alpha }_i-{\mu }_i\text{\hspace{1em}(9)}$$
  我们分别令上面的三个偏导数为0，于是我们可以得到：
$$\omega =\sum _{i=1}^m{\alpha }_i{y}_i{x}_i\text{\hspace{1em}(10)}$$
$$\sum _{i=1}^m{\alpha }_i{y}_i=0\text{\hspace{1em}(11)}$$
$$C={\alpha }_i+{\mu }_i\text{\hspace{1em}(12)}$$
  我们将上面的结果代入到 $L(w,b,\xi ,\alpha ,\mu )$ 中：
$$\begin{array}{rl}L(w,b,\xi ,\alpha ,\mu )& =\frac{1}{2}||\omega |{|}^2+C\sum _{i=1}^m{\xi }_i-\omega \cdot \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-\sum _{i=1}^m{\alpha }_i{\xi }_i-\sum _{i=1}^m{\mu }_i{\xi }_i\\ & =\frac{1}{2}||\omega |{|}^2+C\sum _{i=1}^m{\xi }_i-\omega \cdot \omega +\sum _{i=1}^m{\alpha }_i-\sum _{i=1}^m{\alpha }_i{\xi }_i-\sum _{i=1}^m{\mu }_i{\xi }_i\\ & =\frac{1}{2}||\omega |{|}^2+C\sum _{i=1}^m{\xi }_i-\omega \cdot \omega +\sum _{i=1}^m{\alpha }_i-\sum _{i=1}^m({\alpha }_i{\xi }_i+{\mu }_i{\xi }_i)\\ & =\frac{1}{2}||\omega |{|}^2+C\sum _{i=1}^m{\xi }_i-\omega \cdot \omega +\sum _{i=1}^m{\alpha }_i-\sum _{i=1}^m{\xi }_i({\alpha }_i+{\mu }_i)\\ & =\frac{1}{2}||\omega |{|}^2+C\sum _{i=1}^m{\xi }_i-\omega \cdot \omega +\sum _{i=1}^m{\alpha }_i-C\sum _{i=1}^m{\xi }_i\\ & =-\frac{1}{2}||\omega |{|}^2+\sum _{i=1}^m{\alpha }_i\\ & =-\frac{1}{2}(\sum _{i=1}^m{\alpha }_i{y}_i{x}_i)(\sum _{j=1}^m{\alpha }_j{y}_j{x}_j)+\sum _{i=1}^m{\alpha }_i\\ & =-\frac{1}{2}\sum _{i=1}^m\sum _{j=1}^m{\alpha }_i{\alpha }_j{y}_i{y}_j(x_i\cdot x_j)+\sum _{i=1}^m{\alpha }_i\end{array}\text{\hspace{1em}(13)}$$
  所以问题就变成了求解下面式子的极大值：
$$\underset{\alpha ,\mu }{\max }-\frac{1}{2}\sum _{i=1}^m\sum _{j=1}^m{\alpha }_i{\alpha }_j{y}_i{y}_j(x_i\cdot x_j)+\sum _{i=1}^m{\alpha }_i$$
  反转一下正负号，就会有：
$$\underset{\alpha ,\mu }{\min }\frac{1}{2}\sum _{i=1}^m\sum _{j=1}^m{\alpha }_i{\alpha }_j{y}_i{y}_j(x_i\cdot x_j)-\sum _{i=1}^m{\alpha }_i\text{\hspace{1em}(14)}$$
  所以我们所需要的原问题的对偶问题可以表示如下：
$$\begin{array}{rl}\underset{\alpha ,\mu }{\min }& \frac{1}{2}\sum _{i=1}^m\sum _{j=1}^m{\alpha }_i{\alpha }_j{y}_i{y}_j(x_i\cdot x_j)-\sum _{i=1}^m{\alpha }_i\\ s.t.& \sum _{i=1}^m{\alpha }_i{y}_i=0\\ & C-{\alpha }_i-{\mu }_i=0\\ & {\alpha }_i\ge 0\\ & {\mu }_i\ge 0,1\le i\le m\end{array}\text{\hspace{1em}(15)}$$
  观察到上式中存在一个等式关系，于是，我们可以利用这个等式关系消去 ${\mu }_i$ 这个变量：
$$\begin{gathered} \because C-{\alpha }_i-{\mu }_i=0 \\ \therefore {\mu }_i=C-{\alpha }_i \\ \therefore {\mu }_i=C-{\alpha }_i\ge 0 \\ \therefore C\ge {\alpha }_i \end{gathered}$$
  于是，我们的问题可以有如下的表达：
$$\begin{array}{rl}\underset{\alpha }{\min }& \frac{1}{2}\sum _{i=1}^m\sum _{j=1}^m{\alpha }_i{\alpha }_j{y}_i{y}_j(x_i\cdot x_j)-\sum _{i=1}^m{\alpha }_i\\ s.t.& \sum _{i=1}^m{\alpha }_i{y}_i=0\\ & 0\le {\alpha }_i\le C,1\le i\le m\end{array}\text{\hspace{1em}(16)}$$
  以上就是我们需要求出的线性不可分SVM原问题的对偶形式。

求解 $\omega$ 和 $b$ 参数

  我们求出最后满足公式(16)所设定的条件的最优解的 $\alpha$ 的数值之后，接下来就需要计算超平面的 $\omega$ 和 $b$ 参数了。根据公式(10)，再结合我们计算出的 $\alpha$ 的最优解，就可以计算出超平面的 $\omega$ 参数，即：
$${\omega }^{\ast }=\sum _{i=1}^m{\alpha }_i{y}_i{x}_i\text{\hspace{1em}(17)}$$
  由于超平面方程只和某些特定的数据样本有关，和很多数据样本无关，这些无关的样本点对应的 $\alpha$ 参数为0。因此，我们可以选择某一个 $\alpha$ 的分量 ${\alpha }_i$，该分量满足 $0<{\alpha }_i<C$ ，就可以求得超平面的 $b$ 参数。因为对于这类样本数据，我们可以知道它的函数间隔是1（或者-1），于是有：
$$y_i({\omega }^{\ast }\cdot x_i+b^{\ast })=1\text{\hspace{1em}(18)}$$
  由上可以计算出 $b^{\ast }$ 的数值为：
$$\begin{gathered} b^{\ast }=y_i-{\omega }^{\ast }\cdot x_i \\ {b}^{\ast }=y_i-x_i\cdot \sum _{j=1}^m{\alpha }_j{y}_j{x}_j\text{\hspace{1em}(18)} \end{gathered}$$

三、支持向量和决策函数
  当我们根据上面的对偶问题求解出了需要的 ${\omega }^{\ast }$ 和 $b^{\ast }$ 参数后，我们就可以知道该分类模型的决策函数如下：
$$f(x)=sign({\omega }^{\ast }\cdot x+b^{\ast })\text{\hspace{1em}(19)}$$
  其中， $sign()$ 函数表示的是取数值的符号，如果数值是正数，返回1，如果数值是负数，返回-1。
  在线性不可分SVM的求解过程中，因为引入了松弛变量，所以支持向量的判定会较为复杂。会出现以下的情况：

1. 如果 ${\alpha }_i<C$，这时 ${\xi }_i=0$，我们可以知道此时该样本位于间隔的边界上，正好属于一个支持向量。
2. 如果 ${\alpha }_i=C$，这时我们知道该样本的数据已经偏离了正确的间隔边界，所以又存在以下的几种情况：
   1. 如果此时 $0<{\xi }_i<1$，那我们根据函数间隔的性质可以知道，虽然样本数据已经越过了间隔的边界，但是仍然属于分类正确的样本，此时样本数据位于超平面和间隔边界之间。
   2. 如果此时 ${\xi }_i=1$，那么我们根据函数间隔的性质，可以知道的是样本正好位于分类的超平面上。
   3. 如果此时 ${\xi }_i>1$，那么这时候，样本数据已经越过了分类的超平面，已经属于误分类的样本数据了，且 ${\xi }_i$ 的数值越大，距离超平面越远。

  以上的样本数据统称为支持向量。

  综上就是线性不可分SVM的全部求解过程。