研讨会系列 · KF L17 · 第一部分

条件高斯定理与 MMSE 估计

用概率法重新推导 KF —— 条件高斯分布、控制输入的融入、MMSE = LMMSE 等价性证明与分离原理

📅 2026-05-30 🧩 MMSE 📐 条件高斯 🔄 分离原理

🎯 学完本节,你应能

理解条件高斯定理与MMSE=LMMSE的等价关系

📖 前置知识

建议先完成前面的课程,或具备线性代数与概率论基础。

🎯 本讲目标

从概率角度重新推导 KF:复习条件高斯定理,证明在线性高斯系统中 MMSE = LMMSE,并理解控制输入与分离原理。

🔑 核心思想

在线性高斯系统中,条件均值(MMSE)恰好是测量的线性函数——恰好等于 LMMSE。因此 KF 不仅是线性最优的,更是全局最优的。

📖 目录

1. 回顾:几何方法的局限

前几节课我们用几何方法(投影定理)推导了 KF 的五个核心等式,得到了 LMMSE(线性最小均方误差估计)。但有一个关键问题尚未回答:

❓ 核心问题

在线性高斯系统中,LMMSE 是否等于 MMSE(最小均方误差估计,即条件均值 $E[x_{k} | y_{0}, \dots, y_{k}]$)?

回答这个问题至关重要:

本讲将用概率方法(条件高斯分布)来回答这个问题。顺便,我们也会在系统模型中引入控制输入 $B u_{k}$。

2. 系统模型与符号体系(含控制输入)

现在我们将控制输入 $u_{k}$ 加入到系统模型中:

状态方程

$$ x_{k+1} = A x_{k} + B u_{k} + w_{k} $$

$w_{k} \sim N(0, Q)$,$x_{0} \sim N(\bar x_{0}, \Pi_{0})$

测量方程

$$ y_{k} = C x_{k} + v_{k} $$

$v_{k} \sim N(0, R)$,与 $w_{k}$、$x_{0}$ 独立

符号体系

定义信息集合:

$$ Y_{k} = \{y_{0}, y_{1}, \dots, y_{k}\}, \quad U_{k} = \{u_{0}, u_{1}, \dots, u_{k}\} $$

MMSE 估计定义为条件均值:

估计

$$ \hat x_{k|k} = E[x_{k} \mid Y_{k}, U_{k-1}] $$ $$ P_{k|k} = \text{Cov}[x_{k} \mid Y_{k}, U_{k-1}] $$

预测

$$ \hat x_{k|k-1} = E[x_{k} \mid Y_{k-1}, U_{k-1}] $$ $$ P_{k|k-1} = \text{Cov}[x_{k} \mid Y_{k-1}, U_{k-1}] $$

💡 关于 $U_{k-1}$ 的说明

为什么是 $U_{k-1}$ 而不是 $U_{k}$?因为 $u_{k}$ 是在估计 $x_{k}$ 之后才施加的控制,它影响的是 $x_{k+1}$,而不是 $x_{k}$。所以估计 $x_{k}$ 时只能用过去的控制输入。

3. 条件高斯定理回顾

这是整个概率方法推导 KF 的核心定理:

定理:设 $X$ 和 $Y$ 联合高斯:

$$ \begin{bmatrix} X \\\\ Y \end{bmatrix} \sim N\left( \begin{bmatrix} \bar x \\\\ \bar y \end{bmatrix},\\; \begin{bmatrix} \Sigma_{x} & \Sigma_{xy} \\\\ \Sigma_{yx} & \Sigma_{y} \end{bmatrix} \right) $$

则在给定 $Y = y$ 条件下,$X$ 仍然是高斯分布:

$$ \boxed{\mu_{x|y=y} = \bar x + \Sigma_{xy} \Sigma_{y}^{-1} (y - \bar y)} $$ $$ \boxed{\Sigma_{x|y=y} = \Sigma_{x} - \Sigma_{xy} \Sigma_{y}^{-1} \Sigma_{yx}} $$

特例验证

① $X = Y$ 时

$$ \mu_{x|y=y} = y,\quad \Sigma_{x|y=y} = 0 $$

完全确定,不确定性降为零。

② $X \perp Y$ 时

$$ \mu_{x|y=y} = \bar x,\quad \Sigma_{x|y=y} = \Sigma_{x} $$

Y 不提供 X 的任何信息。

🎯 重要意义

注意条件均值 $\bar x + \Sigma_{xy} \Sigma_{y}^{-1} (y - \bar y)$ 是 $y$ 的线性函数。而 LMMSE 是在所有线性估计器中最小化 MSE 的解——它就是那个线性函数!所以在线性高斯系统中,MMSE = LMMSE。

4. Time Update(概率法)

给定 $\hat x_{k-1|k-1}$ 和 $P_{k-1|k-1}$,计算预测:

预测均值

$$ \begin{aligned} \hat x_{k|k-1} &= E[x_{k} \mid Y_{k-1}, U_{k-1}] \\ &= E[A x_{k-1} + B u_{k-1} + w_{k-1} \mid Y_{k-1}, U_{k-1}] \\ &= A E[x_{k-1} \mid Y_{k-1}, U_{k-1}] + B u_{k-1} + E[w_{k-1} \mid Y_{k-1}, U_{k-1}] \\ &= A \hat x_{k-1|k-1} + B u_{k-1} \end{aligned} $$

其中 $E[w_{k-1} \mid Y_{k-1}, U_{k-1}] = 0$ 因为 $w_{k-1}$ 与历史信息不相关(零均值高斯,且独立于所有历史信息)。

预测协方差

$$ \begin{aligned} P_{k|k-1} &= E\big[(x_{k} - \hat x_{k|k-1})(x_{k} - \hat x_{k|k-1})' \mid Y_{k-1}, U_{k-1}\big] \\ &= E\big[(A \tilde x_{k-1|k-1} + w_{k-1})(A \tilde x_{k-1|k-1} + w_{k-1})' \mid Y_{k-1}, U_{k-1}\big] \\ &= A P_{k-1|k-1} A' + Q \end{aligned} $$

交叉项 $E[\tilde x_{k-1|k-1} w_{k-1}'] = 0$ 因为 $w_{k-1}$ 与历史信息独立。

📌 关键观察

控制输入 $B u_{k}$ 只影响预测均值,不影响预测协方差。这意味着控制与估计可以分离(Separation Principle)——设计控制器时可以暂时忽略估计误差。

5. Measurement Update(概率法)

给定 $\hat x_{k|k-1}$ 和 $P_{k|k-1}$,新的测量 $y_{k}$ 到来。我们需要找到 $E[x_{k} \mid Y_{k}, U_{k-1}]$。

步骤 1:构造联合高斯分布

在给定 $Y_{k-1}, U_{k-1}$ 的条件下,$x_{k}$ 和 $y_{k}$ 是联合高斯的:

$$ \begin{aligned} E[x_{k} \mid \cdot] &= \hat x_{k|k-1} \\ E[y_{k} \mid \cdot] &= C \hat x_{k|k-1} \\ \text{Cov}[x_{k} \mid \cdot] &= P_{k|k-1} \\ \text{Cov}[y_{k} \mid \cdot] &= C P_{k|k-1} C' + R \\ \text{Cov}[x_{k}, y_{k} \mid \cdot] &= P_{k|k-1} C' \end{aligned} $$

步骤 2:代入条件高斯公式

$$ \begin{aligned} \hat x_{k|k} &= \hat x_{k|k-1} + P_{k|k-1} C' [C P_{k|k-1} C' + R]^{-1} (y_{k} - C \hat x_{k|k-1}) \\ P_{k|k} &= P_{k|k-1} - P_{k|k-1} C' [C P_{k|k-1} C' + R]^{-1} C P_{k|k-1} \end{aligned} $$

与几何方法得到的表达式完全一致!

💡 概率法 vs 几何法

概率法自动得到了正确结果,不需要投影定理的几何直觉。它只需要知道"在高斯分布下,条件均值就是 MMSE"。

6. MMSE = LMMSE:线性高斯下的等价性

现在我们终于可以回答核心问题了:

✅ 结论

在线性高斯系统中,MMSE 估计等于 LMMSE 估计。

证明思路

  1. MMSE 定义:$\hat x_{k|k}^{\text{MMSE}} = E[x_{k} \mid Y_{k}, U_{k-1}]$,即条件均值
  2. 条件高斯公式告诉我们,条件均值是测量的线性函数
  3. LMMSE 定义:在所有的线性估计器中找均方误差最小的那个
  4. 既然 MMSE 本身(条件均值)已经是线性函数,那么它恰好就是 LMMSE 的最优解

📐 几何方法

投影 → LMMSE

只利用了内积空间结构,不需要概率假设

优点:直观、可视化

局限:只得到 LMMSE

📊 概率方法

条件均值 → MMSE

利用了高斯分布的性质

优点:证明 MMSE = LMMSE

局限:需要高斯假设

🎯 所以结论是

在线性高斯情况下,$$\hat x_{k|k}^{\text{LMMSE}} = \hat x_{k|k}^{\text{MMSE}}$$两者完美重合。这在非高斯系统中一般不成立——非高斯条件下 MMSE 通常严格小于 LMMSE(因为 MMSE 在所有估计器中最优,而 LMMSE 只在线性估计器中最优)。这就是为什么 KF 在高斯假设下如此强大的原因。

7. 分离原理的含义

从上面的推导中,我们还发现了一个重要的性质——分离原理(Separation Principle):

🔄 分离原理

在线性高斯系统中,最优控制器的设计可以分为两个独立的任务:

  1. 估计器设计:设计 KF 估计状态(不依赖于控制目标)
  2. 控制器设计:假设状态完全已知,设计反馈控制器(不依赖于噪声统计)

两者结合即为最优:$u_{k} = -F \hat x_{k|k}$

为什么控制输入不影响协方差?因为在给定 $U_{k-1}$ 条件下,$B u_{k-1}$ 是确定性的确定量——它不引入随机性,不影响误差和协方差。

🧠 估计器

KF 处理不确定性,提供最优状态估计 $\hat x_{k|k}$

🎮 控制器

LQR 等控制器假定状态已知,设计最优反馈 $u_{k} = -F x_{k}$

结合时,将估计值代入控制器:$u_{k} = -F \hat x_{k|k}$。这就是 LQG 控制(Linear Quadratic Gaussian control)的核心结构。

📝 自测题

问题 1

为什么在线性高斯系统中 MMSE = LMMSE?

点击查看答案

在线性高斯系统中,给定测量后的条件均值(MMSE)是测量的线性函数(由条件高斯公式保证)。而 LMMSE 是所有线性估计器中 MSE 最小的。既然 MMSE 本身也是线性的,它必然就是那个最优的线性估计器。

问题 2

控制输入 $B u_{k}$ 对 KF 的协方差 $P_{k|k}$ 有影响吗?为什么?

点击查看答案

没有影响。控制输入 $B u_{k}$ 在给定 $U_{k-1}$ 条件下是确定量,它只影响预测均值,不影响误差和协方差。这是 Separation Principle 的基础——控制律设计和估计器设计可以独立进行。

问题 3

在非高斯系统中,MMSE 和 LMMSE 哪个更小?

施老师大课堂 · 研讨会系列 KF L17 · 条件高斯定理与 MMSE 估计