Last updated on October 24, 2024 11:01 PM
引例
课上通过如下几个例子引出表示定理相关内容。
SVM
在带松弛变量的 SVM 表示中,
wmini∈[n]∑max(0,1−yi(wTxi+b))+λ∥w∥2
接下来将记号进行改写(本节内容均如此约定):换成 f(xi)=wTφ(xi),优化形式变为
wmini∈[n]∑max(0,1−yi(wTφ(xi)))+λ∥w∥2
(注意这里的 w 包含之前的 b,多进行了对 b 的 norm,但实际上不太影响)
为了推导对偶形式,将其再变回带约束的优化问题,此时重新引回 ξi=max(0,1−yiwTφ(xi)),问题变为
w,ξmini∈[n]∑ξi+2λ∥w∥2s. t. {ξi≥0ξi≥1−yiwTφ(xi),∀i∈[n]
此处 λ 除以二是为了方便求导。
引入两套拉格朗日乘子,写出拉格朗日函数:
L(w,ξ,α,β)=i∈[n]∑ξi+2λwTw+i∈[n]∑αi(1−yiwTφ(xi)−ξi)−i∈[n]∑βiξi
对偶问题则为:
α≥0,β≥0maxw,ξminL(w,ξ,α,β)
根据 K.K.T. 条件,里层要达到最小,则偏导为 0:
∂w∂L∂ξi∂L=λw−i∈[n]∑αiyiφ(xi)=0=1−αi−βi=0
这推出:
wαi+βi=λ1i∈[n]∑αiyiφ(xi)=1,∀i∈[n]
这说明 w 可以写成训练集数据经过某种训练后的线性组合。
带入外层并化简,发现 ξi 的项都没了,注意把 w 换掉:
α,βmax⎩⎪⎨⎪⎧−2λ1⎝⎛i∈[n]∑αiyiφ(xi)T⎠⎞⎝⎛j∈[n]∑αjyjφ(xj)⎠⎞+i∈[n]∑αi⎭⎪⎬⎪⎫s. t. αi≥0,βi≥0,αi+βi=1
βi 多余,
αmaxi∈[n]∑αi−2λ1i∈[n]∑j∈[n]∑αiαjyiyjφ(xi)Tφ(xj)s. t. 0≤αi≤1,∀i
与原先的区别:λ 与 αi≤1 的上界,且没有 ∑αiyi=0。从直觉上来说,αi 表示支持的“强度”,hard constraint 的时候,outlier 的 αi 就会很大,引入松弛变量后,就不会有某一个单独的向量(outlier)贡献特别大的 αi。
岭回归
目标函数:
wmin21i∈[n]∑(yi−wTφ(xi))2+2λ∥w∥2
对 w 求偏导:
∂w∂L=−i∈[n]∑(yi−wTφ(xi))φ(xi)+λw=0
所以
w=λ1i∈[n]∑(yi−wTφ(xi))φ(xi)
定义 αi=λ1(yi−wTφ(xi)),则 w=i∈[n]∑αiφ(xi),注意这里只是说明 w 满足此形式,即可以被表示为 φ(xi) 的线性组合。
逻辑回归 + L2-Norm
假设 y∈{−1,1}
wmini∈[n]∑log(1+exp(−yiwTφ(xi)))+2λ∥w∥2
一样,求偏导:
∂w∂L=−i∈[n]∑1+exp(−zi)exp(−zi)yiφ(xi)+λw=0
所以
w=λ1i∈[n]∑σ(−zi)yiφ(xi)
令 αi=λ1σ(−zi)yi,则 w=i∈[n]∑αiφ(xi)。
如果要预测一个点的话,f(x)=wTφ(x)=i∈[n]∑αiφ(xi)Tφ(x)=i∈[n]∑αik(xi,x),这不是偶然的。
一般性推测
根据上面的例子,我们可以猜想,对于任意线性模型 + ERM + 正则化都可以有
f(x)=i∈[n]∑αik(xi,x)
此为表示定理。为了进一步说明,引入一些数学工具。
再生核希尔伯特空间(RKHS)
在课上听这个的时候有点懵,这里进行重新整理。
由于本节内容不做考察要求,且作者的数学并不太好,所以下面的定义与推导都不严谨,如有错误敬请指出
Hilbert 空间
一个希尔伯特空间(Hilbert Space)是对欧氏空间的一个扩展,其还是向量空间,有内积 ⟨f,g⟩H。由内积可以定义范数 ∥f∥:=⟨f,f⟩H,且完备(no holes,每个柯西列都收敛到某个点)。
内积要满足的性质:
- 对称性:⟨f,g⟩H=⟨g,f⟩H;
- 线性性:⟨af1+bf2,g⟩H=a⟨f1,g⟩H+b⟨f2,g⟩H;
- 正定性:⟨f,f⟩H≥0 且 ⟨f,f⟩H=0⟺f=0
事实上,定义在集合 X 上的全体函数 f:X→R 便可以构成一个希尔伯特空间(函数可以看作为无穷维度的向量,回忆高代书上的内容(x))
RKHS
定义(from Wikipedia):考虑一个希尔伯特空间,其元素为定义在 X 上的函数,若 ∀x∈X,存在函数 Kx∈H 使得 ∀f∈H 有 f(x)=⟨f,Kx⟩H,则 H 为再生核希尔伯特空间(RKHS)
有点抽象,我们换个角度重新理解,因为理论上这个所谓的 RKHS 就是由核“再生”出来的,所以我们自然从核函数考虑起会自然一些。
考虑一个合法的核函数 K(x,y),由于我们可以将函数看成无穷维向量,此二元函数自然就可以看成无穷维矩阵。其正定性可以被解读为 ∫∫f(x)K(x,y)f(y)dxdy≥0,对称性可解读为 K(x,y)=K(y,x)。
那么同样地,类似对普通矩阵特征值分解,我们也可以对核函数进行类似操作,即存在特征值 λ 与特征函数 ψ(x) 使得
∫K(x,y)ψ(x)dx=λψ(x)
特征向量自然是正交的(证明略过),即
⟨ψ1,ψ2⟩H=∫ψ1(x)ψ2(x)dx=0
这说明一个核函数对应着无穷个特征值与无穷个特征方程,与矩阵类似地对其特征值分解:
K(x,y)=i=1∑∞λiψi(x)ψi(y)
这就是 SVM 一文中所说的 Mercer 定理,这说明所有的 ψi 构成了函数空间的一组正交基。
{λiψi}i=1∞ 自然为空间中的一组正交基。用这组基可以张成希尔伯特空间 H,该空间的任意一个函数都可以表示为这组基的线性组合:
f=i=1∑∞fiλiψi
自然,f 就可以被这组基表示为 (f1,f2,⋯)。若另一个函数 g 可以被表示为 (g1,g2,⋯),则 f 与 g 的内积就可以表示为
⟨f,g⟩H=i=1∑∞figi
对于核函数 K,固定一个变量后的 K(x,⋅) 为一元函数,用“特征值分解”的等式拆开:
K(x,⋅)=i=1∑∞λiψi(x)λiψi=i=1∑∞λiψi(x)ψi
用上述基可以表示为 (λ1ψ1(x),λ2ψ2(x),⋯)。
同理,K(⋅,y) 可以表示为 (λ1ψ1(y),λ2ψ2(y),⋯),所以,这二者的内积
⟨K(x,⋅),K(⋅,y)⟩H=i=1∑∞λiψi(x)ψi(y)=K(x,y)
这就是核的可再生性,H 即为 RKHS。
那么对于之前的 f=i=1∑∞fiλiψi,显然就有 f(x)=i=1∑∞fiλiψi(x),
⟨K(x,⋅),f⟩H=i=1∑∞λiψi(x)⋅fi=f(x)
感性理解一下,就是这个 kernel 决定了这个希尔伯特空间。这使得函数的求值可以被等价为在原空间 X 中的某种 similarity。k 从 X “reproduces” 每个 f∈H。
很多函数空间都是 RKHS:
-
考虑 f(x)=wTx(先忽略偏置项),所有的 f 就组成了一个希尔伯特空间 H。考虑 k(z,x)=zTx,f(x)=⟨k(w,⋅),k(⋅,x)⟩H=k(w,x)=wTx。这个核函数将 w 映射到 f。考虑 ∥f∥2,其为 ⟨k(w,⋅),k(⋅,w)⟩H=k(w,w)=∥w∥2。
-
考虑 f(x)=wTφ(x)(更为普遍的形式),定义核函数 k(z,x)=φ(z)Tφ(x)。因为 f(x)=⟨k(φ−1(w),⋅),k(⋅,x)⟩H=k(φ−1(w),x)=wTφ(x)。考虑 ∥f∥2,⟨k(φ−1(w),⋅),k(⋅,φ−1(w))⟩H=wTw。
所以其实对 w 正则化和对 f 正则化是等价的
表示定理
定理内容
考虑 RKHS H 与其再生核 k:X×X→R,给定训练数据 {(x1,y1),⋯,(xn,yn)}⊆X×R,损失函数 L:R×R→R,正则化项 R:[0,+∞)→R,且其严格递增。
对于任意使得 i∈[n]∑L(f(xi),yi)+R(∥f∥) 最小的 f∈H,都可以被表示成
f=i∈[n]∑αik(xi,⋅)
证明
先将 f 分解为两部分:一部分在 $\langle k(x_1,\cdot ),k(x_2,\cdot ), \cdots k(x_n,\cdot ) \rangle $ 上(此处记号为 span),另一部分 v 与前一部分正交。于是 ∀f 可以写成 f=i∈[n]∑αik(xi,⋅)+v,满足 $\langle v,k(x_i,\cdot ) \rangle_{\mathcal{H}} = 0 $。接下来证明 v 一定为 0 即可。
考虑一个训练样本点 xj,有
f(xj)=⟨i∈[n]∑αik(xi,⋅)+v,k(⋅,xj)⟩H=i∈[n]∑αi⟨k(xi,⋅),k(⋅,xj)⟩H+i∈[n]∑⟨v,k(⋅,xj)⟩H=i∈[n]∑αik(xi,xj)
这进而说明损失函数是跟 v 无关的。
考虑 R(∥f∥)=R(∥f0+v∥)。∥f0+v∥2=⟨f0,f0⟩H+⟨v,v⟩H+2⟨f0,v⟩H≥∥f0∥2。而且等于号成立当且仅当 v=0。又 R 严格递增,所以为了使得整体的值最小(损失函数项与 v 无关),我们自然是需要 v=0 的。□
直觉上理解起来应该挺自然的,模型的输出结果必然是与训练数据有关的,而正则化项把无关项去掉了。
优势
- 将一个对 f 的潜在的无穷维优化变成只对 n 个变量的搜索。
考虑 RBF 核,φ(x) 是无穷维的,wT 也应无穷维,这样就无法优化了,但利用表示定理我们就可以对 n 个 αi 进行优化了。
- 表明相当一部分机器学习算法的解能被有限训练数据上的核函数的线性组合表出(搭建参数化与非参数化的桥梁)。
线性回归的对偶形式
下节课的东西。