支持向量机(support vector machines,SVM)是一种二类分类模型。它的基本模型是定义在特征空间上的间隔最大的线性分类器,间隔最大使它有别于感知机;支持向量机还包括核技巧,这使它成为实质上的非线性分类器。支持向量机的学习策略就是间隔最大化,可形式化为一个求解凸二次规划(convex quadratic programming)的问题,也等价于正则化的合页损失函数的最小化问题。支持向量机的学习算法是求解凸二次规划的最优化算法。
支持向量机学习方法包含构建由简至繁的模型:线性可分支持向量机( linear support vector machine in linearly separable case )、线性支持向量机( linear support vector machine)及非线性支持向量机(non-linear support vector machine)。简单模型是复杂模型的基础,也是复杂模型的特殊情况。当训练数据线性可分时,通过硬间隔最大化( hard margin maximization),学习一个线性的分类器,即线性可分支持向量机,又称为硬间隔支持向量机;当训练数据近似线性可分时,通过软间隔最大化( soft margin maximization),也学习一个线性的分类器,即线性支持向量机,又称为软间隔支持向量机;当训练数据线性不可分时,通过使用核技巧(kemel trick)及软间隔最大化,学习非线性支持向量机。
当输入空间为欧氏空间或离散集合、特征空间为希尔伯特空间时,核函数(kernel function)表示将输入从输入空间映射到特征空间得到的特征向量之间的内积。通过使用核函数可以学习非线性支持向量机,等价于隐式地在高维的特征空间中学习线性支持向量机。这样的方法称为核技巧。核方法( kernel method)是比支持向量机更为一般的机器学习方法。
Cortes与Vapnik提出线性支持向童机,Boser、Guyon与Vapnik又引入核技巧,提出非线性支持向量机。
本章按照上述思路介绍3类支持向量机、核函数及一种快速学习算法——序列最小最优化算法(SMO)。
线性可分支持向量机
考虑一个二分类问题,假设输入空间与特征空间为两个不同的空间。输入空间为欧式空间或离散集合,特征空间为欧式空间或希尔伯特空间。
线性可分支持向量机、线性支持向量机假设这两个空间的元素一一对应,并将输入空间中的输入映射为特征空间中的特征向量。非线性支持向量机利用一个从输入空间到特征空间的非线性映射,将输入映射为特征向量,使得在特征空间中能够使用线性支持向量机来学习。所以,输入都由输入空间转换到特征空间,支持向量机的学习是在特征空间上进行的。
假设给定一个特征空间上的训练数据集
%2C(x_2%2Cy_2)%2C…%2C(x_N%2Cy_N)%20%5C%7D%0A#card=math&code=T%3D%5C%7B%20%28x_1%2Cy_1%29%2C%28x_2%2Cy_2%29%2C…%2C%28x_N%2Cy_N%29%20%5C%7D%0A&height=20&width=266)
其中,$x_i\in X=R^n$$,$y_i\in Y={+1,-1}, i=1,2,…,N,x_i为第i个特征向量,也称为实例,y_i为x_i的类标记,当y_i=+1时,称x_i为正例;当y_i=-1时,称x_i为负例,(x_i,y_i)$$称为样本点。再假设训练数据集是线性可分的。
学习的目标是在特征空间中找到一个分离超平面,能将实例分到不同的类。分离超平面对应于方程w\cdot x+b=0,它由法向量w和截距b决定,可用(w,b)来表示。分离超平面将特征空间划分为两部分,一部分是正类,一部分是负类。法向量指向的一侧为正类,另一侧为负类。
一般地,当训练数据集线性可分时,存在无穷个分离超平面可将两类数据正确分开。感知机利用误分类最小的策略,求得分离超平面,不过这时的解有无穷多个。线性可分支持向量机利用间隔最大化求最优分离超平面,这时,解是唯一的。
线性可分支持向量机的定义:
给定线性可分训练数据集,通过间隔最大化或等价地求解相应的凸二次规划问题学习得到的分离超平面为
以及相应的分类决策函数
%3D%5Ctext%7Bsign%7D(w%5E%5Ccdot%20x%2Bb%5E)%0A#card=math&code=f%28x%29%3D%5Ctext%7Bsign%7D%28w%5E%2A%5Ccdot%20x%2Bb%5E%2A%29%0A&height=20&width=173)
称为线性可分支持向量机。
考虑如下图所示的二维特征空间中的分类问题。图中的“。”表示正例,“x”表示负例。训练数据集线性可分,这时有许多直线能将两类数据正确划分。线性可分支持向量机对应将两类数据正确划分并且间隔最大的直线,如下图所示。

间隔最大及相应的约束最优化问题将在下面叙述。这里先介绍函数间隔和集合间隔的概念。
函数间隔和几何间隔
在上图中,有A,B,C三个点,表示三个实例,均在分离超平面的正类一侧,预测他们的类。点A距分离超平面较远,若预测该点为正类,就比较确信预测是正确的;点C距分离超平面较近,若预测该点为正类就不那么确信;点B介于点A与C之间,预测其为正类的确信度也在A与C之间。
一般来说,一个点距离分离超平面的远近可以表示分类预测的确信程度。在超平面w\cdot x+b=0确定的情况下,|w\cdot x+b|能够相对地表示点x距离超平面的远近。而w\cdot x+b的符号与类标记y的符号是否一致能够表示分类是否正确。所以可用量y(w\cdot x+b)来表示分类的正确性及确信度,这就是函数间隔(functional margin)的概念。
疑问:为什么|w\cdot x+b|能够相对地表示点x距离超平面的远近?
因为:根据点到直线的距离计算公式为:设直线L的方程为Ax+By+C=0,点P的坐标为(x_0, y_0),则点P到直线L的距离为:
函数间隔的定义
对于给定的训练数据集T和超平面(w,b),定义超平面(w,b)关于样本点(x,y)的函数间隔为
%0A#card=math&code=%5Chat%7B%5Cgamma%7D_i%3Dy_i%28w%5Ccdot%20x_i%2Bb%29%0A&height=20&width=131)
定义超平面关于训练数据集T的函数间隔为超平面(w,b)关于T中所有样本点(x_i,y_i)的函数间隔之最小值,即
函数间隔可以表示分类预测的正确性及确信度。但是选择分离超平面时,只有函数间隔还不够。因为只要成比例地改变w和b,例如将它们改为2w和2b,超平面并没有改变,但函数间隔却成为原来的2倍。这一事实启示我们,可以对分离超平面的法向量w加某些束,如规范化,||w|| = 1,使得间隔是确定的。这时函数间隔成为几何间隔(geometric margin)。
下图给出了超平面(w,b)及其法向量w。点A表示某一实例其类标记为y_i=+1。点A与超平面(w,b)的距离由线段AB给出,记作\gamma_i
其中,||w||为w的L_2范数。这是点A在超平面正的一侧的情形。如果点Z在超平面负的一侧,即y_i=-1,那么点与超平面的距离为
%0A#card=math&code=%5Cgamma_i%3D-%5Cleft%28%5Cfrac%7Bw%7D%7B%7C%7Cw%7C%7C%7D%5Ccdot%20x_i%2B%5Cfrac%7Bb%7D%7B%7C%7Cw%7C%7C%7D%5Cright%29%0A&height=45&width=199)
一般地,当样本点(x_i,y_i)被超平面(w,b)正确分类时,点x_i与超平面(w,b)的距离是
%0A#card=math&code=%5Cgamma_i%3Dy_i%5Cleft%28%5Cfrac%7Bw%7D%7B%7C%7Cw%7C%7C%7D%5Ccdot%20x_i%2B%5Cfrac%7Bb%7D%7B%7C%7Cw%7C%7C%7D%5Cright%29%0A&height=45&width=200)
由这一事实导出几何间隔的概念。

几何间隔的定义:
对于给定的训练数据集T和超平面(w,b),定义超平面(w,b)关于样本点(x_i,y_i)的几何间隔为
%0A#card=math&code=%5Cgamma_i%3Dy_i%5Cleft%28%5Cfrac%7Bw%7D%7B%7C%7Cw%7C%7C%7D%5Ccdot%20x_i%2B%5Cfrac%7Bb%7D%7B%7C%7Cw%7C%7C%7D%5Cright%29%0A&height=45&width=200)
定义超平面(w,b)关于训练数据集T的几何间隔为超平面(w,b)关于T所有样本点(x_i,y_i)的几何间隔之最小值,即
超平面(w,b)关于样本点(x_i,y_i)的几何间隔一般是实例点到超平面的带符号的距离(signed distance),当样本点被超平面正确分类时就是实例点到超平面的距离。
从函数间隔和几何间隔的定义可知,函数间隔和几何间隔有下面的关系:
如果||w||=1,那么函数间隔和几何间隔相等。如果超平面参数w和b成比例地改变(超平面没有改变),函数间隔也按此比例改变,而几何间隔不变。
间隔最大化
支持向量机学习的基本想法是求解能够正确划分训练数据集并且几何间隔最大的分离超平面。对线性可分的训练数据集而言,线性可分分离超平面有无穷多个(等价于感知机),但是几何间隔最大的分离超平面是唯一的。这里的间隔最大化又称为硬间隔最大化(与将要讨论的训练数据集近似线性可分时的软间隔最大化相对应)。
间隔最大化的直观解释是:对训练数据集找到几何间隔最大的超平面意味着以充分大的确信度对训练数据进行分类。也就是说,不仅将正负实例点分开,而且对最难分的实例点(离超平面最近的点)也有足够大的确信度将它们分开。这样的超平面应该对未知的新实例有很好的分类预测能力。
最大间隔分离超平面
下面考虑如何求得一个几何间隔最大的分离超平面,即最大间隔分离超平面。具体地,这个问题可以表示为下面的约束最优化问题:
%5Cgeqslant%20%5Cgamma%2C%5C%20i%3D1%2C2%2C…%2CN%5C%5C%0A%5Cend%7Baligned%7D%0A#card=math&code=%5Cbegin%7Baligned%7D%0A%26%5Cmathop%7B%5Ctext%7Bmax%7D%7D_%7Bw%2Cb%7D%5Cquad%20%5Cgamma%5C%5C%0A%26%5Ctext%7Bs.t.%7D%5Cquad%20y_i%5Cleft%28%5Cfrac%7Bw%7D%7B%7C%7Cw%7C%7C%7D%5Ccdot%20x_i%2B%5Cfrac%7Bb%7D%7B%7C%7Cw%7C%7C%7D%5Cright%29%5Cgeqslant%20%5Cgamma%2C%5C%20i%3D1%2C2%2C…%2CN%5C%5C%0A%5Cend%7Baligned%7D%0A&height=78&width=356)
即我们希望最大化超平面(w,b)关于训练数据集的几何间隔\gamma,约束条件表示的是超平面(w,b)关于每个训练样本点的几何间隔至少是\gamma。
考虑几何间隔和函数间隔的关系式,可将这个问题改写为
%5Cgeqslant%20%5Chat%7B%5Cgamma%7D%2C%5C%20i%3D1%2C2%2C…%2CN%5C%5C%0A%5Cend%7Baligned%7D%0A#card=math&code=%5Cbegin%7Baligned%7D%0A%26%5Cmathop%7B%5Ctext%7Bmax%7D%7D_%7Bw%2Cb%7D%5Cquad%20%5Cfrac%7B%5Chat%7B%5Cgamma%7D%7D%7B%7C%7Cw%7C%7C%7D%5C%5C%0A%26%5Ctext%7Bs.t.%7D%5Cquad%20y_i%5Cleft%28w%5Ccdot%20x_i%2Bb%5Cright%29%5Cgeqslant%20%5Chat%7B%5Cgamma%7D%2C%5C%20i%3D1%2C2%2C…%2CN%5C%5C%0A%5Cend%7Baligned%7D%0A&height=69&width=290)
函数间隔\hat{\gamma}的取值并不影响最优化问题的解。事实上,假设将w和b按比例改变为\lambda w和\lambda b,这时函数间隔为λ\hat{\gamma}。函数间隔的这一改变对上面最优化问题的不等式约束没有影响,对目标函数的优化也没有影响,也就是说,它产生一个等价的最优化问题。这样,就可以取函数间隔\hat{\gamma}=1。
将\hat{\gamma}=1带入上面的最优化问题,注意到最大化\frac{1}{||w||}和最小化\frac{1}{2}||w||^2是等价的,于是就得到下面的线性可分支持向量机学习的最优化问题
-1%5Cgeqslant%200%2C%5C%20i%3D1%2C2%2C…%2CN%5C%5C%0A%5Cend%7Baligned%7D%0A#card=math&code=%5Cbegin%7Baligned%7D%0A%26%5Cmathop%7B%5Ctext%7Bmin%7D%7D_%7Bw%2Cb%7D%5Cquad%20%5Cfrac%7B1%7D%7B2%7D%7C%7Cw%7C%7C%5E2%5C%5C%0A%26%5Ctext%7Bs.t.%7D%5Cquad%20y_i%5Cleft%28w%5Ccdot%20x_i%2Bb%5Cright%29-1%5Cgeqslant%200%2C%5C%20i%3D1%2C2%2C…%2CN%5C%5C%0A%5Cend%7Baligned%7D%0A&height=64&width=318)
由于$\mathop{\text{min}}_{w,b}\frac{1}{2}||w||^2$是一个凸函数,所以这是一个凸二次规划问题。
凸优化是指约束优化问题
%26%5C%5C%0A%26%5Ctext%7Bs.t.%7D%5C%20%5Cquad%20gi(w)%5Cleqslant%200%2C%5C%20i%3D1%2C2%2C…%2Ck%5C%5C%0A%26%5C%20%5C%20%5Cquad%20%5Cquad%20h_i(w)%3D0%2C%5C%20i%3D1%2C2%2C…%2Cl%5C%5C%0A%5Cend%7Baligned%7D%0A#card=math&code=%5Cbegin%7Baligned%7D%0A%26%5Cmathop%7B%5Ctext%7Bmin%7D%7D%7Bw%7D%5Cquad%20f%28w%29%26%5C%5C%0A%26%5Ctext%7Bs.t.%7D%5C%20%5Cquad%20g_i%28w%29%5Cleqslant%200%2C%5C%20i%3D1%2C2%2C…%2Ck%5C%5C%0A%26%5C%20%5C%20%5Cquad%20%5Cquad%20h_i%28w%29%3D0%2C%5C%20i%3D1%2C2%2C…%2Cl%5C%5C%0A%5Cend%7Baligned%7D%0A&height=74&width=263)
其中,目标函数$f(w)和约束函数g_i(x)都是Rn$$上的仿射函数。
当目标函数f(w)是二次函数且约束函数g_i(w)是仿射函数时,上述凸最优化问题成了凸二次规划问题。
如果求出了约束最优化问题的解w^{}, b\cdot x+b\cdot x+b^),即线性可分支持向量机模型。
综上所述,就有下面的线性可分支持向量机的学习算法——最大间隔法(maximum margin method)。
线性可分支持向量机学习算法:最大间隔法
输入:线性可分训练数据
%2C(x_2%2Cy_2)%2C…%2C(x_N%2Cy_N)%20%5C%7D%0A#card=math&code=T%3D%5C%7B%20%28x_1%2Cy_1%29%2C%28x_2%2Cy_2%29%2C…%2C%28x_N%2Cy_N%29%20%5C%7D%0A&height=20&width=266)
其中,x_i\in X=R^n,y_i\in Y={+1,-1}, i=1,2,…,N;
输出:最大间隔分离超平面和分类决策函数
(1)构造并求解约束最优化问题:
-1%5Cgeqslant%200%2C%5C%20i%3D1%2C2%2C…%2CN%5C%5C%0A%5Cend%7Baligned%7D%0A#card=math&code=%5Cbegin%7Baligned%7D%0A%26%5Cmathop%7B%5Ctext%7Bmin%7D%7D_%7Bw%2Cb%7D%5Cquad%20%5Cfrac%7B1%7D%7B2%7D%7C%7Cw%7C%7C%5E2%5C%5C%0A%26%5Ctext%7Bs.t.%7D%5Cquad%20y_i%5Cleft%28w%5Ccdot%20x_i%2Bb%5Cright%29-1%5Cgeqslant%200%2C%5C%20i%3D1%2C2%2C…%2CN%5C%5C%0A%5Cend%7Baligned%7D%0A&height=64&width=318)
求得最优解w^{}, b^{}。
(2)由此得到分离超平面:
分类决策函数
%3D%5Ctext%7Bsign%7D(x%5E%5Ccdot%20x%2Bb%5E)%0A#card=math&code=f%28x%29%3D%5Ctext%7Bsign%7D%28x%5E%2A%5Ccdot%20x%2Bb%5E%2A%29%0A&height=20&width=171)
最大间隔分离超平面的存在唯一性
线性可分训练数据集的最大间隔分离超平面是存在且唯一的。
证明此处略,详见李航的《统计学习基础》对应章节。
支持向量和间隔边界
在线性可分情况下,训练数据的样本点中与分离超平面最近的样本点的实例称为支持向量。支持向量是使不等式约束条件的等号成立的点,即
-1%3D0%0A#card=math&code=y_i%28w%5Ccdot%20x_i%2Bb%29-1%3D0%0A&height=20&width=154)
对y_i=+1的正例点,支持向量在超平面H_1:w\cdot x+b=1上,对y_i=-1的负例点,支持向量在超平面H_2: w\cdot x+b=-1上。如下图所示,在H_1和H_2上的点就是支持向量。

注意到H_1和H_2平行,并旦没有实例点落在它们中间。在H_1与H_2之间形成一条长带,分离超平面与它们平行且位于它们中央。长带的宽度,即H_1与H_2之间的距离称为间隔(margin)。间隔依赖于分离超平面的法向量w,等于\frac{2}{||w||},H_1和H_2称为间隔边界。
在决定分离超平面时只有支持向量起作用,而其他实例点并不起作用。如果移动支持向量将改变所求的解:但是如果在间隔边界以外移动其他实例点,甚至去掉这些点,则解是不会改变的。由于支持向量在确定分离超平面中起着决定性作用,所以将这种分类模型称为支持向量机。支持向量的个数一般很少,所以支持向量机由很少的“重要的”训练样本确定。
线性可分支持向量机学习的对偶算法
对偶优化问题
为了求解线性可分支持向量机的最优化问题,将它作为原始最优化问题,应用拉格朗日对偶性,通过求解对偶问题,得到原始问题的最优解,这就是线性可分支持向量机的对偶算法。这样做的优点,一是对偶问题往往更容易求解,二是自然引入核函数,进而推广到非线性分类问题。
首先构建拉格朗日函数。为此,对每一个不等式约束引入拉格朗日乘子\alpha_i\geqslant 0, i=1,2,…,N,定义拉格朗日函数:
%3D%5Cfrac%7B1%7D%7B2%7D%7C%7Cw%7C%7C%5E2-%5Csum%7Bi%3D1%7D%5EN%5Calpha_iy_i(w%5Ccdot%20x_i%2Bb)%2B%5Csum%7Bi%3D1%7D%5EN%5Calphai%0A#card=math&code=L%28w%2Cb%2C%5Calpha%29%3D%5Cfrac%7B1%7D%7B2%7D%7C%7Cw%7C%7C%5E2-%5Csum%7Bi%3D1%7D%5EN%5Calphaiy_i%28w%5Ccdot%20x_i%2Bb%29%2B%5Csum%7Bi%3D1%7D%5EN%5Calpha_i%0A&height=53&width=367)
其中,\alpha=(\alpha_1,\alpha_2,…,\alpha_N)为拉格朗日乘子向量。
所以,通过拉格朗日函数我们可以将原始问题
-1%5Cgeqslant%200%2C%5C%20i%3D1%2C2%2C…%2CN%5C%5C%0A%5Cend%7Baligned%7D%0A#card=math&code=%5Cbegin%7Baligned%7D%0A%26%5Cmathop%7B%5Ctext%7Bmin%7D%7D_%7Bw%2Cb%7D%5Cquad%20%5Cfrac%7B1%7D%7B2%7D%7C%7Cw%7C%7C%5E2%5C%5C%0A%26%5Ctext%7Bs.t.%7D%5Cquad%20y_i%5Cleft%28w%5Ccdot%20x_i%2Bb%5Cright%29-1%5Cgeqslant%200%2C%5C%20i%3D1%2C2%2C…%2CN%5C%5C%0A%5Cend%7Baligned%7D%0A&height=64&width=318)
转化为:
%5C%5C%0A%26%5Ctext%7Bs.t.%7D%5Cquad%20%5Calphai%5Cgeqslant%200%5C%5C%0A%5Cend%7Baligned%7D%0A#card=math&code=%5Cbegin%7Baligned%7D%0A%26%5Cmathop%7B%5Ctext%7Bmin%7D%7D%7Bw%2Cb%7D%5Cmathop%7B%5Ctext%7Bmax%7D%7D_%7B%5Clambda%7D%5Cquad%20L%28w%2Cb%2C%5Calpha%29%5C%5C%0A%26%5Ctext%7Bs.t.%7D%5Cquad%20%5Calpha_i%5Cgeqslant%200%5C%5C%0A%5Cend%7Baligned%7D%0A&height=54&width=153)
有的小伙伴这里不理解,为什么这里突然变成了拉格朗日函数最大值的最小化?原因是
%3E0%2C%5Cquad%20%5Cmathop%7B%5Ctext%7Bmax%7D%7D%7B%5Clambda%7DL(w%2Cb%2C%5Calpha)%3D%5Cfrac%7B1%7D%7B2%7D%7C%7Cw%7C%7C%5E2%2B%5Cinfty%3D%5Cinfty%5C%5C%0A%26%5Ctext%7Bif%7D%5Cquad%201-y_i(w%5ETx_i%2Bb)%5Cleqslant%200%2C%5Cquad%20%5Cmathop%7B%5Ctext%7Bmax%7D%7D%7B%5Clambda%7DL(w%2Cb%2C%5Calpha)%3D%5Cfrac%7B1%7D%7B2%7D%7C%7Cw%7C%7C%5E2%2B0%3D%5Cfrac%7B1%7D%7B2%7D%7C%7Cw%7C%7C%5E2%5C%5C%0A%5Cend%7Baligned%7D%0A#card=math&code=%5Cbegin%7Baligned%7D%0A%26%5Ctext%7Bif%7D%5Cquad%201-yi%28w%5ETx_i%2Bb%29%3E0%2C%5Cquad%20%5Cmathop%7B%5Ctext%7Bmax%7D%7D%7B%5Clambda%7DL%28w%2Cb%2C%5Calpha%29%3D%5Cfrac%7B1%7D%7B2%7D%7C%7Cw%7C%7C%5E2%2B%5Cinfty%3D%5Cinfty%5C%5C%0A%26%5Ctext%7Bif%7D%5Cquad%201-yi%28w%5ETx_i%2Bb%29%5Cleqslant%200%2C%5Cquad%20%5Cmathop%7B%5Ctext%7Bmax%7D%7D%7B%5Clambda%7DL%28w%2Cb%2C%5Calpha%29%3D%5Cfrac%7B1%7D%7B2%7D%7C%7Cw%7C%7C%5E2%2B0%3D%5Cfrac%7B1%7D%7B2%7D%7C%7Cw%7C%7C%5E2%5C%5C%0A%5Cend%7Baligned%7D%0A&height=81&width=492)
所以,由于满足1-y_i(w^Tx_i+b)\leqslant 0,所以是等价的。
根据拉格朗日对偶性,原始问题的对偶问题是极大极小问题:
%0A#card=math&code=%5Cmathop%7B%5Ctext%7Bmax%7D%7D%7B%5Calpha%7D%5Cmathop%7B%5Ctext%7Bmin%7D%7D%7Bw%2Cb%7DL%28w%2Cb%2C%5Calpha%29%0A&height=31&width=134)
通俗的讲,这个过程的主要操作就是将min和max互掉位置,并且二者之间有一个性质,即前者>=后者,这就好比在高个子人群中挑一个身高较矮的(min max)要高于在矮个子人群中挑一个身高较高的(max min)。默认情况下二者是呈弱对偶关系的,但在此目标函数和约束条件下是呈强对偶关系(等价关系)的。
**等等,“根据拉格朗日对偶性”是什么意思?你需要看看凸优化一章。。如果拉格朗日函数具有强对偶性,那么原函数(极小极大)的最优值p的值是一样大的。
那怎么判断是否满足强对偶性呢?其中一个就是,在原函数和不等式约束是凸函数的基础上,如果满足Slater条件,那么强对偶性就成立,从而原函数和对偶函数的最优值就相等,就可以通过求对偶函数的最优值来求原函数的最优值了。
那Slater条件又是什么呢?Slater条件是说:存在x,使不等式约束中的“小于等于号”要严格取到“小于号”(也就是弹簧振子系统要能正常运转,等于零就是有约束力)。这个条件一般都能满足。。所以,我们的这个函数显然很容易就满足了Slater条件,而且它又是凸函数,那么就是满足强对偶性质喽,那么原函数(最小最大)和对偶函数(极大极小)的最优值就是一样的啦。所以,既然原函数不好求最优值,那我们就求对偶函数的最优值,因为对偶函数好求最优值嘛。。求出来的对偶函数的最优值,就是原函数的最优值了。。
所以,为了求得对偶问题的解,需要先求L(w,b,\alpha)对w,b的极小,再求对\alpha的极大。
(1)求\mathop{\text{min}}_{w,b}L(w,b,\alpha)
将拉格朗日函数L(w,b,\alpha)分别对w,b求偏导数,并令其等于0。
%26%3Dw-%5Csum%7Bi%3D1%7D%5EN%5Calpha_iy_ix_i%3D0%5C%5C%0A%5Cbigtriangledown_bL(w%2Cb%2C%5Calpha)%26%3D-%5Csum%7Bi%3D1%7D%5EN%5Calphaiy_i%3D0%5C%5C%0A%5Cend%7Baligned%7D%0A#card=math&code=%5Cbegin%7Baligned%7D%0A%5Cbigtriangledown_wL%28w%2Cb%2C%5Calpha%29%26%3Dw-%5Csum%7Bi%3D1%7D%5EN%5Calphaiy_ix_i%3D0%5C%5C%0A%5Cbigtriangledown_bL%28w%2Cb%2C%5Calpha%29%26%3D-%5Csum%7Bi%3D1%7D%5EN%5Calpha_iy_i%3D0%5C%5C%0A%5Cend%7Baligned%7D%0A&height=107&width=259)
得
将上式中第一行式子带入拉格朗日函数,并利用上式中的第二行式子,即得:
%26%3D%5Cfrac%7B1%7D%7B2%7D%7C%7Cw%7C%7C%5E2-%5Csum%7Bi%3D1%7D%5EN%5Calpha_iy_i(w%5Ccdot%20x_i%2Bb)%2B%5Csum%7Bi%3D1%7D%5EN%5Calphai%5C%5C%0A%26%3D%5Cfrac%7B1%7D%7B2%7D%5Csum%7Bi%3D1%7D%5EN%5Csum%7Bj%3D1%7D%5EN%5Calpha_i%5Calpha_jy_iy_j(x_i%5Ccdot%20x_j)-%5Csum%7Bi%3D1%7D%5EN%5Calphaiy_i%5Cleft(%5Cleft(%20%5Csum%7Bj%3D1%7D%5EN%5Calphajy_jx_j%20%5Cright)%5Ccdot%20x_i%2Bb%5Cright)%2B%5Csum%7Bi%3D1%7D%5EN%5Calphai%5C%5C%0A%26%3D-%5Cfrac%7B1%7D%7B2%7D%5Csum%7Bi%3D1%7D%5EN%5Csum%7Bj%3D1%7D%5EN%5Calpha_i%5Calpha_jy_iy_j(x_i%5Ccdot%20x_j)%2B%5Csum%7Bi%3D1%7D%5EN%5Calphai%5C%5C%0A%5Cend%7Baligned%7D%0A#card=math&code=%5Cbegin%7Baligned%7D%0AL%28w%2Cb%2C%5Calpha%29%26%3D%5Cfrac%7B1%7D%7B2%7D%7C%7Cw%7C%7C%5E2-%5Csum%7Bi%3D1%7D%5EN%5Calphaiy_i%28w%5Ccdot%20x_i%2Bb%29%2B%5Csum%7Bi%3D1%7D%5EN%5Calphai%5C%5C%0A%26%3D%5Cfrac%7B1%7D%7B2%7D%5Csum%7Bi%3D1%7D%5EN%5Csum%7Bj%3D1%7D%5EN%5Calpha_i%5Calpha_jy_iy_j%28x_i%5Ccdot%20x_j%29-%5Csum%7Bi%3D1%7D%5EN%5Calphaiy_i%5Cleft%28%5Cleft%28%20%5Csum%7Bj%3D1%7D%5EN%5Calphajy_jx_j%20%5Cright%29%5Ccdot%20x_i%2Bb%5Cright%29%2B%5Csum%7Bi%3D1%7D%5EN%5Calphai%5C%5C%0A%26%3D-%5Cfrac%7B1%7D%7B2%7D%5Csum%7Bi%3D1%7D%5EN%5Csum%7Bj%3D1%7D%5EN%5Calpha_i%5Calpha_jy_iy_j%28x_i%5Ccdot%20x_j%29%2B%5Csum%7Bi%3D1%7D%5EN%5Calpha_i%5C%5C%0A%5Cend%7Baligned%7D%0A&height=166&width=617)
即
%3D-%5Cfrac%7B1%7D%7B2%7D%5Csum%7Bi%3D1%7D%5EN%5Csum%7Bj%3D1%7D%5EN%5Calphai%5Calpha_jy_iy_j(x_i%5Ccdot%20x_j)%2B%5Csum%7Bi%3D1%7D%5EN%5Calphai%0A#card=math&code=%5Cmathop%7B%5Ctext%7Bmin%7D%7D%7Bw%2Cb%7DL%28w%2Cb%2C%5Calpha%29%3D-%5Cfrac%7B1%7D%7B2%7D%5Csum%7Bi%3D1%7D%5EN%5Csum%7Bj%3D1%7D%5EN%5Calphai%5Calpha_jy_iy_j%28x_i%5Ccdot%20x_j%29%2B%5Csum%7Bi%3D1%7D%5EN%5Calpha_i%0A&height=55&width=391)
拉格朗日对偶的重要作用是将w的计算提前并消除w,使得优化函数变为拉格朗日乘子的单一参数优化问题。
(2)求\mathop{\text{min}}_{w,b}L(w,b,\alpha)对\alpha的极大,即是对偶问题
%2B%5Csum%7Bi%3D1%7D%5EN%5Calpha_i%5C%5C%0A%26s.t.%5Cquad%20%5Csum%7Bi%3D1%7D%5EN%5Calphaiy_i%3D0%5C%5C%0A%26%5Cquad%5Cquad%5Cquad%5Calpha_i%5Cgeqslant%200%2C%5Cquad%20i%3D1%2C2%2C…%2CN%0A%5Cend%7Baligned%7D%0A#card=math&code=%5Cbegin%7Baligned%7D%0A%26%5Cmathop%7B%5Ctext%7Bmax%7D%7D%7B%5Calpha%7D%5C%20-%5Cfrac%7B1%7D%7B2%7D%5Csum%7Bi%3D1%7D%5EN%5Csum%7Bj%3D1%7D%5EN%5Calphai%5Calpha_jy_iy_j%28x_i%5Ccdot%20x_j%29%2B%5Csum%7Bi%3D1%7D%5EN%5Calphai%5C%5C%0A%26s.t.%5Cquad%20%5Csum%7Bi%3D1%7D%5EN%5Calpha_iy_i%3D0%5C%5C%0A%26%5Cquad%5Cquad%5Cquad%5Calpha_i%5Cgeqslant%200%2C%5Cquad%20i%3D1%2C2%2C…%2CN%0A%5Cend%7Baligned%7D%0A&height=131&width=319)
将上式从求极大换成求极小,就得到下面与之等价的对偶优化问题:
-%5Csum%7Bi%3D1%7D%5EN%5Calpha_i%5C%5C%0A%26s.t.%5Cquad%20%5Csum%7Bi%3D1%7D%5EN%5Calphaiy_i%3D0%5C%5C%0A%26%5Cquad%5Cquad%5Cquad%5Calpha_i%5Cgeqslant%200%2C%5Cquad%20i%3D1%2C2%2C…%2CN%0A%5Cend%7Baligned%7D%0A#card=math&code=%5Cbegin%7Baligned%7D%0A%26%5Cmathop%7B%5Ctext%7Bmin%7D%7D%7B%5Calpha%7D%5C%20%5C%20%5Cfrac%7B1%7D%7B2%7D%5Csum%7Bi%3D1%7D%5EN%5Csum%7Bj%3D1%7D%5EN%5Calphai%5Calpha_jy_iy_j%28x_i%5Ccdot%20x_j%29-%5Csum%7Bi%3D1%7D%5EN%5Calphai%5C%5C%0A%26s.t.%5Cquad%20%5Csum%7Bi%3D1%7D%5EN%5Calpha_iy_i%3D0%5C%5C%0A%26%5Cquad%5Cquad%5Cquad%5Calpha_i%5Cgeqslant%200%2C%5Cquad%20i%3D1%2C2%2C…%2CN%0A%5Cend%7Baligned%7D%0A&height=131&width=300)
考虑原始最优化问题和对偶最优化问题,原始问题满足Slater条件,使得对偶性为强对偶,所以存在w^, \alpha^, \beta是原始问题的解,\alpha^, \beta^*是对偶问题的解。这意味着求解原始问题可以转换为求解对偶问题。
以KKT条件由对偶解α求原始解w和b
对线性可分训练数据集,假设对偶最优化问题对\alpha的解为\alpha,\alpha_2)求得原始最优化问题对(w,b)的解w。有下面的定理。
定理:
设\alpha,\alpha_2)>0,并可按下式求得原始最优化问题的解w,b*:
%5C%5C%0A%5Cend%7Baligned%7D%0A#card=math&code=%5Cbegin%7Baligned%7D%0A%26w%5E%2A%3D%5Csum%7Bi%3D1%7D%5EN%5Calpha_i%5E%2Ay_ix_i%5C%5C%0A%26b%5E%2A%3Dy_j-%5Csum%7Bi%3D1%7D%5EN%5Calpha_i%5E%2Ay_i%28x_i%5Ccdot%20x_j%29%5C%5C%0A%5Cend%7Baligned%7D%0A&height=107&width=194)
上式的通俗理解:把\alpha理解为接触力,那接触力大于零了,不就说明正好到达约束边界了嘛,所以,此时的约束边界就是真正起作用的支持向量了。
证明:
这里先说一下KKT条件(不知道KKT条件的,建议先看凸优化一章):
在原函数和不等式约束是凸函数的基础上,如果满足Slater条件(存在x可以让所有不等式约束均小于零),那么强对偶性成立。则x,\beta,\alpha*满足下面的KKT条件:
%3D0%5C%5C%0A%26%5Cbigtriangledown%7B%5Calpha%7DL(x%5E%2C%5Calpha%5E%2C%5Cbeta%5E*)%3D0%5C%5C%0A%26%5Cbigtriangledown%7B%5Cbeta%7DL(x%5E%2C%5Calpha%5E%2C%5Cbeta%5E)%3D0%5C%5C%0A%26%5Calpha_i%5Eci(x_i%5E)%3D0%2C%20%5Cquad%20i%3D1%2C2%2C…%2Ck%5Cquad%20%5Ctext%7BKKT%7D%5C%5C%0A%26c_i(x%5E)%5Cleqslant%200%2C%5Cquad%20i%3D1%2C2%2C…%2Ck%5C%5C%0A%26%5Calpha_i%5E%5Cgeqslant%200%5Cquad%20i%3D1%2C2%2C…%2Ck%5C%5C%0A%26h_j(x%5E)%3D0%5Cquad%20i%3D1%2C2%2C…%2Cl%5C%5C%0A%5Cend%7Baligned%7D%0A#card=math&code=%5Cbegin%7Baligned%7D%0A%26%5Cbigtriangledown_xL%28x%5E%2A%2C%5Calpha%5E%2A%2C%5Cbeta%5E%2A%29%3D0%5C%5C%0A%26%5Cbigtriangledown%7B%5Calpha%7DL%28x%5E%2A%2C%5Calpha%5E%2A%2C%5Cbeta%5E%2A%29%3D0%5C%5C%0A%26%5Cbigtriangledown_%7B%5Cbeta%7DL%28x%5E%2A%2C%5Calpha%5E%2A%2C%5Cbeta%5E%2A%29%3D0%5C%5C%0A%26%5Calpha_i%5E%2Ac_i%28x_i%5E%2A%29%3D0%2C%20%5Cquad%20i%3D1%2C2%2C…%2Ck%5Cquad%20%5Ctext%7BKKT%7D%5C%5C%0A%26c_i%28x%5E%2A%29%5Cleqslant%200%2C%5Cquad%20i%3D1%2C2%2C…%2Ck%5C%5C%0A%26%5Calpha_i%5E%2A%5Cgeqslant%200%5Cquad%20i%3D1%2C2%2C…%2Ck%5C%5C%0A%26h_j%28x%5E%2A%29%3D0%5Cquad%20i%3D1%2C2%2C…%2Cl%5C%5C%0A%5Cend%7Baligned%7D%0A&height=160&width=276)
特别指出,上式中的第四行称为KKT的互补松弛条件,粗略地讲,互补松弛条件意味着在最优点处,除了第i个约束起作用的情况(ci(x都为零,也就是说,若\alphai)=0(即若接触力不为零,则弹簧振子接触边界)。
下面开始正式证明:
根据上面所述,KKT条件成立,即得
%3Dw%5E-%5Csum_%7Bi%3D1%7D%5EN%5Calpha_i%5Eyix_i%3D0%5C%5C%0A%26%5Cbigtriangledown%7Bb%7DL(w%5E%2Cb%5E%2C%5Calpha%5E)%3D-%5Csum_%7Bi%3D1%7D%5EN%5Calpha_iy_i%3D0%5C%5C%0A%26%5Calpha_i%5E(yi(w%5E%5Ccdot%20x_i%2Bb%5E)-1)%3D0%2C%20%5Cquad%20i%3D1%2C2%2C…%2CN%5C%5C%0A%26y_i(w%5Ex_i%2Bb%5E)-1%5Cgeqslant%200%2C%5Cquad%20i%3D1%2C2%2C…%2CN%5C%5C%0A%26%5Calpha_i%5E*%5Cgeqslant%200%2C%5Cquad%20i%3D1%2C2%2C…%2CN%5C%5C%0A%5Cend%7Baligned%7D%0A#card=math&code=%5Cbegin%7Baligned%7D%0A%26%5Cbigtriangledown_wL%28w%5E%2A%2Cb%5E%2A%2C%5Calpha%5E%2A%29%3Dw%5E%2A-%5Csum%7Bi%3D1%7D%5EN%5Calphai%5E%2Ay_ix_i%3D0%5C%5C%0A%26%5Cbigtriangledown%7Bb%7DL%28w%5E%2A%2Cb%5E%2A%2C%5Calpha%5E%2A%29%3D-%5Csum_%7Bi%3D1%7D%5EN%5Calpha_iy_i%3D0%5C%5C%0A%26%5Calpha_i%5E%2A%28y_i%28w%5E%2A%5Ccdot%20x_i%2Bb%5E%2A%29-1%29%3D0%2C%20%5Cquad%20i%3D1%2C2%2C…%2CN%5C%5C%0A%26y_i%28w%5E%2Ax_i%2Bb%5E%2A%29-1%5Cgeqslant%200%2C%5Cquad%20i%3D1%2C2%2C…%2CN%5C%5C%0A%26%5Calpha_i%5E%2A%5Cgeqslant%200%2C%5Cquad%20i%3D1%2C2%2C…%2CN%5C%5C%0A%5Cend%7Baligned%7D%0A&height=177&width=335)
由此得
其中至少有一个\alphaj均为0,则w^_为0,不是原始最优化问题的解,则产生矛盾),因此对这个j或这几个j有(根据上面的互补松弛条件,不等式变为等式,这说明处于边界上,所以其实下标j代表的点就位于边界上,也就是所谓的支持向量。所以注意观察上式,最优分界面的w的值其实并不是每一个点都在贡献的,而仅仅是由a不为零的点,也就是边界上的支撑点贡献的):
-1%3D0%0A#card=math&code=y_j%28w%5E%2A%5Ccdot%20x_j%2Bb%5E%2A%29-1%3D0%0A&height=21&width=171)
将上上式带入上式,并注意到y_j^2=1,即得
%5C%5C%0A%5Cend%7Baligned%7D%0A#card=math&code=%5Cbegin%7Baligned%7D%0Ab%5E%2A%26%3Dyj-w%5E%2Ax_j%5C%5C%0A%26%3Dy_j-%5Csum%7Bi%3D1%7D%5EN%5Calpha_i%5E%2Ay_i%28x_i%5Ccdot%20x_j%29%5C%5C%0A%5Cend%7Baligned%7D%0A&height=76&width=194)
故此定理得证,即有
%5C%5C%0A%5Cend%7Baligned%7D%0A#card=math&code=%5Cbegin%7Baligned%7D%0A%26w%5E%2A%3D%5Csum%7Bi%3D1%7D%5EN%5Calpha_i%5E%2Ay_ix_i%5C%5C%0A%26b%5E%2A%3Dy_j-%5Csum%7Bi%3D1%7D%5EN%5Calpha_i%5E%2Ay_i%28x_i%5Ccdot%20x_j%29%5C%5C%0A%5Cend%7Baligned%7D%0A&height=107&width=194)
由此定理可知,分离超平面可以写成
%2Bb%5E*%3D0%0A#card=math&code=%5Csum_%7Bi%3D1%7D%5EN%5Calpha_i%5E%2Ay_i%28x%5Ccdot%20x_i%29%2Bb%5E%2A%3D0%0A&height=53&width=176)
分类决策函数可以写成
%3D%5Ctext%7Bsign%7D%5Cleft(%20%5Csum%7Bi%3D1%7D%5EN%5Calpha_i%5Ey_i(x%5Ccdot%20x_i)%2Bb%5E%20%5Cright)%0A#card=math&code=f%28x%29%3D%5Ctext%7Bsign%7D%5Cleft%28%20%5Csum%7Bi%3D1%7D%5EN%5Calpha_i%5E%2Ay_i%28x%5Ccdot%20x_i%29%2Bb%5E%2A%20%5Cright%29%0A&height=54&width=258)
这就是说,分类决策函数只依赖与输入x和训练样本输入的内积。上式称为线性可分支持向量机的对偶形式。
线性可分支持向量机学习算法
综上所述,对于给定的线性可分训练数据集,可以首先求对偶问题(强对偶性,必须满足Slater条件,即弹簧振子系统能正常运转)
-%5Csum%7Bi%3D1%7D%5EN%5Calpha_i%5C%5C%0A%26s.t.%5Cquad%20%5Csum%7Bi%3D1%7D%5EN%5Calphaiy_i%3D0%5C%5C%0A%26%5Cquad%5Cquad%5Cquad%5Calpha_i%5Cgeqslant%200%2C%5Cquad%20i%3D1%2C2%2C…%2CN%0A%5Cend%7Baligned%7D%0A#card=math&code=%5Cbegin%7Baligned%7D%0A%26%5Cmathop%7B%5Ctext%7Bmin%7D%7D%7B%5Calpha%7D%5C%20%5C%20%5Cfrac%7B1%7D%7B2%7D%5Csum%7Bi%3D1%7D%5EN%5Csum%7Bj%3D1%7D%5EN%5Calphai%5Calpha_jy_iy_j%28x_i%5Ccdot%20x_j%29-%5Csum%7Bi%3D1%7D%5EN%5Calphai%5C%5C%0A%26s.t.%5Cquad%20%5Csum%7Bi%3D1%7D%5EN%5Calpha_iy_i%3D0%5C%5C%0A%26%5Cquad%5Cquad%5Cquad%5Calpha_i%5Cgeqslant%200%2C%5Cquad%20i%3D1%2C2%2C…%2CN%0A%5Cend%7Baligned%7D%0A&height=131&width=300)
的解\alpha^;再利用(*KKT条件)
%5C%5C%0A%5Cend%7Baligned%7D%0A#card=math&code=%5Cbegin%7Baligned%7D%0A%26w%5E%2A%3D%5Csum%7Bi%3D1%7D%5EN%5Calpha_i%5E%2Ay_ix_i%5C%5C%0A%26b%5E%2A%3Dy_j-%5Csum%7Bi%3D1%7D%5EN%5Calpha_i%5E%2Ay_i%28x_i%5Ccdot%20x_j%29%5C%5C%0A%5Cend%7Baligned%7D%0A&height=107&width=194)
求得原始问题的解w;从而得到分离超平面及分类决策函数。这种算法称为*线性可分支持向量机的对偶学习算法,是线性可分支持向量机的基本算法。
线性可分支持向量机学习算法:
输入:线性可分训练集T={(x_1,y_1),(x_2,y_2),…,(x_N,y_N)},其中,x_i\in X=R^n,y_i\in Y={-1,+1},i=1,2,…,N;
输出:分离超平面和分类决策函数。
(1)构造并求解约束最优化问题
-%5Csum%7Bi%3D1%7D%5EN%5Calpha_i%5C%5C%0A%26s.t.%5Cquad%20%5Csum%7Bi%3D1%7D%5EN%5Calphaiy_i%3D0%5C%5C%0A%26%5Cquad%5Cquad%5Cquad%5Calpha_i%5Cgeqslant%200%2C%5Cquad%20i%3D1%2C2%2C…%2CN%0A%5Cend%7Baligned%7D%0A#card=math&code=%5Cbegin%7Baligned%7D%0A%26%5Cmathop%7B%5Ctext%7Bmin%7D%7D%7B%5Calpha%7D%5C%20%5C%20%5Cfrac%7B1%7D%7B2%7D%5Csum%7Bi%3D1%7D%5EN%5Csum%7Bj%3D1%7D%5EN%5Calphai%5Calpha_jy_iy_j%28x_i%5Ccdot%20x_j%29-%5Csum%7Bi%3D1%7D%5EN%5Calphai%5C%5C%0A%26s.t.%5Cquad%20%5Csum%7Bi%3D1%7D%5EN%5Calpha_iy_i%3D0%5C%5C%0A%26%5Cquad%5Cquad%5Cquad%5Calpha_i%5Cgeqslant%200%2C%5Cquad%20i%3D1%2C2%2C…%2CN%0A%5Cend%7Baligned%7D%0A&height=131&width=300)
求得最优解\alpha,\alpha_2)^T。
(2)计算
并选择\alpha*>0,计算
%0A#card=math&code=b%5E%2A%3Dyj-%5Csum%7Bi%3D1%7D%5EN%5Calpha_i%5E%2Ay_i%28x_i%5Ccdot%20x_j%29%0A&height=53&width=189)
(3)求得分离超平面
分类决策函数:
%3D%5Ctext%7Bsign%7D(w%5E%5Ccdot%20x%2Bb%5E)%0A#card=math&code=f%28x%29%3D%5Ctext%7Bsign%7D%28w%5E%2A%5Ccdot%20x%2Bb%5E%2A%29%0A&height=20&width=173)
线性可分数据集的支持向量
在线性可分支持向量机中,由上面的(如下式)式子可知,**w只依赖于训练数据中对应于\alpha_i和b>0的实例点x_i称为支持向量。
%5C%5C%0A%5Cend%7Baligned%7D%0A#card=math&code=%5Cbegin%7Baligned%7D%0A%26w%5E%2A%3D%5Csum%7Bi%3D1%7D%5EN%5Calpha_i%5E%2Ay_ix_i%5C%5C%0A%26b%5E%2A%3Dy_j-%5Csum%7Bi%3D1%7D%5EN%5Calpha_i%5E%2Ay_i%28x_i%5Ccdot%20x_j%29%5C%5C%0A%5Cend%7Baligned%7D%0A&height=107&width=194)
根据这一定义,支持向量一定在分隔边界上,为什么?因为由KKT互补松弛条件可知,
-1)%3D0%2C%5Cquad%20i%3D1%2C2%2C…%2CN%0A#card=math&code=%5Calpha_i%5E%2A%28y_i%28w%5E%2A%5Ccdot%20x_i%2Bb%5E%2A%29-1%29%3D0%2C%5Cquad%20i%3D1%2C2%2C…%2CN%0A&height=21&width=330)
对应于\alpha_i^*>0的实例x_i,有
-1%3D0%0A#card=math&code=y_i%28w%5E%2A%5Ccdot%20x_i%2Bb%5E%2A%29-1%3D0%0A&height=20&width=169)
或
即x_i一定在间隔边界上,这里的支持向量的定义与前面给出的支持向量的定义是一致的。
对于线性可分问题,上述线性可分支持向量机的学习(硬间隔最大化)算法是完美的。但是,训练数据集线性可分是理想的清醒。在现实问题中,训练数据集往往是线性不可分的,即在样本中出现的噪点或特异点。此时,有更一般的学习算法。
参考资料
- 《统计学习方法》李航
本章的结构和大部分内容均参考此书对应章节。
