高斯过程(Gaussian Process)
字数 2495 2025-10-27 23:30:49
好的,我们开始学习新的词条:高斯过程(Gaussian Process)。
高斯过程是概率论和统计学中一个非常重要的概念,它为我们提供了一种在函数空间上进行推理的强大框架。我们将从最基础的概念开始,逐步深入。
步骤 1:从随机变量到随机过程
首先,我们需要理解“过程”的含义。
- 随机变量:您已经学过概率论,知道随机变量(例如 X)是一个其数值由随机事件决定的变量。比如,掷一个骰子的结果就是一个随机变量。
- 随机向量:当我们把多个随机变量放在一起(例如 X = [X₁, X₂, ..., Xₙ]),就形成了一个随机向量。它可以表示同一时间点上多个相关随机事件的集合。
- 随机过程:现在,想象一下随机变量不再只是与一个“点”相关,而是与一个连续的“索引”(最常见的是时间 t 或空间位置 x)相关联。随机过程就是一组无限多个、以某种方式相关的随机变量的集合。例如,我们考虑“明天全天的温度变化”,对于每一个时间点 t(如下午2:00,下午2:01,...),温度 T(t) 都是一个随机变量。这整个集合 {T(t) | t ∈ [0, 24]} 就是一个随机过程。
步骤 2:多元高斯分布(多元正态分布)
高斯过程的核心是高斯(正态)分布。我们从熟悉的开始。
- 一元高斯分布:一个随机变量 X 服从高斯分布,由其均值 μ 和方差 σ² 完全描述。它描述了单个随机变量的分布情况。概率密度函数是著名的钟形曲线。
- 多元高斯分布:一个随机向量 X = [X₁, X₂, ..., Xₙ] 服从多元高斯分布。它由两个参数完全描述:
- 均值向量 μ:每个分量随机变量的均值组成的向量,μ = [μ₁, μ₂, ..., μₙ]。
- 协方差矩阵 Σ:一个 n×n 的矩阵,其第 (i, j) 个元素 Σᵢⱼ 是 Xᵢ 和 Xⱼ 之间的协方差。这个矩阵描述了所有随机变量两两之间的线性相关程度。
- 对角线上的元素 Σᵢᵢ 是 Xᵢ 的方差。
- 非对角线上的元素 Σᵢⱼ (i≠j) 决定了 Xᵢ 和 Xⱼ 的协同变化趋势。
关键性质:多元高斯分布有一个极其优美的性质——其任意子集也服从(多元)高斯分布。例如,从 X 中取出前两个变量 [X₁, X₂],它们的联合分布也是一个高斯分布。
步骤 3:高斯过程的定义
现在,我们将步骤 1 和步骤 2 结合起来。
- 核心思想:一个随机过程是无限多个随机变量的集合。如果我们要求这个集合中任意有限个随机变量的联合分布都是多元高斯分布,那么这个随机过程就被称为高斯过程。
- 直观理解:想象一个函数 f(x)。对于任何一个输入点 x,f(x) 的输出值不是一个确定的数,而是一个随机变量。如果我们在一系列输入点 {x₁, x₂, ..., xₙ} 上观察这个函数,得到的输出集合 {f(x₁), f(x₂), ..., f(xₙ)} 是一个随机向量。高斯过程要求,无论你选择哪 n 个点,这个随机向量的分布都是一个 n 维的高斯分布。
步骤 4:高斯过程的参数
既然任意有限个点的集合都服从高斯分布,那么如何描述整个高斯过程呢?类似于一元高斯分布由均值和方差描述,多元高斯分布由均值向量和协方差矩阵描述,高斯过程也由两个要素完全确定:
- 均值函数 m(x):这是一个确定的函数,给出了在任意点 x 处,随机函数 f(x) 的期望值(平均值)。即 m(x) = E[f(x)]。
- 协方差函数 k(x, x’):也称为核函数。这是一个有两个输入的函数,它定义了任意两个点 x 和 x’ 处的函数值 f(x) 和 f(x’) 之间的协方差。即 k(x, x’) = E[(f(x) - m(x))(f(x’) - m(x’))]。
核函数的重要性:核函数 k(x, x’) 是高斯过程的灵魂。它编码了我们关于函数 f 的先验信念,特别是:
- 平滑度:如果两个点 x 和 x’ 很接近,那么 f(x) 和 f(x’) 的值应该也很接近吗?核函数决定了这种相似性随距离衰减的速度,从而控制了函数的平滑程度。
- 周期性、趋势等:通过选择不同的核函数,我们可以让高斯过程优先表示具有特定结构(如周期性变化)的函数。
步骤 5:高斯过程回归(预测)
高斯过程最强大的应用之一是作为非线性回归的工具,也称为克里金法。
场景:我们有一些观测到的数据点(输入 X 和对应的输出 y),我们想预测在新的输入点 x* 上,函数值 f(x*) 会是多少。
高斯过程回归的流程:
- 先验:在没有任何数据之前,我们假设函数 f 服从一个高斯过程,其均值为 0(或某个简单函数),并选择一个合适的核函数(如径向基函数核,即高斯核)来体现我们对函数平滑度的假设。这定义了我们的“先验分布”。
- 后验:当我们观察到训练数据 (X, y) 后,根据贝叶斯定理,我们可以计算出在已知这些数据的条件下,函数 f 的后验分布。这个后验分布仍然是一个高斯过程!
- 预测:利用多元高斯分布的条件分布性质,我们可以直接写出在后验高斯过程中,在任意新测试点 x* 处的预测分布:
- 这个预测不再是一个单一的值,而是一个完整的概率分布(一个一维高斯分布)。
- 它给出了预测的均值(最可能的值)和方差(不确定性的度量)。
优势:
- 自然的不确定性量化:高斯过程回归不仅给出预测值,还给出了该预测的置信区间。在数据稀疏的区域,方差会变大,表明预测不确定性高。
- 灵活性:通过选择不同的核函数,可以捕捉各种复杂的函数关系。
- 贝叶斯性质:它提供了一个完整的贝叶斯框架,将先验知识与观测数据相结合。
总结
高斯过程 是将多元高斯分布的概念推广到函数空间。它通过一个均值函数和一个协方差函数(核函数) 来完全定义,其核心性质是任意有限点集上的函数值服从联合高斯分布。作为强大的非参数贝叶斯模型,它在回归问题中尤为出色,因为它能提供带有不确定性估计的预测,广泛应用于机器学习、空间统计和优化等领域。