A02 · 第 1 章 · 第一编 核方法

特征映射、正定核与再生核空间

从显式特征映射和 Gram 矩阵建立正半定核条件,解释核技巧、线性核、多项式核与 RBF 核的尺度,准确导入再生性质,并讨论中心化、归一化、计算成本和核选择边界。

报告页面错误
预备知识Hilbert 空间、正交投影与对偶正则化正交性

本章目标

  1. 由显式特征映射计算核与 Gram 矩阵,并验证有限样本上的正半定性。
  2. 说明正半定核条件为何保证内积表示,不把任意对称相似度当作核。
  3. 比较线性、多项式与 RBF 核的尺度、表达偏好和参数边界。
  4. 准确陈述再生性质与有限核展开的适用条件,不把核技巧等同于无限维。
  5. 对 Gram 矩阵做中心化与归一化,评估训练、存储和预测成本。
页面阅读位置0% · 仅保存在此浏览器
章节未开始
本册完成进度0/6 章 · 0%
本页目录

显式特征映射把非线性变成线性

设输入空间为 X\mathcal X,特征映射

ϕ:XH\phi:\mathcal X\to\mathcal H

把输入送入带内积的特征空间 H\mathcal H。模型可以在该空间中保持线性形式

f(x)=w,ϕ(x)H+b,f(x)=\langle w,\phi(x)\rangle_{\mathcal H}+b,

却在原始坐标中表现为非线性。例如一维映射 ϕ(x)=(1,x,x2)\phi(x)=(1,x,x^2) 使线性权重组合成为二次函数。特征映射规定模型能轻易表达哪些变化,是归纳偏置的一部分;它不是把信息无条件变多,错误尺度或不相关高阶项仍会增加估计与计算困难。

有限维映射可以直接构建设计矩阵再训练线性模型。若映射维度很高,算法有时只需要样本特征之间的内积,而不需要单独访问每个坐标。定义

k(x,z)=ϕ(x),ϕ(z)H,k(x,z)=\langle\phi(x),\phi(z)\rangle_{\mathcal H},

就能以一个标量核计算替代显式内积。这个替代叫核技巧;它要求原算法能够完全改写为训练样本内积,并且 kk 确实对应某个内积空间。核技巧不自动降低样本数方向的成本,也不表示所有核都来自无限维映射。

例 1:显式展开二次多项式核

在二维输入上取

k(x,z)=(1+xTz)2.k(x,z)=(1+x^\mathsf Tz)^2.

展开后可选特征映射

ϕ(x)=(1,2x1,2x2,x12,2x1x2,x22).\phi(x)= (1,\sqrt2x_1,\sqrt2x_2,x_1^2,\sqrt2x_1x_2,x_2^2).

x=(1,2)x=(1,2)z=(3,1)z=(3,-1),原始内积为 xTz=1x^\mathsf Tz=1,所以核值为 4。显式特征内积也为

1+64+912+4=4.1+6-4+9-12+4=4.

根号二系数不是装饰,它使交叉项在内积展开中得到系数 2。直接计算核只做原始二维内积与一次平方,不必真的存六维向量;但六维映射有限且很小的时候,显式训练也可能更便于批量计算和解释。

Gram 矩阵与正半定条件

给定样本 x1,,xnx_1,\ldots,x_n,核 Gram 矩阵定义为

Kij=k(xi,xj).K_{ij}=k(x_i,x_j).

若核来自内积映射,则 KK 对称,并且任意 cRnc\in\mathbb R^n 都满足

cTKc=i=1nciϕ(xi)H20.c^\mathsf TKc =\left\lVert\sum_{i=1}^n c_i\phi(x_i)\right\rVert_{\mathcal H}^2 \ge0.

因此每个有限样本集的 Gram 矩阵都必须正半定。反过来,一个实值对称函数若对任意有限输入集合都产生正半定 Gram 矩阵,就能建立相应 Hilbert 特征空间表示。注意要求是“任意集合”,在一份训练样本上偶然没有负特征值不足以证明函数对所有输入有效。

正半定允许零特征值,表示某些有限样本特征线性相关。数值计算得到微小负特征值可能来自浮点误差,也可能来自公式、近似或数据依赖核不再正半定。应结合矩阵尺度、对称误差和构造方法判断,不能不加说明地截断所有负值。

例 2:用二次型识别无效相似函数

一维线性核 k(x,z)=xzk(x,z)=xz 对样本 x1=1,x2=2x_1=1,x_2=2 给出

K=[1224].K=\begin{bmatrix}1&2\\2&4\end{bmatrix}.

它等于向量 (1,2)T(1,2)(1,2)^\mathsf T(1,2),所以任意 cc 都有 cTKc=(c1+2c2)20c^\mathsf TKc=(c_1+2c_2)^2\ge0。矩阵行列式为零,说明它正半定但不正定。

若有人把 s(x,z)=xzs(x,z)=-xz 当作核,只取单点 x1=1x_1=1 就得到 Gram 矩阵 [1][-1]。选 c=1c=1cTKc=1c^\mathsf TKc=-1,立即违反正半定条件。该函数虽然对称,也能给某些输入“相似度”数值,却不能作为普通内积核直接放入依赖凸性的核算法。

有效核可以按保持正半定的规则组合。若 k1,k2k_1,k_2 都是核且 a,b0a,b\ge0,则

k(x,z)=ak1(x,z)+bk2(x,z)k(x,z)=a k_1(x,z)+b k_2(x,z)

仍是核,它相当于拼接并缩放两组特征。逐点乘积 k1(x,z)k2(x,z)k_1(x,z)k_2(x,z) 也保持有效,可由特征张量积或 Gram 矩阵的 Schur 乘积理解。对任意输入变换 ggk1(g(x),g(z))k_1(g(x),g(z)) 仍有效,因此可先用领域规则提取具有不变性的表示再计算核。

相反,两个核相减、给核值逐项套任意非线性函数,或按样本标签修改部分矩阵元素,都不自动保持正半定。只修补当前训练矩阵的特征值也未必定义了可供新样本使用的一致核函数。若核的权重、输入变换或组合方式需要从数据学习,它们应只在训练流程内拟合,并为验证点与测试点提供同一外推公式;否则即使最终矩阵正半定,评估仍可能泄漏。

三种常用核及其尺度

线性核

k(x,z)=xTzk(x,z)=x^\mathsf Tz

对应原始特征本身。它适合维度已经很高、关系接近线性或需要低预测成本的情形。若还需要截距,通常由学习算法单独处理;把常数加到核中则等价于增加常数特征,两种参数化的正则化含义要保持一致。

多项式核常写为

k(x,z)=(γxTz+r)p,k(x,z)=(\gamma x^\mathsf Tz+r)^p,

其中整数 p1p\ge1 控制最高次数,rr 决定是否包含低阶项,γ\gamma 调整输入尺度。它表达固定次数的全局交互,维数随原始维度和次数快速增长。并非任意实数幂或负参数组合都自动产生实值正半定核,使用时应遵循已证明的参数条件。

RBF 核写为

k(x,z)=exp ⁣(xz22σ2)=exp(γxz2),γ=12σ2.k(x,z)=\exp\!\left(-\frac{\lVert x-z\rVert^2}{2\sigma^2}\right) =\exp(-\gamma\lVert x-z\rVert^2), \qquad \gamma=\frac1{2\sigma^2}.

σ\sigma 很大时,多数样本彼此相似,函数变化较平缓;σ\sigma 很小时,非完全相同样本核值迅速接近零,模型可形成局部而曲折的变化。它对应无限维特征空间,但实际算法仍只保存有限训练样本系数。无限维不等于无限容量不受控制,正则化和核尺度共同决定可用函数复杂度。

例 3:同一个 RBF 参数不能跨单位照搬

一维两点 x=0,z=2x=0,z=2,取 σ=1\sigma=1,则

k(x,z)=e20.1353.k(x,z)=e^{-2}\approx0.1353.

若原变量由米改成厘米,数值变为 0 与 200,而仍把 σ\sigma 留为 1 厘米,核值变成 e20000e^{-20000},数值上几乎为零。物理对象没变,模型相似结构却完全改变。

一种做法是把 σ\sigma 同步换成 100 厘米;另一种做法是只用训练集统计量标准化每个特征,再在验证协议内选择 σ\sigmaγ\gamma。对混合单位和高维特征,还要判断欧氏距离是否具有任务意义;标准化不能修复错误的距离语义。

再生性质的准确导论

再生核 Hilbert 空间 Hk\mathcal H_k 是由函数构成的 Hilbert 空间,点值评价必须连续。对每个 xx,存在核截面 kx()=k(x,)k_x(\cdot)=k(x,\cdot) 属于该空间,并满足再生性质

f(x)=f,kxHk.f(x)=\langle f,k_x\rangle_{\mathcal H_k}.

f=kzf=k_z 就得到

k(z,x)=kz,kxHk.k(z,x)=\langle k_z,k_x\rangle_{\mathcal H_k}.

“再生”指函数值由与核截面的内积恢复,不是说核会从数据中自动生成新样本。每个正半定核对应一个唯一到等距同构意义的再生核空间;同一输入域上的不同核定义不同函数集合或不同范数,因此也定义不同平滑性与复杂度偏好。

表示定理的常见形式说明:若目标只通过有限训练点的函数值依赖 ff,并加上关于 fHk\lVert f\rVert_{\mathcal H_k} 的严格递增惩罚,则至少一个最优解具有

f()=i=1nαik(xi,)f(\cdot)=\sum_{i=1}^n\alpha_i k(x_i,\cdot)

的有限展开。结论依赖目标与惩罚条件,不表示再生核空间中每个函数本来就只用当前 nn 个样本展开,也不保证系数唯一。它把可能无限维的优化缩减为有限系数问题,是许多核岭回归和支持向量机算法的基础。

中心化与归一化必须匹配训练空间

若要在特征空间中减去训练均值 ϕˉ=n1iϕ(xi)\bar\phi=n^{-1}\sum_i\phi(x_i),不能只把输入坐标均值减掉后假设等价。令

H=I1n11T,H=I-\frac1n\mathbf1\mathbf1^\mathsf T,

训练 Gram 矩阵中心化为

Kc=HKH.K_c=HKH.

用于新样本时还要使用训练集的行均值、总均值和该新样本对训练点的核均值应用同一中心化公式,不能把测试样本并入均值。中心化常用于核主成分等依赖协方差结构的方法;带独立截距的监督模型是否需要中心化,应由其目标推导。

例 4:不用显式坐标完成特征中心化

两个一维特征值为 1 与 3,线性核的 Gram 矩阵是

K=[1339].K=\begin{bmatrix}1&3\\3&9\end{bmatrix}.

特征均值为 2,中心化后显式值是 1,1-1,1,所以预期中心 Gram 矩阵为

[1111].\begin{bmatrix}1&-1\\-1&1\end{bmatrix}.

使用 H=12[1111]H=\tfrac12\begin{bmatrix}1&-1\\-1&1\end{bmatrix} 直接计算 HKHHKH 得到同一结果。中心矩阵每行每列和均为零,可作为实现检查。若用测试点重新计算均值,训练表示会随预测请求改变,造成数据泄漏和不可复现。

核归一化常写为

k~(x,z)=k(x,z)k(x,x)k(z,z).\widetilde k(x,z)=\frac{k(x,z)}{\sqrt{k(x,x)k(z,z)}}.

它让非零特征向量的自相似为 1,相当于比较特征空间夹角。若某点 k(x,x)=0k(x,x)=0,分母为零,必须明确拒绝或另行处理。归一化会改变模型与 RKHS 范数,不是纯数值美化;中心化和归一化的先后次序也会产生不同核,应作为管线参数固定。

计算成本与核选择边界

nn 个样本,完整 Gram 矩阵有 n2n^2 个元素,存储为 Θ(n2)\Theta(n^2)。若一次核计算成本为 O(d)O(d),构建矩阵通常为 O(n2d)O(n^2d);朴素线性代数求解可达 O(n3)O(n^3),具体成本随算法、稀疏性和迭代方法变化。预测还要对若干训练点求核,支持向量或展开系数很多时延迟明显。

当显式特征维度可控时,直接线性学习可能更省内存并适合批量硬件。大样本核方法可使用 Nyström、随机 Fourier 特征或分块迭代等近似,但近似秩、随机种子、误差和验证协议都要报告。核矩阵缓存能节省重复计算,却不能消除二次存储上限。

核应结合不变性、距离语义、样本规模和部署预算选择,再在训练内部验证集或交叉验证中调参数。不得查看最终测试集后挑核。高验证分数也不能证明核表达了真实机制;分布偏移、异常尺度和相似度捷径仍需分组诊断。非负加权和、乘积及输入变换下的已知核可构造新核,但每一步都要保持正半定条件。

常见误区

常见误区

“相似度函数只要对称就能当核。”有效内积核还必须对任意有限样本产生正半定 Gram 矩阵。

常见误区

“核技巧一定更快,因为不构造特征。”它可能省去高维坐标,却引入样本数平方的 Gram 存储和大量核评估。

常见误区

“RBF 核没有参数化假设,因此是通用默认答案。”欧氏距离、特征尺度、宽度和正则化共同施加强归纳偏置。

练习:从内积到可用核

练习

为二维核 (1+xTz)2(1+x^\mathsf Tz)^2 写出一个显式特征映射并核对内积。

查看提示
展开平方,并把交叉项系数二分到两个特征坐标的乘积中。
查看解答
可取 ϕ(x)=(1,2x1,2x2,x12,2x1x2,x22)\phi(x)=(1,\sqrt{2}x1,\sqrt{2}x2,x1^{2},\sqrt{2}x1x2,x2^{2})。逐项内积得到 1+2x1z1+2x2z2+x12z12+2x1x2z1z2+x22z221+2x1z1+2x2z2+x1^{2}z1^{2}+2x1x2z1z2+x2^{2}z2^{2},正好等于 (1+xTz)2(1+x^{\mathsf T}z)^{2}
练习

证明任何由内积特征映射构造的 Gram 矩阵都是正半定矩阵。

查看提示
把二次型重新组合成特征向量线性组合的范数平方。
查看解答
Kij=ϕ(xi),ϕ(xj)Kij=\langle \phi(xi),\phi(xj)\rangle,有 cTKc=Σijcicjϕ(xi),ϕ(xj)=Σiciϕ(xi)20c^{\mathsf T}Kc=\Sigma ij ci cj\langle \phi(xi),\phi(xj)\rangle=||\Sigma i ci\phi(xi)||^{2}\ge 0,因此任意有限样本的 Gram 矩阵正半定。
练习

用最小反例说明 s(x,z)=xzs(x,z)=-xz 不能作为实内积核。

查看提示
只需寻找一个有限样本和一个系数向量使二次型为负。
查看解答
s(x,z)=xzs(x,z)=-xz,取单点 x=1,Gram 矩阵为 [-1];系数 c=1 时 cTKc=1c^{\mathsf T}Kc=-1,所以它不满足核的正半定条件。对称性本身不充分。
练习

推导训练 Gram 矩阵的中心化形式,并解释新样本为何不能重新估计均值。

查看提示
使用 H=I11T/nH=I-11^{\mathsf T}/n,并检查结果行列和为零。
查看解答
训练特征减均值对应左、右各乘一次 H,因此 Kc=HKHKc=HKH。新样本必须减训练均值,公式使用训练 Gram 的行均值、总均值及新点对训练集核均值;若把测试点纳入均值就改变训练表示并泄漏。
练习

某特征数值因单位转换扩大 aa 倍,RBF 核参数应怎样变化才能保持该维相似度?

查看提示
距离扩大 a 倍时,平方距离扩大 a2a^{2} 倍;让宽度同步缩放。
查看解答
保持同一相似度需把 σ\sigma 同步扩大 a 倍,或把 γ\gamma 缩小到原来的 1/a21/a^{2}。更稳妥的是用训练统计量标准化并在训练内验证宽度;测试集不能参与尺度估计。
练习

说明为什么一个核计算很便宜仍不代表完整核学习适合百万样本,并提出可验证替代方案。

查看提示
分别计算 Gram 构造、存储、求解和单次预测所依赖的样本数。
查看解答
完整 Gram 存储 Θ(n2)\Theta(n^{2}),若核评估 O(d)O(d),构造 O(n2d)O(n^{2}d),朴素求解可能 O(n3)O(n^{3});预测需对展开中非零训练点求核。n 很大时可考虑显式有限特征、低秩近似或随机特征,并用独立验证量化近似误差和延迟。

知识连接与资源

课程 · 年份待核

Stanford CS229 Course Materials

Andrew Ng

用于核对经典机器学习模型的目标函数、推导和适用前提。

打开官方来源

Stanford CS229 官方材料可用于核对特征映射、有效核、核技巧和支持向量展开。再生核空间的完整构造与表示定理证明还需遵守其函数空间和正则条件。