研讨会系列 4 · Kalman Filter L1

矩阵特征值、SVD与正定性

如果把矩阵看作线性变换,特征值就是拉伸系数,特征向量就是拉伸方向——这是理解 Kalman Filter 一切数学的起点。

📅 2025-11-22 ⏱ ~60 分钟 📊 中级 🎯 KF 数学基础

🎯 学完本节,你应能

  • 解释特征值和特征向量的几何意义
  • 理解 SVD 为何是"将矩阵拆成旋转×拉伸×旋转"
  • 用特征值判据判断矩阵的正定性
  • 理解谱分解和 Cholesky 分解的基本思想

📖 前置知识

  • 大学线性代数基础(矩阵乘法、行列式、逆矩阵)
  • 了解向量的线性变换概念

💡 建议先看 0000–0003 建立 KF 学习的动机,再进入本讲。

1. 特征值的几何直觉:矩阵 = 线性变换

矩阵 $A \in \mathbb{R}^{n \times n}$ 可以看作一个线性变换:它将一个向量 $v$ 映射为 $Av$。变换的效果可能是旋转、拉伸、剪切或者它们的组合。

💡 核心直觉:特征向量 $v_{i}$ 是"被 $A$ 变换后方向不变的向量",特征值 $\lambda_{i}$ 是对应的"拉伸倍数"。

$$ A v_{i} = \lambda_{i} v_{i}, \quad i = 1,\dots,n $$

也就是说:在 $v_{i}$ 这个方向上,矩阵 $A$ 的效果退化为标量乘法。

【课堂例子】旋转矩阵没有(实)特征值

考虑 $A = \begin{bmatrix}0 & -1 \\ 1 & 0\end{bmatrix}$(旋转 90°)。任何非零向量被旋转后方向都变了,所以没有实特征值。

求法

$$ \det(\lambda I - A) = 0 \quad\Rightarrow\quad \lambda_{1}, \dots, \lambda_{n} $$

解特征多项式得到 $n$ 个特征值(可能有重复、可能有复数)。

2. 特征值的性质与迹-行列式关系

🔢 迹 = 特征值之和

$$ \sum_{i=1}^n \lambda_{i} = \operatorname{Tr}(A) $$

📊 行列式 = 特征值之积

$$ \prod_{i=1}^n \lambda_{i} = \det(A) $$

矩阵变换对特征值的影响

变换特征值变化
$A^{-1}$$\{\lambda_{i}^{-1}\}$
$A^{T}$$\{\lambda_{i}\}$(相同)
$A^{2}$$\{\lambda_{i}^{2}\}$
$A^{k}$$\{\lambda_{i}^{k}\}$
$cI + A$$\{c + \lambda_{i}\}$

这些性质后面在分析系统稳定性时非常有用——KF L3 会用到。

3. 谱分解/相似变换:对角化的条件

如果 $A$ 有 $n$ 个线性无关的特征向量,则 $A$ 可以对角化:

$$ P^{-1} A P = \begin{bmatrix} \lambda_{1} & & \\ & \ddots & \\ & & \lambda_{n} \end{bmatrix} \triangleq D, \quad P = [v_{1} \; \cdots \; v_{n}] $$

其中 $P$ 的列是特征向量。

⚠️ 注意:并非所有矩阵都可对角化。

反例:Jordan 块

$$ A = \begin{bmatrix} 1 & 1 \\ 0 & 1 \end{bmatrix} $$

这个矩阵只有一个特征向量 $(1,0)^T$,无法张成二维空间。如果能对角化,则 $P^{-1}AP = I \Rightarrow A = I$,矛盾!

🔗 连接 KF:谱分解是理解系统稳定性的关键——$A^k = P D^k P^{-1}$,因此 $x_{k} \to 0$ 当且仅当所有 $|\lambda_{i}| < 1$。

4. SVD:旋转·拉伸·旋转

SVD 是线性代数中最优美的分解之一,适用于任何矩阵(不要求方阵、不要求可对角化):

$$ A = U \Sigma V^T = \sum_{i=1}^{r} \sigma_{i} \, u_{i} v_{i}^T $$

其中 $\sigma_{i}$ 是奇异值($A A^T$ 的特征值的平方根),$U$ 和 $V$ 是正交矩阵。

💡 三层理解:

  1. 符号:$A = U\Sigma V^T$,$U$ 和 $V$ 是正交矩阵,$\Sigma$ 是对角阵
  2. 含义:任何线性变换可以分解为"旋转($V^T$) → 拉伸($\Sigma$) → 旋转($U$)"
  3. 直觉:SVD 告诉你"矩阵在每个方向上的放大倍数"——奇异值越大,该方向越重要

【课堂例子】低秩近似

施老师在课堂上用 SVD 做低秩近似:只保留前 $k$ 个最大的奇异值:

$$ A \approx \sum_{i=1}^{k} \sigma_{i} \, u_{i} v_{i}^T $$

这在图像压缩、PCA、数据降维中都是核心思想。

5. 正定矩阵:特征值 > 0 的几何意义

设 $A \in \mathbb{R}^{n \times n}$ 是对称矩阵。

性质定义特征值条件
正定 (p.d.)$x^T A x > 0,\; \forall x \neq 0$$\lambda_{i} > 0$
半正定 (p.s.d.)$x^T A x \ge 0,\; \forall x \neq 0$$\lambda_{i} \ge 0$

🌰 通俗例子:$A = \begin{bmatrix}1 & 0 \\ 0 & 1\end{bmatrix}$ 是正定的。因为 $x^T A x = x_{1}^2 + x_{2}^2 > 0$ 对任意非零 $x$ 成立。而 $A = \begin{bmatrix}1 & 1 \\ 0 & 1\end{bmatrix}$ 不是正定的——它不对称!

正定矩阵的性质

  1. 若 $A$ 正定,$\alpha > 0$,则 $\alpha A$ 正定
  2. 若 $A,B$ 半正定,则 $A+B$ 半正定
  3. 若 $A$ 正定,$B$ 半正定,则 $A+B$ 正定
  4. 若 $A$ 半正定,$C \in \mathbb{R}^{n \times m}$,则 $C^T A C$ 半正定
  5. 若 $A$ 正定,则 $A^{-1}$ 正定(课后可证)

🔗 连接 KF:协方差矩阵一定是半正定的($P \ge 0$)。正定性保证了二次型优化有唯一最小值——KF L3 中你会看到 $P$ 的正定性如何影响 Kalman Gain 的计算。

6. Cholesky 分解

如果 $A$ 是半正定矩阵,则存在下三角矩阵 $L$(实对角元非负)使得:

$$ A = L L^T $$

应用:解线性方程 $A x = y$

Cholesky 分解将求解转化为两步回代(无需计算逆矩阵):

$$ A x = y \;\Longrightarrow\; \begin{cases} L z = y \\ L^T x = z \end{cases} $$

递归计算过程(板书核心)

$$ A = \begin{bmatrix} A_{11} & A_{12} \\ A_{21} & A_{22} \end{bmatrix},\quad L = \begin{bmatrix} L_{11} & 0 \\ L_{21} & L_{22} \end{bmatrix} $$ $$ L_{11} = \sqrt{A_{11}},\quad L_{21} = A_{21} / L_{11} $$ $$ L_{22} L_{22}^T = A_{22} - L_{21} L_{21}^T $$

不断对右下角子矩阵递归,直到 $n \to 1$。这是数值线性代数的经典算法。

7. 【补充引理】$AB$ 与 $BA$ 的非零特征值相同

设 $A \in \mathbb{R}^{n \times m}$,$B \in \mathbb{R}^{m \times n}$,$n > m$。则:

$$ \lambda(AB) = \{\lambda(BA)\} \cup \{0\} $$

说明 $AB$ 有 $n$ 个特征值,其中 $n-m$ 个为零;$BA$ 有 $m$ 个非零特征值,与 $AB$ 的相同。

证明思路(板书)

构造相似变换 $P^{-1}MP$,其中 $P = \begin{bmatrix} I_n & A \\ 0 & I_m \end{bmatrix}$,$M = \begin{bmatrix} AB & 0 \\ B & 0 \end{bmatrix}$:

$$ P^{-1} M P = \begin{bmatrix} I_n & -A \\ 0 & I_m \end{bmatrix} \begin{bmatrix} AB & 0 \\ B & 0 \end{bmatrix} \begin{bmatrix} I_n & A \\ 0 & I_m \end{bmatrix} = \begin{bmatrix} 0 & 0 \\ B & BA \end{bmatrix} $$

这里 $P$ 可逆,故 $M$ 与 $\begin{bmatrix} 0 & 0 \\ B & BA \end{bmatrix}$ 相似,特征值相同。

两个矩阵相似,故特征值相同。于是 $\lambda(AB) \cup \{0\} = \lambda(BA) \cup \{0\}$。

推论:$\operatorname{Tr}(AB) = \operatorname{Tr}(BA)$。更一般地,$\operatorname{Tr}(ABCD) = \operatorname{Tr}(DABC)$(循环不变性)。

💡 直觉:这个引理在后面矩阵求逆引理(KF L2)的推导中非常关键。

📝 课后练习

  1. 【数值计算】给定 $A = \begin{bmatrix} 2 & 1 \\ 1 & 2 \end{bmatrix}$,求其特征值和特征向量,并判断 $A$ 是否正定。

    查看答案

    特征方程: $\det(\lambda I - A) = (\lambda-2)^2 - 1 = 0$ → $\lambda^2 - 4\lambda + 3 = 0$

    $\lambda_{1} = 3$, $v_{1} = (1,1)^T$;$\lambda_{2} = 1$, $v_{2} = (1,-1)^T$

    两个特征值都 > 0,且 $A$ 对称,所以 $A$ 正定。

  2. 【正定性判断】为什么协方差矩阵 $\Sigma = E[(X-\mu)(X-\mu)^T]$ 一定是半正定的?

    查看答案提示

    对任意 $y \neq 0$,$y^T \Sigma y = E[y^T(X-\mu)(X-\mu)^T y] = E[(y^T(X-\mu))^2] \ge 0$。平方的期望非负,所以半正定。

  3. 【推导】若 $A$ 正定,证明 $A^{-1}$ 也正定。(板书留作 homework)

    查看提示

    $A$ 正定 $\Rightarrow$ $\lambda_{i}(A) > 0$ $\Rightarrow$ $\lambda_{i}(A^{-1}) = 1/\lambda_{i}(A) > 0$。对称性继承,所以 $A^{-1}$ 正定。

  4. 【编程思考】用 Python/NumPy 对随机生成的数据矩阵做 SVD,观察奇异值的衰减速度。保留前 $k$ 个奇异值重构原矩阵,看看重构误差如何随 $k$ 变化。

📌 小结

🔑 特征值

特征向量是"不变方向",特征值是"拉伸倍数"。

🌀 SVD

$A = U\Sigma V^T$:旋转→拉伸→旋转。任何矩阵都能分解。

✅ 正定性

对称矩阵 + 所有特征值 > 0 = 正定。协方差矩阵天然半正定。

👉 下节预告:KF L2 将学习矩阵求逆引理(MIL)——这是后面推导 Kalman Gain 的关键工具。你会看到 MIL 如何把"大矩阵求逆"转化为"小矩阵求逆",大幅降低计算量。

施老师大课堂 · 研讨会系列 4 · KF L1 矩阵特征值、SVD 与正定性

字幕来源:subtitles/20251115_施凌_研讨会系列4_KF_L1_audio_34eaee802b23f1a1.srt

板书来源:pdf-md/20251122_施凌_研讨会系列4_KF_L1 md文档/