← 返回目录

Chapter 15: 截断型奇异值分解

Truncated SVD

📖 本章概要

截断型奇异值分解(Truncated SVD)是完成主成分分析的重要数学工具。 本章从几何视角回顾SVD分解,介绍完全型、经济型、紧缩型、截断型四种SVD形式, 详解秩一矩阵叠加原理,从投影、格拉姆矩阵、瑞利商、F-范数等优化视角理解SVD, 最后用鸢尾花图片演示截断SVD的图像压缩应用。

15.1

几何视角看奇异值分解

矩阵形式 • 线性变换 • 旋转→缩放→旋转

15.1节从几何视角回顾奇异值分解(SVD)的基本概念。图1展示了3×2矩阵A的完全型SVD分解结构:A = U @ S @ Vᵀ, 其中U是3×3正交方阵(左奇异向量),S是3×2对角矩阵(奇异值),V是2×2正交方阵(右奇异向量)。图2说明线性变换Ax = y的几何意义,将向量x从ℝ²映射到ℝ³。图3揭示SVD的核心几何直觉:任意矩阵变换可分解为旋转(Vᵀ)→缩放(S)→旋转(U)三步, 这是理解PCA、降维、图像压缩等应用的数学基础。

✦
15.2

四种SVD分解形式

完全型 • 经济型 • 紧缩型 • 截断型

15.2节介绍奇异值分解的四种主要形式及其演化关系。 从完全型SVD(U和V都是正交方阵)开始, 逐步演化为经济型(去除U中的零列), 再到紧缩型(只保留非零奇异值), 最后得到截断型(只保留前p个最大奇异值用于降维)。

✦
15.3

截断SVD的几何视角

秩一矩阵叠加 • 投影视角 • 秩的增长

15.3节从几何角度理解截断型SVD分解。 原始数据矩阵X可以写成D个秩一矩阵的叠加:X = Σsⱼuⱼvⱼᵀ。 每叠加一个秩一矩阵,重建矩阵X̂的秩就增加1。 从投影视角看,X朝vⱼ投影的结果为sⱼuⱼ,即XV = US。 用前p个主成分重建数据,相当于只保留前p个最大奇异值对应的秩一矩阵。

✦
15.4

优化视角看截断型SVD

最大投影 • 格拉姆矩阵 • 瑞利商 • F-范数

15.4节从优化角度理解截断型SVD分解。 寻找第一主成分v₁等价于求解最大L2范数投影:arg max ||Xv||(约束||v||=1)。 这等价于求格拉姆矩阵G=XᵀX的最大特征值,也等价于求瑞利商的最大值。 从F-范数角度,截断SVD给出秩不超过p的最优近似:min ||X - X̂||_F。 最大奇异值s₁与最大特征值λ₁的关系为λ₁ = s₁²。

✦
15.5

鸢尾花照片分析

图像压缩 • 特征脸 • 降维应用

15.5节用截断SVD分析鸢尾花照片。 黑白化后的图像每个像素都是[0,1]范围内的数字,整幅图可看作数据矩阵。 用1个主成分重建时图像难以辨识,16个主成分已能看到鸢尾花样子,64个主成分图像清晰且数据压缩至2.5%。 这种方法在人脸识别的特征脸(Eigenface)技术中应用广泛: 将人脸图像转换为低维特征向量进行识别。

✦

章节总结

🔀

四种SVD形式

  • • 完全型:U(n×n)和V(D×D)都是正交方阵
  • • 经济型:U(n×D)与X形状相同
  • • 紧缩型:只保留非零奇异值(r ≤ min(n,D))
  • • 截断型:只保留前p个最大奇异值(p < r)
🎯

几何与优化视角

  • • 秩一叠加:X = Σsⱼuⱼvⱼᵀ,每项贡献一个秩
  • • 投影视角:X朝vⱼ投影结果为sⱼuⱼ
  • • 优化问题:arg max ||Xv|| (||v||=1)
  • • 等价关系:λ₁ = s₁²,特征值=奇异值平方
⚡

优化等价形式

  • • 最大投影:max ||Xv||,v为单位向量
  • • 格拉姆矩阵:max vᵀGv,G = XᵀX
  • • 瑞利商:max (xᵀGx)/(xᵀx)
  • • F-范数:min ||X - X̂||_F,秩约束
🖼️

图像压缩应用

  • • 1个PC:图像难以辨识
  • • 16个PC:可以看到鸢尾花
  • • 64个PC:清晰,压缩至2.5%
  • • 特征脸:人脸识别的PCA特征提取

🔑 关键公式汇总

经济型SVD(公式3):
X = U @ S @ Vᵀ
n×D = n×D @ D×D @ D×D
秩一叠加(公式4):
X = Σsⱼuⱼvⱼᵀ
j = 1, 2, ..., D
截断重建(公式5):
X̂ ≈ Σ(j=1 to p) sⱼuⱼvⱼᵀ
p < D,近似重建
投影关系(公式9):
Xvⱼ = sⱼuⱼ
||uⱼ|| = 1,||Xvⱼ|| = sⱼ