如果你平時的工作離不開矩陣分解那你一定繞不開Eckart-Young-Mirsky定理。這個定理還有一個更常用的名字——截斷SVD最優(yōu)性定理。說白了就是任意一個矩陣你想用一個低秩矩陣去近似它只要用譜范數(shù)或Frobenius范數(shù)來衡量誤差最優(yōu)答案都是同一個——保留最大的k個奇異值其余全部置零。這句話聽起來輕飄飄的但整個數(shù)據(jù)降維、主成分分析、推薦系統(tǒng)里的矩陣分解、圖像壓縮的誤差分析全都建立在這個定理之上。這篇文章目標(biāo)很明確把E-Y-M定理的數(shù)學(xué)含義、證明思路、代碼驗(yàn)證和工程中的坑一次講透適合已經(jīng)會調(diào)用SVD、但還沒認(rèn)真想過為什么截斷SVD是最優(yōu)的讀者。1. E-Y-M定理到底在說什么1.1 一句話版本與數(shù)學(xué)表述先說標(biāo)準(zhǔn)表述。設(shè)A是一個m×n的實(shí)矩陣它的奇異值分解為A UΣV^T其中U和V是正交矩陣Σ是對角矩陣對角元滿足σ? ≥ σ? ≥ ... ≥ σ_r 0r是A的秩。截斷SVD的做法是只保留前k個奇異值以及對應(yīng)的左右奇異向量A_k U_k Σ_k V_k^T其中U_k是U的前k列V_k是V的前k列Σ_k是Σ前k個奇異值組成的k×k對角矩陣。因?yàn)楹竺鎚-k或n-k個奇異值都被丟棄了A_k的秩最多為k但它仍然是一個m×n矩陣形狀和A完全一樣。E-Y-M定理說的是在所有秩不超過k的m×n矩陣B里A_k同時是譜范數(shù)和Frobenius范數(shù)意義下的最優(yōu)逼近。換句話說||A - A_k||? ≤ ||A - B||?||A - A_k||_F ≤ ||A - B||_F對任意秩≤k的B都成立。而且誤差有精確表達(dá)式||A - A_k||? σ_{k1}第k1個奇異值||A - A_k||F sqrt(σ{k1}2 σ_{k2}2 ... σ_r2)注意k可以從0取到r。k0時A_0就是零矩陣譜范數(shù)誤差就是σ?k≥r時誤差為0因?yàn)榈椭染仃嚤旧硪呀?jīng)能完全表示A了。這里經(jīng)常有人混淆譜范數(shù)和Frobenius范數(shù)的定位。譜范數(shù)衡量的是最壞情況下的放大倍數(shù)它只關(guān)心最大奇異值Frobenius范數(shù)則是把所有元素的平方差加起來再開方像歐氏距離的矩陣版本衡量的是整體平均能量。E-Y-M定理巧妙的地方在于這兩個性質(zhì)不同、用途不同的范數(shù)最優(yōu)解居然都是同一個截斷SVD。1.2 直覺為什么截斷SVD是最優(yōu)的很多新手會有一個直覺想用一個低秩矩陣近似A直接從A里抽前k列不行嗎實(shí)際試一下就知道了這樣做的誤差通常大得離譜。原因是矩陣的列與列之間高度相關(guān)只拿前k列意味著完全放棄了后面所有列里的信息而后面列未必是不重要的列。SVD之所以比抽列抽行強(qiáng)是因?yàn)樗茸隽艘淮巫鴺?biāo)旋轉(zhuǎn)。V^T把A的行空間旋轉(zhuǎn)到了主軸坐標(biāo)系在這個坐標(biāo)系里每個方向?qū)仃嚹芰康呢暙I(xiàn)正好是對應(yīng)的奇異值平方。奇異值越大這個方向上的能量越集中。截斷SVD相當(dāng)于保留能量最大的k個主軸丟掉其余能量較小的主軸。打個比方矩陣像一張建筑圖紙奇異值就是不同比例尺下的輪廓強(qiáng)度。前幾個大奇異值決定了整體骨架和主要結(jié)構(gòu)后面一堆小奇異值只是門窗細(xì)節(jié)和裝飾線。做低秩逼近就像是先畫大框架再摳細(xì)節(jié)——如果只允許你保留k筆最優(yōu)策略肯定是先畫影響最大的幾筆而不是隨便截掉圖紙的一半。E-Y-M定理就是把這句直覺話變成了嚴(yán)格的數(shù)學(xué)保證。2. 為什么這個定理值得你花時間2.1 PCA是它的直接推論主成分分析PCA是E-Y-M定理最經(jīng)典的應(yīng)用。PCA做的事情是給一個數(shù)據(jù)矩陣Xn個樣本p個特征先對每一列做中心化也就是減去該列的均值然后找到一組新的正交方向讓數(shù)據(jù)投影到這組方向上后方差最大或者等價地投影后重構(gòu)原始數(shù)據(jù)的誤差最小。如果把中心化后的數(shù)據(jù)矩陣記為X_c對它做SVDX_c UΣV^T那么前k個右奇異向量V的前k列就是主成分方向。E-Y-M定理的Frobenius范數(shù)部分保證了在所有的k維子空間里由前k個右奇異向量張成的子空間重建X_c時產(chǎn)生的F范數(shù)誤差最小。這里有一個非常容易被忽略的前提PCA必須先做中心化。E-Y-M定理本身不要求中心化直接對原始X做截斷SVD也能得到一個最優(yōu)低秩近似。但PCA要的是最大方差方向如果數(shù)據(jù)中心不在原點(diǎn)第一主成分會被均值方向帶偏而不是真正的最大方差方向。所以工程上的順序永遠(yuǎn)是先中心化再SVD再取前k個方向。把PCA理解成E-Y-M定理的特例有個額外好處以后看到其他降維方法你可以問一句它的最優(yōu)性來自哪個范數(shù)下的哪個約束思路會清楚很多。2.2 推薦系統(tǒng)與圖像壓縮都在用它推薦系統(tǒng)里的矩陣分解目標(biāo)是把用戶對物品的評分矩陣R分解成兩個低秩矩陣的乘積R ≈ P Q^TP是用戶隱因子矩陣Q是物品隱因子矩陣。如果評分矩陣是完整無缺失的那最優(yōu)的低秩近似就是截斷SVDE-Y-M定理直接給出最優(yōu)解和誤差上界。實(shí)際推薦場景中評分矩陣有大量缺失值不能直接套截斷SVD但低秩假設(shè)本身——即用戶偏好可以被少數(shù)幾個隱因子解釋——正是從E-Y-M定理背后衍生出來的核心思想。圖像壓縮是另一個直觀例子。一張灰度圖就是一個m×n像素矩陣。對它做SVD后只保留前k個奇異值就可以用U_k、Σ_k、V_k^T三塊小矩陣重建圖像。存儲成本從mn個數(shù)值降到k(mn1)個數(shù)值。對人像、風(fēng)景這類結(jié)構(gòu)化的圖像k取30到50時肉眼幾乎看不出和原圖的差別。傳統(tǒng)JPEG壓縮在細(xì)節(jié)保留上當(dāng)然更強(qiáng)但SVD壓縮的教育意義無可替代因?yàn)樗研盘柲芰考性谇皫讉€奇異值這一事實(shí)擺在臺面上。特征臉、潛在語義分析、文本主題模型本質(zhì)也都是同一套思路原始數(shù)據(jù)有冗余低秩結(jié)構(gòu)能解釋大部分方差截斷SVD給出最優(yōu)低秩表示。2.3 低秩假設(shè)的理論基石再往深一層說E-Y-M定理支撐起了整個低秩模型體系。數(shù)據(jù)是低秩加噪聲的低秩部分才是信號噪聲散布在所有方向上。把奇異值中最小的那些置零就是在做最有效的去噪。這些操作看起來像啟發(fā)式技巧背后卻有一個精確的最優(yōu)性定理兜底所以工程上敢放心用。如果有一天你在論文里看到通過低秩近似實(shí)現(xiàn)數(shù)據(jù)壓縮用截斷SVD作為強(qiáng)基線或者低秩矩陣恢復(fù)追根溯源理論起點(diǎn)基本都是E-Y-M定理。理解它等于拿到了一把通用鑰匙。3. 證明思路拆解3.1 譜范數(shù)情況的經(jīng)典論證譜范數(shù)部分的證明不復(fù)雜而且非常漂亮。核心思路分兩步先證明任意秩≤k的B誤差都≥σ_{k1}再證明A_k恰好能達(dá)到這個值。先看下界。B的秩不超過k那么B的零空間null(B)的維度至少是n-k。選擇這個零空間作為候選子空間S那么對任意x∈S且||x||1有Bx 0所以||(A - B)x|| ||Ax||于是||A - B||? max_{||x||1} ||(A-B)x|| ≥ max_{x∈S, ||x||1} ||Ax||接下來需要用到奇異值的一個經(jīng)典變分刻畫第k1個奇異值σ_{k1}等于所有n-k維子空間上最大拉伸的最小值也就是σ_{k1}(A) min_{dim S n-k} max_{x∈S, ||x||1} ||Ax||B的零空間本身就是一個n-k維子空間所以上式的右邊至少是σ_{k1}。連起來就是||A - B||? ≥ σ_{k1}再看構(gòu)造。A_k的秩為k并且A - A_k的奇異值恰好是σ_{k1}, σ_{k2}, ...前k個被消掉了。因此A - A_k的最大奇異值就是σ_{k1}即||A - A_k||? σ_{k1}誤差下界和可達(dá)性同時成立譜范數(shù)部分就證完了。3.2 F范數(shù)誤差為什么是奇異值平方和Frobenius范數(shù)的證明比譜范數(shù)稍微繞一點(diǎn)這里給你一個足夠建立信任的直觀推導(dǎo)。F范數(shù)關(guān)于左右正交變換不變所以可以先把A放進(jìn)SVD坐標(biāo)系里考慮||A - B||_F ||Σ - U^T B V||_F記C U^T B V它的秩仍然不超過k?，F(xiàn)在問題變成在Σ這個對角矩陣旁邊放一個秩不超過k的矩陣C最小化它們差的F范數(shù)。因?yàn)镕范數(shù)的平方等于所有元素差的平方和C的秩最多為k它再厲害也只能同時控制k個方向上接近Σ的對應(yīng)項(xiàng)。第k1個及之后的奇異值方向C并沒有足夠的自由度去覆蓋這些方向上的能量Σ_{k1}, Σ_{k2}, ...就必然作為誤差保留下來。極端理想情況下C精確匹配前k個對角線上的值其他位置全是0這時候誤差就是sqrt(σ_{k1}2 σ_{k2}2 ... σ_r2)而C U^T A_k V恰好能做到這一點(diǎn)所以等號能取到。這個推導(dǎo)雖然偏直覺但和嚴(yán)格證明的結(jié)論完全一致。想要更嚴(yán)謹(jǐn)?shù)臄⑹隹梢匀シ仃嚪治鼋滩睦锏钠娈愔禂_動不等式和交錯定理。對工程應(yīng)用來說記住F范數(shù)誤差平方 被丟棄奇異值的平方和這個結(jié)論比記住繁瑣的證明過程更有用。3.3 這個證明教會我們什么從證明里能讀出三層信息。第一最優(yōu)低秩逼近并不是簡單地把矩陣切成小塊而是要先做正交旋轉(zhuǎn)再截斷。旋轉(zhuǎn)本身不消耗存儲預(yù)算只改變坐標(biāo)表達(dá)。第二第k1個奇異值像一道分水嶺它決定了譜范數(shù)誤差的上限也暗示了一個事實(shí)當(dāng)奇異值快速衰減時低秩逼近效果會好到驚人。第三能取到等號的A_k不是唯一的如果奇異值有重復(fù)對應(yīng)的奇異向量子空間可以旋轉(zhuǎn)但誤差不變。工程上如果發(fā)現(xiàn)多次運(yùn)行結(jié)果的主成分方向不一樣多半是遇到了奇異值簡并這不一定是bug。4. 實(shí)操寫代碼驗(yàn)證并落地4.1 生成含噪聲的低秩矩陣并驗(yàn)證理論說得再好不如跑一段代碼親眼驗(yàn)證。這里用NumPy生成一個真實(shí)的低秩矩陣再加上噪聲然后比較截斷SVD的實(shí)測誤差和定理預(yù)測值。import numpy as np rng np.random.default_rng(42) m, n, true_rank, k 100, 80, 25, 6 # 生成一個低秩矩陣兩個隨機(jī)正交因子 對角奇異值 U0, _ np.linalg.qr(rng.standard_normal((m, true_rank))) V0, _ np.linalg.qr(rng.standard_normal((n, true_rank))) S0 np.logspace(2, 0, true_rank) # 奇異值從100衰減到1 # 真實(shí)信號 高斯噪聲 A U0 np.diag(S0) V0.T 0.1 * rng.standard_normal((m, n)) # 完整SVD U, s, Vt np.linalg.svd(A, full_matricesFalse) # 截斷SVD重建 A_k (U[:, :k] * s[:k]) Vt[:k, :] # 實(shí)測誤差 err_2 np.linalg.norm(A - A_k, 2) err_f np.linalg.norm(A - A_k, fro) # 定理預(yù)測 pred_2 s[k] # 第k1個奇異值索引從0開始 pred_f np.sqrt(np.sum(s[k:] ** 2)) print(f譜范數(shù) 實(shí)測 {err_2:.6f} 預(yù)測 {pred_2:.6f}) print(fF范數(shù) 實(shí)測 {err_f:.6f} 預(yù)測 {pred_f:.6f})跑出來的兩組數(shù)會在浮點(diǎn)精度范圍內(nèi)完全一致。這里有兩個細(xì)節(jié)值得說。一是np.linalg.svd默認(rèn)full_matricesTrue對非方陣會返回兩個大正方陣白白多算很多內(nèi)存。設(shè)成full_matricesFalse后U形狀是m×min(m,n)Vt形狀是min(m,n)×n足夠完成截斷重建。二是s[k]對應(yīng)的是丟棄的第一個奇異值因?yàn)镻ython索引從0開始保留k個意味著從下標(biāo)k開始的奇異值都要丟棄。很多人在這一步對不上數(shù)字基本都是索引沒對齊。4.2 選k的三種實(shí)用方法理論給了誤差公式但k怎么選仍然是個工程問題。我自己常用的方法有三個。第一種是能量占比。奇異值平方表示該方向貢獻(xiàn)的能量定義累計能量占比energy np.cumsum(s ** 2) / np.sum(s ** 2) k_energy np.searchsorted(energy, 0.9) 1 # 找到第一個超過90%的位置想要保留90%的總能量就找累計能量曲線第一次超過0.9的位置。這個方法無監(jiān)督、可解釋、速度快適合探索數(shù)據(jù)時快速估算k的取值范圍。第二種是奇異值缺口。把奇異值從大到小畫成折線圖如果某個位置出現(xiàn)明顯拐點(diǎn)比如σ_k是10σ_{k1}突然掉到1那說明前k個奇異值承載了主要結(jié)構(gòu)后面基本是噪聲。拐點(diǎn)處就是合理的k。這個方法需要看圖帶有主觀成分但往往非常有效。第三種是有監(jiān)督的交叉驗(yàn)證。把矩陣的一部分元素隨機(jī)遮住在剩余元素上做截斷SVD然后預(yù)測被遮住的元素統(tǒng)計預(yù)測誤差。k從1掃到r選驗(yàn)證誤差最小的那一個。這個方法最穩(wěn)健代價是計算量上去了。實(shí)際項(xiàng)目中我會先用能量占比掃一遍范圍再用交叉驗(yàn)證精挑。4.3 大規(guī)模場景下的截斷SVD選擇如果m和n都上千甚至更大直接np.linalg.svd求完整分解會非常慢因?yàn)橥暾鸖VD的復(fù)雜度是O(mn·min(m,n))級別。這時候應(yīng)該用專門的截斷算法。scipy.sparse.linalg.svds支持稀疏矩陣和指定k值只算最大的k個奇異值對。sklearn.decomposition.TruncatedSVD封裝了randomized SVD對大矩陣更友好適合作為降維流水線的一環(huán)。這里有個容易踩的坑svds返回的奇異值默認(rèn)是升序排列的也就是從小到大的順序。如果你沿用np.linalg.svd的慣性直接取s[k:]當(dāng)成丟棄的奇異值結(jié)論會完全反掉。解決方法是拿到結(jié)果后統(tǒng)一用np.sort(svds_s)[::-1]轉(zhuǎn)成降序。當(dāng)k比min(m,n)小很多時randomized SVD的思路是先用一個隨機(jī)矩陣把A投影到一個維度為kp的小空間再對小矩陣做精確SVD。它的理論保證可以理解為先構(gòu)造一個接近最優(yōu)的低秩子空間再用E-Y-M定理保證誤差。實(shí)際中p取5到15就夠速度能比完整SVD快一個數(shù)量級以上。5. 工程上最容易踩的坑5.1 數(shù)據(jù)有沒有中心化和標(biāo)準(zhǔn)化E-Y-M定理本身對輸入矩陣沒有任何均值要求但一旦你把截斷SVD用在PCA語境里就必須先中心化。很多人直接對原始數(shù)據(jù)矩陣做SVD取前k個奇異向量當(dāng)成主成分結(jié)果第一主成分總是和數(shù)據(jù)均值方向攪在一起方差解釋率看起來很高實(shí)際卻跑偏了。標(biāo)準(zhǔn)化是另一回事。比如數(shù)據(jù)里一列是金額取值范圍0到10000另一列是0到1的評分奇異性會被金額列完全支配。這種情況下要先做標(biāo)準(zhǔn)化讓每列同量綱再做SVD。標(biāo)準(zhǔn)化不等于中心化兩者目的不同中心化讓方差方向有意義標(biāo)準(zhǔn)化讓不同特征平等競爭。在實(shí)際項(xiàng)目里我習(xí)慣先看業(yè)務(wù)含義再決定做不做標(biāo)準(zhǔn)化而不是無腦套Z-score。5.2 缺失值場景不能用樸素截斷SVD推薦系統(tǒng)評分矩陣幾乎都有缺失值如果你直接把缺失位置填成0然后對全矩陣做截斷SVDE-Y-M定理依然成立但代價是它會把缺失當(dāng)成真實(shí)評分是0來處理。結(jié)果就是預(yù)測值普遍偏向低分誤差毫無意義。正確的做法是加權(quán)低秩近似只統(tǒng)計有觀測值的位置把缺失位置的誤差權(quán)重設(shè)為0。這不再是標(biāo)準(zhǔn)SVD能解的封閉問題需要用交替最小二乘或隨機(jī)梯度下降迭代求解。正因?yàn)槿绱斯こ讨忻鎸ο∈柙u分矩陣時不會有人真的去調(diào)np.linalg.svd而是會用矩陣補(bǔ)全類方法。低秩假設(shè)還是那個低秩假設(shè)但求解路徑完全不同。5.3 數(shù)值細(xì)節(jié)順序、閾值、shape浮點(diǎn)誤差會讓很小的奇異值變成負(fù)值或者虛數(shù)實(shí)際上NumPy會返回接近0的正數(shù)但有時會帶噪聲在計算能量占比時如果矩陣不是滿秩s[k:]里混入幾個1e-15量級的數(shù)會導(dǎo)致開方后誤差極小但非零這不算bug但會讓自動化測試的斷言失敗。我自己習(xí)慣在截斷前加一個閾值奇異值小于max(m,n) * eps * s[0]的直接視為0。另一個常見問題是特征值分解和SVD混用。對于對稱正定矩陣特征值分解等于SVD對一般矩陣特征分解可能出現(xiàn)負(fù)特征值、復(fù)數(shù)特征值甚至沒有實(shí)特征分解而SVD永遠(yuǎn)存在。低秩逼近一律用SVD不要圖省事去調(diào)np.linalg.eig然后取絕對值。5.4 譜范數(shù)和F范數(shù)怎么選這兩個范數(shù)對應(yīng)不同的工程訴求。譜范數(shù)只關(guān)注最壞情況下誤差有多大適合安全分析、穩(wěn)定性證明或者當(dāng)你對最大異常非常敏感時。F范數(shù)關(guān)注整體誤差能量適合圖像重建、PCA重構(gòu)、推薦系統(tǒng)評估這類看平均效果的任務(wù)。舉個例子圖像壓縮中如果某一塊區(qū)域因?yàn)槠娈愔到財喽霈F(xiàn)劇烈偽影這個偽影在F范數(shù)里占比可能不大卻在譜范數(shù)里非常刺眼。反過來推薦系統(tǒng)評價指標(biāo)通常用均方根誤差它本質(zhì)上和F范數(shù)同源所以F范數(shù)更貼近業(yè)務(wù)指標(biāo)。選哪種取決于你最在乎的是平均好還是最差不離譜。6. 常見問題速查表問題原因解決方案截斷SVD誤差比σ_{k1}小k的索引理解錯誤確認(rèn)保留的是前k個奇異值σ_{k1}對應(yīng)s[k]結(jié)果和理論對不上矩陣非方陣時用了full_matricesTrue重建時形狀對不齊用full_matricesFalse統(tǒng)一U、s、Vt的形狀svds返回的奇異值順序不對SciPy的截斷SVD返回升序用np.sort(s)[::-1]手動降序PCA第一主成分方向不合理漏做中心化每一列減去均值后再SVD列間量綱差異過大大數(shù)值列主導(dǎo)奇異值先做標(biāo)準(zhǔn)化再SVD稀疏矩陣直接截斷SVD預(yù)測值偏低缺失位置被當(dāng)作0參與誤差計算改用矩陣補(bǔ)全或加權(quán)低秩近似小奇異值抖動導(dǎo)致能量占比不穩(wěn)定數(shù)值精度或秩退化設(shè)閾值奇異值小于max(m,n)*eps*s[0]視為0主成分方向每次跑不一樣奇異值簡并這是正?，F(xiàn)象檢查是否有重復(fù)奇異值官方說法是子空間等價這張表基本覆蓋了我這幾年在項(xiàng)目里見過的絕大多數(shù)問題。排查順序建議是先查索引和shape再查是否中心化、標(biāo)準(zhǔn)化最后查奇異值順序。按照這個順序走大多數(shù)神秘bug都能在五秒內(nèi)定位。另外還有一個值得養(yǎng)成的習(xí)慣在寫任何使用SVD的代碼前先手動構(gòu)造一個已知低秩矩陣比如用兩個矩陣相乘生成一個秩為r的矩陣再疊加噪聲然后跑一遍完整流程。這樣既能驗(yàn)證代碼正確性也能讓你對低秩噪聲的數(shù)據(jù)形態(tài)形成直觀手感以后再遇到真實(shí)數(shù)據(jù)時一眼就能判斷低秩假設(shè)是否成立。我個人在實(shí)際操作中最深的體會是E-Y-M定理看起來只是一個小結(jié)論但它把PCA、推薦系統(tǒng)、圖像壓縮、去噪這些看似毫不相干的領(lǐng)域串在了同一條線上。每當(dāng)你遇到一個新問題只要里面出現(xiàn)了低秩近似四個字先問自己三個問題數(shù)據(jù)有沒有缺失要不要中心化和標(biāo)準(zhǔn)化誤差該用譜范數(shù)還是F范數(shù)這三個問題想清楚再用截斷SVD方向基本不會錯。最后一個私房技巧驗(yàn)證定理時一定要用SVD而不是特征值分解因?yàn)镾VD對任何矩陣都存在而且數(shù)值上更穩(wěn)定這一條能幫你避開大量不明不白的報錯。