【机器学习】支持向量机【下】软间隔与核函数

news2024/9/23 11:25:27

有任何的书写错误、排版错误、概念错误等,希望大家包含指正。

在阅读本篇之前建议先学习:
【机器学习】拉格朗日对偶性
【机器学习】核函数

由于字数限制,分成两篇博客。
【机器学习】支持向量机【上】硬间隔
【机器学习】支持向量机【下】软间隔与核函数

线性支持向量机

软间隔最大化

在实际中,线性可分属于比较理想的情况,大多数数据都是线性不可分的,这时可以修改硬间隔最大化,使其成为软间隔最大化,来处理线性不可分问题。这里所谓的“线性不可分”是指,训练数据中有一些特异点(outlier),将这些特异点除去后,剩下大部分的样本点组成的集合是线性可分的。

线性不可分意味着某些样本点 ( x i , y i ) (x_i,y_i) (xi,yi) 不能满足约束条件 y i ( w T x i + b ) ≥ 1 y_i(w^Tx_i+b)\ge 1 yi(wTxi+b)1。为了解决这个问题,可以对每个样本点 ( x i , y i ) (x_i,y_i) (xi,yi) 引进一个松弛变量 ξ ≥ 0 \xi\ge 0 ξ0,使约束条件变为
y i ( w T x i + b ) ≥ 1 − ξ i y_i(w^Tx_i+b)\ge 1-\xi_i yi(wTxi+b)1ξi
对每一个松弛变量 ξ i \xi_i ξi,支付一个代价 ξ i \xi_i ξi。目标函数由原来的 1 2 ∣ ∣ w ∣ ∣ 2 \frac{1}{2} ||w||^2 21w2 变为
1 2 ∣ ∣ w ∣ ∣ 2 + C ∑ i = 1 n ξ i (9) \frac{1}{2} ||w||^2 + C\sum_{i=1}^n \xi_i \tag{9} 21w2+Ci=1nξi(9)
其中, C > 0 C>0 C>0 称为惩罚(超)参数,一般根据应用问题人为决定, C C C 值越大对误分类的惩罚越大。最小化目标函数式 ( 9 ) (9) (9) 包含两层含义:使 1 2 ∣ ∣ w ∣ ∣ 2 \frac{1}{2}||w||^2 21w2 尽量小即间隔尽量大,同时使误分类点的个数尽量少, C C C 是调和二者的系数。

基于上面允许误分类的思想,对于线性不可分的支持向量机学习问题,我们可以采用与线性可分时类似的学习过程。

定义原始问题:
min ⁡ w , b , ξ 1 2 ∣ ∣ w ∣ ∣ 2 + C ∑ i = 1 n ξ i s . t . y i ( w T x i + b ) ≥ 1 − ξ i , i = 1 , 2 , … , n ξ i ≥ 0 , i = 1 , 2 , … , n (10) \begin{matrix} &\min_{w,b,\xi} \frac{1}{2} ||w||^2 + C\sum_{i=1}^n\xi_i& \\\\ s.t.&y_i(w^Tx_i+b)\ge 1-\xi_i,& i=1,2,\dots,n \\ &\xi_i\ge0,& i=1,2,\dots,n \\ \end{matrix} \tag{10} s.t.minw,b,ξ21w2+Ci=1nξiyi(wTxi+b)1ξi,ξi0,i=1,2,,ni=1,2,,n(10)
原始问题式 ( 10 ) (10) (10) 是一个凸二次规划问题,因而关于 ( w , b , ξ ) (w,b,\xi) (w,b,ξ) 的最优解是存在的。而且 w w w 的最优解是唯一的,但 b b b 的最优解可能不唯一,而是存在于一个区间。假设最优解为 w ∗ w^* w b ∗ b^* b,于是可以得到划分超平面 w ∗ T x + b ∗ = 0 {w^*}^Tx+b^*=0 wTx+b=0 及分类决策函数 f ( x ) = s i g n ( w ∗ T x + b ∗ ) f(x)={\rm sign}({w^*}^Tx+b^*) f(x)=sign(wTx+b)。称这样的模型为训练样本线性不可分时的线性支持向量机,简称为线性支持向量机。

对偶问题

构建广义拉格朗日函数
L ( w , b , ξ , α , β ) = 1 2 ∣ ∣ w ∣ ∣ 2 + C ∑ i = 1 n ξ i − ∑ i = 1 n α i ( y i ( w T x i + b ) − 1 + ξ i ) − ∑ i = 1 n β i ξ i L(w, b, \xi,\alpha,\beta) = \frac{1}{2} ||w||^2 + C\sum_{i=1}^n \xi_i -\sum_{i=1}^n \alpha_i \big(y_i(w^Tx_i+b)-1+\xi_i\big) - \sum_{i=1}^n\beta_i\xi_i L(w,b,ξ,α,β)=21w2+Ci=1nξii=1nαi(yi(wTxi+b)1+ξi)i=1nβiξi
其中,拉格朗日乘子 α i ≥ 0 \alpha_i\ge 0 αi0 β i ≥ 0 \beta_i\ge 0 βi0

首先求 L ( w , b , ξ , α , β ) L(w,b,\xi,\alpha,\beta) L(w,b,ξ,α,β) w w w b b b ξ \xi ξ 的极小,由
∇ w L ( w , b , ξ , α , β ) = w − ∑ i = 1 n α i y i x i = 0 ∇ b L ( w , b , ξ , α , β ) = − ∑ i = 1 n α i y i = 0 ∇ ξ i L ( w , b , ξ , α , β ) = C − α i − β i = 0 \nabla_w L(w,b,\xi, \alpha,\beta) = w - \sum_{i=1}^n \alpha_iy_ix_i=0\\ \nabla_b L(w,b,\xi, \alpha,\beta) = - \sum_{i=1}^n \alpha_iy_i=0\\ \nabla_{\xi_i} L(w,b,\xi, \alpha,\beta) = C - \alpha_i - \beta_i=0\\ wL(w,b,ξ,α,β)=wi=1nαiyixi=0bL(w,b,ξ,α,β)=i=1nαiyi=0ξiL(w,b,ξ,α,β)=Cαiβi=0
可得
w = ∑ i = 1 n α i y i x i (11) w = \sum_{i=1}^n \alpha_i y_i x_i\tag{11} w=i=1nαiyixi(11)

∑ i = 1 n α i y i = 0 (12) \sum_{i=1}^n \alpha_i y_i=0\tag{12} i=1nαiyi=0(12)

C − α i − β i = 0 (13) C - \alpha_i - \beta_i = 0\tag{13} Cαiβi=0(13)

将式 ( 11 ) ∼ ( 13 ) (11)\sim(13) (11)(13) 代入拉格朗日函数中得
min ⁡ w , b , ξ L ( w , b , ξ , α , β ) = − 1 2 ∑ i = 1 n ∑ i = 1 n α i α j y u y j ( x i T x j ) + ∑ i = 1 n α i \min_{w,b, \xi} L(w, b, \xi,\alpha, \beta) = -\frac{1}{2} \sum_{i=1}^n\sum_{i=1}^n\alpha_i\alpha_jy_uy_j(x_i^Tx_j)+\sum_{i=1}^n\alpha_i w,b,ξminL(w,b,ξ,α,β)=21i=1ni=1nαiαjyuyj(xiTxj)+i=1nαi
再对 min ⁡ w , b , ξ L ( w , b , ξ , α , β ) \min\limits_{w,b, \xi} L(w, b, \xi,\alpha, \beta) w,b,ξminL(w,b,ξ,α,β) α \alpha α β \beta β 的极大。显然这已经与 β \beta β 无关,得到对偶问题:
max ⁡ α − 1 2 ∑ i = 1 n ∑ i = 1 n α i α j y u y j ( x i T x j ) + ∑ i = 1 n α i s . t . ∑ i = 1 n α i y i = 0 C − α i − β i = 0 α i ≥ 0 β i ≥ 0 ,      i = 1 , 2 , … , n (14) \begin{array}{ll} \max\limits_\alpha & -\frac{1}{2} \sum\limits_{i=1}^n\sum\limits_{i=1}^n\alpha_i\alpha_jy_uy_j(x_i^Tx_j)+\sum\limits_{i=1}^n\alpha_i \\ \\ s.t.& \sum \limits_{i=1}^n\alpha_iy_i = 0\\ &C-\alpha_i - \beta_i = 0 \\ &\alpha_i\ge 0\\ &\beta_i\ge 0,\space\space\space\space i=1,2,\dots, n \end{array} \tag{14} αmaxs.t.21i=1ni=1nαiαjyuyj(xiTxj)+i=1nαii=1nαiyi=0Cαiβi=0αi0βi0,    i=1,2,,n(14)
对对偶问题式 ( 14 ) (14) (14) 进行变形:利用等式约束 C − α i − β i = 0 C-\alpha_i-\beta_i=0 Cαiβi=0 消去 β i \beta_i βi,从而只留下变量 α i \alpha_i αi,并将后三条约束共同表达为
0 ≤ α i ≤ C 0\le \alpha_i\le C 0αiC
最终将问题从求极大转化为求极小,得
min ⁡ α 1 2 ∑ i = 1 n ∑ j = 1 n α i α j y u y j ( x i T x j ) − ∑ i = 1 n α i s . t . ∑ i = 1 n α i y i = 0 0 ≤ α i ≤ C ,      i = 1 , 2 , … , n (15) \begin{matrix} &\min\limits_\alpha \frac{1}{2} \sum\limits_{i=1}^n\sum\limits_{j=1}^n\alpha_i\alpha_jy_uy_j(x_i^Tx_j)-\sum\limits_{i=1}^n\alpha_i \\ \\ s.t.& \sum_{i=1}^n \alpha_iy_i = 0 \\ & 0\le \alpha_i\le C,\space\space\space\space i=1,2,\dots,n \\ \end{matrix} \tag{15} s.t.αmin21i=1nj=1nαiαjyuyj(xiTxj)i=1nαii=1nαiyi=00αiC,    i=1,2,,n(15)
由于问题具有强对偶性,所以 KKT 条件成立,即得
y i ( w ∗ T x i + b ∗ ) − 1 + ξ i ∗ ≥ 0 ξ i ∗ ≥ 0 y_i({w^*}^Tx_i+b^*)-1+\xi_i^* \ge 0\\\\ \xi_i^*\ge0\\ yi(wTxi+b)1+ξi0ξi0

∇ w L ( w ∗ , b ∗ , ξ ∗ , α ∗ , β ∗ ) = w ∗ − ∑ i = 1 n α i ∗ y i x i = 0 (16) \nabla_w L(w^*,b^*,\xi^*, \alpha^*,\beta^*) = w^* - \sum_{i=1}^n \alpha_i^*y_ix_i=0\tag{16} \\ wL(w,b,ξ,α,β)=wi=1nαiyixi=0(16)

∇ b L ( w ∗ , b ∗ , ξ ∗ , α ∗ , β ∗ ) = − ∑ i = 1 n α i ∗ y i = 0 \nabla_b L(w^*,b^*,\xi^*, \alpha^*,\beta^*) = - \sum_{i=1}^n \alpha_i^*y_i=0 bL(w,b,ξ,α,β)=i=1nαiyi=0

∇ ξ L ( w ∗ , b ∗ , ξ ∗ , α ∗ , β ∗ ) = C − α ∗ − β ∗ = 0 (17) \nabla_{\xi} L(w^*,b^*,\xi^*, \alpha^*,\beta^*) = C - \alpha^* - \beta^*=0 \tag{17} ξL(w,b,ξ,α,β)=Cαβ=0(17)

α i ∗ ≥ 0 β i ∗ ≥ 0 \alpha_i^*\ge 0\\\\ \beta_i^*\ge 0 αi0βi0

α i ∗ ( y i ( w ∗ T x i + b ∗ ) − 1 + ξ i ∗ ) = 0 (18) \alpha_i^*\big( y_i({w^*}^Tx_i+b^*)-1+\xi_i^* \big) = 0\tag{18} αi(yi(wTxi+b)1+ξi)=0(18)

β i ∗ ξ i ∗ = 0 (19) \beta^*_i\xi_i^*=0\tag{19} βiξi=0(19)

由式 ( 16 ) (16) (16) 可知
w ∗ = ∑ i = 1 n α i ∗ y i x i (20) w^* = \sum_{i=1}^n \alpha_i^*y_ix_i \tag{20} w=i=1nαiyixi(20)
若存在 α j ∗ \alpha_j^* αj 0 < α j ∗ < C 0<\alpha_j^*<C 0<αj<C,根据式 ( 18 ) (18) (18) 可知 y i ( w ∗ T x i + b ∗ ) − 1 + ξ i ∗ = 0 y_i({w^*}^Tx_i+b^*)-1+\xi_i^* =0 yi(wTxi+b)1+ξi=0;根据式 ( 17 ) (17) (17) 可知 β i ∗ = C − α i ∗ ≠ 0 \beta^*_i = C-\alpha^*_i\ne 0 βi=Cαi=0,所以由式 ( 19 ) (19) (19) 可得 ξ i ∗ = 0 \xi_i^*=0 ξi=0,则 y i ( w ∗ T x i + b ∗ ) − 1 = 0 y_i({w^*}^Tx_i+b^*)-1=0 yi(wTxi+b)1=0,故可得
b ∗ = y j − ∑ i = 1 n y i α i ∗ ( x i T x j ) (21) b^* = y_j - \sum_{i=1}^n y_i\alpha_i^*(x_i^Tx_j)\tag{21} b=yji=1nyiαi(xiTxj)(21)
由此,划分超平面可以写为
∑ i = 1 n α i ∗ y i ( x i T x ) + b ∗ = 0 \sum_{i=1}^n \alpha_i^*y_i(x_i^Tx)+b^*=0 i=1nαiyi(xiTx)+b=0
分类决策函数可以写为
f ( x ) = s i g n ( ∑ i = 1 n α i ∗ y i ( x i T x ) + b ∗ ) f(x) ={\rm sign}\Big( \sum_{i=1}^n \alpha_i^*y_i(x_i^Tx)+b^* \Big) f(x)=sign(i=1nαiyi(xiTx)+b)
上式称为线性支持向量机的对偶形式。线性支持向量机学习算法如下。

输入:   训 练 集   D = { ( x 1 , y 1 ) , ⋅ ⋅ ⋅ , ( x n , y n ) } ,   其 中   x i ∈ R d ,   y i ∈ { + 1 , − 1 } ,   i = 1 , … , n 过程: \begin{array}{ll} \textbf{输入:}&\space训练集\space D = \{(x_1,y_1),···,(x_n,y_n)\},\space 其中\space x_i\in \mathbb R^{d},\space y_i\in \{+1,-1\},\space i=1,\dots,n \\ \textbf{过程:} \end{array} 输入:过程:  D={(x1,y1),,(xn,yn)},  xiRd, yi{+1,1}, i=1,,n

1 : 选 择 惩 罚 参 数   C > 0 , 构 造 并 求 解 凸 二 次 规 划 问 题 min ⁡ α 1 2 ∑ i = 1 n ∑ j = 1 n α i α j y i y j ( x i T x j ) − ∑ i = 1 n α i s . t .      ∑ i = 1 n α i y i = 0 0 ≤ α i ≤ C ,      i = 1 , 2 , … , n 求 得 最 优 解   α ∗ = { α i ∗ } ,   i = 1 , 2 , … , n 2 : 计 算 w ∗ = ∑ i = 1 n α i ∗ y i x i 3 : 选 择   0 < α j ∗ < C ,   计 算 b ∗ = y j − ∑ i = 1 n α i ∗ y i ( x i T x j ) 4 : 求 得 划 分 超 平 面 w ∗ T x + b ∗ = 0 分 类 决 策 函 数 f ( x ) = s i g n ( w ∗ T x + b ∗ ) \begin{array}{rl} 1:& 选择惩罚参数 \space C>0,构造并求解凸二次规划问题\\ \\ &\begin{array}{c} & \min \limits_\alpha \frac{1}{2}\sum \limits_{i=1}^n\sum \limits_{j=1}^n \alpha_i\alpha_jy_iy_j(x_i^Tx_j) -\sum \limits_{i=1}^n \alpha_i &\\ &s.t.\space\space\space\space \sum\limits_{i=1}^n \alpha_iy_i = 0 \\ &0\le\alpha_i\le C,\space\space\space\space i=1,2,\dots, n \\\\ \end{array}\\ & 求得最优解 \space \alpha^* = \{\alpha_i^*\},\space i=1,2,\dots,n \\ 2:& 计算\\ \\ &\begin{array}{c} &&w^* = \sum \limits_{i=1}^n \alpha_i^*y_ix_i \end{array}\\\\ 3:& 选择 \space 0<\alpha_j^*<C,\space 计算\\ \\ &\begin{array}{c} &&b^* = y_j - \sum\limits_{i=1}^n\alpha_i^*y_i(x_i^Tx_j) \end{array}\\\\ 4:& 求得划分超平面\\ \\ &\begin{array}{c} &&{w^*}^Tx+b^*=0 \end{array}\\\\ & 分类决策函数\\ \\ &\begin{array}{c} &&f(x)={\rm sign}({w^*}^Tx+b^*) \end{array} \end{array} 1:2:3:4: C>0,αmin21i=1nj=1nαiαjyiyj(xiTxj)i=1nαis.t.    i=1nαiyi=00αiC,    i=1,2,,n α={αi}, i=1,2,,nw=i=1nαiyixi 0<αj<C, b=yji=1nαiyi(xiTxj)wTx+b=0f(x)=sign(wTx+b)

输出:   划 分 超 平 面 和 分 类 决 策 函 数 \begin{array}{l} \textbf{输出:}\space 划分超平面和分类决策函数 &&&&&&&&&&&&&&&&&& \end{array} 输出: 

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

支持向量

在线性不可分的情况下,将对偶问题式 ( 15 ) (15) (15) 的解 α ∗ = { α i ∗ } α^*=\{\alpha_i^*\} α={αi} 中对应于 α i ∗ > 0 α_i^* > 0 αi>0 的样本点 ( x i , y j ) (x_i, y_j) (xi,yj) 的实例 x i x_i xi 称为支持向量(软间隔的支持向量)。如图 3 3 3 所示,这时的支持向量要比线性可分时的情况复杂一些。图中,划分超平面由实线表示,间隔边界由虚线表示,正样本由“o”表示,负样本由“×”表示。图中还标出了样本 x i x_i xi 到间隔边界的距离 ξ i ∣ ∣ w ∣ ∣ \frac{\xi_i}{||w||} wξi

在这里插入图片描述

图 3    软间隔的支持向量

软间隔的支持向量 x i x_i xi 或者在间隔边界上,或者在间隔边界与划分超平面之间,或者在划分超平面误分一侧。由式 ( 17 ) (17) (17) ( 19 ) (19) (19) 可得,若 α i ∗ < C α_i^* <C αi<C,则 ξ i = 0 \xi_i= 0 ξi=0,支持向量 x i x_i xi 恰好落在间隔边界上;若 α i ∗ = C \alpha_i^*= C αi=C 0 < ξ i < 1 0<\xi_i<1 0<ξi<1,则分类正确, x i x_i xi 在间隔边界与划分超平面之间;若 α i ∗ = C \alpha_i^*=C αi=C ξ i = 1 \xi_i= 1 ξi=1,则 x i x_i xi 在划分超平面上;若 α i ∗ = C \alpha_i^* = C αi=C ξ i > 1 \xi_i> 1 ξi>1,则 x i x_i xi 位于划分超平面误分一侧。

合页损失函数

还可以从另一个角度解释线性支持向量机学习。定义损失函数为划分错误的样本个数,即
J ( w , b ) = ∑ i = 1 n 1 { y i ( w T x i + b ) < 1 } J(w,b) = \sum_{i=1}^n1\{y_i(w^Tx_i+b)<1\} J(w,b)=i=1n1{yi(wTxi+b)<1}
其中, 1 { ⋅ } 1\{·\} 1{} 的作用类似于艾弗森括号,即括号内为真返回 1 1 1,否则返回 0 0 0。显然,函数 J J J 非连续,且存在跳跃,数学性质不好,不便于求导。故将损失函数 J J J 重新定义
J ( w , b ) = ∑ i = 1 n [ 1 − y i ( w T x i + b ) ] + J(w,b) = \sum_{i=1}^n[1-y_i(w^Tx_i+b)]_+ J(w,b)=i=1n[1yi(wTxi+b)]+
其中,函数 [ z ] + [z]_+ [z]+ 为合页函数(hinge function)
[ z ] + = { z , z > 0 0 , z ≤ 0 [z]_+ = \left\{ \begin{matrix} z,&z>0\\ 0,&z\le 0 \end{matrix} \right. [z]+={z,0,z>0z0
也可以等价表示为
[ z ] + = max ⁡ ( 0 , z ) [z]_+ = \max(0, z) [z]+=max(0,z)
函数 J J J 加上正则化项后定义为最终损失函数 L L L
L ( w , b ) = J ( w , b ) + λ ∣ ∣ w ∣ ∣ 2 = ∑ i = 1 n [ 1 − y i ( w T x i + b ) ] + + λ ∣ ∣ w ∣ ∣ 2 \begin{aligned} L(w,b) &= J(w,b) + \lambda ||w||^2 \\ &= \sum_{i=1}^n[1-y_i(w^Tx_i+b)]_+ + \lambda ||w||^2 \end{aligned} L(w,b)=J(w,b)+λw2=i=1n[1yi(wTxi+b)]++λw2
因此我们的优化目标为
min ⁡ w , b ∑ i = 1 n [ 1 − y i ( w T x i + b ) ] + + λ ∣ ∣ w ∣ ∣ 2 (22) \min_{w,b} \sum_{i=1}^n[1-y_i(w^Tx_i+b)]_+ + \lambda ||w||^2 \tag{22} w,bmini=1n[1yi(wTxi+b)]++λw2(22)
观察损失函数 L L L,当样本点 ( x i , y i ) (x_i ,y_i) (xi,yi) 被正确分类且确信度 y i ( w T x i + b ) y_i(w^T x_i+b) yi(wTxi+b) 大于 1 1 1 时,损失为 0 0 0,否则损失为 1 − y i ( w T x i + b ) 1-y_i(w^T x_i+b) 1yi(wTxi+b)。注意到图 3 3 3 中的样本点 x 4 x_4 x4 被正确分类,但损失不是 0 0 0。损失函数的第二项是系数为 λ \lambda λ w w w L 2 L_2 L2 范数,是正则化项。

可以很容易证明优化问题 ( 22 ) (22) (22) 与原始问题 ( 10 ) (10) (10) 等价。令
[ 1 − y i ( w T x i + b ) ] + = ξ i [1-y_i(w^Tx_i+b)]_+ = \xi_i [1yi(wTxi+b)]+=ξi
ξ i ≥ 0 \xi_i\ge0 ξi0,式 ( 10 ) (10) (10) 中的第二个不等式约束成立。由上式,当 1 − y i ( w T x i + b ) > 0 1-y_i(w^Tx_i+b)>0 1yi(wTxi+b)>0 时,有 y i ( w T x i + b ) = 1 − ξ i y_i(w^Tx_i+b)=1-\xi_i yi(wTxi+b)=1ξi;当 1 − y i ( w T x i + b ) ≤ 0 1-y_i(w^Tx_i+b)\le 0 1yi(wTxi+b)0 时, ξ i = 0 \xi_i=0 ξi=0,有 y i ( w T x i + b ) ≥ 1 − ξ i y_i(w^Tx_i+b)\ge 1-\xi_i yi(wTxi+b)1ξi。故式 ( 10 ) (10) (10) 中的第一个不等式约束成立。于是 w w w b b b ξ i \xi_i ξi 满足式 ( 10 ) (10) (10) 的约束条件,所以最优化问题 ( 22 ) (22) (22) 可以写成
min ⁡ w , b ∑ i = 1 n ξ i + λ ∣ ∣ w ∣ ∣ 2 \min_{w,b} \sum_{i=1}^n \xi_i + \lambda||w||^2 w,bmini=1nξi+λw2
若取 λ = 1 2 C \lambda=\frac{1}{2C} λ=2C1,则
min ⁡ w , b 1 C ( 1 2 ∣ ∣ w ∣ ∣ 2 + C ∑ i = 1 n ξ i ) \min_{w,b}\frac{1}{C} \Big(\frac{1}{2} ||w||^2 + C\sum_{i=1}^n\xi_i \Big) w,bminC1(21w2+Ci=1nξi)
与原始问题 ( 10 ) (10) (10) 等价。反之,也可以将优化问题 ( 22 ) (22) (22) 表示成问题 ( 10 ) (10) (10)

合页损失函数 L ( y ( x T x + b ) ) L\big(y(x^Tx+b)\big) L(y(xTx+b)) 的图形如图 4 4 4 所示,横轴为 y ( w T x + b ) y(w^Tx+b) y(wTx+b),纵轴为损失。由于函数形状像一个合页,故名合页损失函数。图中还画出了出 01 01 01 损失函数,可以认为它是二类分类问题的真正的损失函数,而合页损失函数是 0 ~ 1 0\text{\textasciitilde}1 0~1 损失函数的上界。由于 0 ~ 1 0\text{\textasciitilde}1 0~1 损失函数是跳跃的,不利于求导,直接优化由其构成的目标函数比较困难,可以认为线性支持向量机是优化由 0 ~ 1 0\text{\textasciitilde}1 0~1 损失函数的上界(合页损失函数)构成的目标函数。这时的上界损失函数又称为代理损失函数(surrogate loss function)。

在这里插入图片描述

图 4    合页损失函数

4 4 4 中虚线显示的是感知机的损失函数 [ − y i ( w T x i + b ) ] + [-y_i(w^Tx_i+b)]_+ [yi(wTxi+b)]+。这时,当样本点 ( x i , y i ) (x_i, y_i) (xi,yi) 被正确分类时,损失为 0 0 0,否则损失是 − y i ( w T x i + b ) -y_i(w^Tx_i+b) yi(wTxi+b)。相比之下,合页损失函数不仅要分类正确,而且确信度足够高时损失才为 0 0 0。也就是说,合页损失函数对学习有更高的要求。

将式 ( 15 ) (15) (15) 与硬间隔对偶问题对比可以看出,二者唯一的差别就在于对偶变量的约束不同:前者是 0 ≤ α i 0\le \alpha_i 0αi,后者是 0 ≤ α i ≤ C 0\le \alpha_i\le C 0αiC。因此,在引入核函数后,硬间隔和软间隔能够得到同样形式的决策函数。

非线性支持向量机

我们注意到在线性支持向量机的对偶问题中,无论是目标函数还是决策函数(划分超平面)都只涉及输入样本与样本之间的内积。在对偶问题的目标函数 ( 15 ) (15) (15) 中的内积 x i T x j x_i^Tx_j xiTxj 可以用核函数 K ( x i , x j ) = ϕ ( x i ) T ϕ ( x j ) K(x_i,x_j)=\phi(x_i)^T\phi(x_j) K(xi,xj)=ϕ(xi)Tϕ(xj) 来代替。此时对偶问题的目标函数成为
W ( α ) = 1 2 ∑ i = 1 n ∑ j = 1 n α i α j y u y j ( ϕ ( x i ) T ϕ ( x j ) ) − ∑ i = 1 n α i = 1 2 ∑ i = 1 n ∑ j = 1 n α i α j y u y j K ( x i , x j ) − ∑ i = 1 n α i \begin{aligned} W(\alpha) &= \frac{1}{2} \sum\limits_{i=1}^n\sum\limits_{j=1}^n\alpha_i\alpha_jy_uy_j \big(\phi(x_i)^T\phi(x_j)\big)-\sum_{i=1}^n\alpha_i \\ &=\frac{1}{2} \sum\limits_{i=1}^n\sum\limits_{j=1}^n\alpha_i\alpha_jy_uy_jK(x_i,x_j)-\sum_{i=1}^n\alpha_i \end{aligned} W(α)=21i=1nj=1nαiαjyuyj(ϕ(xi)Tϕ(xj))i=1nαi=21i=1nj=1nαiαjyuyjK(xi,xj)i=1nαi
同样,分类决策函数中的内积也可以用核函数替代,而分类决策函数成为
f ( x ) = s i g n ( ∑ i = 1 N α i ∗ y i ( ϕ ( x i ) T ϕ ( x ) ) + b ∗ ) = s i g n ( ∑ i = 1 N α i ∗ y i K ( x i , x ) + b ∗ ) \begin{aligned} f(x) &= {\rm sign}\Big( \sum_{i=1}^N \alpha_i^* y_i \big(\phi(x_i)^T\phi(x)\big) +b^* \Big) \\ &= {\rm sign} \Big( \sum_{i=1}^N \alpha_i^* y_i K(x_i,x) +b^* \Big) \\ \end{aligned} f(x)=sign(i=1Nαiyi(ϕ(xi)Tϕ(x))+b)=sign(i=1NαiyiK(xi,x)+b)
这等价于经过映射函数 ϕ \phi ϕ 将原来的输入空间变换到一个新的特征空间,将输入空间中的内积 x i T x j x_i^Tx_j xiTxj 变换为特征空间中的内积 ϕ ( x i ) T ϕ ( x j ) \phi(x_i)^T\phi(x_j) ϕ(xi)Tϕ(xj),在新的特征空间里从训练样本中学习线性支持线性向量机。当映射函数是非线性函数时,学习到的含有核函数的支持向量机是非线性分类模型。

也就是说,在核函数 K ( x , z ) K(x,z) K(x,z) 给定的条件下,可以利用解线性分类问题的方法求解非线性分类问题的支持向量机。学习是隐式地在特征空间进行的,不需要显式地定义特征空间和映射函数。这样的技巧称为核技巧,它是巧妙地利用线性分类学习方法与核函数解决非线性问题的技术。在实际应用中,往往依赖领域知识直接选择核函数,核函数选择的有效性需要通过实验验证。

非线性支持向量机学习算法如下。

输入:   线 性 可 分 训 练 集   D = { ( x 1 , y 1 ) , ⋅ ⋅ ⋅ , ( x n , y n ) } ,   其 中   x i ∈ R d ,   y i ∈ { + 1 , − 1 } ,   i = 1 , … , n 过程: \begin{array}{ll} \textbf{输入:}&\space线性可分训练集\space D = \{(x_1,y_1),···,(x_n,y_n)\},\space 其中\space x_i\in \mathbb R^{d},\space y_i\in \{+1,-1\},\space i=1,\dots,n \\ \textbf{过程:} \end{array} 输入:过程: 线 D={(x1,y1),,(xn,yn)},  xiRd, yi{+1,1}, i=1,,n

1 : 选 取 适 当 的 核 函 数   K ( x , z )   和 适 当 的 参 数   C , 构 造 并 求 解 最 优 化 问 题 min ⁡ α 1 2 ∑ i = 1 n ∑ j = 1 n α i α j y i y j ( x i T x j ) − ∑ i = 1 n α i s . t .      ∑ i = 1 n α i y i = 0 0 ≤ α i ≤ C ,      i = 1 , 2 , … , n 求 得 最 优 解   α ∗ = { α i ∗ } ,   i = 1 , 2 , … , n 2 : 计 算 w ∗ = ∑ i = 1 n α i ∗ y i x i 3 : 选 择   α j ∗ > 0 ,   计 算 b ∗ = y j − ∑ i = 1 n α i ∗ y i K ( x i , x j ) 4 : 分 类 决 策 函 数 f ( x ) = s i g n ( ∑ i = 1 n α i ∗ y i K ( x i , x ) + b ∗ ) \begin{array}{rl} 1:& 选取适当的核函数 \space K(x,z) \space 和适当的参数\space C,构造并求解最优化问题\\ \\ &\begin{array}{c} & \min \limits_\alpha \frac{1}{2}\sum \limits_{i=1}^n\sum \limits_{j=1}^n \alpha_i\alpha_jy_iy_j(x_i^Tx_j) -\sum \limits_{i=1}^n \alpha_i &\\ &s.t.\space\space\space\space \sum\limits_{i=1}^n \alpha_iy_i = 0 \\ &0\le\alpha_i\le C,\space\space\space\space i=1,2,\dots, n \\\\ \end{array}\\ & 求得最优解 \space \alpha^* = \{\alpha_i^*\},\space i=1,2,\dots,n \\ 2:& 计算\\ \\ &\begin{array}{c} &&w^* = \sum \limits_{i=1}^n \alpha_i^*y_ix_i \end{array}\\\\ 3:& 选择 \space \alpha_j^*>0,\space 计算\\ \\ &\begin{array}{c} &&b^* = y_j - \sum\limits_{i=1}^n\alpha_i^*y_iK(x_i,x_j) \end{array}\\\\ 4: & 分类决策函数\\ \\ &\begin{array}{c} &&f(x)={\rm sign}\left(\sum\limits_{i=1}^n\alpha_i^*y_i K(x_i,x)+b^*\right) \end{array} \end{array} 1:2:3:4: K(x,z)  C,αmin21i=1nj=1nαiαjyiyj(xiTxj)i=1nαis.t.    i=1nαiyi=00αiC,    i=1,2,,n α={αi}, i=1,2,,nw=i=1nαiyixi αj>0, b=yji=1nαiyiK(xi,xj)f(x)=sign(i=1nαiyiK(xi,x)+b)

输出:   分 类 决 策 函 数 \begin{array}{l} \textbf{输出:}\space 分类决策函数 &&&&&&&&&&&&&&&&&& \end{array} 输出: 

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

REF

[1]《统计学习方法(第二版)》李航著

[2]《机器学习》周志华著

[2] 机器学习-白板推导系列(六)-支持向量机SVM(Support Vector Machine)- bilibili

[3] 机器学习-白板推导系列(六)-支持向量机SVM(Support Vector Machine)- 知乎

[4] 【机器学习】拉格朗日对偶性 - CSDN

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/61039.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

研发效能工程实践-利用Superset快速打造大数据BI平台

大数据BI平台自研之殇 随着互联网发展&#xff0c;现在随便哪个公司都手握大量数据。如何利用这些数据为公司商业带来价值&#xff0c;触使各个公司投入大量人力财力去做商业智能。 早期的BI可能就是公司Leader叫开发小哥写几句SQL导出数据&#xff0c;然后导入到Excel里绘制几…

echarts:nuxt项目使用echarts

一、项目环境 nuxt 2.X vue2.X vuex webpack 二、安装 yarn add echarts 三、使用 3.1、plugins目录下创建echarts.js import Vue from vue import * as echarts from echarts // 引入echarts Vue.prototype.$echarts echarts // 引入组件&#xff08;将echarts注册为全…

认证服务-----技术点及亮点

大技术 Nacos做注册中心 把新建的微服务注册到Nacos上去 两个步骤 在配置文件中配置应用名称、nacos的发现注册ip地址&#xff0c;端口号在启动类上用EnableDiscoveryClient注解开启注册功能 使用Redis存验证码信息 加入依赖配置地址和端口号即可 直接注入StringRedisTempla…

HTML静态网页作业——关于我的家乡介绍安庆景点

家乡旅游景点网页作业制作 网页代码运用了DIV盒子的使用方法&#xff0c;如盒子的嵌套、浮动、margin、border、background等属性的使用&#xff0c;外部大盒子设定居中&#xff0c;内部左中右布局&#xff0c;下方横向浮动排列&#xff0c;大学学习的前端知识点和布局方式都有…

Armbian搭建本地Gitea服务器

Armbian搭建本地Gitea服务器 1 安装Docker Docker 是一个用于开发、发布和运行应用程序的开放平台。 Docker 是一个开源的应用容器引擎&#xff0c;Docker 可以让开发者打包他们的应用以及依赖包到一个轻量级、可移植的容器中&#xff0c;然后发布到任何流行的 Linux 机器上&…

R语言中的prophet预测时间序列数据模型

本文 将针对R进行的几次建模练习的结果&#xff0c;以魁北克数据为依据&#xff0c;分为13年的训练和1年的测试。prophet与基本线性模型&#xff08;lm&#xff09;&#xff0c;一般加性模型&#xff08;gam&#xff09;和随机森林&#xff08;randomForest&#xff09;进行了比…

ES6:ES6 的新增语法

什么是 ES6 ? ES 的全称是 ECMAScript , 它是由 ECMA 国际标准化组织,制定的一项脚本语言的标准化规范。 年份 版本 2015年6月 ES2015 2016年6月 ES2016 2017年6月 ES2017 2018年6月 ES2018 … … ES6 实际上是一个泛指&#xff0c;泛指 ES2015 及后续的版本。 …

基于最大熵图像插值Maximum Entropy插值算法的图像超分辨重构研究-附Matlab代码

⭕⭕ 目 录 ⭕⭕✳️ 一、引言✳️ 二、图像复原基本原理✳️ 三、最大熵图像插值原理✳️ 四、实验验证✳️ 五、参考文献✳️ 六、Matlab程序获取与验证✳️ 一、引言 图像是一种表达信息的形式&#xff0c;其中&#xff0c;数字图像反馈的信息更加丰富。 在获取图像的过程中…

基于N32G45的按键驱动

基于N32G45的按键驱动 1.N32G45简介 N32G45系列集成了最新一代嵌入式ARM Cortex™-M4F处理器&#xff0c;在Cortex™-M3内核的基础上强化了运算能力、新增加了浮点运算处理单元&#xff08;FPU&#xff09;、DSP和并行计算指令&#xff0c;提供1.25DMIPS/MHz的优异性能。同时其…

JAVA复习【11】单列集合Collection:ArrayList、 LinkedList、HashSet、TreeSet学习与使用

1.首先思考一个问题&#xff1a;为什么要有集合&#xff1f; 我们也知道&#xff0c;数组可以保存多个对象&#xff0c;但是在某些情况下无法确定到底需要保存多少个对象&#xff0c;此时数组不再适用没因为数组的长度不可变&#xff0c;例如&#xff0c;要保存一个学校的学生信…

移动WEB开发之rem布局--less基础

维护 css 的弊端 CSS 是一门非程序式语言&#xff0c;没有变量、函数、SCOPE&#xff08;作用域&#xff09;等概念。 CSS 需要书写大量看似没有逻辑的代码&#xff0c;CSS 冗余度是比较高的。 不方便维护及扩展&#xff0c;不利于复用。 CSS 没有很好的计算能力 非前端开…

前馈神经网络与支持向量机实战 --- 手写数字识别

前馈神经网络与支持向量机实战 — 手写数字识别 文章目录前馈神经网络与支持向量机实战 --- 手写数字识别一、前馈神经网络介绍二、支持向量机介绍三、数据集说明四、环境准备五、实验要求六、Python代码tutorial_minst_fnn-keras.py&#xff1a;使用TensorFlow的Sequential实现…

Linux开发常用ps命令选项详解

【摘要】本文介绍了在Linux应用/内核开发调试中&#xff0c;经常需要用到的两个选项组合&#xff0c;当然&#xff0c;如果你需要查看更多更详尽的选项说明&#xff0c;可以参考man说明文档&#xff0c;即命令行下输入man ps进行查看。 aux选项组合 使用场景&#xff1a;更多…

算法日常训练12.4(最接近目标价格甜点成本)

只能说回溯实在是诡异&#xff0c;刚看到这题目思路一点不清晰&#xff0c;想着用回溯想到一点写一点&#xff0c;就这样诡异的出来了。 主要回溯思想&#xff0c;由于冰淇淋基料只能选一种&#xff0c;那就对数组遍历&#xff0c;每次对一种冰淇淋基料继续回溯&#xff0c;用r…

Proxmox VE 修改集群名称

作者&#xff1a;田逸&#xff08;formyz) Proxmox VE集群一旦创建&#xff0c;其集群的名称就固定下来。在Proxmox VE Web管理后台&#xff0c;没有相应的菜单或按钮对应与集群名称的修改&#xff08;仅仅发现修改虚拟机选项有修改名称的地方&#xff09;。在宿主系统Debian下…

网课题库接口

网课题库接口 本平台优点&#xff1a; 多题库查题、独立后台、响应速度快、全网平台可查、功能最全&#xff01; 1.想要给自己的公众号获得查题接口&#xff0c;只需要两步&#xff01; 2.题库&#xff1a; 题库&#xff1a;题库后台&#xff08;点击跳转&#xff09; 题库…

[附源码]计算机毕业设计文曦家教预约系统Springboot程序

项目运行 环境配置&#xff1a; Jdk1.8 Tomcat7.0 Mysql HBuilderX&#xff08;Webstorm也行&#xff09; Eclispe&#xff08;IntelliJ IDEA,Eclispe,MyEclispe,Sts都支持&#xff09;。 项目技术&#xff1a; SSM mybatis Maven Vue 等等组成&#xff0c;B/S模式 M…

java计算机毕业设计ssm社团管理系统9e73v(附源码、数据库)

java计算机毕业设计ssm社团管理系统9e73v&#xff08;附源码、数据库&#xff09; 项目运行 环境配置&#xff1a; Jdk1.8 Tomcat8.5 Mysql HBuilderX&#xff08;Webstorm也行&#xff09; Eclispe&#xff08;IntelliJ IDEA,Eclispe,MyEclispe,Sts都支持&#xff09;。…

计及新能源出力不确定性的电气设备综合能源系统协同优化(Matlab代码实现)

目录 第一部分 文献一《计及新能源出力不确定性的电气设备综合能源系统协同优化》 0 引言 &#xff11; 新能源出力不确定性处理 1.1 新能源出力预测误差分布 1.2 新能源出力的时间相关性 1.3 场景生成 &#xff12; 计及温控负荷调节能力的电气综合能源系统协同优化建模…

[附源码]计算机毕业设计文具商城购物系统Springboot程序

项目运行 环境配置&#xff1a; Jdk1.8 Tomcat7.0 Mysql HBuilderX&#xff08;Webstorm也行&#xff09; Eclispe&#xff08;IntelliJ IDEA,Eclispe,MyEclispe,Sts都支持&#xff09;。 项目技术&#xff1a; SSM mybatis Maven Vue 等等组成&#xff0c;B/S模式 M…