支持向量机

机器学习 专栏收录该内容
5 篇文章 1 订阅

支持向量机

点到平面的距离

平面的一般式方程

w1x1+w2x2+w3x3+b=wTx+b=0

其中 w=(w1,w2,w3) 是平面的法向量, b 是将平面平移到坐标原点所需距离(所以 b=0 时,平面过原点)

向量的模

给定一个向量 w=(w1,w2,w3) , 则 |w|=sqrt(w21+w22+w23)

向量的内积

给定两个向量 w=(w1,w2,w3) x=(x1,x2,x3) 则他们的内积是

wx=w1x1+w2x2+w3x3

点到平面的距离

!Alt text

最优间隔分类器与支持向量

函数间隔和几何间隔

令划分超平面的线性方程为

wTx+b=0

假设超平面能正确分类,当 wTx+b>0 的点对应 y=1 的数据点, wTx+b<0 的点对应 y=1 的点。
给定一个训练样本 (x(i),y(i)) 。我们定义函数间隔如下:
γ(i)=y(i)(wTx(i)+b)

实际上
γ(i)=wTx(i)+b

为了使函数间隔最大(更大的信心确定该例是正例还是反例),当 y(i)=1 时, wTx(i)+b 应该是个大正数,反之是个大负数。因此函数间隔代表了我们认为特征是正例还是反例的确信度。

继续考虑 w b,如果同时加大 w b,比如在 wTx(i)+b 前面乘个系数比如2,那么所有点的函数间隔都会增大二倍,这个对求解问题来说不应该有影响,因为我们要求解的是 wTx(i)+b=0 ,同时扩大 w b对结果是无影响的。这样,我们为了限制 w b,可能需要加入归一化条件

现在我们定义全局样本上的函数间隔

γ=mini=1,...,mγ(i)

说白了就是在训练样本上分类正例和负例确信度最小那个函数间隔。

样本空间任意点 x(i) 到超平面的几何间隔为

r(i)=wTx(i)+bw=γ(i)w

函数间隔归一化结果就是几何间隔。同时扩大 w b w 扩大几倍,w 就扩大几倍,结果无影响。
定义全局的几何间隔
r=mini=1,...,mr(i)

分类学习最基本的想法就是基于给定训练集在样本空间中找到一个划分超平面,将不同的类别样本分开。但满足这种要求的划分超平面有很多。
!Alt text

如何确定这个超平面

直观上讲,这个超平面应该是最适合分开两类数据的直线。而判定“最适合”的标准就是寻找一个超平面,使得离超平面比较近的点能有更大的间距。也就是我们不考虑所有的点都必须远离超平面,我们关心求得的超平面能够让所有点中离它最近的点具有最大间距。

形式化表示为

maxr,w,bγws.t.y(i)(wTx(i)+b)γ,i=1,2,...,m

前面说到同时扩大 w b 对结果没有影响,但我们最后要求的仍然是 w b 的确定值。因此,我们需要对 γ 做一些限制,以保证我们解是唯一的。这里为了简便我们取 γ=1 。这样的意义是将全局的函数间隔定义为1。即将离超平面最近的点的距离定义为 1w
即令

{wTxi+b+1,yi=+1wTxi+b1,yi=1

如下图所示
Alt text

上式使得等号成立的几个样本点被称为支持向量

两个异类支持向量到超平面距离之和为:

γ=2w

被称为 间隔

最大间隔划分超平面

欲找到最大间隔的划分超平面,我们得到改写后的目标函数为

maxw,b2ws.t.yi(wTxi+b)1,i=1,2,...,m

我们在优化时喜欢求最小值,将上式转化正等价的求最小值如下

minw,b12w2s.t.yi(wTxi+b)1,i=1,2,...,m

这就是支持向量机的 基本型

对偶问题

对偶问题转化

上式是一个凸二次规划问题,我们可以使用拉格朗日乘数法进行优化,对每个约束添加拉格朗日乘子 αi0 。则问题的拉格朗日函数可写为

L(w,b,α)=12w2+i=1mαi(1yi(wTxi+b))

式中 αi 是拉格朗日乘子.
我们令
θ(w)=maxαi0L(w,b,α)

如果约束条件都满足, θ(w) 的最优值就是 12w2 ,和目标函数一样。

因此我们可以直接求 θ(w) 的最小值,等价于求原目标函数。因此目标函数变成如下

minw,bθ(w)=minw,bmaxαi0L(w,b,α)=p

将求最大值和最小值交换位置,得到原始问题的对偶问题

maxαi0minw,bL(w,b,α)=d

满足

dp

在满足某些条件的情况下(满足KKT条件),这两者相等,这个时候就可以通过求解对偶问题来间接地求解原始问题。

对偶问题求解

分别求 L w,b 的导数,并令其为0,得如下结果

Lw=0w=i=1mαiyixiLb=0i=1mαiyi=0

带入 L(w,b,α) 得到 minw,bL(w,b,α)

minw,bL(w,b,α)=12i,j=1mαiαjyiyjxTixji,j=1mαiαjyiyjxTixjbi=1mαiyi+i=1mαi=i=1mαi12i,j=1mαiαjyiyjxTixj

然后求最大值,从上面的式子得到对偶问题

maxαi0i=1mαi12i,j=1mαiαjyiyjxTixjs.t.i=1mαiyi=0,αi0,i=1,2,...,m

上式是关于 α 的式子,如果能求出 α ,则通过 w=i=1mαiyixi 可以求出 w ,再根据前面函数距离等于1的假设求出 b

得到模型

f(x)=wTx+b=i=1mαiyixTix+b=i=1mαiyixi,x+b

对于新点 x 的预测,只需计算它与训练数据点的内积即可。且所有非支持向量所对应的系数 α 都等于零(这里为满足 maxαi0L(w,b,α) 最大化必须等于0),因此对于内积计算只针对支持向量即可。

拉格朗日乘子法和KKT条件

拉格朗日乘子法和KKT(Karush-Kuhn-Tucker)条件是求解约束优化问题的重要方法,在有等式约束时使用拉格朗日乘子法,在有不等约束时使用KKT条件。

前提:只有当目标函数为凸函数时,使用这两种方法才保证求得的是最优解。

最优化问题分类:
1). 如果目标函数和约束条件都为变量的线性函数, 称该最优化问题为线性规划;
2). 如果目标函数为变量的二次函数, 约束条件为变量的线性函数, 称该最优化问题为二次规划;
3). 如果目标函数或者约束条件为变量的非线性函数, 称该最优化问题为非线性规划.

拉格朗日乘子法–等式约束问题

对于等式约束问题

x=argminxf(x)subject to hi(x)=0,i=1,..,m

利用拉格朗日乘子法,原问题转换为
x=argmin xL(x,λ)=argmin xf(x)+mi=1λihi(x),

其中 λi0 ,称为拉格朗日乘子,这样我们有了 k+m 个变量。

通过对 x 求偏导,找最优解,得到最优的条件

xL(x,λ)=xf(x)+iλixhi(x)=0

我们得到了 k 个方程式。进一步对拉格朗日乘子求偏导,得到另外 m 个方程式

L(x,λ)λi=hi(x)=0,

这里通过强制偏导数为零,我们同样可以限制可能的解满足原始约束条件。现在变量数与方程数相同,使得方程解唯一。

等式约束问题的合理性

举个二维最优化的例子

minf(x,y)s.t.g(x,y)=c

我们可以画图来辅助思考
!Alt text这里写图片描述

蓝线是 f(x,y) 的等高线。箭头表示斜率,和等高线的法线平行。从梯度的方向上来看,显然有 d1>d2
绿线标出的是约束 g(x,y)=c 的点的轨迹,也就是说,只要正好落在这条绿线上的点才可能是满足要求的点。

最小值点应该在哪里呢?显然应该是在 f(x,y) 的等高线正好和约束线相切的位置,因为如果只是相交意味着肯定还存在其它的等高线在该条等高线的内部或者外部,使得新的等高线与目标函数的交点的值更大或者更小,只有到等高线与目标函数的曲线相切的时候,可能取得最优值。

很容易看出来,要想让目标函数 f(x,y) 的等高线和约束相切,则他们切点的梯度一定在一条直线上。即在最优解处, f g 的斜率平行。
因此,我们通过观察可以得到优化取到最小值时满足

[f(x,y)+λ(g(x,y)c)]=0λ0

这个结果正好等价于方程 F(x,y)=f(x,y)+λ(g(x,y)c) x 求偏导。方程 F(x,y) 在达到极值时与 f(x,y) 相等,因为 F(x,y) 达到极值时 g(x,y)c 总等于零。

拉格朗日乘子法–不等式约束问题

当同时存在不等式约束时,新的问题描述为

x=argminxf(x)

subject to hi(x)=0,i=1,..,msubject to gi(x)0,i=1,..,n

此时不等式约束下的最优化问题描述为

x=argmin xL(x,λ,μ)=argmin xf(x)+mi=1λihi(x)+ni=1μigi(x)

其中 L(x,λ,μ) 为拉格朗日函数, λi 对应等式约束系数, μk 对应不等式约束系数。现在总共有 k+m+n 个变量。

跟等式约束一样,我们想通过对 x λ 求偏导,得到 k+m 个方程式。接下来的问题是怎么通过不等式约束得到更多的 n 个方程式。

对于不等式约束,当极值点存在 gi(x)<0 区域,那么这个约束不影响极值的取值,等同于没有约束。因此它的系数 μi 可以设置为0。另一方面,如果极值在约束边界上,此时, gi(x)=0 。如下图所示
Alt text
综合可得到

μigi(x)=0

为问题的必要条件。至此,我们得到了 n 个方程式通过不等式约束。系数 λi 可以为任意值,没有符号的限制。然而 μk0 ,下图解释说明。
Alt text

假设 x gi(x)=0 区域内,那么 μi 可以不为零。

{x^*} = \underset{{x}}{\textrm{argmin }} f({x})  + \mu_i g_i({x})
\\0 = \nabla f({x})  + \mu_i \nabla g_i({x})
\\\begin{equation}\mu_i = -\frac{\nabla f({x})}{\nabla g_i({x})}\label{eq:signofmu}\end{equation}

x 点处, f(x) gi(x) 的梯度方向像相反。所以 μi0

KKT条件

对于一个标准形式的最优化数学模型

x=argminxf(x)

subject to hi(x)=0,i=1,..,msubject to gi(x)0,i=1,..,n

KKT条件是一个非线性规划问题能有最优化解法的必要和充分条件。KKT条件就是指上面最优化问题中的最小点 x 必须满足下面的条件:
Stationarity:Primalfeasibility:Dualfeasibility:Complementaryslackness:0=f(x)+iλihi(x)+iμigi(x)hi(x)=0,i=1,..,mgi(x)0,i=1,..,nμi0,i=1,..,nμigi(x)=0,i=1,..,m

不等式约束问题的合理性

为什么可对偶问题求解

L(x,μ)=f(x)+k=1qμkgk(x)μk0gk(x)0

μk0gk(x)0}=>μg(x)0

\therefore \begin{equation}\max_{\mu}L(x,\mu)=f(x)\label{a}\end{equation}

\therefore\begin{equation}\min_{x}f(x)=\min_{x}\max_{\mu}L(x,\mu)\label{firsthalf}\end{equation}

maxμminxL(x,μ)=maxμ[minxf(x)+minxμg(x)]=maxμminxf(x)+maxμminxμg(x)=minxf(x)+maxμminxμg(x)

μk0gk(x)0}=>minxμg(x)={0ifμ=0org(x)=0ifμ>0andg(x)<0

maxμminxμg(x)=0

此时 μ=0org(x)=0
\begin{equation}\therefore \max_{\mu}\min_{x}L(x,\mu)=\min_{x}f(x)+\max_{\mu}\min_{x}\mu{g(x)}=\min_{x}f(x)\label{secondhalf}\end{equation}

此时 μ=0org(x)=0
联合 (???),(???) 我们得到 >
minxmaxμL(x,μ)=maxμminxL(x,μ)

L(x,μ)=f(x)+k=1qμkgk(x)μk0gk(x)0
=>
minxmaxμL(x,μ)=maxμminxL(x,μ)=minxf(x)

我们把 maxμminxL(x,μ) 称为原问题 minxmaxμL(x,μ) 的对偶问题,上式表明当满足一定条件时原问题、对偶的解、以及 minxf(x) 是相同的,且在最优解 x μ=0org(x)=0 。把 x 代入 (???) maxμL(x,μ)=f(x) ,由 (???) maxμminxL(x,μ)=f(x) ,所以 L(x,μ)=minxL(x,μ) ,这说明 x 也是 L(x,μ) 的极值点,即 L(x,μ)x|x=x=0

最后总结一下:

L(x,μ)=f(x)+k=1qμkgk(x)μk0gk(x)0
=>
minxmaxμL(x,μ)=maxμminxL(x,μ)=minxf(x)=f(x)μkgk(x)=0L(x,μ)x|x=x=0

KKT条件是拉格朗日乘子法的泛化,如果我们把等式约束和不等式约束一并纳入进来则表现为:

L(x,λ,μ)=f(x)+i=1nλihi(x)+k=1qμkgk(x)λi0hi(x)=0μk0gk(x)0

=>
minxmaxμL(x,λ,μ)=maxμminxL(x,λ,μ)=minxf(x)=f(x)μkgk(x)=0L(x,λ,μ)x|x=x=0

L(x,λ,μ)x|x=x=0 表明 f(x) 在极值点 x 处的梯度是各个 hi(x) gk(x) 梯度的线性组合。

核函数

支持向量机是建立在统计学习理论基础之上的新一代机器学习算法,支持向量机的优势主要体现在解决线性不可分问题,它通过引入核函数,巧妙地解决了在高维空间中的内积运算,从而很好地解决了非线性分类问题。

特征空间的隐式映射

前面的讨论中,我们假设训练样本是线性可分的,即存在一个划分超平面能将训练样本正确分类。事实上大部分时候数据并不是线性可分的,这个时候满足这样条件的超平面就根本不存在。

对于非线性的情况,SVM 的处理方法是选择一个核函数 κ(,) ,通过将数据从原始空间映射到一个更高维的特征空间,使得样本在特征空间内线性可分,来解决在原始空间中线性不可分的问题。

具体来说,在线性不可分的情况下,支持向量机首先在低维空间中完成计算,然后通过核函数将输入空间映射到高维特征空间,最终在高维特征空间中构造出最优分离超平面,从而把平面上本身不好分的非线性数据分开。

ϕ(x) 表示将 x 映射后的特征向量,于是,在特征空间中的划分超平面所对应的模型可表示为

f(x)=wTϕ(x)+b

这意味着建立非线性学习器分为两步:首先使用一个非线性映射将数据变换到一个特征空间。然后在特征空间使用线性学习器分类

类似于原始空间线性可分的问题,可以得到特征空间的SVM基本型

minw,b12w2s.t.yi(wTϕ(xi)+b)1,i=1,2,...,m

其对偶问题是

maxαi0i=1mαi12i,j=1mαiαjyiyjϕ(xi)Tϕ(xj)s.t.i=1mαiyi=0,αi0,i=1,2,...,m

求解上式涉及到计算 ϕ(xi)Tϕ(xj) ,这是样本在特征空间的内积。

ϕ(xi)Tϕ(xj) 可以有两种方法
1、先找到这种映射,然后将输入空间中的样本映射到新的空间中,最后在新空间中去求内积
2、或者是找到某种方法,它不需要显式的将输入空间中的样本映射到新的空间中,而能够在输入空间中直接计算出内积,这就是传说中的核函数方法

由于特征空间维数可能很高,直接计算通常困难。为了避开障碍,可以设想一个这样的函数

k(xi,xj)=ϕ(xi),ϕ(xj)=ϕ(xi)Tϕ(xj)

即特征空间的内积等于在原始空间中通过函数 k(.,.) 计算而得。这样就不必直接去计算高维特征空间的内积。于是上式可重写为

maxαi0i=1mαi12i,j=1mαiαjyiyjk(xi,xj)s.t.i=1mαiyi=0,αi0,i=1,2,...,m

求解后可得到

f(x)=wTϕ(x)+b=i=1mαiyiϕ(xi)Tϕ(x)+b=i=1mαiyik(xi,x)+b

这里的函数 k(.,.) 就是核函数。

定义:核是一个函数 K ,对于所有的 x1,x2X 满足, K<x1,x2>=<Φ(x1),Φ(x2)> ,这里的 Φ 为从 X 到内积特征空间 F 的映射。

核函数有效性判定

先看一个例子,假设 x z 都是 n维的 K(x,z)=(xTz)2 展开后,得

k(x,z)=(xTz)2=(i=1nxizi)j=1nxjzj=i=1nj=1n(xixj)(zizj)=ϕ(x)Tϕ(z)

现在看一下映射函数(n=3时),根据上面的公式,得到
ϕ(x)=[x1x1x1x2x1x3x2x1x2x2x2x3x3x1x3x2x3x3]T

这个时候发现我们可以只计算原始特征 x z 内积的平方,就等价与计算映射后特征的内积。也就是说核函数 K(x,z)=(xTz)2 只能在选择这样的 ϕ 作为映射函数时才能够等价于映射后特征的内积。

我们再看一个核函数

k(x,z)=(xTz+c)2=(i=1nxizi+c)j=1nxjzj+c=i=1nj=1n(xixj)(zizj)+(i=1n(2cxi)(2czi))+c2

对应的映射函数(n=3时)是

ϕ(x)=[x1x1x1x2x1x3x2x1x2x2x2x3x3x1x3x2x3x32cx12cx22cx3c]T

半正定矩阵

实对称矩阵:如果有 n 阶矩阵 A,其各个元素都为实数,矩阵A的转置等于其本身,则称 A 为实对称矩阵。
主要性质
1.实对称矩阵 A 的不同特征值对应的特征向量是正交的。
2.实对称矩阵 A 的特征值都是实数,特征向量都是实向量。
3.n 阶实对称矩阵 A 必可对角化,且相似对角阵上的元素即为矩阵本身特征值。

埃尔米特矩阵n 阶复方阵 A 的对称单元互为共轭,即 A 的共轭转置矩阵等于它本身,则 A 是埃尔米特矩阵。显然它是实对称阵的推广。

半正定矩阵:一个 n×n 的埃尔米特矩阵 M 是正定的当且仅当对于每个非零的复向量 z,都有 zMz>0 ,则称 M 为正定矩阵,其中z 表示 z 的共轭转置。当zMz>0弱化为 zMz0 时,称 M 是半正定矩阵。
判定:所有主子式大于或等于零

问题描述与分析

给定一个函数 K,我们能否使用 K 来替代计算 ϕ(x)Tϕ(z),也就说,是否能够找出一个 ϕ ,使得对于所有的 x z,都有 k(x,z)=ϕ(x)Tϕ(z) 。 比如给出了 k(x,z)=(xTz)2 ,是否能够认为 K 是一个有效的核函数。

下面来解决这个问题,给定 m 个训练样本 {x(1),x(2),...,x(m)} ,每一个 x(i) 对应一个特征向量。那么,我们可以将任意两个 x(i) x(j) 带入 K 中,计算得到 Kij=K(x(i),x(j)) i 可以从1到 m j 可以从1到 m,这样可以计算出 mm 的核函数矩阵。我们用 K 来表示。

有效地核函数定义

如果假设 K 是有效地核函数,那么根据核函数定义

Kij=K(x(i),x(j))=ϕ(xi)Tϕ(xj)=ϕ(xj)Tϕ(xi)=K(x(j),x(i))=Kji

可见,矩阵K应该是个对称阵。
再者,对于任意向量 z ,得

zTKz=i=1nj=1n(ziKijzj)=i=1nj=1n(ziϕ