1408 字
7 分钟
线性代数 10|正定、矩阵微分与数值稳定性

数学总导航 · 本课程目录 · ← 上一章

两次称重几乎一样,为什么单件重量会乱跳?#

你想知道两个零件各重多少。第一次称它们各一个,得到 2 克;第二次把第二个零件数量的系数只增加一点点,模型是 x+1.001y=2.001x+1.001y=2.001。联立 x+y=2x+y=2,算出 x=y=1x=y=1。

但第二次读数若从 2.001 抖到 2.002,只差 0.001 克,结果就变成 y=2,x=0y=2,x=0。观测误差很小,推断误差却很大,因为两条测量几乎在问同一个问题。

矩阵仍可逆,但这是一种难以稳定求解的测量设计。条件数是在量化它有多敏感;正则化是在承认信息不足后,增加偏好来压住不合理的大幅波动。

测量为什么会把误差放大一千倍#

两条零件测量相减,得到 0.001y=b2−b10.001y=b_2-b_1。因此 y=1000(b2−b1)y=1000(b_2-b_1)。你最终依赖的是两个很接近读数的微小差值,而差值里也带着测量噪声。

如果 b2b_2 有 10−310^{-3} 的误差,yy 就有 1 的误差。问题不在消元写得不够漂亮,而在测量本来就几乎不能区分两个零件。重新设计差别更大的测量,可能比更高精度算术更有效。

一个碗为什么会又长又窄#

平方损失的等高线若呈细长椭圆,沿窄方向变化一点损失就很大,沿长方向移动很多损失却几乎不变。Hessian 的大、小特征值分别衡量这两种弯曲强度。

梯度下降若步长照顾窄方向,就可能在长方向进展很慢;若步长太大,又可能跨过窄谷不断振荡。这把正定矩阵的几何与优化算法的实际表现连在一起。

正则化具体往碗里加了什么#

在原损失上加 λ∥x∥2\lambda\|x\|^2,相当于给每个方向额外增加一份向原点拉回的二次代价。原来几乎平的方向也开始有了曲率。

在奇异方向上,未正则化解使用 1/σ1/\sigma 放大读数;加惩罚后使用 σ/(σ2+λ)\sigma/(\sigma^2+\lambda)。当 σ\sigma 很小时,后者约为 σ/λ\sigma/\lambda,不会像倒数那样爆炸。

但如果真相本来就在一个很弱的方向上,它也会被压小。因此正则化不是免费去噪,而是在信息不足时明确接受偏差,以换取稳定。以后看到某种范数惩罚,应继续追问:它偏爱哪类答案,可能压掉哪类真实信息?

带你做:比较正则化前后的噪声放大#

在一个奇异方向上,机器倍率为 σ=0.001\sigma=0.001。直接反解要乘 1000,若观测误差为 0.0001,恢复误差就达到 0.1。

取 λ=0.01\lambda=0.01,正则化倍率变成 0.001/(0.000001+0.01)≈0.099990.001/(0.000001+0.01)\approx0.09999,同样噪声只贡献约 10−510^{-5}。稳定性改善巨大。

但若真实输入在这个方向上为 1,精确观测为 0.001,正则化恢复值也只有约 0.0001,真信号几乎被全部压掉。这说明如此大的 λ\lambda 对这个弱方向非常保守。

调参时不能只看重建是否平滑,还要看哪些真实结构被压低。一个有用的练习是分别在强奇异方向和弱奇异方向放同样大小的真信号,观察同一个 λ\lambda 对它们的偏差为何不同。

回头整理:定义、工具与练习#

二次型描述一个碗#

实对称矩阵 AA 的二次型为 xTAxx^TAx。正定指对所有 x≠0x\ne0 都严格为正,半正定允许等于零。通过谱分解,它等于 ∑iλiyi2\sum_i\lambda_i y_i^2。

正定矩阵的等高线是椭球;有正有负的特征值产生鞍形。零特征值对应平坦方向。对非对称矩阵,二次型只看到其对称部分:xTAx=xT(A+AT)x/2x^TAx=x^T(A+A^T)x/2。

梯度与 Hessian#

对 f(x)=12xTAx−bTxf(x)=\tfrac12x^TAx-b^Tx,若 AA 对称,则

∇f=Ax−b,Hf=A.\nabla f=Ax-b,\qquad H_f=A.

若 AA 正定,唯一最小点为 A−1bA^{-1}b。一般矩阵的 ∇(xTAx)=(A+AT)x\nabla(x^TAx)=(A+A^T)x,不能无条件写成 2Ax2Ax。

若 AA 正定,存在 Cholesky 分解 A=LLTA=LL^T,所以 xTAx=∥LTx∥2x^TAx=\|L^Tx\|^2。这说明正定二次型也可以理解为另一种坐标几何中的长度平方。

条件数区分可逆与好求解#

对可逆矩阵,κ2(A)=∥A∥2∥A−1∥2=σmax⁡/σmin⁡\kappa_2(A)=\|A\|_2\|A^{-1}\|_2=\sigma_{\max}/\sigma_{\min}。若 Ax=bAx=b 且仅右侧有误差,

∥δx∥∥x∥≤κ2(A)∥δb∥∥b∥(b≠0).\frac{\|\delta x\|}{\|x\|} \le\kappa_2(A)\frac{\|\delta b\|}{\|b\|}\quad(b\ne0).

例:A=diag⁡(1,10−6)A=\operatorname{diag}(1,10^{-6}) 可逆,但第二个方向的观测噪声会被放大一百万倍。数值误差来自问题本身的敏感性,也来自算法;选择好算法不能彻底消除问题本身的病态。

正则化是一种有偏的稳定化#

求解

min⁡x(∥Ax−b∥2+λ∥x∥2),λ>0\min_x\bigl(\|Ax-b\|^2+\lambda\|x\|^2\bigr),\qquad\lambda>0

得到 (ATA+λI)x=ATb(A^TA+\lambda I)x=A^Tb。SVD 坐标中的系数由 1/σi1/\sigma_i 改为 σi/(σi2+λ)\sigma_i/(\sigma_i^2+\lambda),抑制了小奇异值的噪声放大。代价是引入偏差,并需要选择 λ\lambda。

通向泛函分析的门槛#

到这里应能回答:核与像是什么?正交投影为何最优?特征值与奇异值差在哪?有限维中所有范数等价,而无限维中长度的选择会改变连续性和完备性。下一步进入 为什么要把函数当作点。

综合练习与解答#

  1. A=(2112)A=\begin{pmatrix}2&1\\1&2\end{pmatrix} 是否正定,条件数多少?
  2. 为什么 ATAA^TA 总半正定?何时正定?
  3. 序列空间中 T(x)n=xn/nT(x)_n=x_n/n 每个坐标都可逆,为何整体逆可能不稳定?

解答:1. 特征值 3、1,故正定且二范数条件数为 3。2. xTATAx=∥Ax∥2x^TA^TAx=\|Ax\|^2;正定当且仅当 AA 列满秩。3. 逆把第 nn 个方向放大 nn 倍,不存在统一的放大上界。这将是紧算子反问题中的核心例子。


数学总导航 · 本课程目录 · ← 上一章

线性代数 10|正定、矩阵微分与数值稳定性
https://f1yingwhite.github.io/posts/数学/线性代数教程/10-正定、矩阵微分与数值稳定性/
作者
F1yingWhite
发布于
2026-11-03
许可协议
CC BY-NC-SA 4.0