跳到正文
格致开物MATHWIKI

高斯过程

高斯过程(Gaussian process,GP)是一族随机变量,其任意有限子集都具有联合正态分布。索引可以是时间、空间位置或一般输入;一次抽样得到的是整个索引集合上的一套随机函数值。它由均值函数和协方差函数确定有限维分布,常作为未知函数或随机场的模型。

有限个点上的联合约束

设索引集合为 𝒯。若对任意 n≥1 和 t1,…,tn∈𝒯,都有 (Yt1,…,Ytn)𝖳∼Nn(mT,KT), 其中 (mT)i=m(ti)、(KT)ij=k(ti,tj),就写 Y∼GP(m,k)。联合正态允许退化:某些函数值之间可以存在精确线性关系。[1]

协方差核必须对称,且对任意有限点和实数权重 a1,…,an 满足 ∑i,jaiajk(ti,tj)≥0. 这就是随机线性组合 ∑iaiYti 的方差非负。单独要求 k(t,t)≥0 不够;不同位置间的协方差也受整个矩阵约束。

两个正态数生成的过程

令独立标准正态变量为 Z0,Z1,定义 Yt=Z0+tZ1,t∈ℝ。任意有限线性组合可写为 ∑iaiYti=(∑iai)Z0+(∑iaiti)Z1, 所以仍是正态变量。这证明所有有限向量联合正态。均值为零,核为 k(s,t)=1+st;其半正定性还可直接核对为 ∑i,jaiaj(1+titj)=(∑iai)2+(∑iaiti)2≥0.

固定三对潜在正态数值所生成的三条仿射过程教学路径
每条路径只由截距和斜率两个随机数控制;图中三组数值用于说明结构,不是经验观测。一般高斯过程不一定只由两个随机数决定。

在 t=0,1,2,协方差矩阵为 KT=(111123135). 它秩为二,因为每条路径都满足 Y0−2Y1+Y2=0。矩阵奇异并不使这个高斯过程非法;只是不能用普通三维体积密度描述它。

变式。 把核改为 k(s,t)=1+2st,可用 Yt=Z0+2tZ1 构造。t=0,1 的方差分别为一与三,协方差仍为一。相关系数因此从原来的 1/2 变为 1/3,协方差不变并不表示相关系数不变。

正态边缘不等于联合高斯

令 Z∼N(0,1),另取独立等概率符号 R∈{−1,1}。两个变量 Y1=Z,Y2=RZ 各自都是标准正态,且协方差为零;但它们满足 Y12=Y22,明显不独立。联合分布集中在 y2=±y1 两条直线上,不能是具有单位协方差的二维正态。这说明逐点画出正态直方图,仍不足以验证高斯过程定义。

另一个非法例子是试图令两点方差都为一、协方差为 3/2。此时差值方差为 1+1−2(3/2)=−1,不可能由任何随机变量实现。核的整体半正定条件在建模时承担了这类一致性检查。

核函数与函数的规则性

径向基核、指数核等不同选择,会给函数变化快慢及平滑性带来不同假设。但“每个有限向量正态”本身不要求路径处处光滑,也不自动保证所选版本连续。讨论路径连续性或可微性,需要核的附加条件及相应概率论结果。确定一个合法核后,在有限网格上可通过协方差分解生成函数值;这只是网格上的有限维抽样。

高斯过程回归再加入观测模型并条件化,研究观测以后未知函数的分布;本文的过程定义和核合法性,是那一步之前的先验结构。

参考文献

  1. ↑ C. E. Rasmussen、C. K. I. Williams:Gaussian Processes for Machine Learning,Chapter 2,§2.2高斯过程定义;Chapter 4,§4.1关于协方差核的半正定条件。