定理:每个m×n实矩阵A都可以写成#
A=PSQT
其中P,Q是m阶与n阶正交矩阵,
S=σ1...σr0m×nr=rank(A),σ1≥σ2≥...≥...σr≥0
其中σi为ATA正特征值的算术平方根
S为A的奇异矩阵,σi为A的奇异值,对应的向量为奇异向量
设A为m×n实矩阵,由于AAT实对称
存在正交矩阵P使得
AAT=[β1,β2,...,βm]λ1...λr0PT且AAT的特征值非负:
λ1≥λ2≥...≥λr>λr+1=...=0
由
(ATP)TATP=PTAATP=λ1...λr0知ATP的列向量ATβ,...,ATβm
两两正交且长度分别为
λ1,...,λr,0,...,0
令γi=λi1ATβi
则γ1,...,γr可扩充成Rn的标准正交基
γ1,...,γr,...,γn
⎩⎨⎧γr+1γ1=0γr+1γ2=0...γr+1γr=0解方程即可
记正交矩阵Q=[γ1,...,γn]我们有
ATP=Qλ1...λr0n×m故每个m×n实对称A都能写成
A=PSQT=[β1,β2,...,βm]λ1...λr0m×n[γ1,γ2,...,γn]T=λ1β1γ1T+λ2β2γ2T+...+λrβrγrT空间r+rn+rm−>O((m+n)r)
后面的奇异值很小接近零,
在秩k的限制下矩阵A的最佳逼近:#
Ak=λ1β1γ1T+λ2β2γ2T+...+λkβkγkT→A
由
AATP=[β1,β2,...,βm]λ1...λr0m×n得
A[ATβ1,ATβ2,...,ATβr]=[λ1β1,λ2β2,...,λrβr]
A[γ1,γ2,...,γr]=[λ1β1,λ2β2,...,λrβr]