登录
原创

SVM 支持向量机

专栏苏州谷歌开发者社区
发布于 2020-11-20 阅读 231
  • 人工智能
  • 机器学习
原创

支持向量机(support vector machines, SVM)是一种二分分类模型。他的基本模型是定义在特征空间上的间隔最大的线性分类器,间隔最大使它有别与别的感知机;支持向量机还包括核技巧,这使它成为实质上的非线性分类器。支持向量机的学习策略就是间隔最大化,可形式化为一个求解凸二次规划(convex quadratic programming)的问题,也等价于正则化的合页损失函数最小化问题。支持向量机的学习算法是求解凸二次规划的最优算法。

注:凸二次规划中,局部最优解就是全局最优解。

支持向量机学习方法包含构建由简至繁的模型:线性可分支持向量机(linear support vector machine in linear separable case)、线性支持向量机(linear support vector machine)及非线性支持向量机(non-liner support vector machine).

简单模型是复杂模型的基础和特殊情况。当训练数据线性可分时,通过硬间隔最大化(hard margin maximization),学习一个线性分类器,即线性可分支持向量机,又成为硬间隔支持向量机;

当训练数据近似线性可分时,通过软间隔最大化(soft margin maximization),也学习一个线性分类,即线性支持向量机,又称为软间隔支持向量机;

当训练模型线性不可分是,通过核技巧(kernel trick)及软间隔最大化,学习分线性支持向量机。

当输入空间为欧式空间或离散集合、特征空间为希伯特空间时,核函数(kernel function)表示将输入从输入空间映射到特征空间得到的 特征向量之间的内积 。通过核函数可以学习非线性支持向量机,等价于隐式地在高维特征空间中学习线性支持向量机。这样的方法成为核技巧,核技巧(kernel method)是比支持向量机更为一般的机器学习方法。

线性可分支持向量机与硬间隔最大

线性可分支持向量机

一般地,当训练数据集线性可分时,存在无穷个分离超平面可将两类数据正确分开。感知机利用误分类最小策略,求得分离超平面,不过这时的解有无穷多个。线性可分支持向量机利用间隔最大化求最优分离超平面,这时,解是唯一的。

定义 (线性可分支持向量机) 给定线性可分训练数据集,通过间隔最大化或等价地求解相应的凸二次规划问题学习得到的分离超平面为

<math><semantics><mrow><msup><mi>w</mi><mo>∗</mo></msup><mi>x</mi><mo>+</mo><msup><mi>b</mi><mo>∗</mo></msup><mo>=</mo><mn>0</mn><mspace width="2em"></mspace><mo>(</mo><mn>7</mn><mi mathvariant="normal">.</mi><mn>1</mn><mo>)</mo></mrow><annotation encoding="application/x-tex">w^*x+b^*=0 \qquad (7.1) </annotation></semantics></math>wx+b=0(7.1)

以及相应的分类决策函数

<math><semantics><mrow><mi>f</mi><mo>(</mo><mi>x</mi><mo>)</mo><mo>=</mo><mi>s</mi><mi>i</mi><mi>g</mi><mi>n</mi><mo>(</mo><msup><mi>w</mi><mo>∗</mo></msup><mi>x</mi><mo>+</mo><msup><mi>b</mi><mo>∗</mo></msup><mo>)</mo><mspace width="2em"></mspace><mo>(</mo><mn>7</mn><mi mathvariant="normal">.</mi><mn>2</mn><mo>)</mo></mrow><annotation encoding="application/x-tex">f(x)=sign(w^*x+b^*) \qquad (7.2) </annotation></semantics></math>f(x)=sign(wx+b)(7.2)

称为线性可分支持向量机

函数间隔和几何间隔

一般来说,一个点距离分的远近可以表示分类预测的确信程度。在超平面$w \centerdot x +b =0 确定的情况下,|w \centerdot x +b|能够相对地表示点x距离超平面的远近。而w \centerdot x +b 的符号与类标记y的符号是否一致能够表示分类是否正确。所以可以用量y(w \centerdot x+ b)$来表示分类的正确性与确信度,这就是函数间隔(functional margin)的概念。

函数间隔

对于给定的训练数据集T和超平面(w,b)定义超平面(w,b)关于样本点<math><semantics><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo separator="true">,</mo><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>)</mo></mrow><annotation encoding="application/x-tex">(x_{i},y_{i})</annotation></semantics></math>(xi,yi)的函数间隔为

<math><semantics><mrow><mover accent="true"><mrow><msub><mi>γ</mi><mrow><mi>i</mi></mrow></msub></mrow><mo>^</mo></mover><mo>=</mo><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>(</mo><mi>w</mi><mo>⋅</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo>+</mo><mi>b</mi><mo>)</mo></mrow><annotation encoding="application/x-tex">\hat{ \gamma_{i}}=y_{i}(w \cdot x_{i} +b) </annotation></semantics></math>γi^=yi(wxi+b)

定义超平面(w,b)关于训练集T的函数间隔为超平面(w,b),定义超平面(w,b)关于T中所有样本点<math><semantics><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo separator="true">,</mo><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>)</mo></mrow><annotation encoding="application/x-tex">(x_{i},y_{i})</annotation></semantics></math>(xi,yi)的函数间隔之最小值,即

<math><semantics><mrow><mover accent="true"><mrow><mi>γ</mi><mo>=</mo><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>(</mo><mi>w</mi><mo>⋅</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo>+</mo><mi>b</mi><mo>)</mo></mrow><mo>^</mo></mover></mrow><annotation encoding="application/x-tex">\hat{\gamma=y_{i}(w \cdot x_{i}+b)} </annotation></semantics></math>γ=yi(wxi+b)^

定义超平面<math><semantics><mrow><mo>(</mo><mi>w</mi><mo separator="true">,</mo><mi>b</mi><mo>)</mo></mrow><annotation encoding="application/x-tex">(w,b)</annotation></semantics></math>(w,b)关于样本点<math><semantics><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo separator="true">,</mo><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>)</mo></mrow><annotation encoding="application/x-tex">(x_{i},y_{i})</annotation></semantics></math>(xi,yi)的函数间隔为

\hat{\gamma}=\min_{i=1,\dots,N} \hat{\gamma_{i}}

函数间隔可以表示为分类预测的正确性及确信度。但是选择分离平面时,只有函数间隔还不够。因为只要成比例地改变w和b,超平面没有变化,函数间隔却城北变化。所以我们可以对分离超平面的法向量w加某些约束,如规范化,<math><semantics><mrow><mi mathvariant="normal">∣</mi><mi mathvariant="normal">∣</mi><mi>w</mi><mi mathvariant="normal">∣</mi><mi mathvariant="normal">∣</mi><mo>=</mo><mn>1</mn></mrow><annotation encoding="application/x-tex">||w||=1</annotation></semantics></math>w=1,使得间隔是确定的。这时函数间隔成为几何间隔(geometric margin)。

几何间隔

对于给定的训练集T和超平面(w,b),定义超平面(w,b)关于样本点<math><semantics><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo separator="true">,</mo><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>)</mo></mrow><annotation encoding="application/x-tex">(x_{i},y_{i})</annotation></semantics></math>(xi,yi)的几何间隔为

<math><semantics><mrow><msub><mi>γ</mi><mrow><mi>i</mi></mrow></msub><mo>=</mo><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>(</mo><mfrac><mrow><mi>w</mi></mrow><mrow><mi mathvariant="normal">∣</mi><mi mathvariant="normal">∣</mi><mi>w</mi><mi mathvariant="normal">∣</mi><mi mathvariant="normal">∣</mi></mrow></mfrac><mo>⋅</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo>+</mo><mfrac><mrow><mi>b</mi></mrow><mrow><mi mathvariant="normal">∣</mi><mi mathvariant="normal">∣</mi><mi>w</mi><mi mathvariant="normal">∣</mi><mi mathvariant="normal">∣</mi></mrow></mfrac><mo>)</mo></mrow><annotation encoding="application/x-tex">\gamma_{i} = y_{i}(\frac{w}{||w|| }\cdot x_{i}+\frac{b}{||w||}) </annotation></semantics></math>γi=yi(wwxi+wb)

定义超平面(w,b)关于训练集T的几何间隔为超平面(w,b)关于T中所有样本点<math><semantics><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo separator="true">,</mo><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>)</mo></mrow><annotation encoding="application/x-tex">(x_{i},y_{i})</annotation></semantics></math>(xi,yi)的几何间隔最小值

\gamma=\min_{i=1,\dots,N} \gamma_{i}

间隔最大化

支持向量机学习的基本思想就是求解能够正确划分训练集并且几何间隔最大的分离超平面。

间隔最大化意味着不仅将正负实例点分开,而且对最难分的实例点(离超平面最近的点)也有足够的确信度将它们分开。这样的超平面应该对未知的新实例有很好的分类预测能力。

间隔最大分离超平面

可以将求解几何间隔最大化问题表示为以下的约束最优化问题:

<math><semantics><mrow><msub><mi>max</mi><mrow><mi>w</mi><mo separator="true">,</mo><mi>b</mi></mrow></msub><mspace width="2em"></mspace><mi>γ</mi></mrow><annotation encoding="application/x-tex">\max_{w,b} \qquad \gamma </annotation></semantics></math>w,bmaxγ

s.t. \qquad y_{i}\left(\frac{w}{||w||} \cdot x_{i} + \frac{b}{||w||}\right) \ge \gamma , \quad i = 1,2,\dots,N

即我们希望最大化超平面(w,b)关于训练集的几何间隔<math><semantics><mrow><mi>γ</mi></mrow><annotation encoding="application/x-tex">\gamma</annotation></semantics></math>γ,约束条件表示的是超平面(w,b)关于每个训练样本点的几何间隔至少是<math><semantics><mrow><mi>γ</mi></mrow><annotation encoding="application/x-tex">\gamma</annotation></semantics></math>γ

考虑几何间隔和函数间隔的关系式,可将这个问题改写为

<math><semantics><mrow><msub><mi>max</mi><mrow><mi>w</mi><mo separator="true">,</mo><mi>b</mi></mrow></msub><mspace width="2em"></mspace><mfrac><mrow><mover accent="true"><mrow><mi>γ</mi></mrow><mo>^</mo></mover></mrow><mrow><mi mathvariant="normal">∣</mi><mi mathvariant="normal">∣</mi><mi>w</mi><mi mathvariant="normal">∣</mi><mi mathvariant="normal">∣</mi></mrow></mfrac></mrow><annotation encoding="application/x-tex">\max_{w,b} \qquad \frac{\hat{\gamma}}{||w||} </annotation></semantics></math>w,bmaxwγ^

s.t. \qquad y_{i}(w \cdot x_{i} + b) \ge \hat{\gamma} , \quad i = 1,2,\dots,N

函数间隔<math><semantics><mrow><mi>γ</mi></mrow><annotation encoding="application/x-tex">\gamma</annotation></semantics></math>γ的取值并不影响最优化问题的解,也就是说,可以产生一个等价的最优化问题。这样,可以取<math><semantics><mrow><mover accent="true"><mrow><mi>γ</mi></mrow><mo>^</mo></mover><mo>=</mo><mn>1</mn></mrow><annotation encoding="application/x-tex">\hat{\gamma}=1</annotation></semantics></math>γ^=1,于是就得到线性可分支持向量机学习的最优化问题(最大化<math><semantics><mrow><mfrac><mrow><mn>1</mn></mrow><mrow><mi mathvariant="normal">∣</mi><mi mathvariant="normal">∣</mi><mi>w</mi><mi mathvariant="normal">∣</mi><mi mathvariant="normal">∣</mi></mrow></mfrac></mrow><annotation encoding="application/x-tex">\frac{1}{||w||}</annotation></semantics></math>w1和最小化<math><semantics><mrow><mfrac><mrow><mn>1</mn></mrow><mrow><mn>2</mn></mrow></mfrac><mi mathvariant="normal">∣</mi><mi mathvariant="normal">∣</mi><mi>w</mi><mi mathvariant="normal">∣</mi><msup><mi mathvariant="normal">∣</mi><mn>2</mn></msup></mrow><annotation encoding="application/x-tex">\frac{1}{2}||w||^2</annotation></semantics></math>21w2是等价的):

<math><semantics><mrow><msub><mi>min</mi><mrow><mi>w</mi><mo separator="true">,</mo><mi>b</mi></mrow></msub><mspace width="2em"></mspace><mfrac><mrow><mn>1</mn></mrow><mrow><mn>2</mn></mrow></mfrac><mi mathvariant="normal">∣</mi><mi mathvariant="normal">∣</mi><mi>w</mi><mi mathvariant="normal">∣</mi><msup><mi mathvariant="normal">∣</mi><mn>2</mn></msup><mspace width="2em"></mspace></mrow><annotation encoding="application/x-tex">\min_{w,b} \qquad \frac{1}{2}||w||^2 \qquad </annotation></semantics></math>w,bmin21w2

s.t. \qquad y_{i}(w \cdot x_{i} + b)-1 \ge 0 , \quad i = 1,2,\dots,N \qquad (7.14)

这是一个凸二次规划(convex quadratic programming)问题

凸优化问题是指优化问题

<math><semantics><mrow><msub><mi>min</mi><mrow><mi>w</mi></mrow></msub><mi>f</mi><mo>(</mo><mi>w</mi><mo>)</mo></mrow><annotation encoding="application/x-tex">\min_{w} f(w) </annotation></semantics></math>wminf(w)

<math><semantics><mrow><mi>s</mi><mi mathvariant="normal">.</mi><mi>t</mi><mi mathvariant="normal">.</mi><mspace width="2em"></mspace><msub><mi>g</mi><mrow><mi>i</mi></mrow></msub><mo>≤</mo><mn>0</mn><mo separator="true">,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo separator="true">,</mo><mn>2</mn><mo separator="true">,</mo><mo>⋯</mo><mo separator="true">,</mo><mi>k</mi></mrow><annotation encoding="application/x-tex">s.t. \qquad g_{i} \le 0, i = 1,2,\cdots , k </annotation></semantics></math>s.t.gi0,i=1,2,,k

<math><semantics><mrow><mspace width="2em"></mspace><msub><mi>h</mi><mrow><mi>i</mi></mrow></msub><mo>(</mo><mi>w</mi><mo>)</mo><mo>=</mo><mn>0</mn><mo separator="true">,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo separator="true">,</mo><mn>2</mn><mo separator="true">,</mo><mo>⋯</mo><mo separator="true">,</mo><mi>l</mi></mrow><annotation encoding="application/x-tex">\qquad h_{i}(w)=0 , i=1,2,\cdots,l </annotation></semantics></math>hi(w)=0,i=1,2,,l

其中,目标函数<math><semantics><mrow><mi>f</mi><mo>(</mo><mi>w</mi><mo>)</mo></mrow><annotation encoding="application/x-tex">f(w)</annotation></semantics></math>f(w)和约束函数<math><semantics><mrow><msub><mi>g</mi><mrow><mi>i</mi></mrow></msub><mo>(</mo><mi>w</mi><mo>)</mo></mrow><annotation encoding="application/x-tex">g_{i}(w)</annotation></semantics></math>gi(w)都是<math><semantics><mrow><msup><mi>R</mi><mrow><mi>n</mi></mrow></msup></mrow><annotation encoding="application/x-tex">R^{n}</annotation></semantics></math>Rn上的连续可微的凸函数,约束函数<math><semantics><mrow><msub><mi>h</mi><mrow><mi>i</mi></mrow></msub><mo>(</mo><mi>w</mi><mo>)</mo></mrow><annotation encoding="application/x-tex">h_{i}(w)</annotation></semantics></math>hi(w)<math><semantics><mrow><msup><mi>R</mi><mi>n</mi></msup></mrow><annotation encoding="application/x-tex">R^n</annotation></semantics></math>Rn上的仿射函数。

<math><semantics><mrow><mi>f</mi><mo>(</mo><mi>x</mi><mo>)</mo></mrow><annotation encoding="application/x-tex">f(x)</annotation></semantics></math>f(x)称为仿射函数,如果它满足<math><semantics><mrow><mi>f</mi><mo>(</mo><mi>x</mi><mo>)</mo><mo>=</mo><mi>a</mi><mo>⋅</mo><mi>x</mi><mo>+</mo><mi>b</mi><mo separator="true">,</mo><mi>a</mi><mo>∈</mo><msup><mi>R</mi><mi>n</mi></msup><mo separator="true">,</mo><mi>b</mi><mo>∈</mo><mi>R</mi><mo separator="true">,</mo><mi>x</mi><mo>∈</mo><msup><mi>R</mi><mi>n</mi></msup><mi mathvariant="normal">.</mi></mrow><annotation encoding="application/x-tex">f(x)=a \cdot x+b ,a \in R^n, b \in R, x \in R^n.</annotation></semantics></math>f(x)=ax+b,aRn,bR,xRn.

当目标函数<math><semantics><mrow><mi>f</mi><mo>(</mo><mi>x</mi><mo>)</mo></mrow><annotation encoding="application/x-tex">f(x)</annotation></semantics></math>f(x)是二次函数且约束函数<math><semantics><mrow><msub><mi>g</mi><mrow><mi>i</mi></mrow></msub><mo>(</mo><mi>w</mi><mo>)</mo></mrow><annotation encoding="application/x-tex">g_{i}(w)</annotation></semantics></math>gi(w)是仿射函数时,上述凸最优化问题成为凸二次规划问题。

最大间隔分离超平面存在的唯一性

线性可分数据集的最大间隔分离超平面是存在且唯一的。
(1)存在性
(2)唯一性

证明略,详见《统计学习方法》

支持向量和间隔边界

在线性可分情况下,训练数据集的样本点中与分离超平面距离最近的样本点的实例称为支持向量。支持向量是使约束条件式(7.14)等号成立的点。

s.t. \qquad y_{i}(w \cdot x_{i} + b)-1 \ge 0 , \quad i = 1,2,\dots,N \qquad (7.14)

<math><semantics><mrow><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>(</mo><mi>w</mi><mo>⋅</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo>+</mo><mi>b</mi><mo>)</mo><mo>−</mo><mn>1</mn><mo>=</mo><mn>0</mn></mrow><annotation encoding="application/x-tex">y_{i}(w \cdot x_{i} + b)-1 = 0 </annotation></semantics></math>yi(wxi+b)1=0

<math><semantics><mrow><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>=</mo><mo>+</mo><mn>1</mn></mrow><annotation encoding="application/x-tex">y_{i}=+1</annotation></semantics></math>yi=+1的正例点,支持向量在超平面

<math><semantics><mrow><msub><mi>H</mi><mrow><mn>1</mn></mrow></msub><mo>:</mo><mi>w</mi><mo>⋅</mo><mi>x</mi><mo>+</mo><mi>b</mi><mo>=</mo><mn>1</mn></mrow><annotation encoding="application/x-tex">H_{1}:w \cdot x + b =1 </annotation></semantics></math>H1:wx+b=1

上,对<math><semantics><mrow><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>=</mo><mo>−</mo><mn>1</mn></mrow><annotation encoding="application/x-tex">y_{i}=-1</annotation></semantics></math>yi=1的负实例点,支持向量在超平面

<math><semantics><mrow><msub><mi>H</mi><mrow><mn>2</mn></mrow></msub><mo>:</mo><mi>w</mi><mo>⋅</mo><mi>x</mi><mo>+</mo><mi>b</mi><mo>=</mo><mo>−</mo><mn>1</mn></mrow><annotation encoding="application/x-tex">H_{2}:w \cdot x + b =-1 </annotation></semantics></math>H2:wx+b=1

上。

注意到<math><semantics><mrow><msub><mi>H</mi><mrow><mn>1</mn></mrow></msub></mrow><annotation encoding="application/x-tex">H_{1}</annotation></semantics></math>H1<math><semantics><mrow><msub><mi>H</mi><mrow><mn>2</mn></mrow></msub></mrow><annotation encoding="application/x-tex">H_{2}</annotation></semantics></math>H2平行,并且没有实例点落在它们中间。在<math><semantics><mrow><msub><mi>H</mi><mrow><mn>1</mn></mrow></msub></mrow><annotation encoding="application/x-tex">H_{1}</annotation></semantics></math>H1<math><semantics><mrow><msub><mi>H</mi><mrow><mn>2</mn></mrow></msub></mrow><annotation encoding="application/x-tex">H_{2}</annotation></semantics></math>H2之间形成一条长带,分离超平面与它们平行且位于它们中央。长带的宽度,即<math><semantics><mrow><msub><mi>H</mi><mrow><mn>1</mn></mrow></msub></mrow><annotation encoding="application/x-tex">H_{1}</annotation></semantics></math>H1<math><semantics><mrow><msub><mi>H</mi><mrow><mn>2</mn></mrow></msub></mrow><annotation encoding="application/x-tex">H_{2}</annotation></semantics></math>H2之间的距离成为间隔(margin). 间隔依赖于分离超平面的法向量w,等于<math><semantics><mrow><mfrac><mrow><mn>2</mn></mrow><mrow><mi mathvariant="normal">∣</mi><mi mathvariant="normal">∣</mi><mi>w</mi><mi mathvariant="normal">∣</mi><mi mathvariant="normal">∣</mi></mrow></mfrac></mrow><annotation encoding="application/x-tex">\frac{2}{||w||}</annotation></semantics></math>w2.<math><semantics><mrow><msub><mi>H</mi><mrow><mn>1</mn></mrow></msub></mrow><annotation encoding="application/x-tex">H_{1}</annotation></semantics></math>H1<math><semantics><mrow><msub><mi>H</mi><mrow><mn>2</mn></mrow></msub></mrow><annotation encoding="application/x-tex">H_{2}</annotation></semantics></math>H2成为间隔边界。

在决定分离超平面时只有支持向量起作用,而其他实例点并不起作用。如果移动支持向量将改变所求的解;

由于支持向量在确定分离超平面中起着决定性作用,所以将这种分类模型称为支持向量机。支持向量一般很少,所以支持向量机由很少的“重要的”的训练样本确定。## 学习的对偶算法

为了求解线性可分支持向量机的最优化问题(7.13)-(7.14),将它作为原始最优化问题,应用拉格朗日对偶性,通过求解对偶(dual problem)问题得到原始问题(primal problem)的最优解,这就是线性可分支持向量机的对偶算法(dual algorithm)。这样做的优点,一是对偶问题往往更容易求解;而是自然引入核函数,进而推广到非线性分类问题。

KKT条件

<math><semantics><mrow><msub><mo>▽</mo><mrow><mi>w</mi></mrow></msub><mi>L</mi><mo>(</mo><msup><mi>w</mi><mo>∗</mo></msup><mo separator="true">,</mo><msup><mi>b</mi><mo>∗</mo></msup><mo separator="true">,</mo><msup><mi>α</mi><mo>∗</mo></msup><mo>)</mo><mo>=</mo><msup><mi>w</mi><mo>∗</mo></msup><mo>−</mo><msubsup><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi></mrow></msubsup><msubsup><mi>α</mi><mrow><mi>i</mi></mrow><mo>∗</mo></msubsup><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo>=</mo><mn>0</mn></mrow><annotation encoding="application/x-tex">\bigtriangledown_{w} L(w^*,b^*,\alpha^*)=w^*-\sum_{i=1}^{N}\alpha_{i}^*y_{i}x_{i}=0 </annotation></semantics></math>wL(w,b,α)=wi=1Nαiyixi=0

<math><semantics><mrow><msub><mo>▽</mo><mrow><mi>b</mi></mrow></msub><mi>L</mi><mo>(</mo><msup><mi>w</mi><mo>∗</mo></msup><mo separator="true">,</mo><msup><mi>b</mi><mo>∗</mo></msup><mo separator="true">,</mo><msup><mi>α</mi><mo>∗</mo></msup><mo>)</mo><mo>=</mo><mo>−</mo><msubsup><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi></mrow></msubsup><msubsup><mi>α</mi><mrow><mi>i</mi></mrow><mo>∗</mo></msubsup><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>=</mo><mn>0</mn></mrow><annotation encoding="application/x-tex">\bigtriangledown_{b}L(w^*,b^*,\alpha^*)=-\sum_{i=1}^{N}\alpha_{i}^*y_{i}=0 </annotation></semantics></math>bL(w,b,α)=i=1Nαiyi=0

<math><semantics><mrow><msubsup><mi>α</mi><mrow><mi>i</mi></mrow><mo>∗</mo></msubsup><mo>(</mo><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>(</mo><msup><mi>w</mi><mo>∗</mo></msup><mo>⋅</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo>+</mo><msup><mi>b</mi><mo>∗</mo></msup><mo>)</mo><mo>−</mo><mn>1</mn><mo>=</mo><mn>0</mn><mo separator="true">,</mo><mspace width="2em"></mspace><mi>i</mi><mo>=</mo><mn>1</mn><mo separator="true">,</mo><mn>2</mn><mo separator="true">,</mo><mo>⋯</mo><mo separator="true">,</mo><mi>N</mi></mrow><annotation encoding="application/x-tex">\alpha_{i}^*(y_{i}(w^*\cdot x_{i}+b^*)-1=0, \qquad i=1,2,\cdots,N </annotation></semantics></math>αi(yi(wxi+b)1=0,i=1,2,,N

<math><semantics><mrow><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>(</mo><msup><mi>w</mi><mo>∗</mo></msup><mo>⋅</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo>+</mo><msup><mi>b</mi><mo>∗</mo></msup><mo>)</mo><mo>−</mo><mn>1</mn><mo>≥</mo><mn>0</mn><mo separator="true">,</mo><mspace width="2em"></mspace><mi>i</mi><mo>=</mo><mn>1</mn><mo separator="true">,</mo><mn>2</mn><mo separator="true">,</mo><mo>⋯</mo><mo separator="true">,</mo><mi>N</mi></mrow><annotation encoding="application/x-tex">y_{i}(w^*\cdot x_{i}+b^*)-1 \ge 0, \qquad i=1,2,\cdots,N </annotation></semantics></math>yi(wxi+b)10,i=1,2,,N

<math><semantics><mrow><msubsup><mi>α</mi><mrow><mi>i</mi></mrow><mo>∗</mo></msubsup><mo>≥</mo><mn>0</mn><mo separator="true">,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo separator="true">,</mo><mn>2</mn><mo separator="true">,</mo><mo>⋯</mo><mo separator="true">,</mo><mi>N</mi></mrow><annotation encoding="application/x-tex">\alpha^*_{i}\ge 0, i=1,2,\cdots,N </annotation></semantics></math>αi0,i=1,2,,N

算法 线性可分支持向量机学习算法

输入: 线性可分训练集T,<math><semantics><mrow><mi>y</mi><mo>=</mo><mrow><mo>−</mo><mn>1</mn><mo separator="true">,</mo><mo>+</mo><mn>1</mn></mrow><mo separator="true">,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo separator="true">,</mo><mn>2</mn><mo separator="true">,</mo><mo>⋯</mo><mo separator="true">,</mo><mi>N</mi></mrow><annotation encoding="application/x-tex">y={-1,+1},i=1,2,\cdots,N</annotation></semantics></math>y=1,+1,i=1,2,,N
输出: 分离超平面和分类决策函数

(1)构造并求解约束最优化问题

<math><semantics><mrow><msub><mi>min</mi><mrow><mi>α</mi></mrow></msub><mfrac><mrow><mn>1</mn></mrow><mrow><mn>2</mn></mrow></mfrac><msubsup><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi></mrow></msubsup><msubsup><mo>∑</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi></mrow></msubsup><msub><mi>α</mi><mrow><mi>i</mi></mrow></msub><msub><mi>α</mi><mrow><mi>j</mi></mrow></msub><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><msub><mi>y</mi><mrow><mi>j</mi></mrow></msub><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo>⋅</mo><msub><mi>x</mi><mrow><mi>j</mi></mrow></msub><mo>)</mo><mo>−</mo><msubsup><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi></mrow></msubsup><msub><mi>α</mi><mrow><mi>i</mi></mrow></msub></mrow><annotation encoding="application/x-tex">\min_{\alpha} \frac{1}{2} \sum_{i=1}^{N} \sum_{j=1}^{N} \alpha_{i} \alpha_{j} y_{i} y_{j}(x_{i}\cdot x_{j})-\sum_{i=1}^{N}\alpha_{i} </annotation></semantics></math>αmin21i=1Nj=1Nαiαjyiyj(xixj)i=1Nαi

<math><semantics><mrow><msubsup><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi></mrow></msubsup><msub><mi>α</mi><mrow><mi>i</mi></mrow></msub><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>=</mo><mn>0</mn></mrow><annotation encoding="application/x-tex">\sum_{i=1}^{N}\alpha_{i}y_{i}=0 </annotation></semantics></math>i=1Nαiyi=0

<math><semantics><mrow><msub><mi>α</mi><mrow><mi>i</mi></mrow></msub><mo>≥</mo><mn>0</mn><mo separator="true">,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo separator="true">,</mo><mn>2</mn><mo separator="true">,</mo><mo>⋯</mo><mo separator="true">,</mo><mi>N</mi></mrow><annotation encoding="application/x-tex">\alpha_{i} \ge 0,i=1,2,\cdots,N </annotation></semantics></math>αi0,i=1,2,,N

求得最优解$$\alpha*=(\alpha_{1},\alpha_{2}*,\cdots,\alpha_{N})^T$$
<math><semantics><mrow><mi>α</mi></mrow><annotation encoding="application/x-tex">\alpha</annotation></semantics></math>α是拉格朗日乘子向量。

(2)计算

<math><semantics><mrow><msup><mi>w</mi><mo>∗</mo></msup><mo>=</mo><msubsup><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi></mrow></msubsup><msubsup><mi>α</mi><mrow><mi>i</mi></mrow><mo>∗</mo></msubsup><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub></mrow><annotation encoding="application/x-tex">w^*=\sum_{i=1}^{N}\alpha_{i}^*y_{i}x_{i} </annotation></semantics></math>w=i=1Nαiyixi

并选择<math><semantics><mrow><msup><mi>α</mi><mo>∗</mo></msup></mrow><annotation encoding="application/x-tex">\alpha^*</annotation></semantics></math>α的一个正分量<math><semantics><mrow><msubsup><mi>α</mi><mrow><mi>j</mi></mrow><mo>∗</mo></msubsup><mo>></mo><mn>0</mn></mrow><annotation encoding="application/x-tex">\alpha_{j}^*>0</annotation></semantics></math>αj>0,计算

<math><semantics><mrow><msup><mi>b</mi><mo>∗</mo></msup><mo>=</mo><msub><mi>y</mi><mrow><mi>j</mi></mrow></msub><mo>−</mo><msubsup><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi></mrow></msubsup><msubsup><mi>α</mi><mrow><mi>i</mi></mrow><mo>∗</mo></msubsup><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo>⋅</mo><msub><mi>x</mi><mrow><mi>j</mi></mrow></msub><mo>)</mo></mrow><annotation encoding="application/x-tex">b^*=y_{j}-\sum_{i=1}^{N}\alpha_{i}^*y_{i}(x_{i}\cdot x_{j}) </annotation></semantics></math>b=yji=1Nαiyi(xixj)

(3)求得分离超平面

<math><semantics><mrow><msup><mi>w</mi><mo>∗</mo></msup><mo>⋅</mo><mi>x</mi><mo>+</mo><msup><mi>b</mi><mo>∗</mo></msup><mo>=</mo><mn>0</mn></mrow><annotation encoding="application/x-tex">w^*\cdot x + b^*=0 </annotation></semantics></math>wx+b=0

分类决策函数:

<math><semantics><mrow><mi>f</mi><mo>(</mo><mi>x</mi><mo>)</mo><mo>=</mo><mi>s</mi><mi>i</mi><mi>g</mi><mi>n</mi><mo>(</mo><msup><mi>w</mi><mo>∗</mo></msup><mo>⋅</mo><mi>x</mi><mo>+</mo><msup><mi>b</mi><mo>∗</mo></msup><mo>)</mo></mrow><annotation encoding="application/x-tex">f(x)=sign(w^* \cdot x +b ^*) </annotation></semantics></math>f(x)=sign(wx+b)

<math><semantics><mrow><msup><mi>w</mi><mo>∗</mo></msup></mrow><annotation encoding="application/x-tex">w^*</annotation></semantics></math>w<math><semantics><mrow><msup><mi>b</mi><mo>∗</mo></msup></mrow><annotation encoding="application/x-tex">b^*</annotation></semantics></math>b只依赖于训练数据集中对应<math><semantics><mrow><msubsup><mi>α</mi><mrow><mi>i</mi></mrow><mo>∗</mo></msubsup><mo>></mo><mn>0</mn></mrow><annotation encoding="application/x-tex">\alpha_{i}^*>0</annotation></semantics></math>αi>0的样本点<math><semantics><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo separator="true">,</mo><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>)</mo></mrow><annotation encoding="application/x-tex">(x_{i},y_{i})</annotation></semantics></math>(xi,yi),而其他样本点对<math><semantics><mrow><msup><mi>w</mi><mo>∗</mo></msup></mrow><annotation encoding="application/x-tex">w^*</annotation></semantics></math>w<math><semantics><mrow><msup><mi>b</mi><mo>∗</mo></msup></mrow><annotation encoding="application/x-tex">b^*</annotation></semantics></math>b没有影响。我们将训练数据中对应于<math><semantics><mrow><msubsup><mi>α</mi><mrow><mi>i</mi></mrow><mo>∗</mo></msubsup><mo>></mo><mn>0</mn></mrow><annotation encoding="application/x-tex">\alpha_{i}^*>0</annotation></semantics></math>αi>0的实例点<math><semantics><mrow><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo>∈</mo><msup><mi>R</mi><mi>n</mi></msup></mrow><annotation encoding="application/x-tex">x_{i}\in R^n</annotation></semantics></math>xiRn称为支持向量。

支持向量一定在间隔边界上。

<math><semantics><mrow><msubsup><mi>α</mi><mrow><mi>i</mi></mrow><mo>∗</mo></msubsup><mo>(</mo><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>(</mo><msup><mi>w</mi><mo>∗</mo></msup><mo>⋅</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo>+</mo><msup><mi>b</mi><mo>∗</mo></msup><mo>)</mo><mo>−</mo><mn>1</mn><mo>)</mo><mo>=</mo><mn>0</mn><mo separator="true">,</mo><mspace width="2em"></mspace><mi>i</mi><mo>=</mo><mn>1</mn><mo separator="true">,</mo><mn>2</mn><mo separator="true">,</mo><mo>⋯</mo><mo separator="true">,</mo><mi>N</mi></mrow><annotation encoding="application/x-tex">\alpha_{i}^*(y_{i}(w^* \cdot x_{i} +b^*)-1)=0, \qquad i=1,2,\cdots,N </annotation></semantics></math>αi(yi(wxi+b)1)=0,i=1,2,,N

<math><semantics><mrow><msup><mi>w</mi><mo>∗</mo></msup><mo>⋅</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo>+</mo><msup><mi>b</mi><mo>∗</mo></msup><mo>=</mo><mo>±</mo><mn>1</mn></mrow><annotation encoding="application/x-tex">w^* \cdot x_{i} +b ^*= \pm 1 </annotation></semantics></math>wxi+b=±1

<math><semantics><mrow><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub></mrow><annotation encoding="application/x-tex">x_{i}</annotation></semantics></math>xi一定在间隔边界上.这里的支持向量的定义与前面给出的支持向量的定义是一致的。

对于线性可分问题,上述线性可分支持向量机的学习(硬间隔最大化)算法是完美的。但是,训练集线性可分是理想的情形。在现实问题中,训练数据集往往是线性不可分的。即在样本中出现噪声或离群点。此时,有更一般的学习算法。

讨论:为什么要用对偶算法?

  1. 不等式约束一直是优化问题中的难题,求解对偶问题可以将支持向量机原问题约束中的不等式约束转化为等式约束;参看KKT条件
  1. 支持向量机中用到了高维映射,但是映射函数的具体形式几乎完全不可确定,而求解对偶问题之后,可以使用核函数来解决这个问题。

线性支持向量机与软间隔最大化

线性可分问题的支持向量学习方法,对线性不可分数据是不适用的,因为这时上述方法中的不等式约束并不都能成立。这就需要修改硬间隔最大化,使其成为软间隔最大化。

线性不可分意味着某些样本点<math><semantics><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo separator="true">,</mo><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>)</mo></mrow><annotation encoding="application/x-tex">(x_{i},y_{i})</annotation></semantics></math>(xi,yi)不能满足函数间隔大于等于1的约束条件。为了解决这个问题,可以对每个样本点<math><semantics><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo separator="true">,</mo><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>)</mo></mrow><annotation encoding="application/x-tex">(x_{i},y_{i})</annotation></semantics></math>(xi,yi)引进一个松弛变量<math><semantics><mrow><msub><mi>ξ</mi><mrow><mi>i</mi></mrow></msub><mo>≥</mo><mn>0</mn></mrow><annotation encoding="application/x-tex">\xi_{i} \ge 0</annotation></semantics></math>ξi0,使函数间隔加上松弛变量大于等于1.这样约束条件变为

<math><semantics><mrow><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>(</mo><mi>w</mi><mo>⋅</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo>+</mo><mi>b</mi><mo>)</mo><mo>≥</mo><mn>1</mn><mo>−</mo><msub><mi>ξ</mi><mrow><mi>i</mi></mrow></msub></mrow><annotation encoding="application/x-tex">y_{i}(w \cdot x_{i} +b) \ge 1- \xi_{i} </annotation></semantics></math>yi(wxi+b)1ξi

同时,对每个松弛变量<math><semantics><mrow><msub><mi>ξ</mi><mrow><mi>i</mi></mrow></msub></mrow><annotation encoding="application/x-tex">\xi_{i}</annotation></semantics></math>ξi,支付一个代价<math><semantics><mrow><msub><mi>ξ</mi><mrow><mi>i</mi></mrow></msub></mrow><annotation encoding="application/x-tex">\xi_{i}</annotation></semantics></math>ξi,目标函数由原来的<math><semantics><mrow><mfrac><mrow><mn>1</mn></mrow><mrow><mn>2</mn></mrow></mfrac><mi mathvariant="normal">∣</mi><mi mathvariant="normal">∣</mi><mi>w</mi><mi mathvariant="normal">∣</mi><msup><mi mathvariant="normal">∣</mi><mn>2</mn></msup></mrow><annotation encoding="application/x-tex">\frac{1}{2}||w||^2</annotation></semantics></math>21w2变成

<math><semantics><mrow><mfrac><mrow><mn>1</mn></mrow><mrow><mn>2</mn></mrow></mfrac><mi mathvariant="normal">∣</mi><mi mathvariant="normal">∣</mi><mi>w</mi><mi mathvariant="normal">∣</mi><msup><mi mathvariant="normal">∣</mi><mn>2</mn></msup><mo>+</mo><mi>C</mi><msubsup><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi></mrow></msubsup><msub><mi>ξ</mi><mrow><mi>i</mi></mrow></msub></mrow><annotation encoding="application/x-tex">\frac{1}{2}||w||^2+C \sum_{i=1}^{N} \xi_{i} </annotation></semantics></math>21w2+Ci=1Nξi

这里,C>0称为惩罚参数,一般由应用问题决定,C值大时对误分类的惩罚增大,C值小时对误分类的惩罚减小。最小化目标函数包含两层含义:使<math><semantics><mrow><mfrac><mrow><mn>1</mn></mrow><mrow><mn>2</mn></mrow></mfrac><mi mathvariant="normal">∣</mi><mi mathvariant="normal">∣</mi><mi>w</mi><mi mathvariant="normal">∣</mi><msup><mi mathvariant="normal">∣</mi><mn>2</mn></msup></mrow><annotation encoding="application/x-tex">\frac{1}{2}||w||^2</annotation></semantics></math>21w2尽量小即间隔尽量大,同时使误分类点的个数尽量小,C是调和二者的系数。

有了松弛变量和惩罚,我们就可以用处理线性可分数据的思路来处理线性不可分数据。但相对于硬间隔可分的绝对分割,我们称这种方法为软间隔最大化。

线性不可分的支持向量机的学习问题变成凸二次规划(convex quadratic programming)问题(原始问题):

算法 线性支持向量机学习算法

输入: 线性可分训练集T,<math><semantics><mrow><mi>y</mi><mo>=</mo><mrow><mo>−</mo><mn>1</mn><mo separator="true">,</mo><mo>+</mo><mn>1</mn></mrow><mo separator="true">,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo separator="true">,</mo><mn>2</mn><mo separator="true">,</mo><mo>⋯</mo><mo separator="true">,</mo><mi>N</mi></mrow><annotation encoding="application/x-tex">y={-1,+1},i=1,2,\cdots,N</annotation></semantics></math>y=1,+1,i=1,2,,N
输出: 分离超平面和分类决策函数

(1)选择惩罚参数C>0,构造并求解凸二次规划问题

<math><semantics><mrow><msub><mi>min</mi><mrow><mi>α</mi></mrow></msub><mfrac><mrow><mn>1</mn></mrow><mrow><mn>2</mn></mrow></mfrac><msubsup><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi></mrow></msubsup><msubsup><mo>∑</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi></mrow></msubsup><msub><mi>α</mi><mrow><mi>i</mi></mrow></msub><msub><mi>α</mi><mrow><mi>j</mi></mrow></msub><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><msub><mi>y</mi><mrow><mi>j</mi></mrow></msub><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo>⋅</mo><msub><mi>x</mi><mrow><mi>j</mi></mrow></msub><mo>)</mo><mo>−</mo><msubsup><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi></mrow></msubsup><msub><mi>α</mi><mrow><mi>i</mi></mrow></msub></mrow><annotation encoding="application/x-tex">\min_{\alpha} \frac{1}{2} \sum_{i=1}^{N} \sum_{j=1}^{N} \alpha_{i} \alpha_{j} y_{i} y_{j}(x_{i}\cdot x_{j})-\sum_{i=1}^{N}\alpha_{i} </annotation></semantics></math>αmin21i=1Nj=1Nαiαjyiyj(xixj)i=1Nαi

<math><semantics><mrow><msubsup><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi></mrow></msubsup><msub><mi>α</mi><mrow><mi>i</mi></mrow></msub><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>=</mo><mn>0</mn></mrow><annotation encoding="application/x-tex">\sum_{i=1}^{N}\alpha_{i}y_{i}=0 </annotation></semantics></math>i=1Nαiyi=0

<math><semantics><mrow><mn>0</mn><mo>≤</mo><msub><mi>α</mi><mrow><mi>i</mi></mrow></msub><mo>≤</mo><mi>C</mi><mo separator="true">,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo separator="true">,</mo><mn>2</mn><mo separator="true">,</mo><mspace width="2em"></mspace><mo>⋯</mo><mo separator="true">,</mo><mi>N</mi></mrow><annotation encoding="application/x-tex">0 \le \alpha_{i} \le C,i=1,2, \qquad \cdots,N </annotation></semantics></math>0αiC,i=1,2,,N

求得最优解$$\alpha*=(\alpha_{1},\alpha_{2}*,\cdots,\alpha_{N})^T$$
<math><semantics><mrow><mi>α</mi></mrow><annotation encoding="application/x-tex">\alpha</annotation></semantics></math>α是拉格朗日乘子向量。

(2)计算

<math><semantics><mrow><msup><mi>w</mi><mo>∗</mo></msup><mo>=</mo><msubsup><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi></mrow></msubsup><msubsup><mi>α</mi><mrow><mi>i</mi></mrow><mo>∗</mo></msubsup><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub></mrow><annotation encoding="application/x-tex">w^*=\sum_{i=1}^{N}\alpha_{i}^*y_{i}x_{i} </annotation></semantics></math>w=i=1Nαiyixi

并选择<math><semantics><mrow><msup><mi>α</mi><mo>∗</mo></msup></mrow><annotation encoding="application/x-tex">\alpha^*</annotation></semantics></math>α的一个正分量<math><semantics><mrow><mn>0</mn><mo>≤</mo><msubsup><mi>α</mi><mrow><mi>j</mi></mrow><mo>∗</mo></msubsup><mo>≤</mo><mi>C</mi></mrow><annotation encoding="application/x-tex">0 \le \alpha_{j}^* \le C</annotation></semantics></math>0αjC,计算

<math><semantics><mrow><msup><mi>b</mi><mo>∗</mo></msup><mo>=</mo><msub><mi>y</mi><mrow><mi>j</mi></mrow></msub><mo>−</mo><msubsup><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi></mrow></msubsup><msubsup><mi>α</mi><mrow><mi>i</mi></mrow><mo>∗</mo></msubsup><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo>⋅</mo><msub><mi>x</mi><mrow><mi>j</mi></mrow></msub><mo>)</mo></mrow><annotation encoding="application/x-tex">b^*=y_{j}-\sum_{i=1}^{N}\alpha_{i}^*y_{i}(x_{i}\cdot x_{j}) </annotation></semantics></math>b=yji=1Nαiyi(xixj)

(3)求得分离超平面

<math><semantics><mrow><msup><mi>w</mi><mo>∗</mo></msup><mo>⋅</mo><mi>x</mi><mo>+</mo><msup><mi>b</mi><mo>∗</mo></msup><mo>=</mo><mn>0</mn></mrow><annotation encoding="application/x-tex">w^*\cdot x + b^*=0 </annotation></semantics></math>wx+b=0

分类决策函数:

<math><semantics><mrow><mi>f</mi><mo>(</mo><mi>x</mi><mo>)</mo><mo>=</mo><mi>s</mi><mi>i</mi><mi>g</mi><mi>n</mi><mo>(</mo><msup><mi>w</mi><mo>∗</mo></msup><mo>⋅</mo><mi>x</mi><mo>+</mo><msup><mi>b</mi><mo>∗</mo></msup><mo>)</mo></mrow><annotation encoding="application/x-tex">f(x)=sign(w^* \cdot x +b ^*) </annotation></semantics></math>f(x)=sign(wx+b)

支持向量

软间隔的支持向量<math><semantics><mrow><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub></mrow><annotation encoding="application/x-tex">x_{i}</annotation></semantics></math>xi在间隔边界上(<math><semantics><mrow><msubsup><mi>α</mi><mrow><mi>i</mi></mrow><mo>∗</mo></msubsup><mo><</mo><mi>C</mi><mo separator="true">,</mo><msub><mi>ξ</mi><mrow><mi>i</mi></mrow></msub><mo>=</mo><mn>0</mn></mrow><annotation encoding="application/x-tex">\alpha_{i}^*<C,\xi_{i}=0</annotation></semantics></math>αi<C,ξi=0)或者在间隔边界与超平面之间(<math><semantics><mrow><msubsup><mi>α</mi><mrow><mi>i</mi></mrow><mo>∗</mo></msubsup><mo>=</mo><mi>C</mi><mo separator="true">,</mo><mn>0</mn><mo><</mo><msub><mi>ξ</mi><mrow><mi>i</mi></mrow></msub><mo><</mo><mn>1</mn></mrow><annotation encoding="application/x-tex">\alpha_{i}^*=C,0 < \xi_{i}<1</annotation></semantics></math>αi=C,0<ξi<1),在超平面上((<math><semantics><mrow><msubsup><mi>α</mi><mrow><mi>i</mi></mrow><mo>∗</mo></msubsup><mo>=</mo><mi>C</mi><mo separator="true">,</mo><msub><mi>ξ</mi><mrow><mi>i</mi></mrow></msub><mo>=</mo><mn>1</mn></mrow><annotation encoding="application/x-tex">\alpha_{i}^*=C,\xi_{i}=1</annotation></semantics></math>αi=C,ξi=1)),甚至在分离超平面误分的一侧(<math><semantics><mrow><msubsup><mi>α</mi><mrow><mi>i</mi></mrow><mo>∗</mo></msubsup><mo>=</mo><mi>C</mi><mo separator="true">,</mo><msub><mi>ξ</mi><mrow><mi>i</mi></mrow></msub><mo>></mo><mn>1</mn></mrow><annotation encoding="application/x-tex">\alpha_{i}^*=C,\xi_{i}>1</annotation></semantics></math>αi=C,ξi>1)。

合页损失函数

对于线性支持向量机来说,还有另外一种解释,就是最小化以下目标函数:

<math><semantics><mrow><msubsup><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi></mrow></msubsup><mo>[</mo><mn>1</mn><mo>−</mo><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>(</mo><mi>w</mi><mo>⋅</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo>+</mo><mi>b</mi><msub><mo>)</mo><mrow><mo>+</mo></mrow></msub><mo>+</mo><mi>λ</mi><mi mathvariant="normal">∣</mi><mi mathvariant="normal">∣</mi><mi>w</mi><mi mathvariant="normal">∣</mi><msup><mi mathvariant="normal">∣</mi><mn>2</mn></msup><mo>]</mo></mrow><annotation encoding="application/x-tex">\sum_{i=1}^{N}[1-y_{i}(w \cdot x_{i} + b )_{+} + \lambda ||w||^2] </annotation></semantics></math>i=1N[1yi(wxi+b)++λw2]

目标函数的第一项是经验损失或经验风险函数,第二项是正则项。

<math><semantics><mrow><mi>L</mi><mo>(</mo><mi>y</mi><mo>(</mo><mi>w</mi><mo>⋅</mo><mi>x</mi><mo>+</mo><mi>b</mi><mo>)</mo><mo>)</mo><mo>=</mo><mo>[</mo><mn>1</mn><mo>−</mo><mi>y</mi><mo>(</mo><mi>w</mi><mo>⋅</mo><mi>w</mi><mo>+</mo><mi>b</mi><mo>)</mo><msub><mo>]</mo><mrow><mo>+</mo></mrow></msub></mrow><annotation encoding="application/x-tex">L(y(w \cdot x + b))=[1-y(w \cdot w +b)]_{+} </annotation></semantics></math>L(y(wx+b))=[1y(ww+b)]+

称为合页损失函数(hinge loss function)。下标"+"表示以下取正值的函数。

<math><semantics><mrow><mo>[</mo><mi>z</mi><msub><mo>]</mo><mrow><mo>+</mo></mrow></msub><mo>=</mo><msubsup><mo>{</mo><mrow><mi>z</mi><mo separator="true">,</mo><mspace width="1em"></mspace><mi>z</mi><mo>≤</mo><mn>0</mn></mrow><mrow><mn>0</mn><mo separator="true">,</mo><mspace width="1em"></mspace><mi>z</mi><mo>></mo><mn>0</mn></mrow></msubsup></mrow><annotation encoding="application/x-tex">[z]_{+}= \lbrace_{z, \quad z \le 0}^{0, \quad z >0} </annotation></semantics></math>[z]+={z,z00,z>0

就是说,当样本点<math><semantics><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo separator="true">,</mo><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>)</mo></mrow><annotation encoding="application/x-tex">(x_{i},y_{i})</annotation></semantics></math>(xi,yi)被正确分类且函数间隔(确信度)$$y_{i}(w \cdot x_{i} +b)$$大于1时,损失是0,否者损失是<math><semantics><mrow><mn>1</mn><mo>−</mo><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>(</mo><mi>w</mi><mo>⋅</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mo>+</mo><mi>b</mi><mo>)</mo></mrow><annotation encoding="application/x-tex">1-y_{i}(w \cdot x_{i} +b)</annotation></semantics></math>1yi(wxi+b)

合页损失函数横轴是函数间隔,纵轴是损失。由于其函数形状像一个合页,故名合页损失函数。

可以认为他是二分类问题真正的损失函数,而合页损失函数是0-1损失函数的上界,但0-1损失函数不是连续可导的,直接优化尤其构成的目标函数比较困难,可以认为线性支持向量机优化0-1损失函数的上界构成的目标函数。这时的上界损失函数又称为代理损失函数(surrogate loss function)。# 非线性支持向量机与核函数

核技巧

非线性问题往往不好求解,所以希望能用解线性分类问题的方法解决这个问题。所采取的方法是进行一个非线性变换,将非线性问题变换为线性问题,通过解变换后的线性问题的方法求解原来的非线性问题。

核技巧应用到支持向量机,其基本想法就是通过一个非线性变换将输入空间(欧式空间<math><semantics><mrow><msup><mi>R</mi><mi>n</mi></msup></mrow><annotation encoding="application/x-tex">R^n</annotation></semantics></math>Rn或离散集合)对应一个特征空间(希尔伯特空间<math><semantics><mrow><mrow><mi mathvariant="script">H</mi></mrow></mrow><annotation encoding="application/x-tex">\mathcal{H}</annotation></semantics></math>H),使得在输入空间<math><semantics><mrow><msup><mi>R</mi><mi>n</mi></msup></mrow><annotation encoding="application/x-tex">R^n</annotation></semantics></math>Rn中的超平面模型对应于特征空间<math><semantics><mrow><mrow><mi mathvariant="script">H</mi></mrow></mrow><annotation encoding="application/x-tex">\mathcal{H}</annotation></semantics></math>H中的超平面模型(支持向量)。这样,分类问题的学习任务通过在特征空间中求解线性支持向量机的任务就可以完成。

核函数

<math><semantics><mrow><mrow><mi mathvariant="script">X</mi></mrow></mrow><annotation encoding="application/x-tex">\mathcal{X}</annotation></semantics></math>X是输入空间(欧式空间<math><semantics><mrow><msup><mi>R</mi><mi>n</mi></msup></mrow><annotation encoding="application/x-tex">R^n</annotation></semantics></math>Rn的子集或离散集合),又设<math><semantics><mrow><mrow><mi mathvariant="script">H</mi></mrow></mrow><annotation encoding="application/x-tex">\mathcal{H}</annotation></semantics></math>H是特征空间(希尔伯特空间),如果存在一个从<math><semantics><mrow><mrow><mi mathvariant="script">X</mi></mrow></mrow><annotation encoding="application/x-tex">\mathcal{X}</annotation></semantics></math>X<math><semantics><mrow><mrow><mi mathvariant="script">H</mi></mrow></mrow><annotation encoding="application/x-tex">\mathcal{H}</annotation></semantics></math>H的映射<math><semantics><mrow><mi>ϕ</mi><mo>(</mo><mi>x</mi><mo>)</mo><mo>:</mo><mrow><mi mathvariant="script">X</mi></mrow><mo>→</mo><mrow><mi mathvariant="script">H</mi></mrow></mrow><annotation encoding="application/x-tex">\phi (x): \mathcal{X} \rightarrow \mathcal{H}</annotation></semantics></math>ϕ(x):XH
使得对所有<math><semantics><mrow><mi>x</mi><mo separator="true">,</mo><mi>z</mi><mo>∈</mo><mrow><mi mathvariant="script">X</mi></mrow></mrow><annotation encoding="application/x-tex">x,z \in \mathcal{X}</annotation></semantics></math>x,zX,函数$$K(x,z)=\phi (x) \cdot \phi(z)$$
则成K(x,z)为核函数,<math><semantics><mrow><mi>ϕ</mi><mo>(</mo><mi>x</mi><mo>)</mo></mrow><annotation encoding="application/x-tex">\phi(x)</annotation></semantics></math>ϕ(x)为映射函数,式中<math><semantics><mrow><mi>ϕ</mi><mo>(</mo><mi>x</mi><mo>)</mo><mo>⋅</mo><mi>ϕ</mi><mo>(</mo><mi>z</mi><mo>)</mo></mrow><annotation encoding="application/x-tex">\phi(x) \cdot \phi(z)</annotation></semantics></math>ϕ(x)ϕ(z)<math><semantics><mrow><mi>ϕ</mi><mo>(</mo><mi>x</mi><mo>)</mo></mrow><annotation encoding="application/x-tex">\phi(x)</annotation></semantics></math>ϕ(x)<math><semantics><mrow><mi>ϕ</mi><mo>(</mo><mi>z</mi><mo>)</mo></mrow><annotation encoding="application/x-tex">\phi(z)</annotation></semantics></math>ϕ(z)的内积。

核技巧的想法是,在学习与预测中定义核函数K(x,z),而不是显示地定义映射函数。因为计算核函数更容易。

常用核函数

  1. 多项式核函数 (polynomial kernl function)

<math><semantics><mrow><mi>K</mi><mo>(</mo><mi>x</mi><mo separator="true">,</mo><mi>z</mi><mo>)</mo><mo>=</mo><mo>(</mo><mi>x</mi><mo>⋅</mo><mi>z</mi><mo>+</mo><mn>1</mn><msup><mo>)</mo><mi>p</mi></msup></mrow><annotation encoding="application/x-tex">K(x, z)=(x \cdot z +1)^p </annotation></semantics></math>K(x,z)=(xz+1)p

对应的支持向量机是一个p次多项式分类器。在此情形下,分类决策函数成为 <math><semantics><mrow><mi>f</mi><mo>(</mo><mi>x</mi><mo>)</mo><mo>=</mo><mi>s</mi><mi>i</mi><mi>g</mi><mi>n</mi><mo>(</mo><msub><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow></msub><mrow><msub><mi>N</mi><mrow><mi>s</mi></mrow></msub></mrow><msubsup><mi>a</mi><mrow><mi>i</mi></mrow><mo>∗</mo></msubsup><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi></mrow></msub><mover accent="true"><mi>x</mi><mo>˙</mo></mover><mo>+</mo><mn>1</mn><msup><mo>)</mo><mi>p</mi></msup><mo>+</mo><msup><mi>b</mi><mo>∗</mo></msup><mo>)</mo></mrow><annotation encoding="application/x-tex">f(x)=sign( \sum_{i=1}{N_{s}} a_{i}^* y_{i}(x_{i} \dot x +1)^p + b^*)</annotation></semantics></math>f(x)=sign(i=1Nsaiyi(xix˙+1)p+b)

  1. 高斯核函数(Gaussian kernel function)

<math><semantics><mrow><mi>K</mi><mo>(</mo><mi>x</mi><mo separator="true">,</mo><mi>z</mi><mo>)</mo><mo>=</mo><mi>e</mi><mi>x</mi><mi>p</mi><mo>(</mo><mo>−</mo><mfrac><mrow><mi mathvariant="normal">∣</mi><mi mathvariant="normal">∣</mi><mi>x</mi><mo>−</mo><mi>z</mi><mi mathvariant="normal">∣</mi><msup><mi mathvariant="normal">∣</mi><mn>2</mn></msup></mrow><mrow><mn>2</mn><mi>σ</mi></mrow></mfrac><mo>)</mo></mrow><annotation encoding="application/x-tex">K(x,z)=exp(- \frac{||x-z||^2}{2 \sigma}) </annotation></semantics></math>K(x,z)=exp(2σxz2)

对应的支持向量机是高斯径向基核函数(radial basis function)分类器。在此情形下,分类决策函数成为

<math><semantics><mrow><mi>f</mi><mo>(</mo><mi>x</mi><mo>)</mo><mo>=</mo><mi>s</mi><mi>i</mi><mi>g</mi><mi>n</mi><mo>(</mo><msubsup><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><msub><mi>N</mi><mrow><mi>s</mi></mrow></msub></mrow></msubsup><msubsup><mi>a</mi><mrow><mi>i</mi></mrow><mo>∗</mo></msubsup><msub><mi>y</mi><mrow><mi>i</mi></mrow></msub><mi>e</mi><mi>x</mi><mi>p</mi><mo>(</mo><mo>−</mo><mfrac><mrow><mi mathvariant="normal">∣</mi><mi mathvariant="normal">∣</mi><mi>x</mi><mo>−</mo><mi>z</mi><mi mathvariant="normal">∣</mi><msup><mi mathvariant="normal">∣</mi><mn>2</mn></msup></mrow><mrow><mn>2</mn><msup><mi>σ</mi><mn>2</mn></msup></mrow></mfrac><mo>)</mo><mo>+</mo><mi>b</mi><mo>)</mo></mrow><annotation encoding="application/x-tex">f(x)=sign(\sum_{i=1}^{N_{s}}a_{i}^*y_{i}exp(-\frac{||x-z||^2}{2 \sigma ^2})+b) </annotation></semantics></math>f(x)=sign(i=1Nsaiyiexp(2σ2xz2)+b)

讨论:如何选择核函数

n=特征数,m=训练样本数目

如果n相对m比较大,使用逻辑回归或者线性SVM,例如n=10000,m=10-1000

如果n比较小,m数值适中,使用高斯核的SVM,例如n=1-1000,m=10-10000

如果n比较小,m很大,构建更多特征,然后使用逻辑回归或者线性SVM参考阅读:
支持向量机SVM http://www.cnblogs.com/jerrylead/archive/2011/03/13/1982639.html
机器学习与数据挖掘-支持向量机(SVM) https://wizardforcel.gitbooks.io/dm-algo-top10/content/svm.html

评论区

PhD Candidate in Machine Learning

0

0

0

举报