陈述
设 A∈Mm×n(R)(或 C),秩为 r。则存在正交(酉)矩阵 U∈Mm、V∈Mn 与「对角」矩阵 Σ(只有前 r 个对角元非零)使
A=UΣV†,Σ=diag(σ1,…,σr,0,…),σ1≥⋯≥σr>0σi 叫 A 的奇异值,它们是 A†A(或 AA†)的特征值的平方根。
说明
⭐ SVD 的几何读法:任何线性映射 = 旋转/反射 → 沿坐标轴拉伸 → 再旋转/反射。于是它把任意矩阵化成了「U 与 V 取好基之后」的对角形 —— 这是等价(A↦PAQ)分类的完整答案(对角元即奇异值,非零个数即秩)。
⭐ 与谱定理的分工:谱定理处理方形且自伴(相似分类);SVD 处理任意形状(等价分类)。前者要求 m=n 且 A∗=A,后者没有任何要求。
⭐ 用途(数值线性代数里最重要的一条定理):降秩逼近(Eckart–Young:截断到前 k 个奇异值给出最优低秩逼近)、主成分分析、伪逆、条件数。
谁指向它(入边)
推出 1
-
设 A∈Mm×n(R),秩 r。
- ATA∈Mn 是对称半正定的(xTATAx=∥Ax∥2≥0),故由谱定理存在标准正交基 v1,…,vn 与 λ1≥⋯≥λn≥0 使 ATAvi=λivi。
- 令 σi:=λi,设 σ1≥⋯≥σr>0=σr+1=⋯。对 i≤r 令 ui:=Avi/σi,则 ∥ui∥=1 且 ⟨ui,uj⟩=δij(用 ⟨Avi,Avj⟩=λjδij)。
- 把 u1,…,ur 扩充成 Rm 的标准正交基 u1,…,um;把 v1,…,vn 补全成标准正交基。则 Avi=σiui(i≤r)且 Avi=0(i>r,因为 ∥Avi∥2=λi=0)。
- 以 vi 为定义域基、ui 为值域基写出矩阵:A 在这两组基下是 Σ=diag(σ1,…,σr,0,…),即 A=UΣVT。∎
关键观察只有一条:ATA 半正定 —— 于是谱定理给出正交的特征基,而 A 把这组基送到一组正交向量上(差伸缩 σi)。SVD 的全部内容就是这个观察的坐标化。
参考
- Horn–Johnson, Matrix Analysis, Ch. 7
- Roman, Advanced Linear Algebra, Ch. 11