本文转载至 跳转
支持向量机,因其英文名为 Support Vector Machine,故一般简称SVM,通俗来讲,它是一种二类分类模型,其基本模型定义为特征空间上的间隔最大的线性分类器,其学习策略便是间隔最大化,最终可转化为一个凸二次规划问题的求解。
理解SVM,咱们必须先弄清楚一个概念:线性分类器
给定一些数据点,它们分别属于两个不同的类,现在要找到一个线性分类器把这些数据分成两类。如果用 x x x 表示数据点,用 y y y 表示类别( y y y 可以取 1 1 1 或者 − 1 -1 −1,分别代表两个不同的类),一个线性分类器的学习目标便是要在 n n n 维的数据空间中找到一个 超平面(hyper plane),这个超平面的方程可以表示为( w T w^T wT 中的 T T T 代表转置):
(1) w T + b = 0 w^T + b =0 \tag1 wT+b=0(1)
可能有读者对类别取 1 1 1 或 − 1 -1 −1 有疑问,事实上,这个 1 1 1 或 − 1 -1 −1 的分类标准起源于 Logistic回归。
Logistic回归 目的是从特征学习出一个 0 / 1 0/1 0/1 分类模型,而这个模型是将特性的线性组合作为自变量,由于自变量的取值范围是负无穷到正无穷。因此,使用 Logistic回归(或称作Sigmoid函数)将自变量映射到(0,1)上,映射后的值被认为是属于 y = 1 y=1 y=1 的概率。
假设函数:
(2) h θ ( x ) = g ( θ T x ) = 1 1 + e − θ T x h_{\theta}(x) = g({\theta^T}x) = \frac{1}{1+e^{-{\theta^T}x}} \tag2 hθ(x)=g(θTx)=1+e−θTx1(2)
其中 x x x 是 n n n 维特征向量,函数 g g g 就是 Logistic 函数,而 g ( z ) = 1 1 + e − z g(z)=\frac{1}{1+e^{-z}} g(z)=1+e−z1 即 Sigmoid函数 的图像是
可以看到,将无穷映射到了(0, 1),而假设函数就是特征属于 y = 1 y=1 y=1 的概率。
p ( y = 1 ∣ x ; θ ) = h θ ( x ) p(y = 1 \ | \ x \ ; \ \theta ) = h_{\theta}(x) p(y=1 ∣ x ; θ)=hθ(x)
p ( y = 0 ∣ x ; θ ) = 1 − h θ ( x ) p(y = 0 \ | \ x \ ; \ \theta ) = 1 - h_{\theta}(x) p(y=0 ∣ x ; θ)=1−hθ(x)
从而,当我们要判别一个新来的特征属于哪个类时,只需求 h θ ( x ) h_{\theta}(x) hθ(x) 即可
若 h θ ( x ) > 0.5 h_{\theta}(x) > 0.5 hθ(x)>0.5 就是 y = 1 y=1 y=1 的类
若 h θ ( x ) ⩽ 0.5 h_{\theta}(x) \leqslant 0.5 hθ(x)⩽0.5 就是 y = 0 y=0 y=0 的类
此外, h θ ( x ) h_{\theta}(x) hθ(x) 只和 θ T x \theta^Tx θTx 有关, 而 g ( z ) g(z) g(z) 只是用来映射,真实的类别决定权还是在于 θ T x \theta^Tx θTx。
若 θ T x > 0 \theta^Tx > 0 θTx>0 ,那么 h θ ( x ) > 0.5 h_{\theta}(x)>0.5 hθ(x)>0.5;
若 θ T x ≫ 0 \theta^Tx \gg 0 θTx≫0,那么 h θ ( x ) = 1 h_{\theta}(x) =1 hθ(x)=1;
若 θ T x ≪ 0 \theta^Tx \ll 0 θTx≪0,那么 h θ ( x ) = 0 h_{\theta}(x) =0 hθ(x)=0
如果我们只从 θ T x \theta^Tx θTx 出发,希望模型达到的目标就是让训练数据中 y = 1 y=1 y=1 的特征 θ T x ≫ 0 \theta^Tx \gg 0 θTx≫0,而是 y = 0 y=0 y=0 的特征 θ T x ≪ 0 \theta^Tx \ll 0 θTx≪0。Logistic回归就是要学习得到 θ \theta θ,使得正例的特征远大于0,负例的特征远小于0,而且要在全部训练实例上达到这个目标。
接下来,尝试把 Logistic回归 做个变形。首先,将使用的结果标签 y = 0 y=0 y=0 和 y = 1 y=1 y=1 替换为 y = − 1 y=-1 y=−1 和 y = 1 y=1 y=1,然后将 θ T x = θ 0 x 0 + θ 1 x 1 + θ 2 x 2 + ⋯ + θ n x n ( x 0 = 1 ) \theta^Tx = \theta_0x_0 + \theta_1x_1 + \theta_2x_2 + \cdots + \theta_nx_n(x_0 = 1) θTx=θ0x0+θ1x1+θ2x2+⋯+θnxn(x0=1) 中的 θ 0 \theta_0 θ0 替换为 b b b,最后将后面的 θ 1 x 1 + θ 2 x 2 + ⋯ + θ n x n \theta_1x_1 + \theta_2x_2 + \cdots + \theta_nx_n θ1x1+θ2x2+⋯+θnxn 替换为 w T x w^Tx wTx 。如此,则有了 θ T x = w T x + b \theta^Tx = w^Tx +b θTx=wTx+b 。也就是说除了 y y y 由 y = 0 y=0 y=0 变为 y = − 1 y=-1 y=−1 外,线性分类函数跟Logistic回归的形式化表示没区别。
进一步,可以将假设函数 h w , b ( x ) = g ( w T x + b ) h_{w,b}(x) = g(w^Tx +b) hw,b(x)=g(wTx+b) 中的 g ( z ) g(z) g(z) 做一个简化,将其简单映射到 y = − 1 y=-1 y=−1 和 y = 1 y=1 y=1 上。映射关系如下:
(3) g ( z ) = { 1 z ≥ 0 − 1 z < 0 g(z)= \left\{\begin{matrix} 1 & z \geq 0\\ -1 & z < 0 \end{matrix}\right. \tag3 g(z)={1−1z≥0z<0(3)
下面举个简单的例子。如下图所示,现在有一个二维平面,平面上有两种不同的数据,分别用圈和叉表示。由于这些数据是线性可分的,所以可以用一条直线将这两类数据分开,这条直线就相当于一个超平面,超平面一边的数据点所对应的 y y y 全是 − 1 -1 −1 ,另一边所对应的 y y y 全是 1 1 1 。
这个超平面可以用分类函数 f ( x ) = w T + b f(x) = w^T + b f(x)=wT+b 表示,
当 f ( x ) = 0 f(x)=0 f(x)=0 的时候, x x x 便是位于超平面上的点,
而 f ( x ) > 0 f(x)>0 f(x)>0 的时候, x x x 对应 y = 1 y=1 y=1 的点,
而 f ( x ) < 0 f(x)<0 f(x)<0 的时候, x x x 对应 y = − 1 y=-1 y=−1 的点,
如下图所示:
注:有的资料上定义特征到结果的输出函数 u = w ˉ x ˉ − b u =\bar{w}\bar{x} -b u=wˉxˉ−b,与这里定义的 f ( x ) = w T + b f(x) = w^T + b f(x)=wT+b 实质是一样的。为什么?因为无论是 u = w ˉ x ˉ − b u =\bar{w}\bar{x} -b u=wˉxˉ−b,还是 f ( x ) = w T + b f(x) = w^T + b f(x)=wT+b,不影响最终优化结果。下文你将看到,当我们转化到优化 max 1 ∣ ∣ w ∣ ∣ s . t . , y i ( w T + b ) ≥ 1 , i = 1 , ⋯   , n \max \frac{1}{||w||} \ s.t.,\ y_i(w^T + b) \geq 1, \ i=1,\cdots ,n max∣∣w∣∣1 s.t., yi(wT+b)≥1, i=1,⋯,n 的时候,为了求解方便,会把 y f ( x ) yf(x) yf(x) 令为 1 1 1,即 y f ( x ) yf(x) yf(x) 是 y ( w T x + b ) y(w^Tx + b) y(wTx+b),还是 y ( w T x − b ) y(w^Tx - b) y(wTx−b),对我们要优化的式子 max 1 ∣ ∣ w ∣ ∣ \max\frac{1}{||w||} max∣∣w∣∣1 已无影响。
(有一朋友飞狗来自Mare_Desiderii,看了上面的定义之后,问道:请教一下SVM functional margin 为 γ ^ = y f ( x ) = y ( w T x + b ) \hat{\gamma}=yf(x)=y(w^Tx+b) γ^=yf(x)=y(wTx+b) 中的 y y y 是只取 1 1 1 和 − 1 -1 −1 吗? y y y 的唯一作用就是确保functional margin的非负性?真是这样的么?当然不是,详情请见本文评论下第43楼)
换言之,在进行分类的时候,遇到一个新的数据点 x x x,将 x x x 代入 f ( x ) f(x) f(x) 中,如果 f ( x ) < 0 f(x) < 0 f(x)<0 则将 x x x 的类别赋为 − 1 -1 −1,如果 f ( x ) > 0 f(x) >0 f(x)>0 则将 x x x 的类别赋为 1 1 1 。
接下来的问题是,如何确定这个超平面呢?从直观上而言,这个超平面应该是最适合分开两类数据的直线。而判定“最适合”的标准就是这条直线离直线两边的数据的间隔最大。所以,得寻找有着最大间隔的超平面。
在超平面 w ∗ x + b = 0 w*x+b=0 w∗x+b=0 确定的情况下, ∣ w ∗ x + b ∣ |w*x+b| ∣w∗x+b∣ 能够表示点 x x x 到距离超平面的远近,而通过观察 w ∗ x + b w*x+b w∗x+b 的符号与类标记 y y y 的符号是否一致可判断分类是否正确,所以,可以用 y ∗ ( w ∗ x + b ) y*(w*x+b) y∗(w∗x+b) 的正负性来判定或表示分类的正确性 。于此,我们便引出了函数间隔(functional margin)的概念。
定义 函数间隔(用 γ ^ \hat{\gamma} γ^ 表示) 为: γ ^ = y f ( x ) = y ( w T x + b ) \hat{\gamma} = yf(x)= y(w^Tx + b) γ^=yf(x)=y(wTx+b)
而超平面 ( w , b ) (w,b) (w,b) 关于 T T T 中所有样本点 ( x i , y i ) (x_i,y_i) (xi,yi) 的函数间隔最小值(其中, x x x 是特征, y y y 是结果标签, i i i 表示第 i i i 个样本),便为超平面 ( w , b ) (w, b) (w,b) 关于训练数据集 T T T 的 函数间隔: γ ^ = min γ ^ i ( i = 1 , 2 , ⋯   , n ) \hat{\gamma} = \min\hat{\gamma} \ i \ (i=1,2,\cdots,n) γ^=minγ^ i (i=1,2,⋯,n)
但这样定义的函数间隔有问题,即如果成比例的改变 w w w 和 b b b(如将它们改成 2 w 2w 2w 和 2 b 2b 2b),则函数间隔的值 f ( x ) f(x) f(x) 却变成了原来的 2 倍(虽然此时超平面没有改变),所以只有函数间隔还远远不够。
事实上,我们可以对法向量 w w w 加些约束条件,从而引出真正定义点到超平面的距离 – 几何间隔(Geometrical margin) 的概念。
假定对于一个点 x x x ,令其垂直投影到超平面上的对应点为 x 0 x_0 x0 , w w w 是垂直于超平面的一个向量,为样本 x x x 到超平面的距离,如下图所示:
根据平面几何知识,有 x = x 0 + γ w ∣ ∣ w ∣ ∣ x=x_0+\gamma\frac{w}{||w||} x=x0+γ∣∣w∣∣w,其中 ∣ ∣ w ∣ ∣ ||w|| ∣∣w∣∣ 为 w w w 的二阶范数(范数是一个类似于模的表示长度的概念), w ∣ ∣ w ∣ ∣ \frac{w}{||w||} ∣∣w∣∣w 是单位向量(一个向量除以它的模称之为单位向量)。又由于 x 0 x_0 x0 是超平面上的点,满足 f ( x 0 ) = 0 f(x_0)=0 f(x0)=0,代入超平面的方程 w T x + b = 0 w^Tx + b=0 wTx+b=0,可得 w T x 0 + b = 0 w^Tx_0 + b=0 wTx0+b=0,即 w T x = − b w^Tx = -b wTx=−b。随即让此式 x = x 0 + γ w ∣ ∣ w ∣ ∣ x=x_0+\gamma\frac{w}{||w||} x=x0+γ∣∣w∣∣w 的两边同时乘以 w T w^T wT,再根据 w T x 0 + b = 0 w^Tx_0 + b=0 wTx0+b=0 和 w T w = ∣ ∣ w ∣ ∣ 2 w^Tw = {||w||}^2 wTw=∣∣w∣∣2,即可算出 γ \gamma γ: γ = w T x + b ∣ ∣ w ∣ ∣ = f ( x ) ∣ ∣ w ∣ ∣ \gamma = \frac{w^Tx + b}{||w||}=\frac{f(x)}{||w||} γ=∣∣w∣∣wTx+b=∣∣w∣∣f(x),为了得到 γ \gamma γ 的绝对值,令 γ \gamma γ 乘上对应的类别 y y y,即可得出几何间隔(用 γ ~ \tilde{\gamma} γ~ 表示)的定义: γ ~ = y γ = γ ^ ∣ ∣ w ∣ ∣ \tilde{\gamma} = y\gamma = \frac{\hat\gamma}{||w||} γ~=yγ=∣∣w∣∣γ^
函数间隔和几何间隔详细的转换过程如下:
x = x 0 + γ w ∣ ∣ w ∣ ∣ w T x = w T x 0 + γ w T w ∣ ∣ w ∣ ∣ ( w T w = ∣ ∣ w ∣ ∣ 2 ) w T x = w T x 0 + γ ∣ ∣ w ∣ ∣ 2 ∣ ∣ w ∣ ∣ ( w T x 0 = − b ) − b = w T x 0 + γ ∣ ∣ w ∣ ∣ γ = − b + w T x 0 ∣ ∣ w ∣ ∣ γ = b + w T x 0 ∣ ∣ w ∣ ∣ ( γ ^ = y f ( x ) = y ( w T x + b ) ) y γ = y ( b + w T x 0 ) ∣ ∣ w ∣ ∣ γ ~ = y γ = γ ^ ∣ ∣ w ∣ ∣ \begin{aligned} x & = x_0+\gamma\frac{w}{||w||} \\ w^Tx & = w^Tx_0+\gamma\frac{w^Tw}{||w||} \left ( w^Tw = {||w||}^2 \right ) \\ w^Tx & = w^Tx_0+\gamma\frac{{||w||}^2}{||w||} \left (w^Tx_0 = -b \right ) \\ -b & = w^Tx_0+\gamma{||w||} \\ \gamma & = -\frac{b+w^Tx_0}{||w||} \\ \gamma & = \frac{b+w^Tx_0}{||w||} \left (\hat{\gamma} = yf(x)= y(w^Tx + b) \right )\\ y\gamma & = \frac{y(b+w^Tx_0)}{||w||} \\ \tilde{\gamma} & = y\gamma = \frac{\hat\gamma}{||w||} \end{aligned} xwTxwTx−bγγyγγ~=x0+γ∣∣w∣∣w=wTx0+γ∣∣w∣∣wTw(wTw=∣∣w∣∣2)=wTx0+γ∣∣w∣∣∣∣w∣∣2(wTx0=−b)=wTx0+γ∣∣w∣∣=−∣∣w∣∣b+wTx0=∣∣w∣∣b+wTx0(γ^=yf(x)=y(wTx+b))=∣∣w∣∣y(b+wTx0)=yγ=∣∣w∣∣γ^
其中, γ \gamma γ 表示样本到超平面的距离, γ ^ \hat{\gamma} γ^ 表示函数间隔, γ ~ \tilde{\gamma} γ~ 表示几何间隔
从上述函数间隔和几何间隔的定义可以看出:几何间隔就是函数间隔除以 ∣ ∣ w ∣ ∣ ||w|| ∣∣w∣∣,而且函数间隔 y ∗ ( w x + b ) = y ∗ f ( x ) y*(wx+b) = y*f(x) y∗(wx+b)=y∗f(x) 实际上就是 ∣ f ( x ) ∣ |f(x)| ∣f(x)∣,只是人为定义的一个间隔度量,而几何间隔 ∣ f ( x ) ∣ ∣ ∣ w ∣ ∣ \frac{|f(x)|}{||w||} ∣∣w∣∣∣f(x)∣才是直观上的点到超平面的距离。
对一个数据点进行分类,当超平面离数据点的“间隔”越大,分类的确信度(confidence)也越大。所以,为了使得分类的确信度尽量高,需要让所选择的超平面能够最大化这个“间隔”值。这个间隔就是下图中的 Gap 的一半。
通过由前面的分析可知:函数间隔不适合用来最大化间隔值,因为在超平面固定以后,可以等比例地缩放 w w w 的长度和 b b b 的值,这样可以使得 f ( x ) = w T x + b f(x) = w^Tx + b f(x)=wTx+b 的值任意大,亦即函数间隔 γ \gamma γ 可以在超平面保持不变的情况下被取得任意大。但几何间隔因为除上了 ∣ ∣ w ∣ ∣ ||w|| ∣∣w∣∣,使得在缩放 w w w 和 b b b 的时候几何间隔 γ \gamma γ 的值是不会改变的,它只随着超平面的变动而变动,因此,这是更加合适的一个间隔。换言之,这里要找的最大间隔分类超平面中的“间隔”指的是几何间隔。
于是最大间隔分类器(maximum margin classifier)的目标函数可以定义为: max γ ^ \max \hat\gamma maxγ^
同时需满足一些条件,根据间隔的定义,有 y i ( w T x i + b ) = γ ^ i ≥ γ , i = 1 , 2 , ⋯   , n y_i(w^Tx_i+b)=\hat{\gamma}_i \geq\gamma, \ i=1,2,\cdots ,n yi(wTxi+b)=γ^i≥γ, i=1,2,⋯,n
其中, s . t . s.t. s.t. 即 subject to 的意思,它导出的是约束条件。
回顾下几何间隔的定义 γ = y γ = γ ^ ∣ ∣ w ∣ ∣ \gamma = y\gamma = \frac{\hat\gamma}{||w||} γ=yγ=∣∣w∣∣γ^,可知:如果令函数间隔 γ ^ = 1 \hat{\gamma} = 1 γ^=1(之所以令等于1,是为了方便推导和优化,且这样做对目标函数的优化没有影响,至于为什么,请见本文评论下第42楼回复),则有 γ ^ = 1 ∣ ∣ w ∣ ∣ \hat{\gamma} = \frac{1}{||w||} γ^=∣∣w∣∣1 且 y i ( w T + b ) ≥ 1 , i = 1 , 2 , ⋯   , n y_i(w^T + b) \geq 1, \ i=1,2,\cdots ,n yi(wT+b)≥1, i=1,2,⋯,n,从而上述目标函数转化成了 max 1 ∣ ∣ w ∣ ∣ s . t . , y i ( w T + b ) ≥ 1 , i = 1 , 2 , ⋯   , n \max \frac{1}{||w||} \ s.t.,\ y_i(w^T + b) \geq 1, \ i=1,2,\cdots ,n max∣∣w∣∣1 s.t., yi(wT+b)≥1, i=1,2,⋯,n
相当于在相应的约束条件 y i ( w T + b ) ≥ 1 , i = 1 , 2 , ⋯   , n y_i(w^T + b) \geq 1, \ i=1,2,\cdots ,n yi(wT+b)≥1, i=1,2,⋯,n 下,最大化这个 1 ∣ ∣ w ∣ ∣ \frac{1}{||w||} ∣∣w∣∣1 值,而 1 ∣ ∣ w ∣ ∣ \frac{1}{||w||} ∣∣w∣∣1 便是几何间隔。
如下图所示,中间的实线便是寻找到的最优超平面(Optimal Hyper Plane),其到两条虚线边界的距离相等,这个距离便是几何间隔 γ ^ \hat\gamma γ^,两条虚线间隔边界之间的距离等于 2 γ ^ 2\hat\gamma 2γ^,而虚线间隔边界上的点则是支持向量。由于这些支持向量刚好在虚线间隔边界上,所以它们满足 y ( w T x + b ) = 1 y(w^Tx +b) =1 y(wTx+b)=1(还记得我们把 functional margin 定为 1 1 1 了吗?上节中:处于方便推导和优化的目的,我们可以令 γ ^ = 1 \hat\gamma =1 γ^=1),而对于所有不是支持向量的点,则显然有 y ( w T x + b ) > 1 y(w^Tx +b) >1 y(wTx+b)>1。
OK,到此为止,算是了解到了SVM的第一层,对于那些只关心怎么用SVM的朋友便已足够,不必再更进一层深究其更深的原理。
接着考虑之前得到的目标函数:
max 1 ∣ ∣ w ∣ ∣ s . t . , y i ( w T + b ) ≥ 1 , i = 1 , 2 , ⋯   , n \max \frac{1}{||w||} \ s.t.,\ y_i(w^T + b) \geq 1, \ i=1,2,\cdots ,n max∣∣w∣∣1 s.t., yi(wT+b)≥1, i=1,2,⋯,n
由于求 1 ∣ ∣ w ∣ ∣ \frac{1}{||w||} ∣∣w∣∣1 的最大值相当于求 1 2 ∣ ∣ w ∣ ∣ 2 \frac{1}{2}{||w||}^2 21∣∣w∣∣2 的最小值,所以上述目标函数等价于( w w w 由分母变成分子,从而也有原来的 max \max max 问题变为 min \min min 问题,很明显,两者问题等价):
min 1 2 ∣ ∣ w ∣ ∣ 2 s . t . , y i ( w T + b ) ≥ 1 , i = 1 , 2 , ⋯   , n \min \frac{1}{2}{||w||}^2 \ s.t.,\ y_i(w^T + b) \geq 1, \ i=1,2,\cdots ,n min21∣∣w∣∣2 s.t., yi(wT+b)≥1, i=1,2,⋯,n
因为现在的目标函数是二次的,约束条件是线性的,所以它是一个凸二次规划问题。这个问题可以用现成的QP (Quadratic Programming) 优化包进行求解。一言以蔽之:在一定的约束条件下,目标最优,损失最小。
此外,由于这个问题的特殊结构,还可以通过拉格朗日对偶性(Lagrange Duality)变换到对偶变量 (Dual Variable) 的优化问题,即通过求解与原问题等价的对偶问题(Dual Problem)得到原始问题的最优解,这就是线性可分条件下支持向量机的对偶算法,这样做的优点在于:一者对偶问题往往更容易求解;二者可以自然的引入核函数,进而推广到非线性分类问题。
那什么是拉格朗日对偶性呢?简单来讲,通过给每一个约束条件加上一个拉格朗日乘子(Lagrange Multiplier),定义拉格朗日函数(通过拉格朗日函数将约束条件融合到目标函数里去,从而只用一个函数表达式便能清楚的表达出我们的问题):
L ( w , b , α ) = 1 2 ∣ ∣ w ∣ ∣ 2 − ∑ i = 1 n α ( y i ( w T x i + b ) − 1 ) \mathcal{L}(w,b,\alpha) = \frac{1}{2}{||w||}^2 - \sum_{i=1}^{n}\alpha \left(y_i \left( w^Tx_i +b\right ) -1\right ) L(w,b,α)=21∣∣w∣∣2−∑i=1nα(yi(wTxi+b)−1)
然后令
θ ( w ) = max α i ≥ 0 L ( w , b , α ) \theta (w) = \max \limits_{\alpha_i \geq \ 0} \mathcal{L}(w,b,\alpha) θ(w)=αi≥ 0maxL(w,b,α)
容易验证,当某个约束条件不满足时,例如 y i ( w T x i + b ) < 1 y_i \left( w^Tx_i +b\right ) < 1 yi(wTxi+b)<1,那么显然有 θ ( w ) = ∞ \theta (w) = \infty θ(w)=∞(只要令即可 α i = ∞ \alpha_i = \infty αi=∞)。而当所有约束条件都满足时,则最优值为 θ ( w ) = 1 2 ∣ ∣ w ∣ ∣ 2 \theta (w) = \frac{1}{2}{||w||}^2 θ(w)=21∣∣w∣∣2 ,亦即最初要最小化的量。
因此,在要求约束条件得到满足的情况下最小化 θ ( w ) = 1 2 ∣ ∣ w ∣ ∣ 2 \theta (w) = \frac{1}{2}{||w||}^2 θ(w)=21∣∣w∣∣2,实际上等价于直接最小化 θ ( w ) \theta (w) θ(w) (当然,这里也有约束条件,就是≥0,i=1,…,n) ,因为如果约束条件没有得到满足, θ ( w ) \theta (w) θ(w) 会等于无穷大,自然不会是我们所要求的最小值。
具体写出来,目标函数变成了:
min w , b θ ( w ) = min w , b max α i ≥ 0 L ( w , b , α ) = p ∗ \min \limits_{w,\ b} \theta(w) = \min \limits_{w,\ b}\max \limits_{\alpha_i \geq \ 0}\mathcal{L}(w,b,\alpha) = p^* w, bminθ(w)=w, bminαi≥ 0maxL(w,b,α)=p∗
这里用 p ∗ p^* p∗ 表示这个问题的最优值,且和最初的问题是等价的。如果直接求解,那么一上来便得面对 w w w 和 b b b 两个参数,而 α i \alpha_i αi 又是不等式约束,这个求解过程不好做。不妨把最小和最大的位置交换一下,变成:
max α i ≥ 0 min w , b L ( w , b , α ) = d ∗ \max \limits_{\alpha_i \geq \ 0}\min \limits_{w,\ b}\mathcal{L}(w,b,\alpha) = d^* αi≥ 0maxw, bminL(w,b,α)=d∗
交换以后的新问题是原始问题的对偶问题,这个新问题的最优值用 d ∗ d^* d∗ 来表示。而且有 d ∗ ≤ p ∗ d^* \leq p^* d∗≤p∗,在满足某些条件的情况下,这两者相等,这个时候就可以通过求解对偶问题来间接地求解原始问题。
换言之,之所以从 min max \min\max minmax 的原始问题 p ∗ p^* p∗,转化为 min max \min\max minmax 的对偶问题 d ∗ d^* d∗,一者因为 d ∗ d^* d∗ 是 p ∗ p^* p∗ 的近似解,二者,转化为对偶问题后,更容易求解。
下面可以先求 L L L 对 w w w、 b b b 的极小,再求 L \mathcal{L} L 对 α \alpha α 的极大。
上文中提到 d ∗ ≤ p ∗ d^* \leq p^* d∗≤p∗ 在满足某些条件的情况下,两者等价”,这所谓的“满足某些条件”就是要满足KKT条件。
勘误:经读者 qq_28543029 指出,这里的条件不应该是KKT条件,要让两者等价需满足strong duality (强对偶),而后有学者在强对偶下提出了KKT条件,且KKT条件的成立要满足constraint qualifications,而constraint qualifications之一就是Slater条件。所谓Slater 条件,即指:凸优化问题,如果存在一个点 x x x,使得所有等式约束都成立,并且所有不等式约束都严格成立(即取严格不等号,而非等号),则满足Slater 条件。对于此处,Slater 条件成立,所以 d ∗ ≤ p ∗ d^* \leq p^* d∗≤p∗ 可以取等号。
一般地,一个最优化数学模型能够表示成下列标准形式:
min . f ( x ) s . t . h j ( x ) = 0 , j = 1 , ⋯   , p , g k ( x ) ≤ 0 , k = 1 , ⋯   , p , x ∈ X ⊂ R n \begin{aligned} \min.f(x) & \\ s.t. \ \ \ h_j(x) & = 0, j =1, \cdots, p, \\ g_k(x) & \leq 0, k =1, \cdots, p, \\ x \in X & \subset \mathbb{R}^n \end{aligned} min.f(x)s.t. hj(x)gk(x)x∈X=0,j=1,⋯,p,≤0,k=1,⋯,p,⊂Rn
其中, f ( x ) f(x) f(x) 是需要最小化的函数, h ( x ) h(x) h(x) 是等式约束, g ( x ) g(x) g(x) 是不等式约束, p p p 和 q q q 分别为等式约束和不等式约束的数量。
同时,得明白以下两点:
而KKT条件就是指上面最优化数学模型的标准形式中的最小点 x* 必须满足下面的条件:
经过论证,我们这里的问题是满足 KKT 条件的(首先已经满足Slater条件,再者f和gi也都是可微的,即L对w和b都可导),因此现在我们便转化为求解第二个问题。
也就是说,原始问题通过满足KKT条件,已经转化成了对偶问题。而求解这个对偶学习问题,分为3个步骤:
(1)、首先固定 α \alpha α,要让 L L L 关于 w 和 b 最小化,我们分别对w,b求偏导数,即令 ∂L/∂w 和 ∂L/∂b 等于零(对w求导结果的解释请看本文评论下第45楼回复):
∂ L ∂ w = 0 ⇒ w = ∑ i = 1 n α i y i x i \frac{\partial \mathcal{L}}{\partial w} = 0 \Rightarrow w = \sum_{i=1}^{n}\alpha_iy_ix_i ∂w∂L=0⇒w=∑i=1nαiyixi
∂ L ∂ b = 0 ⇒ ∑ i = 1 n α i y i \frac{\partial \mathcal{L}}{\partial b} = 0 \Rightarrow \sum_{i=1}^{n}\alpha_iy_i ∂b∂L=0⇒∑i=1nαiyi
将以上结果代入之前的L:
L ( w , b , α ) = 1 2 ∣ ∣ w ∣ ∣ 2 − ∑ i = 1 n α ( y i ( w T x i + b ) − 1 ) \mathcal{L}(w,b,\alpha) = \frac{1}{2}{||w||}^2 - \sum_{i=1}^{n}\alpha \left(y_i \left( w^Tx_i +b\right ) -1\right ) L(w,b,α)=21∣∣w∣∣2−∑i=1nα(yi(wTxi+b)−1)
得到:
L ( w , b , α ) = 1 2 ∑ i , j = 1 n α i α j y i y j x i T x j − ∑ i , j = 1 n α i α j y i y j x i T x j − b ∑ i = 1 n α i y i + ∑ i = 1 n α i = ∑ i = 1 n α i − 1 2 ∑ i , j = 1 n α i α j y i y j x i T x j \begin{aligned} \mathcal{L}(w,b,\alpha) & = \frac{1}{2}\sum_{i,j=1}^{n}\alpha_i\alpha_jy_iy_jx_i^Tx_j - \sum_{i,j=1}^{n}\alpha_i\alpha_jy_iy_jx_i^Tx_j -b\sum_{i=1}^{n}\alpha_iy_i + \sum_{i=1}^{n}\alpha_i \\ & = \sum_{i=1}^{n}\alpha_i - \frac{1}{2}\sum_{i,j=1}^{n}\alpha_i\alpha_jy_iy_jx_i^Tx_j \end{aligned} L(w,b,α)=21i,j=1∑nαiαjyiyjxiTxj−i,j=1∑nαiαjyiyjxiTxj−bi=1∑nαiyi+i=1∑nαi=i=1∑nαi−21i,j=1∑nαiαjyiyjxiTxj
提醒:有读者可能会问上述推导过程如何而来?说实话,其具体推导过程是比较复杂的,如下图所示:
L ( w , b , α ) = 1 2 ∣ ∣ w ∣ ∣ 2 − ∑ i = 1 m α i [ y ( i ) ( w T x ( i ) + b ) − 1 ] = 1 2 w T w − ∑ i = 1 m α i y ( i ) w T x ( i ) − ∑ i = 1 m α i y ( i ) b + ∑ i = 1 m α i = 1 2 w T ∑ i = 1 m α i y ( i ) x ( i ) − ∑ i = 1 m α i y ( i ) w T x ( i ) − ∑ i = 1 m α i y ( i ) b + ∑ i = 1 m α i = 1 2 w T ∑ i = 1 m α i y ( i ) x ( i ) − w T ∑ i = 1 m α i y ( i ) x ( i ) − ∑ i = 1 m α i y ( i ) b + ∑ i = 1 m α i = − 1 2 w T ∑ i = 1 m α i y ( i ) x ( i ) − ∑ i = 1 m α i y ( i ) b + ∑ i = 1 m α i = − 1 2 w T ∑ i = 1 m α i y ( i ) x ( i ) − b ∑ i = 1 m α i y ( i ) + ∑ i = 1 m α i = − 1 2 ( ∑ i = 1 m α i y ( i ) x ( i ) ) T ∑ i = 1 m α i y ( i ) x ( i ) − b ∑ i = 1 m α i y ( i ) + ∑ i = 1 m α i = − 1 2 ∑ i = 1 m α i y ( i ) ( x ( i ) ) T ∑ i = 1 m α i y ( i ) x ( i ) − b ∑ i = 1 m α i y ( i ) + ∑ i = 1 m α i = − 1 2 ∑ i = 1 , j = 1 m α i y ( i ) ( x ( i ) ) T α j y ( j ) x ( j ) − b ∑ i = 1 m α i y ( i ) + ∑ i = 1 m α i \begin{aligned} \mathcal{L}(w,b,\alpha) &= \frac{1}{2}{||w||}^2 - \sum_{i=1}^{m}\alpha_i \left[y^{(i)} \left( w^Tx^{(i)} +b\right ) -1\right ] \\ & = \frac{1}{2}w^Tw - \sum_{i=1}^{m}\alpha_iy^{(i)}w^Tx^{(i)}- \sum_{i=1}^{m}\alpha_iy^{(i)}b + \sum_{i=1}^{m}\alpha_i \\ & = \frac{1}{2}w^T \sum_{i=1}^{m}\alpha_iy^{(i)}x^{(i)}- \sum_{i=1}^{m}\alpha_iy^{(i)}w^Tx^{(i)} - \sum_{i=1}^{m}\alpha_iy^{(i)}b + \sum_{i=1}^{m}\alpha_i\\ & = \frac{1}{2}w^T \sum_{i=1}^{m}\alpha_iy^{(i)}x^{(i)}- w^T\sum_{i=1}^{m}\alpha_iy^{(i)}x^{(i)} - \sum_{i=1}^{m}\alpha_iy^{(i)}b + \sum_{i=1}^{m}\alpha_i \\ & = - \frac{1}{2}w^T \sum_{i=1}^{m}\alpha_iy^{(i)}x^{(i)} - \sum_{i=1}^{m}\alpha_iy^{(i)}b + \sum_{i=1}^{m}\alpha_i \\ & = - \frac{1}{2}w^T \sum_{i=1}^{m}\alpha_iy^{(i)}x^{(i)} - b\sum_{i=1}^{m}\alpha_iy^{(i)} + \sum_{i=1}^{m}\alpha_i\\ & = - \frac{1}{2} \left( \sum_{i=1}^{m}\alpha_iy^{(i)}x^{(i)}\right )^T\sum_{i=1}^{m}\alpha_iy^{(i)}x^{(i)}- b\sum_{i=1}^{m}\alpha_iy^{(i)} + \sum_{i=1}^{m}\alpha_i\\ & = - \frac{1}{2}\sum_{i=1}^{m}\alpha_iy^{(i)}(x^{(i)})^T\sum_{i=1}^{m}\alpha_iy^{(i)}x^{(i)} - b\sum_{i=1}^{m}\alpha_iy^{(i)} + \sum_{i=1}^{m}\alpha_i \\ & = - \frac{1}{2}\sum_{i=1,j=1}^{m}\alpha_iy^{(i)}(x^{(i)})^T\alpha_jy^{(j)}x^{(j)} - b\sum_{i=1}^{m}\alpha_iy^{(i)} + \sum_{i=1}^{m}\alpha_i \\ \end{aligned} L(w,b,α)=21∣∣w∣∣2−i=1∑mαi[y(i)(wTx(i)+b)−1]=21wTw−i=1∑mαiy(i)wTx(i)−i=1∑mαiy(i)b+i=1∑mαi=21wTi=1∑mαiy(i)x(i)−i=1∑mαiy(i)wTx(i)−i=1∑mαiy(i)b+i=1∑mαi=21wTi=1∑mαiy(i)x(i)−wTi=1∑mαiy(i)x(i)−i=1∑mαiy(i)b+i=1∑mαi=−21wTi=1∑mαiy(i)x(i)−i=1∑mαiy(i)b+i=1∑mαi=−21wTi=1∑mαiy(i)x(i)−bi=1∑mαiy(i)+i=1∑mαi=−21(i=1∑mαiy(i)x(i))Ti=1∑mαiy(i)x(i)−bi=1∑mαiy(i)+i=1∑mαi=−21i=1∑mαiy(i)(x(i))Ti=1∑mαiy(i)x(i)−bi=1∑mαiy(i)+i=1∑mαi=−21i=1,j=1∑mαiy(i)(x(i))Tαjy(j)x(j)−bi=1∑mαiy(i)+i=1∑mαi
最后,得到:
L ( w , b , α ) = 1 2 ∑ i , j = 1 n α i α j y i y j x i T x j − ∑ i , j = 1 n α i α j y i y j x i T x j − b ∑ i = 1 n α i y i + ∑ i = 1 n α i = ∑ i = 1 n α i − 1 2 ∑ i , j = 1 n α i α j y i y j x i T x j \begin{aligned} \mathcal{L}(w,b,\alpha) & = \frac{1}{2}\sum_{i,j=1}^{n}\alpha_i\alpha_jy_iy_jx_i^Tx_j - \sum_{i,j=1}^{n}\alpha_i\alpha_jy_iy_jx_i^Tx_j -b\sum_{i=1}^{n}\alpha_iy_i + \sum_{i=1}^{n}\alpha_i \\ & = \sum_{i=1}^{n}\alpha_i - \frac{1}{2}\sum_{i,j=1}^{n}\alpha_i\alpha_jy_iy_jx_i^Tx_j \end{aligned} L(w,b,α)=21i,j=1∑nαiαjyiyjxiTxj−i,j=1∑nαiαjyiyjxiTxj−bi=1∑nαiyi+i=1∑nαi=i=1∑nαi−21i,j=1∑nαiαjyiyjxiTxj
如 jerrylead所说:“倒数第4步”推导到“倒数第3步”使用了线性代数的转置运算,由于 a i a_i ai 和 y i y_i yi 都是实数,因此转置后与自身一样。“倒数第3步”推导到“倒数第2步”使用了 ( a + b + c + ⋯   ) ( a + b + c + ⋯   ) = a a + a b + a c + b a + b b + b c + ⋯ (a+b+c+\cdots)(a+b+c+\cdots)=aa+ab+ac+ba+bb+bc+ \cdots (a+b+c+⋯)(a+b+c+⋯)=aa+ab+ac+ba+bb+bc+⋯ 的乘法运算法则。最后一步是上一步的顺序调整。
从上面的最后一个式子,我们可以看出,此时的拉格朗日函数只包含了一个变量,那就是 α i \alpha_i αi(求出了便能求出 w w w 和 b b b,由此可见,上文第1.2节提出来的核心问题:分类函数 f ( x ) = w T + b f(x)=w^T+b f(x)=wT+b 也就可以轻而易举的求出来了)。
(2)、求对 α \alpha α 的极大,即是关于对偶问题的最优化问题。经过上面第一个步骤的求 w w w 和 b b b,得到的拉格朗日函数式子已经没有了变量 w w w , b b b,只有 α \alpha α。从上面的式子得到:
max α ∑ i = 1 n α i − 1 2 ∑ i , j = 1 n α i α j y i y j x i T x j s . t . , a i ≥ 0 , i = 1 , ⋯   , n ∑ i = 1 n α i y i = 0 \max_{\alpha}\sum_{i=1}^{n}\alpha_i - \frac{1} {2}\sum_{i,j=1}^{n}\alpha_i\alpha_jy_iy_jx_i^Tx_j \\ s.t.,a_i\geq0,i=1, \cdots,n \\ \sum_{i=1}^{n}\alpha_iy_i=0 αmaxi=1∑nαi−21i,j=1∑nαiαjyiyjxiTxjs.t.,ai≥0,i=1,⋯,ni=1∑nαiyi=0
这样,求出了 α i \alpha_i αi,根据 w = ∑ i = 1 m α i y ( i ) x ( i ) w=\sum_{i=1}^{m}\alpha_iy^{(i)}x^{(i)} w=∑i=1mαiy(i)x(i), 即可求出 w w w,然后通过 b ∗ = m a x i : y ( i ) = − 1 w ∗ T x ( i ) + m i n i : y ( i ) = 1 w ∗ T x ( i ) 2 b^*=\frac{max_{i:y^{(i)}=-1}w^{*T}x^{(i)}+min_{i:y^{(i)}=1}w^{*T}x^{(i)}}{2} b∗=2maxi:y(i)=−1w∗Tx(i)+mini:y(i)=1w∗Tx(i),即可求出 b b b,最终得出分离超平面和分类决策函数。
(3)在求得 L ( w , b , α ) \mathcal{L}(w,b,\alpha) L(w,b,α) 关于 w w w 和 b b b 最小化,以及对 α \alpha α 的极大之后,最后一步则可以利用SMO算法求解对偶问题中的拉格朗日乘子 α \alpha α 。
上述式子要解决的是在参数 { α 1 , α 2 , ⋯   , α n } \left \{\alpha_1,\alpha_2, \cdots,\alpha_n\right \} {α1,α2,⋯,αn} 上求最大值W的问题,至于 x ( i ) x^{(i)} x(i) 和 y ( i ) y^{(i)} y(i) 都是已知数。要了解这个SMO算法是如何推导的,请跳到下文第3.5节、SMO算法。
到目前为止,我们的 SVM 还比较弱,只能处理线性的情况,下面我们将引入核函数,进而推广到非线性分类问题。
OK,为过渡到下节2.2节所介绍的核函数,让我们再来看看上述推导过程中得到的一些有趣的形式。首先就是关于我们的 hyper plane ,对于一个数据点 x x x 进行分类,实际上是通过把 x x x 带入到 f ( x ) = w T x + b f(x) = w^Tx + b f(x)=wTx+b 算出结果然后根据其正负号来进行类别划分的。而前面的推导中我们得到 w = ∑ i = 1 n α i y i x i w=\sum_{i=1}^n \alpha_iy_ix_i w=∑i=1nαiyixi
因此分类函数为:
f ( x ) = ( ∑ i = 1 n α i y i x i ) T x + b = ∑ i = 1 n α i y i ⟨ x i , x ⟩ + b \begin{aligned} f(x) & = \left( \sum_{i=1}^n \alpha_iy_ix_i\right)^Tx+b \\ & = \sum_{i=1}^n \alpha_iy_i\left \langle x_i,x\right \rangle +b\\ \end{aligned} f(x)=(i=1∑nαiyixi)Tx+b=i=1∑nαiyi⟨xi,x⟩+b
这里的形式的有趣之处在于,对于新点 x x x 的预测,只需要计算它与训练数据点的内积即可(表示向量内积),这一点至关重要,是之后使用 Kernel 进行非线性推广的基本前提。此外,所谓 Supporting Vector 也在这里显示出来——事实上,所有非Supporting Vector 所对应的系数都是等于零的,因此对于新点的内积计算实际上只要针对少量的“支持向量”而不是所有的训练数据即可。
为什么非支持向量对应的等于零呢?直观上来理解的话,就是这些“后方”的点——正如我们之前分析过的一样,对超平面是没有影响的,由于分类完全有超平面决定,所以这些无关的点并不会参与分类问题的计算,因而也就不会产生任何影响了。
回忆一下我们2.1.1节中通过 Lagrange multiplier得到的目标函数:
max a i ≥ 0 L ( w , b , α ) = max a i ≥ 0 1 2 ∣ ∣ w ∣ ∣ 2 − ∑ i = 1 n α ( y i ( w T x i + b ) − 1 ) \max_{a_i \geq 0}\mathcal{L}(w,b,\alpha) = \max_{a_i \geq 0}\frac{1}{2}{||w||}^2 - \sum_{i=1}^{n}\alpha \left(y_i \left( w^Tx_i +b\right ) -1\right ) ai≥0maxL(w,b,α)=ai≥0max21∣∣w∣∣2−i=1∑nα(yi(wTxi+b)−1)
注意到如果 x i x_i xi 是支持向量的话,上式中红颜色的部分是等于 0 的(因为支持向量的 functional margin 等于 1 ),而对于非支持向量来说,functional margin 会大于 1 ,因此红颜色部分是大于零的,而 a i a_i ai 又是非负的,为了满足最大化, a i a_i ai 必须等于 0 。这也就是这些非Supporting Vector 的点的局限性。
至此,我们便得到了一个maximum margin hyper plane classifier,这就是所谓的支持向量机(Support Vector Machine)。当然,到目前为止,我们的 SVM 还比较弱,只能处理线性的情况,不过,在得到了对偶dual 形式之后,通过 Kernel 推广到非线性的情况就变成了一件非常容易的事情了(相信,你还记得本节开头所说的:“通过求解对偶问题得到最优解,这就是线性可分条件下支持向量机的对偶算法,这样做的优点在于:一者对偶问题往往更容易求解;二者可以自然的引入核函数,进而推广到非线性分类问题”)。
事实上,大部分时候数据并不是线性可分的,这个时候满足这样条件的超平面就根本不存在。在上文中,我们已经了解到了SVM处理线性可分的情况,那对于非线性的数据SVM咋处理呢?对于非线性的情况,SVM 的处理方法是选择一个核函数 K ⟨ ⋅ , ⋅ ⟩ K\left \langle ·,·\right \rangle K⟨⋅,⋅⟩,通过将数据映射到高维空间,来解决在原始空间中线性不可分的问题。
具体来说,在线性不可分的情况下,支持向量机首先在低维空间中完成计算,然后通过核函数将输入空间映射到高维特征空间,最终在高维特征空间中构造出最优分离超平面,从而把平面上本身不好分的非线性数据分开。如图所示,一堆数据在二维空间无法划分,从而映射到三维空间里划分:
而在我们遇到核函数之前,如果用原始的方法,那么在用线性学习器学习一个非线性关系,需要选择一个非线性特征集,并且将数据写成新的表达形式,这等价于应用一个固定的非线性映射,将数据映射到特征空间,在特征空间中使用线性学习器,因此,考虑的假设集是这种类型的函数:
f ( x ) = ∑ i = 1 N w i ϕ i ( x ) + b f(x) = \sum_{i=1}^{N}w_i\phi_i(x)+b f(x)=i=1∑Nwiϕi(x)+b
这里 ϕ : X → F \phi:X \rightarrow F ϕ:X→F是从输入空间到某个特征空间的映射,这意味着建立非线性学习器分为两步:
而由于对偶形式就是线性学习器的一个重要性质,这意味着假设可以表达为训练点的线性组合,因此决策规则可以用测试点和训练点的内积来表示:
f ( x ) = ∑ i = 1 ι α i y i ⟨ ϕ ( x i ) , ϕ ( x ) ⟩ + b f(x) = \sum_{i=1}^{\iota}\alpha_iy_i \left \langle \phi(x_i),\phi(x)\right \rangle +b f(x)=i=1∑ιαiyi⟨ϕ(xi),ϕ(x)⟩+b
如果有一种方式可以在特征空间中直接计算内积 ⟨ ϕ ( x i ) , ϕ ( x ) ⟩ \left \langle \phi(x_i),\phi(x)\right \rangle ⟨ϕ(xi),ϕ(x)⟩,就像在原始输入点的函数中一样,就有可能将两个步骤融合到一起建立一个非线性的学习器,这样直接计算法的方法称为核函数方法:
核是一个函数K,对所有 x , z ∈ X x,z \in X x,z∈X,满足,这里 ϕ \phi ϕ 是从 X X X 到内积特征空间 F F F 的映射。
来看个核函数的例子。如下图所示的两类数据,分别分布为两个圆圈的形状,这样的数据本身就是线性不可分的,此时咱们该如何把这两类数据分开呢(下文将会有一个相应的三维空间图)?
事实上,上图所述的这个数据集,是用两个半径不同的圆圈加上了少量的噪音生成得到的,所以,一个理想的分界应该是一个“圆圈”而不是一条线(超平面)。如果用 X 1 X_1 X1 和 X 2 X_2 X2 来表示这个二维平面的两个坐标的话,我们知道一条二次曲线(圆圈是二次曲线的一种特殊情况)的方程可以写作这样的形式:
a 1 X 1 + a 2 X 1 2 + a 3 X 2 + a 4 X 2 2 + a 5 X 1 X 2 + a 6 = 0 a_1X_1+a_2X_1^2+a_3X_2+a_4X_2^2+a_5X_1X_2+a_6=0 a1X1+a2X12+a3X2+a4X22+a5X1X2+a6=0
注意上面的形式,如果我们构造另外一个五维的空间,其中五个坐标的值分别为 Z 1 = X 1 , Z 2 = X 1 2 , Z 3 = X 2 , Z 4 = X 2 2 , Z 5 = X 1 X 2 Z_1=X_1,Z_2=X_1^2,Z_3=X_2,Z_4=X_2^2,Z_5=X_1X_2 Z1=X1,Z2=X12,Z3=X2,Z4=X22,Z5=X1X2,那么显然,上面的方程在新的坐标系下可以写作:
∑ i = 1 5 α i Z i + a 6 = 0 \sum_{i=1}^{5}\alpha_iZ_i + a_6 = 0 i=1∑5αiZi+a6=0
关于新的坐标 Z Z Z,这正是一个 hyper plane 的方程!也就是说,如果我们做一个映射 ϕ : R 2 → R 5 \phi:R_2 \rightarrow R_5 ϕ:R2→R5,将 X X X 按照上面的规则映射为 Z Z Z,那么在新的空间中原来的数据将变成线性可分的,从而使用之前我们推导的线性分类算法就可以进行处理了。这正是 Kernel 方法处理非线性问题的基本思想。
再进一步描述 Kernel 的细节之前,不妨再来看看上述例子在映射过后的直观形态。当然,你我可能无法把 5 维空间画出来,不过由于我这里生成数据的时候用了特殊的情形,所以这里的超平面实际的方程是这个样子的(圆心在 X 2 X_2 X2 轴上的一个正圆):
∑ i = 1 5 α i Z i + a 6 = 0 \sum_{i=1}^{5}\alpha_iZ_i + a_6 = 0 i=1∑5αiZi+a6=0
因此我只需要把它映射到 Z 1 = X 1 2 , Z 2 = X 2 2 , Z 3 = X 2 Z_1=X_1^2,Z_2=X_2^2,Z_3=X_2 Z1=X12,Z2=X22,Z3=X2,这样一个三维空间中即可,下图即是映射之后的结果,将坐标轴经过适当的旋转,就可以很明显地看出,数据是可以通过一个平面来分开的( pluskid:下面的 gif 动画,先用 Matlab 画出一张张图片,再用 Imagemagick 拼贴成):
核函数相当于把原来的分类函数:
f ( x ) = ∑ i = 1 n α i y i ⟨ x i , x ⟩ + b f(x) = \sum_{i=1}^{n}\alpha_iy_i \left \langle x_i,x\right \rangle +b f(x)=i=1∑nαiyi⟨xi,x⟩+b
映射成:
f ( x ) = ∑ i = 1 ι α i y i ⟨ ϕ ( x i ) , ϕ ( x ) ⟩ + b f(x) = \sum_{i=1}^{\iota}\alpha_iy_i \left \langle \phi(x_i),\phi(x)\right \rangle +b f(x)=i=1∑ιαiyi⟨ϕ(xi),ϕ(x)⟩+b
而其中的 α \alpha α 可以通过求解如下 dual 问题而得到的:
max α ∑ i = 1 n α i − 1 2 ∑ i , j = 1 n α i α j y i y j x i T x j s . t . , a i ≥ 0 , i = 1 , ⋯   , n ∑ i = 1 n α i y i = 0 \max_{\alpha}\sum_{i=1}^{n}\alpha_i - \frac{1} {2}\sum_{i,j=1}^{n}\alpha_i\alpha_jy_iy_jx_i^Tx_j \\ s.t.,a_i\geq0,i=1, \cdots,n \\ \sum_{i=1}^{n}\alpha_iy_i=0 αmaxi=1∑nαi−21i,j=1∑nαiαjyiyjxiTxjs.t.,ai≥0,i=1,⋯,ni=1∑nαiyi=0
这样一来问题就解决了吗?似乎是的:拿到非线性数据,就找一个映射 ϕ ( ⋅ ) \phi(·) ϕ(⋅),然后一股脑把原来的数据映射到新空间中,再做线性 SVM 即可。不过事实上好像并没有这么简单。
细想一下,刚才的方法是不是有问题?
这个时候,可能就需要 Kernel 出马了。
不妨还是从最开始的简单例子出发,设两个向量 x 1 = ( η 1 , η 2 ) T x_1 = (\eta_1,\eta_2)^T x1=(η1,η2)T 和 x 1 = ( ξ 1 , ξ 2 ) T x_1 = (\xi _1,\xi _2)^T x1=(ξ1,ξ2)T,而 ϕ ( ⋅ ) \phi(·) ϕ(⋅) 即是到前面说的五维空间的映射,因此映射过后的内积为:
⟨ ϕ ( x 1 ) , ϕ ( x 2 ) ⟩ = η 1 ξ 1 + η 1 2 ξ 1 2 + η 2 ξ 2 + η 2 2 ξ 2 2 + η 1 η 2 ξ 1 ξ 2 \left \langle \phi(x_1),\phi(x_2)\right \rangle = \eta_1\xi_1+\eta_1^2\xi_1^2+\eta_2\xi_2+\eta_2^2\xi_2^2+\eta_1\eta_2\xi_1\xi_2 ⟨ϕ(x1),ϕ(x2)⟩=η1ξ1+η12ξ12+η2ξ2+η22ξ22+η1η2ξ1ξ2
(公式说明:上面的这两个推导过程中,所说的前面的五维空间的映射,这里说的前面便是文中2.2.1节的所述的映射方式,回顾下之前的映射规则,再看那第一个推导,其实就是计算 x 1 x1 x1, x 2 x2 x2 各自的内积,然后相乘相加即可,第二个推导则是直接平方,去掉括号,也很容易推出来)
另外,我们又注意到:
( ⟨ x 1 , x 2 ⟩ + 1 ) 2 = 2 η 1 ξ 1 + η 1 2 ξ 1 2 + 2 η 2 ξ 2 + η 2 2 ξ 2 2 + 2 η 1 η 2 ξ 1 ξ 2 + 1 (\left \langle x_1,x_2\right \rangle +1)^2 = 2\eta_1\xi_1+\eta_1^2\xi_1^2+2\eta_2\xi_2+\eta_2^2\xi_2^2+2\eta_1\eta_2\xi_1\xi_2+1 (⟨x1,x2⟩+1)2=2η1ξ1+η12ξ12+2η2ξ2+η22ξ22+2η1η2ξ1ξ2+1
二者有很多相似的地方,实际上,我们只要把某几个维度线性缩放一下,然后再加上一个常数维度,具体来说,上面这个式子的计算结果实际上和映射
φ ( X 1 , X 2 ) = ( 2 X 1 , X 1 2 , 2 X 2 , X 2 2 , 2 X 1 X 2 , 1 ) T \varphi (X_1,X_2)=(\sqrt{2}X_1,X_1^2,\sqrt{2}X_2,X_2^2,\sqrt{2}X_1X_2,1)^T φ(X1,X2)=(2X1,X12,2X2,X22,2X1X2,1)T
之后的内积 ⟨ φ ( x 1 ) , φ ( x 2 ) ⟩ \left \langle \varphi(x_1),\varphi(x_2)\right \rangle ⟨φ(x1),φ(x2)⟩的结果是相等的,那么区别在于什么地方呢?
(公式说明:上面之中,最后的两个式子,第一个算式,是带内积的完全平方式,可以拆开,然后,通过凑一个得到,第二个算式,也是根据第一个算式凑出来的)
回忆刚才提到的映射的维度爆炸,在前一种方法已经无法计算的情况下,后一种方法却依旧能从容处理,甚至是无穷维度的情况也没有问题。
我们把这里的计算两个向量在隐式映射过后的空间中的内积的函数叫做核函数 (Kernel Function) ,例如,在刚才的例子中,我们的核函数为:
κ ( x 1 , x 2 ) = ( ⟨ x 1 , x 2 ⟩ + 1 ) 2 \kappa (x_1,x_2)=(\left \langle x_1,x_2\right \rangle +1)^2 κ(x1,x2)=(⟨x1,x2⟩+1)2
核函数能简化映射空间中的内积运算——刚好“碰巧”的是,在我们的 SVM 里需要计算的地方数据向量总是以内积的形式出现的。对比刚才我们上面写出来的式子,现在我们的分类函数为:
∑ i = 1 n α i y i κ ( x i , x j ) \sum_{i=1}^{n}\alpha_iy_i\kappa(x_i,x_j) i=1∑nαiyiκ(xi,xj)
其中 由如下 dual 问题计算而得:
max α ∑ i = 1 n α i − 1 2 ∑ i , j = 1 n α i α j y i y j x i κ ( x i , x j ) s . t . , a i ≥ 0 , i = 1 , ⋯   , n ∑ i = 1 n α i y i = 0 \max_{\alpha}\sum_{i=1}^{n}\alpha_i - \frac{1} {2}\sum_{i,j=1}^{n}\alpha_i\alpha_jy_iy_jx_i\kappa(x_i,x_j) \\ s.t.,a_i\geq0,i=1, \cdots,n \\ \sum_{i=1}^{n}\alpha_iy_i=0 αmaxi=1∑nαi−21i,j=1∑nαiαjyiyjxiκ(xi,xj)s.t.,ai≥0,i=1,⋯,ni=1∑nαiyi=0
这样一来计算的问题就算解决了,避开了直接在高维空间中进行计算,而结果却是等价的!当然,因为我们这里的例子非常简单,所以我可以手工构造出对应于 ϕ ( ⋅ ) \phi(·) ϕ(⋅) 的核函数出来,如果对于任意一个映射,想要构造出对应的核函数就很困难了。
通常人们会从一些常用的核函数中选择(根据问题和数据的不同,选择不同的参数,实际上就是得到了不同的核函数),例如:
上面说了这么一大堆,读者可能还是没明白核函数到底是个什么东西?我再简要概括下,即以下三点:
最后引用这里的一个例子举例说明下核函数解决非线性问题的直观效果。
假设现在你是一个农场主,圈养了一批羊群,但为预防狼群袭击羊群,你需要搭建一个篱笆来把羊群围起来。但是篱笆应该建在哪里呢?你很可能需要依据牛群和狼群的位置建立一个“分类器”,比较下图这几种不同的分类器,我们可以看到SVM完成了一个很完美的解决方案。
OK,不再做过多介绍了,对核函数有进一步兴趣的,还可以看看此文。
在本文第一节最开始讨论支持向量机的时候,我们就假定,数据是线性可分的,亦即我们可以找到一个可行的超平面将数据完全分开。后来为了处理非线性数据,在上文2.2节使用 Kernel 方法对原来的线性 SVM 进行了推广,使得非线性的的情况也能处理。虽然通过映射 ϕ ( ⋅ ) \phi(·) ϕ(⋅) 将原始数据映射到高维空间之后,能够线性分隔的概率大大增加,但是对于某些情况还是很难处理。
例如可能并不是因为数据本身是非线性结构的,而只是因为数据有噪音。对于这种偏离正常位置很远的数据点,我们称之为 outlier ,在我们原来的 SVM 模型里,outlier 的存在有可能造成很大的影响,因为超平面本身就是只有少数几个 support vector 组成的,如果这些 support vector 里又存在 outlier 的话,其影响就很大了。例如下图:
用黑圈圈起来的那个蓝点是一个 outlier ,它偏离了自己原本所应该在的那个半空间,如果直接忽略掉它的话,原来的分隔超平面还是挺好的,但是由于这个 outlier 的出现,导致分隔超平面不得不被挤歪了,变成途中黑色虚线所示(这只是一个示意图,并没有严格计算精确坐标),同时 margin 也相应变小了。当然,更严重的情况是,如果这个 outlier 再往右上移动一些距离的话,我们将无法构造出能将数据分开的超平面来。
为了处理这种情况,SVM 允许数据点在一定程度上偏离一下超平面。例如上图中,黑色实线所对应的距离,就是该 outlier 偏离的距离,如果把它移动回来,就刚好落在原来的 超平面 蓝色间隔边界上,而不会使得超平面发生变形了。
插播下一位读者@Copper_PKU的理解:“换言之,在有松弛的情况下outline点也属于支持向量SV,同时,对于不同的支持向量,拉格朗日参数的值也不同,如此篇论文《Large Scale Machine Learning》中的下图所示:
对于远离分类平面的点值为0;对于边缘上的点值在[0, 1/L]之间,其中,L为训练数据集个数,即数据集大小;对于outline数据和内部的数据值为1/L。更多请参看本文文末参考条目第51条。”
OK,继续回到咱们的问题。我们,原来的约束条件为:
y i ( w T x i + b ) ≥ 1 , i = 1 , ⋯   , n y_i(w^Tx_i +b) \geq 1, i=1,\cdots ,n yi(wTxi+b)≥1,i=1,⋯,n
现在考虑到outlier问题,约束条件变成了:
y i ( w T x i + b ) ≥ 1 − ξ i , i = 1 , ⋯   , n y_i(w^Tx_i +b) \geq 1 - \xi_i , i=1,\cdots ,n yi(wTxi+b)≥1−ξi,i=1,⋯,n
其中 ξ i ≥ 0 \xi_i \geq 0 ξi≥0 称为松弛变量 (slack variable) ,对应数据点 x i x_i xi 允许偏离的 functional margin 的量。当然,如果我们运行 ξ i \xi_i ξi 任意大的话,那任意的超平面都是符合条件的了。所以,我们在原来的目标函数后面加上一项,使得这些 ξ i \xi_i ξi 的总和也要最小:
min 1 2 ∣ ∣ w ∣ ∣ 2 + C ∑ i = 1 n ξ i \min\frac{1}{2}{||w||}^2 +\mathbb{C}\sum_{i=1}^{n}\xi_i min21∣∣w∣∣2+Ci=1∑nξi
其中 C \mathbb{C} C 是一个参数,用于控制目标函数中两项(“寻找 margin 最大的超平面”和“保证数据点偏差量最小”)之间的权重。注意,其中 ξ \xi ξ 是需要优化的变量(之一),而 C \mathbb{C} C 是一个事先确定好的常量。完整地写出来是这个样子:
min 1 2 ∣ ∣ w ∣ ∣ 2 + C ∑ i = 1 n ξ i s . t . , y i ( w T x i + b ) ≥ 1 − ξ i , i = 1 , ⋯   , n ξ i ≥ 0 , i = 1 , ⋯   , n \min\frac{1}{2}{||w||}^2 +\mathbb{C}\sum_{i=1}^{n}\xi_i \\ s.t.,y_i(w^Tx_i+b) \geq1- \xi_i,i=1, \cdots,n \\ \xi_i \geq0,i=1, \cdots,n min21∣∣w∣∣2+Ci=1∑nξis.t.,yi(wTxi+b)≥1−ξi,i=1,⋯,nξi≥0,i=1,⋯,n
用之前的方法将限制或约束条件加入到目标函数中,得到新的拉格朗日函数,如下所示:
分析方法和前面一样,转换为另一个问题之后,我们先让针对、和最小化:
将 带回 并化简,得到和原来一样的目标函数:
不过,由于我们得到而又有(作为 Lagrange multiplier 的条件),因此有,所以整个 dual 问题现在写作:
把前后的结果对比一下(错误修正:图中的Dual formulation中的Minimize应为maxmize):
可以看到唯一的区别就是现在 dual variable 多了一个上限 。而 Kernel 化的非线性形式也是一样的,只要把换成即可。这样一来,一个完整的,可以处理线性和非线性并能容忍噪音和 outliers 的支持向量机才终于介绍完毕了。
行文至此,可以做个小结,不准确的说,SVM它本质上即是一个分类方法,用wT+b定义分类函数,于是求w、b,为寻最大间隔,引出1/2||w||2,继而引入拉格朗日因子,化为对拉格朗日乘子a的求解(求解过程中会涉及到一系列最优化或凸二次规划等问题),如此,求w.b与求a等价,而a的求解可以用一种快速学习算法SMO,至于核函数,是为处理非线性情况,若直接映射到高维计算恐维度爆炸,故在低维计算,等效高维表现。
OK,理解到这第二层,已经能满足绝大部分人一窥SVM原理的好奇心,然对于那些想在证明层面理解SVM的则还很不够,但进入第三层理解境界之前,你必须要有比较好的数理基础和逻辑证明能力,不然你会跟我一样,吃不少苦头的。
说实话,凡是涉及到要证明的东西.理论,便一般不是怎么好惹的东西。绝大部分时候,看懂一个东西不难,但证明一个东西则需要点数学功底,进一步,证明一个东西也不是特别难,难的是从零开始发明创造这个东西的时候,则显艰难(因为任何时代,大部分人的研究所得都不过是基于前人的研究成果,前人所做的是开创性工作,而这往往是最艰难最有价值的,他们被称为真正的先驱。牛顿也曾说过,他不过是站在巨人的肩上。你,我则更是如此)。
正如陈希孺院士在他的著作《数理统计学简史》的第4章、最小二乘法中所讲:在科研上诸多观念的革新和突破是有着很多的不易的,或许某个定理在某个时期由某个人点破了,现在的我们看来一切都是理所当然,但在一切没有发现之前,可能许许多多的顶级学者毕其功于一役,耗尽一生,努力了几十年最终也是无功而返。
话休絮烦,要证明一个东西先要弄清楚它的根基在哪,即构成它的基础是哪些理论。OK,以下内容基本是上文中未讲到的一些定理的证明,包括其背后的逻辑、来源背景等东西,还是读书笔记。
本部分导述
这个感知机算法是1956年提出的,年代久远,依然影响着当今,当然,可以肯定的是,此算法亦非最优,后续会有更详尽阐述。不过,有一点,你必须清楚,这个算法是为了干嘛的:不断的训练试错以期寻找一个合适的超平面(是的,就这么简单)。
下面,举个例子。如下图所示,凭我们的直觉可以看出,图中的红线是最优超平面,蓝线则是根据感知机算法在不断的训练中,最终,若蓝线能通过不断的训练移动到红线位置上,则代表训练成功。
既然需要通过不断的训练以让蓝线最终成为最优分类超平面,那么,到底需要训练多少次呢?Novikoff定理告诉我们当间隔是正的时候感知机算法会在有限次数的迭代中收敛,也就是说Novikoff定理证明了感知机算法的收敛性,即能得到一个界,不至于无穷循环下去。
Novikoff定理:如果分类超平面存在, 仅需在序列S上迭代几次,在界为的错误次数下就可以找到分类超平面,算法停止。
这里,为扩充间隔。根据误分次数公式可知, 迭代次数与对应于扩充(包括偏置)权重的训练集的间隔有关。
顺便再解释下这个所谓的扩充间隔,即为样本到分类间隔的距离,即从引出的最大分类间隔。OK,还记得上文第1.3节开头的内容么?如下:“
在给出几何间隔的定义之前,咱们首先来看下,如上图所示,对于一个点x,令其垂直投影到超平面上的对应的为x0,由于w是垂直于超平面的一个向量,为样本x到分类间隔的距离,我们有
”
然后后续怎么推导出最大分类间隔请回到本文第一、二部分,此处不重复板书。
同时有一点得注意:感知机算法虽然可以通过简单迭代对线性可分数据生成正确分类的超平面,但不是最优效果,那怎样才能得到最优效果呢,就是上文中第一部分所讲的寻找最大分类间隔超平面。此外,Novikoff定理的证明请见这里。
Mercer定理 :如果函数K是上的映射(也就是从两个n维向量映射到实数域)。那么如果K是一个有效核函数(也称为Mercer核函数),那么当且仅当对于训练样例,其相应的核函数矩阵是对称半正定的。
要理解这个Mercer定理,先要了解什么是半正定矩阵,要了解什么是半正定矩阵,先得知道什么是正定矩阵(矩阵理论博大精深,关于矩阵推荐我正在看的一本《矩阵分析与应用》)。然后这里有一个此定理的证明,可以看下。
正如@Copper_PKU所说:核函数在SVM的分类效果中起了重要的作用,最后这里有个tutorial可以看看。
在本文1.0节有这么一句话“支持向量机(SVM)是90年代中期发展起来的基于统计学习理论的一种机器学习方法,通过寻求结构化风险最小来提高学习机泛化能力,实现经验风险和置信范围的最小化,从而达到在统计样本量较少的情况下,亦能获得良好统计规律的目的。”但初次看到的读者可能并不了解什么是结构化风险,什么又是经验风险。要了解这两个所谓的“风险”,还得又从监督学习说起。
监督学习实际上就是一个经验风险或者结构风险函数的最优化问题。风险函数度量平均意义下模型预测的好坏,模型每一次预测的好坏用损失函数来度量。它从假设空间F中选择模型f作为决策函数,对于给定的输入X,由f(X)给出相应的输出Y,这个输出的预测值f(X)与真实值Y可能一致也可能不一致,用一个损失函数来度量预测错误的程度。损失函数记为L(Y, f(X))。
常用的损失函数有以下几种(基本引用自《统计学习方法》):
如此,SVM有第二种理解,即最优化+损失最小,或如@夏粉_百度所说“可从损失函数和优化算法角度看SVM,boosting,LR等算法,可能会有不同收获”。
OK,关于更多统计学习方法的问题,请参看此文。
关于损失函数,如下文读者评论中所述:可以看看张潼的这篇《Statistical behavior and consistency of classification methods based on convex risk minimization》。各种算法中常用的损失函数基本都具有fisher一致性,优化这些损失函数得到的分类器可以看作是后验概率的“代理”。此外,张潼还有另外一篇论文《Statistical analysis of some multi-category large margin classification methods》,在多分类情况下margin loss的分析,这两篇对Boosting和SVM使用的损失函数分析的很透彻。
既然本节开始之前提到了最小二乘法,那么下面引用《正态分布的前世今生》里的内容稍微简单阐述下。
我们口头中经常说:一般来说,平均来说。如平均来说,不吸烟的健康优于吸烟者,之所以要加“平均”二字,是因为凡事皆有例外,总存在某个特别的人他吸烟但由于经常锻炼所以他的健康状况可能会优于他身边不吸烟的朋友。而最小二乘法的一个最简单的例子便是算术平均。
最小二乘法(又称最小平方法)是一种数学优化技术。它通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据,并使得这些求得的数据与实际数据之间误差的平方和为最小。用函数表示为:
使误差「所谓误差,当然是观察值与实际真实值的差量」平方和达到最小以寻求估计值的方法,就叫做最小二乘法,用最小二乘法得到的估计,叫做最小二乘估计。当然,取平方和作为目标函数只是众多可取的方法之一。
最小二乘法的一般形式可表示为:
有效的最小二乘法是勒让德在 1805 年发表的,基本思想就是认为测量中有误差,所以所有方程的累积误差为
我们求解出导致累积误差最小的参数即可:
勒让德在论文中对最小二乘法的优良性做了几点说明:
最小二乘使得误差平方和最小,并在各个方程的误差之间建立了一种平衡,从而防止某一个极端误差取得支配地位
计算中只要求偏导后求解线性方程组,计算过程明确便捷
最小二乘可以导出算术平均值作为估计值
对于最后一点,从统计学的角度来看是很重要的一个性质。推理如下:假设真值为,x1, … , xn为n次测量值, 每次测量的误差为ei = xi - ,按最小二乘法,误差累积为
求解 使达到最小,正好是算术平均。
由于算术平均是一个历经考验的方法,而以上的推理说明,算术平均是最小二乘的一个特例,所以从另一个角度说明了最小二乘方法的优良性,使我们对最小二乘法更加有信心。
最小二乘法发表之后很快得到了大家的认可接受,并迅速的在数据分析实践中被广泛使用。不过历史上又有人把最小二乘法的发明归功于高斯,这又是怎么一回事呢。高斯在1809年也发表了最小二乘法,并且声称自己已经使用这个方法多年。高斯发明了小行星定位的数学方法,并在数据分析中使用最小二乘方法进行计算,准确的预测了谷神星的位置。
说了这么多,貌似跟本文的主题SVM没啥关系呀,别急,请让我继续阐述。本质上说,最小二乘法即是一种参数估计方法,说到参数估计,咱们得从一元线性模型说起。
什么是一元线性模型呢? 请允许我引用这里的内容,先来梳理下几个基本概念:
监督学习中,如果预测的变量是离散的,我们称其为分类(如决策树,支持向量机等),如果预测的变量是连续的,我们称其为回归。
回归分析中,如果只包括一个自变量和一个因变量,且二者的关系可用一条直线近似表示,这种回归分析称为一元线性回归分析。
如果回归分析中包括两个或两个以上的自变量,且因变量和自变量之间是线性关系,则称为多元线性回归分析。
对于二维空间线性是一条直线;对于三维空间线性是一个平面,对于多维空间线性是一个超平面…
对于一元线性回归模型, 假设从总体中获取了n组观察值(X1,Y1),(X2,Y2), …,(Xn,Yn)。对于平面中的这n个点,可以使用无数条曲线来拟合。要求样本回归函数尽可能好地拟合这组值。综合起来看,这条直线处于样本数据的中心位置最合理。
选择最佳拟合曲线的标准可以确定为:使总的拟合误差(即总残差)达到最小。有以下三个标准可以选择:
用“残差和最小”确定直线位置是一个途径。但很快发现计算“残差和”存在相互抵消的问题。
用“残差绝对值和最小”确定直线位置也是一个途径。但绝对值的计算比较麻烦。
最小二乘法的原则是以“残差平方和最小”确定直线位置。用最小二乘法除了计算比较方便外,得到的估计量还具有优良特性。这种方法对异常值非常敏感。
最常用的是普通最小二乘法( Ordinary Least Square,OLS):所选择的回归模型应该使所有观察值的残差平方和达到最小,即采用平方损失函数。
我们定义样本回归模型为:
其中ei为样本(Xi, Yi)的误差。
接着,定义平方损失函数Q:
则通过Q最小确定这条直线,即确定,以为变量,把它们看作是Q的函数,就变成了一个求极值的问题,可以通过求导数得到。
求Q对两个待估参数的偏导数:
根据数学知识我们知道,函数的极值点为偏导为0的点。
解得:
这就是最小二乘法的解法,就是求得平方损失函数的极值点。自此,你看到求解最小二乘法与求解SVM问题何等相似,尤其是定义损失函数,而后通过偏导求得极值。
OK,更多请参看陈希孺院士的《数理统计学简史》的第4章、最小二乘法。
在上文中,我们提到了求解对偶问题的序列最小最优化SMO算法,但并未提到其具体解法。首先看下最后悬而未决的问题:
等价于求解:
1998年,Microsoft Research的John C. Platt在论文《Sequential Minimal Optimization:A Fast Algorithm for Training Support Vector Machines》中提出针对上述问题的解法:SMO算法,它很快便成为最快的二次规划优化算法,特别是在针对线性SVM和数据稀疏时性能更优。
接下来,咱们便参考John C. Platt的这篇文章来看看SMO的解法是怎样的。
咱们首先来定义特征到结果的输出函数:
注:这个u与我们之前定义的实质是一样的。
接着,重新定义下咱们原始的优化问题,权当重新回顾,如下:
求导得到:
代入中,可得。
通过引入拉格朗日乘子转换为对偶问题后,得:
s.t:
且
注:这里得到的min函数与我们之前的max函数实质也是一样,因为把符号变下,即由min转化为max的问题,且yi也与之前的等价,yj亦如此。
经过加入松弛变量后,模型修改为:
从而最终我们的问题变为:
下面要解决的问题是:在上求上述目标函数的最小值。为了求解这些乘子,每次从中任意抽取两个乘子和,然后固定和以外的其它乘子,使得目标函数只是关于和的函数。这样,不断的从一堆乘子中任意抽取两个求解,不断的迭代求解子问题,最终达到求解原问题的目的。
而原对偶问题的子问题的目标函数可以表达为:
其中
为了解决这个子问题,首要问题便是每次如何选取和。实际上,其中一个乘子是违法KKT条件最严重的,另外一个乘子则由另一个约束条件选取。
根据KKT条件可以得出目标函数中取值的意义:
这里的还是拉格朗日乘子:
对于第1种情况,表明是正常分类,在间隔边界内部(我们知道正确分类的点);
对于第2种情况,表明了是支持向量,在间隔边界上;
对于第3种情况,表明了是在两条间隔边界之间;
而最优解需要满足KKT条件,即上述3个条件都得满足,以下几种情况出现将会出现不满足:
<=1但是 =1但是>0则是不满足的,而原本=0 对于,即第一个乘子,可以通过刚刚说的那3种不满足KKT的条件来找; 第一步选取一对和,选取方法使用启发式方法; 步骤1:先“扫描”所有乘子,把第一个违反KKT条件的作为更新对象,令为a1;
=1但是=0或者=C则表明不满足的,而原本应该是0<因此,如果假设选择的两个乘子和,它们在更新之前分别是、,更新之后分别是、,那么更新前后的值需要满足以下等式才能保证和为0的约束:
其中,是常数。
两个因子不好同时求解,所以可先求第二个乘子的解(),得到的解()之后,再用的解()表示的解()。
为了求解,得先确定的取值范围。假设它的上下边界分别为H和L,那么有:
接下来,综合和这两个约束条件,求取的取值范围。
当y1 != y2时,根据可得,所以有,,如下图所示:
当y1 = y2时,同样根据可得:,所以有,,如下图所示:
如此,根据y1和y2异号或同号,可得出的上下界分别为:
回顾下第二个约束条件,令上式两边乘以y1,可得
其中,。
因此可以用表示,,从而把子问题的目标函数转换为只含的问题:
对求导,可得
化简下:
然后将、、和代入上式可得:
令(表示预测值与真实值之差),,然后上式两边同时除以,得到一个关于单变量的解:
这个解没有考虑其约束条件,即是未经剪辑时的解。
然后考虑约束可得到经过剪辑后的的解析解为:
求出了后,便可以求出,得。
那么如何选择乘子和呢?
而对于第二个乘子可以寻找满足条件:的乘子。
而b在满足下述条件:下更新b:
且每次更新完两个乘子的优化后,都需要再重新计算b,及对应的Ei值。
最后更新所有,y和b,这样模型就出来了,从而即可求出咱们开头提出的分类函数:
此外,这里也有一篇类似的文章,大家可以参考下。
3.5.2 SMO算法的步骤
综上,总结下SMO的主要步骤,如下:
意思是,
第二步,固定除和之外的其他参数,确定W极值条件下的,由表示。
假定在某一次迭代中,需要更新,对应的拉格朗日乘子,,那么这个小规模的二次规划问题写为:那么在每次迭代中,如何更新乘子呢?引用这里的两张PPT说明下:
知道了如何更新乘子,那么选取哪些乘子进行更新呢?具体选择方法有以下两个步骤:
步骤2:在所有不违反KKT条件的乘子中,选择使|E1 −E2|最大的a2进行更新,使得能最大限度增大目标函数的值(类似于梯度下降. 此外,而,求出来的E代表函数ui对输入xi的预测值与真实输出类标记yi之差)。
最后,每次更新完两个乘子的优化后,都需要再重新计算b,及对应的Ei值。综上,SMO算法的基本思想是将Vapnik在1982年提出的Chunking方法推到极致,SMO算法每次迭代只选出两个分量ai和aj进行调整,其它分量则保持固定不变,在得到解ai和aj之后,再用ai和aj改进其它分量。与通常的分解算法比较,尽管它可能需要更多的迭代次数,但每次迭代的计算量比较小,所以该算法表现出较好的快速收敛性,且不需要存储核矩阵,也没有矩阵运算。
3.5.3 SMO算法的实现
行文至此,我相信,SVM理解到了一定程度后,是的确能在脑海里从头至尾推导出相关公式的,最初分类函数,最大化分类间隔,max1/||w||,min1/2||w||^2,凸二次规划,拉格朗日函数,转化为对偶问题,SMO算法,都为寻找一个最优解,一个最优分类平面。一步步梳理下来,为什么这样那样,太多东西可以追究,最后实现。如下图所示:
至于下文中将阐述的核函数则为是为了更好的处理非线性可分的情况,而松弛变量则是为了纠正或约束少量“不安分”或脱离集体不好归类的因子。
台湾的林智仁教授写了一个封装SVM算法的libsvm库,大家可以看看,此外这里还有一份libsvm的注释文档。
除了在这篇论文《fast training of support vector machines using sequential minimal optimization》中platt给出了SMO算法的逻辑代码之外,这里也有一份SMO的实现代码,大家可以看下。
3.6 SVM的应用
或许我们已经听到过,SVM在很多诸如文本分类,图像分类,生物序列分析和生物数据挖掘,手写字符识别等领域有很多的应用,但或许你并没强烈的意识到,SVM可以成功应用的领域远远超出现在已经在开发应用了的领域。
3.6.1 文本分类
一个文本分类系统不仅是一个自然语言处理系统,也是一个典型的模式识别系统,系统的输入是需要进行分类处理的文本,系统的输出则是与文本关联的类别。由于篇幅所限,其它更具体内容本文将不再详述。
OK,本节虽取标题为证明SVM,但聪明的读者们想必早已看出,其实本部分并无多少证明部分(特此致歉),怎么办呢?可以参阅《支持向量机导论》一书,此书精简而有趣。本节完。