📖 本章概要
截断型奇异值分解(Truncated SVD)是完成主成分分析的重要数学工具。 本章从几何视角回顾SVD分解,介绍完全型、经济型、紧缩型、截断型四种SVD形式, 详解秩一矩阵叠加原理,从投影、格拉姆矩阵、瑞利商、F-范数等优化视角理解SVD, 最后用鸢尾花图片演示截断SVD的图像压缩应用。
几何视角看奇异值分解
矩阵形式 • 线性变换 • 旋转→缩放→旋转
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、降维、图像压缩等应用的数学基础。
四种SVD分解形式
完全型 • 经济型 • 紧缩型 • 截断型
15.2节介绍奇异值分解的四种主要形式及其演化关系。 从完全型SVD(U和V都是正交方阵)开始, 逐步演化为经济型(去除U中的零列), 再到紧缩型(只保留非零奇异值), 最后得到截断型(只保留前p个最大奇异值用于降维)。
截断SVD的几何视角
秩一矩阵叠加 • 投影视角 • 秩的增长
15.3节从几何角度理解截断型SVD分解。 原始数据矩阵X可以写成D个秩一矩阵的叠加:X = Σsⱼuⱼvⱼᵀ。 每叠加一个秩一矩阵,重建矩阵X̂的秩就增加1。 从投影视角看,X朝vⱼ投影的结果为sⱼuⱼ,即XV = US。 用前p个主成分重建数据,相当于只保留前p个最大奇异值对应的秩一矩阵。
优化视角看截断型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节用截断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特征提取