研讨会系列 8 · Kalman Filter L5

贝叶斯估计与 Bayes Rule

从"系统能不能被观测"到"观测到数据后如何修正信念"——Bayes Rule 是估计理论的灵魂。

📅 2025-12-27 ⏱ ~60 分钟 📊 中级 🎯 KF 概率基础

🎯 学完本节,你应能

  • 区分 Bayesian 与 Non-Bayesian 估计的本质差异
  • 解释先验、似然、后验的含义与关系
  • 写出 Bayes Rule 的连续与离散形式
  • 用 Bayes Rule 解决简单的参数估计问题
  • 理解"测量 = 缩小不确定性"的核心思想

📖 前置知识

  • KF L4 可观性(0007)或系统理论基础
  • 条件概率基础($P(A|B)$ 的含义)
  • 连续/离散随机变量的 PDF/PMF 概念

💡 本讲是从"系统理论"到"估计理论"的桥梁课。建议先看 L1-L4 再进入。

1. 估计问题是什么

估计问题的基本形式:

有一个未知量 $\theta$,我们想知道它,但只能通过测量 $X$ 间接推断。

场景未知 $\theta$测量 $X$
房间温度真实温度温度计读数
车辆定位真实位置GPS 信号
KF 中的状态$x_{k}$$y_{k}$
硬币正反正面概率 $p$投掷结果序列
Juliet 迟到最大迟到时间 $\theta$实际迟到时间 $x$

💡 核心问题:给定有限的、带噪的测量,我们到底能多"确定"地推断 $\theta$?这个问题的答案,就是整个估计理论要回答的核心。

2. Bayesian vs Non-Bayesian 估计

📍 Bayesian 估计

把 $\theta$ 视为随机变量,给它一个分布 $p(\theta)$。

目标:更新 $\theta$ 的分布(后验)。

"先相信,再看数据修正"

📍 Non-Bayesian 估计

把 $\theta$ 视为固定未知常数。

目标:找一个最合理的数值。

"数据说话,不做先验假设"

视角$\theta$ 是什么随机性来源目标
Bayesian随机变量$\theta$ 和 $X$ 都是随机的求后验分布 $p(\theta|X)$
Non-Bayesian固定未知常数只有 $X$ 是随机的求 $\theta$ 的点估计值

🔗 连接 KF:Kalman Filter 本质上是 Bayesian 的——每一步都在做"用先验(预测)× 似然(测量)→ 后验(估计)"的递归计算。

3. 【课堂例子】房间温度估计

施老师在黑板上用房间温度的例子直观展示两种框架的差异。

设真实温度为 $\theta$,测量为:

$$ Y = \theta + v, \quad v \sim \mathcal{N}(0, 1) $$

Non-Bayesian

$\theta$ 是固定常数:

$$ Y \sim \mathcal{N}(\theta, 1) $$

$\theta$ 的取值范围是 $[0, 30^\circ\text{C}]$

Bayesian

还要给 $\theta$ 一个先验分布:

$$ \theta \sim p(\theta) $$

例如:认为大概率在 20°C–30°C 之间

问题:Bayesian 框架下,$Y$ 的分布是什么?

答案:需要知道 $p(\theta)$ 才能确定 $Y$ 的边缘分布——这就是 Bayesian 与非 Bayesian 的分水岭。

4. 先验、似然、后验

贝叶斯估计的三大支柱:

📦 先验 $p(\theta)$

观测数据之前,我们对 $\theta$ 的信念。

反映已有知识(或无知)

📐 似然 $p(x|\theta)$

如果 $\theta$ 是真值,观测到 $x$ 的概率。

也叫"传感器模型"

🎯 后验 $p(\theta|x)$

观测到 $x$ 后,对 $\theta$ 的更新信念。

这就是我们要估计的

📊 证据 $p(x)$

观测数据的边缘概率。

就是归一化常数

💡 核心关系:后验 $\propto$ 先验 $\times$ 似然。

$$ p(\theta|x) = \frac{p(\theta) \cdot p(x|\theta)}{p(x)} = \frac{\text{旧信念} \times \text{数据支持}}{\text{归一化}} $$

5. Bayes Rule 的数学形式

形式 1:离散 $\theta$ + 离散 $X$

$$ P(\theta|x) = \frac{P(\theta) P(x|\theta)}{\sum_{\theta'} P(\theta') P(x|\theta')} $$

形式 2:连续 $\theta$ + 连续 $X$

$$ f(\theta|x) = \frac{f(\theta) f(x|\theta)}{\int f(\theta') f(x|\theta') \, d\theta'} $$

形式 3:混合(离散 $\theta$ + 连续 $X$ 等)

$$ p(\theta|x) = \frac{p(\theta) f(x|\theta)}{\sum_{\theta'} p(\theta') f(x|\theta')} \quad\text{或}\quad f(\theta|x) = \frac{f(\theta) p(x|\theta)}{\int f(\theta') p(x|\theta') \, d\theta'} $$

💡 三层阅读:

  1. 符号:分母是求和/积分,保证概率归一化
  2. 含义:后验 = 先验 × 似然 × (1/证据)
  3. 直觉:如果数据在某个 $\theta$ 值下很可能是(似然大),那这个 $\theta$ 的后验概率就会调高;同时如果先验就认为这个 $\theta$ 不太可能,那后验也不会太高。

6. 【课堂例子】抽奖号码与测量减熵

施老师用中奖号码为例,直观展示测量如何逐步缩小不确定性。

测量 1:有人告诉你第一位是 6。

$$ P(\theta|x_{1}) = \begin{cases} 1/10^5 & \text{如果 $\theta$ 以 6 开头} \\ 0 & \text{否则} \end{cases} $$

候选号码从 $10^6$ 缩小到 $10^5$。

测量 2:第二位是 2。候选变为 $10^4$。

... 直到 6 位全知道:

$$ P(\theta|x_{1},\dots,x_{6}) \sim \mathcal{N}(627589, 0) $$

不确定性几乎降到零。

💡 核心启发:测量的作用不是"制造真相",而是缩小可能性空间。每次测量都在减少 $\theta$ 的不确定性(方差减小)。Kalman Filter 的协方差更新 $P_{k|k} < P_{k|k-1}$ 正是这个思想的数学体现。

7. 【课堂例题】Romeo & Juliet 迟到问题

这是 KF L5 的经典例题,板书有完整推导。

问题设定

Romeo 和 Juliet 开始约会。Juliet 每次迟到时间 $X$ 均匀分布在 $[0, \theta]$,其中 $\theta$ 是 Juliet 的"最大迟到时间"。$\theta$ 未知,Romeo 的先验是 $\theta \sim U[0,1]$。

第一次约会,Romeo 观测到 Juliet 迟到 $x$ 小时。Romeo 如何更新对 $\theta$ 的信念?

Step 1:写先验

$$ f(\theta) = \begin{cases} 1, & 0 \le \theta \le 1 \\ 0, & \text{otherwise} \end{cases} $$

Step 2:写似然(传感器模型)

$$ f(x|\theta) = \begin{cases} 1/\theta, & 0 \le x \le \theta \\ 0, & \text{otherwise} \end{cases} $$

Step 3:套 Bayes Rule

后验非零的条件:$x \le \theta \le 1$。

$$ f(\theta|x) = \frac{f(\theta) f(x|\theta)}{\int_{x}^1 f(\theta') f(x|\theta') \, d\theta'} = \frac{1/\theta}{\int_{x}^1 (1/\theta') \, d\theta'} = \frac{1/\theta}{|\log x|} = \frac{1}{\theta \, |\log x|}, \quad x \le \theta \le 1 $$

通俗解释

如果 Juliet 已经迟到了 $x$ 小时,那 $\theta < x$ 是不可能的——因为最大迟到时间不能小于实际迟到时间。后验 $\theta$ 的范围被截断为 $[x, 1]$,且密度在 $\theta = x$ 处最大,随 $\theta$ 增大递减。

极值情况:$x=1$

如果 Juliet 迟到了整整 1 小时,则 $\theta = 1$ 唯一可能——后验退化成一个点。

📝 课堂板书练习:画出 $x=0.1,\; x=0.5,\; x=0.9$ 时的后验分布 $f(\theta|x)$。注意 $x$ 越大,$\theta$ 的定义域越窄,后验越"集中"。

📝 课后练习

  1. 【概念题】为什么 Bayes Rule 的分母被称为"归一化常数"?

    查看答案

    分母确保后验 $p(\theta|x)$ 对 $\theta$ 的积分/求和等于 1。它不依赖于 $\theta$,只依赖于数据 $x$,因此只起归一化作用。

  2. 【数值计算】在 Romeo & Juliet 例子中,设 $x = 0.4$。计算 $f(\theta|x)$ 在 $\theta=0.5$ 和 $\theta=0.8$ 处的值。

    查看计算

    $|\log(0.4)| = | -0.9163| = 0.9163$

    $f(0.5|0.4) = 1/(0.5 \times 0.9163) = 2.183$

    $f(0.8|0.4) = 1/(0.8 \times 0.9163) = 1.364$

    $\theta$ 越小,后验密度越大——因为先验均匀 + 似然 $1/\theta$ 随 $\theta$ 增大而减小。

  3. 【推导】设 $\theta$ 连续、$x$ 离散。推导 Bayes Rule 公式 $f(\theta|x) = \frac{f(\theta) P(x|\theta)}{\int f(\theta') P(x|\theta') d\theta'}$。提示:由联合分布 $f(\theta, x) = f(\theta) P(x|\theta)$ 出发。

    查看提示

    $f(\theta|x) = \frac{f(\theta, x)}{P(x)} = \frac{f(\theta) P(x|\theta)}{\int f(\theta') P(x|\theta') d\theta'}$,由条件概率定义得到。

  4. 【开放思考】在抽奖号码例子中,如果只给了你前 3 位数字,你对 $\theta$ 的不确定性减少了多少?用"可能性空间缩小比例"来量化。

    查看答案

    前 3 位确定,候选从 $10^6$ 缩小到 $10^3$,缩小了 1000 倍。不确定性(熵)从 $\log(10^6)$ 降到 $\log(10^3)$,减少了一半。

📌 小结

🎲 Bayesian 估计

$\theta$ 是随机变量。后验 $\propto$ 先验 $\times$ 似然。

📏 Non-Bayesian 估计

$\theta$ 是常数。目标:找一个最佳数值。

🔗 下节预告

KF L6:MAP、MMSE 与假设检验——从后验分布到点估计的三种方法。

📌 本讲在全系列中的位置:L1–L4 打下系统理论基础,L5 引入估计理论的核心思维——Bayesian 更新。接下来 L6(MAP/MMSE)和 L7(最大似然)将基于这讲的概念构建具体的估计器。最终这些全部汇入 Kalman Filter 的递归 Bayesian 更新框架。

施老师大课堂 · 研讨会系列 8 · KF L5 贝叶斯估计与 Bayes Rule

板书来源:pdf-md/20251227_施凌_研讨会系列8_KF_L5 md文档/

课程笔记:course-notes/08-KF L5:贝叶斯估计与 Bayes Rule.md