支持向量机学习(1)——理论推导

支持向量机学习

前言

对于支持向量机，最早是在大学毕业设计时接触到的，那时已经大概把原理看了看，知道是利用最优的超平面将两类数据分开。但好像也仅限于此了，对于其中的核函数以及算法的实现什么的都没有过深的了解，最后也只是调用MATLAB里面以有的支持向量机的函数，然后把输入的特征输入进去，进行训练就可以实现对不同目标的分类了。然而时过两年，现在又要再次用到支持向量机，还是打算把其中的一些原理搞懂，今天不欠账，以后就不用还。

问题描述

也已经无数次的看到过SVM可以找出一个最优超平面将不同的数据分开 ，也无数次的见识过各种文章写道核函数可以使得SVM有可能处理无数维的数据。但是却还是似懂非懂，现在也还是有点似懂非懂，所以还是打算先从简单入手，不涉及核函数和线性不可分。就考虑一个实际的问题：如下图，在直角坐标系上存在两类点，请用支持向量机把这两类点分开。这里，数据的特征就是各个点的横纵坐标了，也就是一个二维的特征。

在直角坐标系下，我们很容易可以理解，就是只找到直线2的解析式，那么我们就可以顺利的将两类数据分开了，至于直线1和直线3，一眼可以看出没有直线2分的效果好。当然，我们肯定需要严格的理论推导来确定直线2，不可以靠眼睛看就断定直线2比直线1和3好。这里的判断方法就是希望有最大的分类间隔。什么是间隔呢？如图，也可以说是离分类的直线距离最近的点（其实这些点就是支持向量）到分类直线的距离。也就是在准确分类的基础上，改变改变直线的斜率什么的，使得这个间隔最大。当然，我们找直线2的这个过程不可以真用笔在纸上画画用数学知识直接把直线2的解析式算出来，肯定利用计算机来找啊，因为一旦数据量变得很大，计算变得很复杂人们很难直接算出来。而通过计算机找直线2的方法也有多种，后面介绍一种叫序列最小优化的算法（smo，Sequential Minimal Optimization）

问题转化

对于我们拿到的训练集$D=(x_1,y_1),(x_2,y_2)，(x_3,y_3)...$，可以包含N个样本点，由$x$和$y$组成，我们也可以把每一个样本点叫做一个实例。$x$作为输入变量，可以是一个数，也可以是一个向量。对于本文中的例子，$x$就是一个二维向量，也就是每个点的横纵坐标。$y$为输出，在二分类问题中，$y$只能取+1或-1（或者0或1）。
例如$x_1=(x_1^1,x_1^2)$,其中$x_1^1,x_1^2$分别对应横纵坐标。这一个点对应一个$y_1=+1(或-1)$类。我们的任务就是通过这些训练集的点来找出一个条直线，然后在测试集时或者有新的数据输入时能够正确的将每一个数据分类。
前面已经说过我们想找一个有最大分类间隔的超平面，我们定义这个超平面，在简单的二维情况下就是一条直线，设为：$0=w^{T}x+b$，我们要求$w$和$b$，使得这个直线到最近点的几何距离最大（可以看看几何距离函数距离的差别，以及训练集到超平面距离的定义）。其中几何距离就是我们高中经常用到的点到直线的距离，写成数学形式：
$$\underset{w,b}{\max }\underset{x_i}{\min }(\frac{|w^T{x}_i+b|}{||w||})\text{\hspace{1em}(*)}$$
稍微理解一下这个式子，其中$||w||$代表求$w$的二范数，也就是相当于常用的$\sqrt{a^2+b^2}$。这个式子中也出现了绝对值，$|w^{T}x+b|$，一般绝对值运算比较麻烦，想着把他去掉，我们知道，y只能是+1或者-1，而一般规定为正样本时，y为+1，此时$w^{T}x+b$也为正的。为负样本时，y为-1，此时$w^{T}x+b$也为负的。（一条直线的一边的点全为正，为正样本；另外一边的点全为负，为负样本。这样规定方便简化计算）也就是$|w^{T}x+b|=y(w^{T}x+b)$,另外$||w||$是$x_i,y_i$无关的，所以可以移出来。最后前面的$(\ast )$式可以改写为：
$$\underset{w,b}{\max }\frac{1}{||w||}\underset{x_i,y_i}{\min }(y_i(w^T{x}_i+b))\text{\hspace{1em}(**)}$$
还是稍微理解一下这个式子，这个式子的意思就是，首先通过尝试不同的$w,b$使得后面的式子能够取得最大值，而后面的式子是指尝试所有已知点，也就是我们的训练数据，找出离这个直线最近距离的点，然后求出来的这个最近的距离。合起来也就是前面所说的求一条直线$0=w^{T}x+b$，使得这个直线到最近点的几何距离最大，我们要求的就是$w,b$。

接着转化

假如支持向量为$(x_v,y_v)$,因为支持向量就是离超平面最近的点，所以对于(**)式，可以简化写成：
$$\underset{w,b}{\max }\frac{1}{||w||}(y_v(w^T{x}_v+b))\text{\hspace{1em}(***)}$$
对于一条直线x+y=0，我们将等式两边同时乘2得到2x+2y=0，很容易知道这两者并没有本质区别。对于$w^{T}x+b=0$也一样，两边同乘一个常数并没有改变这条直线。这里我们就进行这样的操作使得$y_v(w^T{x}_v+b)=1$（直线上的点带入$w^{T}x+b$才等于0，支持向量带入是不等于0的）这样式$(\ast \ast \ast )$就变成了：
$$\underset{w,b}{\max }\frac{1}{||w||}$$
也就是求$w,b$使得$\frac{1}{||w||}$最大。为了方便后面的计算，改写成求$\frac{1}{2}||w|{|}^2$最小。但是是在${\min }_{x_i,y_i}(y_i(w^T{x}_i+b))=1$的情况下（与$y_v(w^T{x}_v+b)=1$等价）求$||w||$最小 。所以现在的问题变成了：
$$\underset{w,b}{\min }\frac{1}{2}||w|{|}^2$$$$s.t.\underset{x_i,y_i}{\min }(y_i(w^T{x}_i+b))=1\text{\hspace{1em}(-)}$$
但是，一般最后约束条件都会写成：
$$s.t.y_i(w^T{x}_i+b)\ge 1$$

其实对这个约束条件是否和前面约束(-)的是否等价,我抱有疑问。我在想会不会出现所有点都离超平面的距离都比1大很多，而没有等于1的，这样会不会造成分类出现问题。最后和同学稍微讨论一下，感觉最后对$w,b$进行优化时，最先是假设一个点为支持向量，然后让这个点到这个超平面的距离为1来改变$w,b$，接着考验其余的点，来一点点优化的。所以开始应该会有一个点离超平面的距离为1，或者说接近1。
对上面的问题稍微有了一点不同的理解，就是所谓约束就是肯定可以取到约束情况下所有点的啊。就比如限制$x^2+y^2\le 1$，这就是一个圆的区域，我们这个范围上的点都是可以取得到的，所以圆上的点当然也可以取到，所以会有等于1的啊

拉格朗日乘数法

这部分是高数中学过的知识，记得还算比较牢固。但是也只是学了最简单的等式约束的情况。例如：
在$xy=1$的情况下，求$x^2+y^2$的最大值。可以写成数学形式：
$$\max (x^2+y^2)$$$$s.t.xy=1$$
该怎么求呢？书中给出了一个比较简单的方法，就是可以构造一个拉格朗日函数：
$$L(x,y,\lambda )=x^2+y^2+\lambda (xy-1)$$
然后分别对$x,y,\lambda$求导联立就可以求出最后的结果。现在回过头来想，这只是给了我们一个既定的解决问题的程序，但却好像没有告诉我们为什么可以这样。下面稍微深入理解一下
考虑更一般的情况，求：
$$\max f(x,y)$$$$s.t.g(x,y)=0$$
对于凸函数来说，极值点就是最值点,设为$(x_0,y_0)$。通过画图就可以容易知道（应该也有严谨的推导）极值点的梯度方向一样,也就是平行，由向量平行的性质，可以得到：
$$\frac{f_x^{{}^{\prime }}(x_0,y_0)}{g_x^{{}^{\prime }}(x_0,y_0)}=\frac{f_y^{{}^{\prime }}(x_0,y_0)}{g_y^{{}^{\prime }}(x_0,y_0)}=-{\lambda }_0$$
上式变形得：
$$f_x^{{}^{\prime }}(x_0,y_0)+{\lambda }_0{g}_x^{{}^{\prime }}(x_0,y_0)=0$$$$f_y^{{}^{\prime }}(x_0,y_0)+{\lambda }_0{g}_y^{{}^{\prime }}(x_0,y_0)=0$$$$g(x_0,y_0)=0$$
通过上面三个式子就可以求出最后我们想要的结果了，这时我们可以回头看看前面的拉格朗日函数$L(x,y,\lambda )$,知道其实我最后通过$L(x,y,\lambda )$分别对$x,y,\lambda$求导得到的就是上面的三个式子。这也就是为什么按照书中规定的程序构造拉格朗日函数然后求导联立就可以求得最后结果的原因。这也只是一种比较简单的理解，还有些通过等高线什么来理解的，还没搞清楚，而且对于如果有两个或者多个约束条件的情况，上面的推导肯定就会更麻烦了，更更复杂的还有不等式约束的情况，慢慢理解。下面还是要理解理解不等式约束的情况，因为我们前面支持向量机的约束就是不等式的约束。

不等式约束下的拉格朗日法

这里加一下对于对偶问题的简单描述，所谓对偶问题，就是在一些条件下，$\min \max f(x)$和$\max \min f(x)$是相等的。具体的看后面的图片里的关于KKT条件的推导，其中会提到。
\quad 参考https://blog.csdn.net/lijil168/article/details/69395023中的推导，感觉自己才大概理解拉格朗日乘数法的由来，也顺带理解了KKT条件和对偶问题。就直接照着里面的图片敲一遍公式啦，加深一下印象。

让$L(x,\mu )=f(x)+{\sum }_{k=1}^q{\mu }_k{g}_k(x)$
其中${\mu }_k\ge 0,g_k(x)\le 0$
$$\because \begin{array}{c}{\mu }_k\ge 0\\ g_k(x)\le 0\end{array}\}=>\mu g(x)\le 0$$
$$\therefore \underset{\mu }{\max }L(x,\mu )=f(x)\text{\hspace{1em}(2)}$$
$$\therefore \min f(x)=\underset{x}{\min }\underset{\mu }{\max }L(x,\mu )\text{\hspace{1em}(3)}$$
$$\underset{\mu }{\max }\underset{x}{\min }L(x,\mu )=\underset{\mu }{\max }[\underset{x}{\min }f(x)+\underset{x}{\min }\mu g(x)]=\underset{\mu }{\max }\underset{x}{\min }f(x)+\underset{\mu }{\max }\underset{x}{\min }\mu g(x)=\min f(x)+\underset{\mu }{\max }\underset{x}{\min }\mu g(x)$$
又$\because \begin{array}{c}{\mu }_k\ge 0\\ g_k(x)\le 0\end{array}\}=>{\min }_x\mu g(x)=\left\{\begin{array}{cc}0& if\mu =0org(x)=0\\ -\infty & if\mu >0org(x)<0\end{array}\right\}$
$\therefore $放弃了，打公式太麻烦了，直接贴个图。打的这么多还先是留着吧，说不定后面可以复制粘贴。直接看图片吧

考虑等式约束和不等式约束都存在的情况，设目标函数 $f(x)$,不等式约束为$g(x)$,等式约束$h(x)$，此时约束优化问题描述如下：
$$\min f(x)$$ $$s.t.h_j(x)=0j=1,2...p$$ $$s.t.g_k(x)\le 0k=1,2...q$$
则定义拉格朗日函数为：
$$L(x,\mu )=f(x)+\sum _{j=1}^p{\lambda }_j{h}_j(x)+\sum _{k=1}^q{\mu }_k{g}_k(x)$$
其中$f(x)$是原目标函数，$h_j(x)$是第j个等式约束条件，${\lambda }_j$对应的约束系数,$g_k(x)$是不等式约束，${\mu }_k$是对应的约束系数想要得到我们最后的结论 ${\min }_x{\max }_{\mu }f(x)={\max }_{\mu }{\min }_{x}L(x,\mu )={\min }_{x}f(x)$, 必须满足下述条件,也就是KKT条件:(推导过程会用到这些条件）
KKT条件：
$$\frac{\partial L}{\partial x}{|}_{x=x^{\ast }}=0\text{\hspace{1em}(1*)}$$ $${\lambda }_j\ne 0\text{\hspace{1em}(2*)}$$ $${\mu }_k\ge 0\text{\hspace{1em}(3*)}$$ $${\mu }_k{g}_k(x^{\ast })=0\text{\hspace{1em}(4*)}$$ $$h_j(x^{\ast })=0j=1,2...,p\text{\hspace{1em}(5*)}$$ $$g_k(x\ast )\le 0k=1,2...,q\text{\hspace{1em}(6*)}$$
这些求解条件就是KKT条件。(1*)是对拉格朗日函数取极值时候带来的一个必要条件，要注意的是，这里的x可能是多维，比如$x_1,x_2$,所以是对$x_1,x_2$的偏导都等于0，(2*)是拉格朗日系数约束（同等式情况），(3*)是不等式约束情况，(4*)是互补松弛条件，(5*)、(6*)是原约束条件。
对于一般的任意问题而言，KKT条件是使一组解成为最优解的必要条件，当原问题是凸问题的时候，KKT条件也是充分条件。
稍微谈一下自己的理解，对于前面图片中的推导，最开始比较疑惑为什么由$\mu g(x)\le 0$就可以推出(2)中的${\max }_{\mu }L(x,\mu )=f(x)$，因为我认为既然$\mu g(x)\le 0$，那么可能不会取到0这个值，而是为一个负数，也就是$L(x,\mu )$的最大值可能比$f(x)$的最大值要大。想了想，感觉应该是这里的$g(x)\le 0$是一个约束范围，就比如限制$x^2+y^2\le 1$，这就是一个圆的区域，我们这个范围上的点都是可以取得到的，所以圆上的点当然也可以取到，而当取圆上的点时就是$g(x)=0$时，所以g(x)肯定是可以取到0这个值的。
还有为什么${\max }_{\mu }{\max }_{x}f(x)={\min }_{x}f(x)$?这是因为后面的f(x)和$\mu$没有关系啊

对支持向量机中条件约束问题求解

前面的问题为：
$$\underset{w,b}{\min }\frac{1}{2}||w|{|}^2$$ $$s.t.y_i(w^T{x}_i+b)\ge 1$$

将约束改写成
$$s.t.1-y_i(w^T{x}_i+b)\le 0$$
由前面图片中的推导得到的结论：${\min }_{x}f(x)={\max }_{\mu }{\min }_{x}L(x,\mu )$，我们可以先求${\min }_{x}L(x,\mu )$，也就是先求${\min }_{w,b}L(w,b,{\alpha }_1,...,{\alpha }_n)$。
\quad 这里是要保证所有的点都满足上面的约束，所以最后约束条件的个数和总共的点数（或者说实例数）是一样的。则拉格朗日函数$L(w,b,{\alpha }_1,...,{\alpha }_n)$可以写成为：
$$L(w,b,{\alpha }_1,{\alpha }_2,...{\alpha }_n)=\frac{1}{2}||w|{|}^2+\sum _{i=1}^n{\alpha }_i(1-y_i(w^T{x}_i+b))$$
由上面的拉格朗日函数和前面说的KKT条件，考虑只有不等式约束的情况，我们可以得到下面的一些式子。${\alpha }_i$就相当于高数中构造拉格朗日函数中的$\lambda$,这里有多个约束，所以会有很多$\alpha$
由(1*)可以得到：
$$\frac{\partial L}{\partial w}=w-\sum _{i=1}^n{\alpha }_i{y}_i{x}_i=0\text{\hspace{1em}(1\#)}$$
$$\frac{\partial L}{\partial b}=\sum _{i=1}^n{\alpha }_i{y}_i=0\text{\hspace{1em}(2\#)}$$
这里面可能涉及到2范数求导的问题，这里就先不研究了，总之上面的式子是对的。还要说明的是$x_1,x_2$…是不同的点，但是每一个不同的点会有许多维数，例如$x_1$可能会有横坐标也会有纵坐标可以记成$(x_1^1,x_1^2)$。我们式子中的$y$是对应的类别额，可以取-1,+1。（最前面说过一遍了已经）
由(3*)可以得到：
$${\alpha }_1\ge 0,{\alpha }_1\ge 0...{\alpha }_n\ge 0\text{\hspace{1em}(4\#)}$$
由(4*)可以得到前面的(2#)
$$\frac{\partial L}{\partial b}=\sum _{i=1}^n{\alpha }_i{y}_i=0\text{\hspace{1em}(2\#)}$$
由(6*)可以得到：
$$s.t.1-y_i(w^T{x}_i+b)\le 0\text{\hspace{1em}(3\#)}$$
说明一下，也就是说要满足KKT条件，我们才可以通过拉格朗日法求最大或者最小值。这里面有的KKT条件最开始已经满足了例如(3#)，有的是为了满足KKT条件而出现的例如(1#)(2#)。
\quad 将(1#)和(2#)带入拉格朗日函数$L(w,b,{\alpha }_1,...{\alpha }_n)$中，特别注意(1#)(2#)在化简式子中的作用，注意$x_i,x_j,w$为列向量，注意求和符号后的变量如果和求和符号的下标i or j没有关系，就可以移到求和符号的外面或者里面 有：

min ⁡ w , b L ( w , b , α 1 , . . . α n ) = min ⁡ w , b { 1 2 ∣ ∣ w ∣ ∣ 2 + ∑ i = 1 n α i ( 1 − y i ( w T x i + b ) ) } = min ⁡ w , b { 1 2 w T w + ∑ i = 1 n α i − ∑ i = 1 n α i y i w T x i − ∑ i = 1 n α i y i b } = min ⁡ w , b { 1 2 ∑ i = 1 n α i y i x i T ∑ i = 1 n α i y i x i + ∑ i = 1 n α i − ∑ i = 1 n α i y i w T x i − b ∑ i = 1 n α i y i = min ⁡ w , b { 1 2 ∑ i = 1 n ∑ j = 1 n α i y i α j y j x i T x j + ∑ i = 1 n α i − ∑ i = 1 n α i y i ∑ j = 1 n α j y j x i T x i } = − 1 2 ∑ i = 1 n ∑ j = 1 n α i y i α j y j x i T x j + ∑ i = 1 n α i \begin {aligned} \min_{w,b}L(w,b,\alpha_1,...\alpha_n)&=\min_{w,b}\{ \frac{1}{2}||w||^2+ \sum_{i=1}^n \alpha_i(1-y_i(w^Tx_i+b) ) \} \\ &=\min_{w,b}\{\frac{1}{2} w^Tw +\sum_{i=1}^n\alpha_i -\sum_{i=1}^n\alpha_iy_iw^Tx_i-\sum_{i=1}^n\alpha_i y_i b \} \\ &=\min_{w,b}\{\frac{1}{2} \sum_{i=1}^n \alpha_i y_i x_i ^T \sum_{i=1}^n \alpha_i y_i x_i +\sum_{i=1}^n\alpha_i -\sum_{i=1}^n\alpha_iy_iw^Tx_i - b\sum_{i=1}^n\alpha_i y_i \\ &=\min_{w,b} \{ \frac{1}{2}\sum_{i=1}^n \sum_{j=1}^n\alpha_iy_i \alpha_jy_j x_i^T x_j + \sum_{i=1}^n\alpha_i - \sum_{i=1}^n \alpha_iy_i \sum_{j=1}^n\alpha_jy_j x_i^T x_i \}\\ &= -\frac{1}{2}\sum_{i=1}^n \sum_{j=1}^n\alpha_iy_i \alpha_jy_j x_i^T x_j + \sum_{i=1}^n\alpha_i \end{aligned} w,bmin​L(w,b,α1​,...αn​)​=w,bmin​{21​∣∣w∣∣2+i=1∑n​αi​(1−yi​(wTxi​+b))}=w,bmin​{21​wTw+i=1∑n​αi​−i=1∑n​αi​yi​wTxi​−i=1∑n​αi​yi​b}=w,bmin​{21​i=1∑n​αi​yi​xiT​i=1∑n​αi​yi​xi​+i=1∑n​αi​−i=1∑n​αi​yi​wTxi​−bi=1∑n​αi​yi​=w,bmin​{21​i=1∑n​j=1∑n​αi​yi​αj​yj​xiT​xj​+i=1∑n​αi​−i=1∑n​αi​yi​j=1∑n​αj​yj​xiT​xi​}=−21​i=1∑n​j=1∑n​αi​yi​αj​yj​xiT​xj​+i=1∑n​αi​​

所以：
min ⁡ w , b 1 2 ∣ ∣ w ∣ ∣ 2 = max ⁡ α 1 , . . . , α n min ⁡ w , b L ( w , b , α 1 , . . . α n ) = max ⁡ α 1 , . . . , α n { − 1 2 ∑ i = 1 n ∑ j = 1 n α i y i α j y j x i T x j + ∑ i = 1 n α i } \begin{aligned} \min_{w,b} \frac{1}{2} ||w||^2 &= \max_{\alpha_1,...,\alpha_n} \min_{w,b}L(w,b,\alpha_1,...\alpha_n) \\ &=\max_{\alpha_1,...,\alpha_n} \{- \frac{1}{2}\sum_{i=1}^n \sum_{j=1}^n\alpha_iy_i \alpha_jy_j x_i ^Tx_j + \sum_{i=1}^n\alpha_i \} \end{aligned} w,bmin​21​∣∣w∣∣2​=α1​,...,αn​max​w,bmin​L(w,b,α1​,...αn​)=α1​,...,αn​max​{−21​i=1∑n​j=1∑n​αi​yi​αj​yj​xiT​xj​+i=1∑n​αi​}​
我们已经将最开始的问题简化成只和$\alpha ,x,y$有关的问题了，其中$x$点的坐标，$y$是其所对应的类别 ，这两个是我们已知的。$\alpha$是未知的，需要我们优化，但是要注意$\alpha$也是有约束条件的，如(4#)(2#) 或所以最终的问题变成了：
max ⁡ α 1 , . . . , α n { ∑ i = 1 n α i − 1 2 ∑ i = 1 n ∑ j = 1 n α i y i α j y j x i T x j } \max_{\alpha_1,...,\alpha_n} \{ \sum_{i=1}^n\alpha_i - \frac{1}{2}\sum_{i=1}^n \sum_{j=1}^n\alpha_iy_i \alpha_jy_j x_i ^Tx_j \} α1​,...,αn​max​{i=1∑n​αi​−21​i=1∑n​j=1∑n​αi​yi​αj​yj​xiT​xj​}$$s.t.{\alpha }_1\ge 0,{\alpha }_1\ge 0...{\alpha }_n\ge 0$$$$\sum _{i=1}^n{\alpha }_i{y}_i=0$$
一般还会引入松弛变量C（现在没搞明白松弛变量是怎么回事），令所有的$\alpha <C$，所以最后可以写成：
max ⁡ α 1 , . . . , α n { ∑ i = 1 n α i − 1 2 ∑ i = 1 n ∑ j = 1 n α i y i α j y j x i T x j } \max_{\alpha_1,...,\alpha_n} \{ \sum_{i=1}^n\alpha_i - \frac{1}{2}\sum_{i=1}^n \sum_{j=1}^n\alpha_iy_i \alpha_jy_j x_i^Tx_j \} α1​,...,αn​max​{i=1∑n​αi​−21​i=1∑n​j=1∑n​αi​yi​αj​yj​xiT​xj​}$$s.t.0\le \alpha \le C$$$$\sum _{i=1}^n{\alpha }_i{y}_i=0$$
这一篇就先写到这吧，~

参考

https://www.jianshu.com/p/e74eb43960a1(这是latex公式的一些用法)
https://zhuanlan.zhihu.com/p/88249368?utm_source=qq(这是latex公式的一些用法)
https://www.cnblogs.com/xxrxxr/p/7535587.html
b站深度之眼官方帐号
https://blog.csdn.net/hao5335156/article/details/82320082
https://www.cnblogs.com/ooon/p/5721119.html
https://zhuanlan.zhihu.com/p/99945521（拉格朗日乘数法推导）
https://blog.csdn.net/lijil168/article/details/69395023（不等式约束拉格朗日乘数法推导）