EMD×SLEEP
Wake

黃鍔院士 · Empirical Mode Decomposition

一夜 之旅

你的腦一整晚都在說話。
傅立葉聽不懂它的口音——EMD 可以。

fs = 128 HzN = 1024IMF × 6
模擬時間 · 向下捲動 = 時間前進 23:00 尚未入睡·清醒
模擬訊號 · 依各睡眠階段的
頻譜特徵合成,非真實病人資料
23:0001:0003:0005:0007:00
01 / WAKE

清醒

閉上眼。後腦浮起 Alpha 波 (8–13 Hz),混著零碎的 Beta。這是傅立葉最擅長的時刻——訊號還算穩定。

EMD 觀察高頻 IMF1–IMF2 能量最高,Alpha 被乾淨地切進 IMF3。

02 / N1

入睡邊界

Alpha 崩解,Theta (4–8 Hz) 湧上來。這幾十秒的過渡是非穩態的教科書:頻率在滑動,而不是跳躍。

EMD 觀察瞬時頻率沿著 IMF3→IMF4 連續下滑——FFT 只會看到兩團模糊。

03 / N2

紡錘與 K 複合波

每隔幾秒,一串 12–14 Hz 的紡錘波閃過,還有一記深沉的 K 複合波。它們只持續 0.5–2 秒——是瞬態,不是「頻帶」。

EMD 觀察紡錘波被單獨隔離在一個 IMF,用 Hilbert 瞬時頻率能精確標出起訖時間 (Causa et al., 2010)。

04 / N3

慢波深海

高振幅 Delta (0.5–4 Hz) 佔據一切。這是身體修復、記憶固化的時段;慢波活動 (SWA) 是「睡眠深度」的黃金指標。

EMD 觀察EEMD 量測的 SWA 比固定頻帶的 FFT 功率譜更靈敏 (Hou et al., 2018)。

05 / REM

夢的鋸齒

腦波突然變得像清醒,卻夾著 鋸齒波與頻率調變的 Theta。眼球快速轉動,你正在做夢。

EMD 觀察各 IMF 能量分布均勻、瞬時頻率變化劇烈——正是「非線性 + 非穩態」的組合。

Why not Fourier?

傅立葉的盲點

給兩種方法同一段「頻率從 3 Hz 滑到 24 Hz」的腦波(底下疊著一條 1.2 Hz 的 delta)。傅立葉把時間全部揉在一起,回答「有哪些頻率」;Hilbert–Huang 回答「每一刻是什麼頻率」。

FFT · 功率譜時間消失了

一片從 3 到 24 Hz 的模糊隆起。它說得沒錯,只是什麼都沒說——你無法知道 20 Hz 出現在第幾秒。

HHT · 瞬時頻率時間被保留

一條清晰的斜線:頻率隨時間線性上升。這就是 EMD + Hilbert 轉換的答案。

輸入訊號 · chirp 3→24 Hz + 1.2 Hz delta

The Sifting Process

篩選:一次一層,剝出 IMF

EMD 沒有預設的基底函數。它只做一件事,反覆做:找極值、畫包絡、減去均值。直到剩下的東西「對稱」——那就是一個固有模態函數。

待命
  1. 01 辨識所有局部極大/極小值
  2. 02 三次樣條連接極大值 → 上包絡
  3. 03 三次樣條連接極小值 → 下包絡
  4. 04 取上下包絡平均 → 均值 m(t)
  5. 05 h(t) = x(t) − m(t),檢查 IMF 條件
x(t) = Σ IMFᵢ(t) + rₙ(t) · IMF 條件:|#極值 − #零交叉| ≤ 1,且局部均值 ≈ 0

Decomposition Lab

分解實驗室

調配一段你自己的睡眠腦波,按下分解,看它被剝成一層層 IMF——然後它們。每個 IMF 會被映射成一個音高(×40 進入可聽範圍),振幅包絡來自 Hilbert 轉換。

Hilbert Spectrum

Hilbert 譜:時間 × 頻率 × 能量

把每個 IMF 的瞬時頻率與瞬時振幅堆疊起來,就是 Hilbert 譜。這裡用一段「從清醒滑入深睡」的 8 秒腦波:能量的山脈從高頻一路往低頻遷移。拖曳旋轉。

時間 → 頻率 (Hz) → ↑ 振幅

IMF ↔ Brain Rhythm

每一層 IMF,都是一種腦的節律

IMF1Gamma 30–60 Hz

高頻皮質活動。REM 與清醒時出現,與認知處理相關。

IMF2Beta 13–30 Hz

活躍清醒與淺睡。隨睡眠加深而衰退,是覺醒程度的指標。

IMF3Alpha 8–13 Hz

閉眼放鬆。在入睡過渡期最突出,用來標記入睡時刻。

IMF4Theta 4–8 Hz

N1/N2 主導。紡錘波 (11–16 Hz) 常落在 IMF3–4 的交界。

IMF5Delta 0.5–4 Hz

N3 深睡主導。其振幅是睡眠深度的黃金指標。

IMF6+Sub-delta <0.5 Hz

超慢振盪,與 90 分鐘的睡眠週期 (ultradian rhythm) 相關。

Limits & Beyond

EMD 的限制,以及它怎麼被修好

EMD 是經驗方法,沒有封閉的數學式。強大,但也帶著三個大家都會踩到的坑——黃鍔團隊後續的工作幾乎都在補這些洞。

01 · MODE MIXING 模態混疊

頻率相近或間歇出現的成分(例如只在 N2 出現的紡錘波)可能被拆到同一個 IMF,或一個成分跨兩個 IMF。你在實驗室裡把雜訊調大就會看到。

EEMD加入多組白雜訊、各做一次 EMD 再取平均,雜訊互相抵消、模態被逼回正確的尺度 (Wu & Huang, 2009)。

02 · END EFFECTS 端點效應

樣條包絡在訊號兩端沒有極值可依靠,會亂飄,誤差往內傳播。篩選動畫裡兩端的包絡線常常是最不像樣的地方。

鏡射延拓把訊號兩端做鏡像延伸再分解,或以極值外推補點;睡眠分析中通常直接丟棄每段前後幾秒。

03 · STOPPING 停止準則與計算量

篩多少次才算 IMF?篩太少均值不為零,篩太多會把振幅調變抹平、失去物理意義。經典做法是 Cauchy 型 SD 準則(SD < 0.2–0.3)。

Holo-Hilbert 譜分析對 IMF 的振幅包絡再做一層 EMD,把「調幅」本身當成訊號,揭露跨頻率的耦合 (Huang et al., 2016)。

References

文獻時間線

從 1998 年的奠基論文到今日的自動偵測。橫向捲動,點擊卡片開啟論文頁面(DOI)。

開啟論文 ↗ 1998

The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis

Huang, N. E., Shen, Z., Long, S. R., Wu, M. C., Shih, H. H., Zheng, Q., Yen, N.-C., Tung, C. C., & Liu, H. H.

Proc. R. Soc. Lond. A, 454(1971), 903–995

奠基論文 doi.org/10.1098/rspa.1998.0193
開啟論文 ↗ 2009

Ensemble Empirical Mode Decomposition: A Noise-Assisted Data Analysis Method

Wu, Z., & Huang, N. E.

Advances in Adaptive Data Analysis, 1(1), 1–41

EEMD doi.org/10.1142/S1793536909000047
開啟論文 ↗ 2009

On instantaneous frequency

Huang, N. E., Wu, Z., Long, S. R., Arnold, K. C., Chen, X., & Blank, K.

Advances in Adaptive Data Analysis, 1(2), 177–229

理論基礎 doi.org/10.1142/S1793536909000096
開啟論文 ↗ 2009

Sleep stage classification based on EEG Hilbert-Huang transform

Li, Y., Yingle, F., Gu, L., & Qinye, T.

4th IEEE Conf. Industrial Electronics and Applications (ICIEA), 3676–3681

睡眠分期 · 81.7% doi.org/10.1109/ICIEA.2009.5138842
開啟論文 ↗ 2010

Automated sleep-spindle detection in healthy children polysomnograms

Causa, L., Held, C. M., Causa, J., Estévez, P. A., Perez, C. A., Chamorro, R., Garrido, M., Algarín, C., & Peirano, P.

IEEE Transactions on Biomedical Engineering, 57(9), 2135–2146

紡錘波偵測 · HHT doi.org/10.1109/TBME.2010.2052924
開啟論文 ↗ 2018

Complexity of Wake Electroencephalography Correlates With Slow Wave Activity After Sleep Onset

Hou, F., Yu, Z., Peng, C.-K., Yang, A., Wu, C., & Ma, Y.

Frontiers in Neuroscience, 12, 809

EEMD · 慢波 doi.org/10.3389/fnins.2018.00809
開啟論文 ↗ 2024

Automated detection of sleep spindles using a single EEG channel based on TEO and EMD

Li, Song, Zhang, & Karray

Expert Systems with Applications, 123661

單通道 · 自動化 doi.org/10.1016/j.eswa.2024.123661

Norden E. Huang

關於黃鍔院士

黃鍔 (Norden E. Huang) 是 Hilbert–Huang Transform 與經驗模態分解的發明者。1998 年,他在 NASA Goddard 太空飛行中心研究海洋波浪時,發現傅立葉分析對非線性、非穩態訊號的根本限制,於是提出一種不預設任何基底函數、完全由資料自身決定分解方式的方法——這就是 EMD。

這篇發表於《Proceedings of the Royal Society A》的論文,如今已是非線性時間序列分析領域被引用最多的文獻之一。方法從海洋學出發,擴散到地震學、結構健康監測、金融、氣候,以及生醫訊號——其中睡眠腦波正是它最自然的舞台,因為腦的節律從來不穩態。

返台後,他在國立中央大學創立數據分析方法研究中心,持續推動 EEMD、瞬時頻率理論與 Holo-Hilbert 譜分析等後續發展。

1998EMD / HHT 奠基論文Proc. R. Soc. Lond. A, 454, 903–995
NASAGoddard 太空飛行中心HHT 於此誕生,最初用於海洋波浪分析
NAE美國國家工程院院士亦為中央研究院院士
NCU中央大學數據分析方法研究中心創辦人,推動 EEMD 與 Holo-Hilbert 譜分析