1986-Daubechies小波#
| 案件檔案 | |
|---|---|
| 案發年份 | 1986（構造完成；理論 1988–1989 發表） |
| 主嫌 | 緊支撐正交小波 \(dbN\) |
| 受害者 | Fourier 分析：全域基無法局部化，Gabor 基不穩定 |
| 關鍵證人 | Ingrid Daubechies、Yves Meyer、Stéphane Mallat（1986/1989） |
| 結案結果 | 緊支撐正交小波問世，MRA 與濾波器組理論，JPEG2000 伏筆 |
案發現場#
1980 年代中期，信號處理與調和分析的共同懸案是：沒有一個基同時擁有「好性質」。已有的候選人各有死角：
| 基 | 局部化 | 頻率局部化 | 正交性 | 快速演算法 |
|---|---|---|---|---|
| Fourier 基 \(e^{inx}\) | 無（全域） | 有 | 有 | FFT \(O(N\log N)\) |
| 短時 Fourier（Gabor） | 有 | 有 | 有（過完備） | 無穩定快速法 |
| Haar 基（1910） | 有 | 差（不連續） | 有 | 有 |
| 理想小波 | 有 | 有 | 有 | ** \(O(N)\) ** |
謎題的兩難：
- 數學上，Balian–Low 定理（1985）證明：Gabor 基若局部化太好，則不可能同時正交且穩定。
- 工程上，地震信號處理（Morlet 1984 的小波直覺）需要「多尺度、局部化、正交、快速」的基。
誰能同時滿足四個條件？這是 Ingrid Daubechies 在 1986 年要接的案子。
偵查過程#
Daubechies 的偵查策略是：從多分辨分析（MRA）的結構出發，把構造化約為解一個濾波器方程。
第一步：MRA 的五條公理#
多分辨分析（Mallat 1986、Meyer 1986）：一列閉子空間 \(\{V_j\}\subset L^2(\mathbb{R})\) 滿足：
- 嵌套： \(V_j\subset V_{j+1}\) ；
- 伸縮： \(f(x)\in V_j\iff f(2x)\in V_{j+1}\) ；
- 平移： \(f(x)\in V_0\iff f(x-k)\in V_0\) ；
- 稠密： \(\bigcup_j V_j\) 在 \(L^2\) 中稠密；
- 平凡交： \(\bigcap_j V_j=\{0\}\) 。
且存在尺度函數 \(\phi\in V_0\) 使 \(\{\phi(x-k)\}_k\) 是 \(V_0\) 的正交基。
第二步：兩尺度方程#
由 \(V_0\subset V_1\) 與伸縮， \(\phi\) 滿足兩尺度方程：
定義小波 \(\psi(x)=\sum_k g_k\,\phi(2x-k)\) （ \(g_k=(-1)^k h_{1-k}\) 型，QMF 條件），則 \(\{\psi_{j,k}(x)=2^{j/2}\psi(2^j x-k)\}_{j,k\in\mathbb{Z}}\) 是 \(L^2(\mathbb{R})\) 的正交基。
第三步：濾波器組的翻譯#
在 Fourier 域，兩尺度方程化為
其中 \(m_0\) 是低通濾波器。正交性化為 QMF 條件：
Mallat 的濾波器組理論（1989）：正交小波基 ⟺ 完全重構濾波器組（分析＋合成），信號處理與數學在此合流。
第四步：緊支撐的關鍵#
緊支撐 \(\phi,\ \psi\) ⟺ 濾波器 \(h_k\) 有限長度（FIR）。Daubechies 的高招：令 \(m_0\) 有 \(\frac12(1+e^{-i\xi})^N\) 型的「平坦度因子」：
其中 \(L\) 選得使 \(|m_0(\xi)|^2\) 在 \(\xi=\pi\) 處有 \(2N\) 階零點（頻率響應的極致平坦），並滿足 QMF。由此 \(h_k\) 有限長、 \(\phi\) 緊支撐、正交性成立。
構造偵查表：
| 條件 | 濾波器語言 | 效果 |
|---|---|---|
| 緊支撐 | \(h_k\) FIR 有限長 | 局部化 |
| 正交性 | QMF： $ | m_0 |
| 光滑性 | \(\xi=\pi\) 處零點階數 \(N\) | \(\phi\in C^{\alpha N}\) |
| 對稱性 | 無（除 Haar） | 線性相位犧牲 |
第五步： \(dbN\) 的誕生#
1986 年 Daubechies 構造出無限族緊支撐正交小波 \(dbN\) （ \(N=1,2,3,\dots\) ）：
- \(db1\) ：Haar 小波（支撐 \([0,1]\) ，不連續）。
- \(db2\) ：支撐 \([0,3]\) ，連續但不可微。
- 一般 \(dbN\) ：支撐 \([0,2N-1]\) ，光滑性 \(\alpha_N\sim 0.2N\) 。
性質隨 \(N\) 的變化：
| 小波 | 支撐長度 | 光滑性 | 消失矩 |
|---|---|---|---|
| \(db1\) （Haar） | \(1\) | 無 | \(1\) |
| \(db2\) | \(3\) | \(C^{0.\sim}\) | \(2\) |
| \(dbN\) | \(2N-1\) | \(C^{\alpha_N}\) | \(N\) |
消失矩（ \(\int x^m\psi(x)dx=0\) ， \(m<N\) ）是壓縮的關鍵：多項式信號在小波係數中只剩少量非零。
第六步：Mallat 的快速演算法#
Mallat（1986/1989）的快速小波變換：把信號逐層過濾波器組（低通＋高通＋降取樣）：
總複雜度 \(O(N)\) ——比 FFT 更快。重構由完全重構濾波器組反向進行。
結案報告#
1986–1989 年，Daubechies 與 Mallat 聯手結案：緊支撐正交小波存在（ \(dbN\) 族），構造化約為解 QMF 濾波器方程，快速演算法 \(O(N)\) 。案件遺產如下：
- 四個條件同時滿足：局部化＋頻率局部化＋正交性＋快速演算法，Balian–Low 的死角被 MRA 繞過（Gabor 用格子，小波用二進伸縮格）。
- 數學與工程的合流：Mallat 的濾波器組理論把信號處理（QMF、子帶編碼）與調和分析（MRA）統一。
- Daubechies 的榮譽：1992 年《小波十講》（Ten Lectures on Wavelets）成為經典；1994 年 MacArthur 獎、2000 年 NAS 獎。
- JPEG2000 伏筆：1990 年代末，小波（特別是 CDF 5/3、9/7 雙正交小波）成為 JPEG2000 標準的核心；消失矩與緊支撐使壓縮效率大幅超越 DCT。
- 後續浪潮：雙正交小波（Cohen–Daubechies–Feauveau 1992）、小波框（Daubechies 1990）、脊小波與曲線小波（Candès 1998）、直到 2006 年壓縮感知（見 2006-壓縮感知）中稀疏表示的角色。
證據與工具#
關鍵公式一覽表：
| 工具 | 公式／陳述 | 用途 |
|---|---|---|
| 兩尺度方程 | \(\phi(x)=\sum_k h_k\phi(2x-k)\) | MRA 的心臟 |
| QMF 條件 | $ | m_0(\xi) |
| 平坦度因子 | \(m_0=(\frac{1+e^{-i\xi}}{2})^N L(\xi)\) | 緊支撐＋光滑 |
| 小波基 | \(\psi_{j,k}=2^{j/2}\psi(2^jx-k)\) | \(L^2\) 的正交基 |
| 消失矩 | \(\int x^m\psi=0\) ， \(m<N\) | 壓縮關鍵 |
| 快速變換 | \(a_{j+1,k}=\sum_n h_{n-2k}a_{j,n}\) | \(O(N)\) |
證明骨架（MRA ⟹ 正交小波基）：
- 兩尺度方程由 \(V_0\subset V_1\) 與平移正交性取得。
- QMF 條件等價於 \(\{\phi(x-k)\}\) 的正交性。
- 定義 \(\psi\) （ \(g_k\) 為 QMF 共軛），由伸縮與稠密性證 \(\{\psi_{j,k}\}\) 是 \(L^2\) 的正交基。 \(\blacksquare\)
典型例子： \(db2\) 的濾波器： \(h_0=\frac{1+\sqrt3}{4\sqrt2},\ h_1=\frac{3+\sqrt3}{4\sqrt2},\ h_2=\frac{3-\sqrt3}{4\sqrt2},\ h_3=\frac{1-\sqrt3}{4\sqrt2}\) 。四個係數、支撐 \([0,3]\) 、二階消失矩。對多項式信號 \(p(x)=ax+b\) ，小波係數 \(d_j=0\) （消失矩吃掉低階多項式），只剩邊界附近的少量非零——這正是壓縮的數學本質，JPEG2000 的伏筆在此埋下。