JOVANA
Explore Library Glossary Getting Started Three Levels Fields How it works Mission
Join the mission
All guides

為何用共變異數,以及它為何活在彎曲空間上

從第一冊向上銜接:空間共變異數矩陣是最自然的 BCI 特徵,但它並不住在平坦的向量空間裡——把它當成向量正是黎曼幾何要修正的錯誤。

從第一冊的特徵到共變異數矩陣

在第一冊裡,你解碼動作想像的做法是:對感覺運動節律做帶通濾波、學習CSP空間濾波器、取濾波後各通道的對數變異數,再送進 LDA。這條流程可行,但它是一連串經過調校、依賴類別的選擇:選哪個頻帶、哪些濾波器、幾個、如何正則化。本路徑改以一個「幾何優先」的單一物件重建動作想像解碼,而這個物件其實默默地包含了上述一切。

這個物件就是一段試次的空間共變異數矩陣。給定一段帶通濾波後的資料片段 X \in \mathbb{R}^{n \times N}n 個通道、N 個時間取樣點),它的樣本共變異數摘要了整段試次中每一對電極如何共同變動。

C = \frac{1}{N-1} X X^{\top}, \qquad C \in \mathbb{R}^{n \times n}, \; C = C^{\top}, \; C \succ 0

一段帶通濾波試次的樣本空間共變異數——對稱,且(在取樣點足夠時)正定。

用通道之間如何共同變動——哪幾對電極一起升降——來概括一段試次,而非用它的原始波形。這份概括就是一個對稱、正定的矩陣。

X
帶通濾波後的試次(通道 × 時間取樣)。
C
概括該試次的空間共變異數矩陣。
N
時間取樣點數;用 N-1 讓估計無偏。
C \succ 0
正定——每個變異數都確實為正。

整個黎曼章節就是在這個共變異數特徵上運作;參見空間共變異數估計

陷阱:把共變異數當成向量

最直覺的做法是把 C 攤平成一個 n^2 維向量(或它 n(n+1)/2 個上三角元素),再像對待任何其他特徵一樣用歐氏距離與算術平均。距離與平均正是每個分類器所需的兩項運算——最近平均、LDA、k-NN 全都要用。於是問題來了:兩個共變異數矩陣之間,直線的歐氏距離真的是那個「對」的距離嗎?

\delta_E(C_1, C_2) = \lVert C_1 - C_2 \rVert_F = \Big( \sum_{i,j} (C_1 - C_2)_{ij}^2 \Big)^{1/2}

歐氏(Frobenius)距離把矩陣看成一袋攤平的數字——它完全無視正定這項約束。

比較兩個共變異數矩陣最天真的方法,就是逐元素相減、量剩下多少。但這把矩陣當成一袋攤平的數字,忽略了共變異數其實住在一個彎曲的曲面上。

\delta_E(C_1, C_2)
兩矩陣之間的歐氏(Frobenius)距離。
C_1, C_2
被比較的兩個共變異數矩陣。
\lVert \cdot \rVert_F
Frobenius 範數:各元素平方和的平方根。

它快速又熟悉,但如下一個區塊所示,它會扭曲平均並使變異數向上偏。

答案是否定的,有兩個具體理由。其一,對稱正定(SPD)矩陣構成一個開放的凸錐,而非完整的向量空間:把兩個相減,你可能落到錐體之外(得到一個不定矩陣,根本不是合法的共變異數)。其二,歐氏平均會出現膨脹效應(swelling)——它的行列式(可視為總變異數/訊號功率的代理量)會被抬高到超過真正的幾何平均,於是歐氏平均後的共變異數看起來會比任何產生它的試次都「更響亮」,這是人為的假象。

\det\!\Big( \tfrac{1}{2}(C_1 + C_2) \Big) \; \ge \; \sqrt{\det C_1 \, \det C_2} \; = \; \det(G)

膨脹效應:歐氏中點的行列式超過幾何平均 G 的行列式——這是表觀變異數上一種系統性的向上偏差。

用攤平的方式平均兩個共變異數矩陣會膨脹表觀的分散度:直線中點的行列式(體積)總是大於真正的幾何平均。這種內建的偏差稱為膨脹效應。

\det(\cdot)
行列式——總變異數或體積的度量。
\tfrac{1}{2}(C_1 + C_2)
兩矩陣攤平的歐氏平均。
G
幾何平均——正確、不膨脹的中點。

在兩個拉長的橢圓之間內插,除非改用幾何平均,否則歐氏中點會鼓成比任一端點都胖。

流形的觀點

解法是別再假裝空間是平的。所有 n \times n 的 SPD 矩陣構成一個光滑、彎曲的黎曼流形:在這個曲面上,距離是沿著停留在曲面上的最短曲線(測地線)來量測,而不是沿著會穿越禁區的直線弦。在這個流形上,平均會順著曲率走,膨脹效應也隨之消失——幾何平均正是最自然的中心。

共變異數矩陣是 SPD 矩陣錐面上的一個個點。距離沿曲面彎曲;某參考點處的平坦切空間,正是我們稍後要利用的線性化。

讓這番工夫值得的唯一關鍵性質,是仿射不變性——下一篇會完整推展:若所有共變異數都被同一個可逆的感測器線性變換所映射,這個自然的流形距離不會改變。重新參考 EEG、重新縮放通道增益、讓電極通過某個未知的導程場(lead field)混合——幾何都渾然不覺。而這種不變性,正是一個埋在「場次之間、受試者之間」擾動變異底下的訊號所迫切需要的,也正是黎曼方法能遷移的原因。

本路徑將建立的內容

接下來四篇要開展的收穫是:一個幾乎沒有超參數的解碼器(第三篇的 MDM)、一條在流形特徵上釋放完整機器學習工具箱的途徑(切空間技巧),以及藉由重新置中流形而幾乎免費的跨場次、跨受試者遷移(第四篇)。第五篇則帶你走到研究前沿——具幾何意識的神經網路、其他度量,以及這套典範誠實的極限。