跳轉到

2.3 Gibbs 現象 1899#

本章導言#

Dirichlet 定理（2.2 節）說：方波的傅立葉級數在連續點收斂到函數值。但如果你真的用電腦畫出部分和 \(S_N(x)\) 的圖，會看到一幅令人不安的景象：在跳躍點兩側，曲線衝過頭，隆起一個約 9% 的「耳廓」，而且無論你加多少項，這個耳廓不會消失——它只是變窄、向跳躍點擠壓。這就是 Gibbs 現象（Gibbs phenomenon）：傅立葉級數在間斷點附近的固定超調。本節將計算這個 9% 的精確數值、追溯它的發現史（一份手稿的遺失、一封信、一次公開更正），並說明它為何是「截斷頻譜」這一行為的內在稅收——任何用有限頻寬表示突變的系統都要付。

發現史：Wilbraham（1848）→ Gibbs（1899）#

第一位發現者其實不是 Gibbs。 1848 年，英國劍橋的Henry Wilbraham（1825–1883） 在《劍橋與都柏林數學雜誌》發表論文，精確計算了方波級數部分和在跳躍點附近的行為，畫出了帶耳廓的圖，並給出了超調值的表達式。這篇論文超前時代近半世紀，卻無人理會——當時連收斂概念都尚未定型，一個「極限不收斂到函數圖形」的現象無處安放。Wilbraham 的貢獻要到 1925 年才被 Bôcher 在其經典論文《Gibbs 現象引論》中追認。

Gibbs 的版本：Josiah Willard Gibbs（1839–1903），耶魯大學教授、美國第一位理論物理博士級的偉大學者（化學熱力學與統計力學的奠基人），在 1898 年 4 月給《Nature》寫了一篇短評，反駁Albert Michelson（邁克生，以光速實驗與 Michelson–Morley 實驗聞名）的抱怨。Michelson 當時設計了一台「諧波分析機」（69 個齒輪的機械式傅立葉分析儀），發現合成方波時端點總有突刺，認為數學（級數理論）與機器出了矛盾。Gibbs 第一封信（1898）的解釋不對；1899 年 4 月的第二封信給出了正確的極限程序：先取 \(N\to\infty\) （固定離跳躍點的距離）、再讓點靠近跳躍點——順序不同，結果不同。1906 年，Maxime Bôcher（1867–1918）（哈佛）把整個現象系統化並命名為「Gibbs 現象」，同時補上 Wilbraham 的優先權。這段歷史教訓：極限交換順序的陷阱（本節核心數學）曾經連 Michelson 與 Gibbs 都要過兩關。

9% 的計算#

考慮方波 \(f(x) = \mathrm{sgn}(\sin x)\) ，級數 \(S_N(x) = \frac{4}{\pi}\sum_{k=0}^{N}\frac{\sin(2k+1)x}{2k+1}\) （奇數項部分和）。在跳躍點 \(x=0\) 右側，超調峰值的位置與值可精確求出。

求峰值位置：對 \(S_N\) 求導：

\[ S_N'(x) = \frac{4}{\pi}\sum_{k=0}^{N}\cos(2k+1)x = \frac{4}{\pi}\cdot\frac{\sin\big(2(N+1)x\big)}{2\sin x} = \frac{2\sin\big(2(N+1)x\big)}{\pi\sin x}, \]

（用 \(\sum_{k=0}^{N}\cos(2k+1)x = \frac{\sin 2(N+1)x}{2\sin x}\) ）。故第一個駐點在 \(\sin\big(2(N+1)x\big) = 0\) 的最小正根： \(x^* = \dfrac{\pi}{2(N+1)}\) ——恰在跳躍點右側一個主瓣寬度處，且 \(x^* \to 0\) （耳廓向跳躍點擠壓）。

求峰值：計算

\[ S_N(x^*) = \frac{4}{\pi}\sum_{k=0}^{N}\frac{\sin\left(\dfrac{(2k+1)\pi}{2(N+1)}\right)}{2k+1}. \]

把它看作函數 \(\frac{\sin t}{t}\) 的 Riemann 和：令 \(t_k = \frac{(2k+1)\pi}{2(N+1)}\) ，則 \(\frac{\sin t_k}{2k+1} = \frac{\pi}{2(N+1)}\cdot\frac{\sin t_k}{t_k}\) ，步長 \(\Delta t = \frac{\pi}{N+1}\) 。因此

\[ S_N(x^*) \to \frac{4}{\pi}\cdot\frac{\pi}{2}\int_0^{\pi}\frac{\sin t}{t}\,dt = 2\,\mathrm{Si}(\pi), \]

其中 \(\mathrm{Si}(\pi) = \int_0^{\pi}\frac{\sin t}{t}dt \approx 1.85194\) 。故極限峰值為

\[ 2(1.85194) = 3.7039\ldots \times \frac{2}{\pi}\cdot\frac{\pi}{2} \;\Rightarrow\; \text{以跳幅 }2\text{ 為基準，超調比} = \frac{3.7039/2 - 1}{2}\cdot 2 = \frac{\mathrm{Si}(\pi)}{\pi} - \frac12 \approx 8.949\%. \]

換個更常用的寫法：超調約為跳躍幅度的 8.95%（方波從 \(-1\) 跳到 \(+1\) ，跳幅 2，峰值約 1.179，超出 0.179 = 8.95% × 2）。同樣有下衝（undershoot） \(-1.179\) 。這個數值與 \(N\) 無關——加項只會讓耳廓更窄、更靠近跳躍點，高度紋絲不動。 \(\blacksquare\)

為何不可消除（概念版）：部分和是有限頻寬的訊號，有限頻寬的連續函數不可能有真正的跳躍；要「最經濟地」用帶限函數貼近跳躍，必然過衝——這是帶寬與陡峭度的根本交換（8.1 節測不準原理的一個化身）。想消滅超調，只有兩條路：換核（Fejér 平均，3.3 節，犧牲銳度）或加窗（6.2 節，犧牲解析度）。

更深一層：極限順序的陷阱#

Gibbs 現象的數學本質是兩個極限不可交換：

\[ \lim_{x \to 0^+}\lim_{N\to\infty} S_N(x) = \lim_{x\to 0^+} 1 = 1 \quad(\text{Dirichlet：每個固定 } x \text{ 上收斂}), \]
\[ \lim_{N\to\infty}\lim_{x\to 0^+} S_N(x) = \lim_{N\to\infty} S_N(x^*_N) \approx 1.179 \quad(\text{先取峰再取極限}). \]

順序一換，答案差 18%。這不是詭辯，而是「部分和函數族不一致收斂」的具體病理：若 \(S_N \to f\) 一致，兩個極限必然相等（3.2 節）。方波級數在跳躍點附近正是非一致收斂的教科書案例——Gibbs 現象就是非一致收斂的照片。

一個具體例題#

例題 2.3.1：對方波級數（ \(N\) 為奇數項個數，記 \(N = 2M+1\) ），（1）用電腦或手算列表給出 \(M = 3, 7, 15\) 時峰值高度與峰位 \(x^*\) ；（2）驗證峰值趨於 1.179、峰位按 \(\frac{\pi}{2(N+1)}\) 縮小；（3）計算「耳廓寬度」（超調大於 1% 的區間）如何隨 \(N\) 縮小，並討論對示波器顯示的影響。

解：

（1）逐項計算 \(S_N(x^*)\) ， \(x^* = \frac{\pi}{2(N+1)}\) ：

\(M\) （ \(N = 2M+1\) ） \(x^*\) （弧度） \(S_N(x^*)\) 超調幅度
3（7 項） \(0.1963\) \(1.1792\) \(0.1792\)
7（15 項） \(0.0982\) \(1.1791\) \(0.1791\)
15（31 項） \(0.0507\) \(1.1790\) \(0.1790\)

（2）峰值極限 \(2\,\mathrm{Si}(\pi)/1 \cdot \frac{1}{\pi}\cdot\pi/2\cdot\frac{2}{\pi}\cdot\pi = 1.17898\) （即 \(1 + \frac{2\mathrm{Si}(\pi) - \pi}{\pi} = 1.17898\) ），表中年年命中 ✓；峰位 \(\frac{\pi}{2(N+1)}\) ： \(N=7\) 時 \(\frac{\pi}{16} = 0.1963\) ✓、 \(N=15\) 時 \(\frac{\pi}{32} = 0.0982\) ✓——峰高不動、峰位減半。

（3）超調量 \(S_N(x) - 1 > 0.01\) 的右側區間近似為 \(0 < x < \frac{c}{N}\) （常數 \(c \approx 2\pi/3\) ，由 \(\mathrm{Si}\) 的性質估出）：寬度 \(O(1/N)\) 。示波器（頻寬有限 \(B\) ）顯示一個快沿時，看到的正是 Gibbs 耳廓：上升時間 \(\approx \frac{1}{2B}\) ，過衝約 9%——工程師叫它「振鈴」（ringing）。兩百年前的數學現象，就是今日每一台示波器螢幕上的固定瑕疵。

歷史啟示#

  • 優先權的延遲歸還：Wilbraham 1848 → Gibbs 1899 → Bôcher 1906 追認。「超前時代的正確結果」若缺少概念框架，會沉睡半世紀。閱讀舊文獻與追認優先權，是數學史的義務。
  • 錯誤也可以分兩步：連 Gibbs 都在 1898 年給過錯誤解釋，一年後才修正。科學的自我糾錯靠公開辯論 + 明確的極限語言，不靠天才的一次性正確。
  • 測量儀器會「說謊」：Michelson 的機器沒有錯、級數理論也沒有錯——錯的是「儀器顯示 = 函數圖形」的默認假設。任何有限頻寬的觀測都攜帶 Gibbs 型偏見；科學儀器的校正清單裡應該有一欄叫「我的帶寬多少」。

與其他章節及領域的關聯#

  • 2.2 節：Dirichlet 定理只談收斂值；本節補上「附近長什麼樣」，兩者合成完整的跳躍點心理學。
  • 3.2 節：非一致收斂的理論根源；一致收斂成立時 Gibbs 現象消失。
  • 3.3 節：Fejér 核消滅超調（代價：邊緣變鈍）——數學給出的工程選項。
  • 4.3 節：理想低通濾波器（sinc 核）對階躍輸入的響應正是 Gibbs 振鈴。
  • 6.1–6.2 節：取樣與截頻的實務中，加窗（Hann、Hamming）是標準防超調手段。
  • 8.1 節：超調的存在 = 帶寬與時域陡峭度的測不準交換。
  • 影像處理：JPEG 銳邊附近的「蚊子噪聲」、MRI 的「Gibbs 振鈴偽影」，都是同一現象的行業名字。

練習#

  1. 用 \(\sum_{k=0}^{N}\cos(2k+1)x = \frac{\sin(2(N+1)x)}{2\sin x}\) 驗證峰值位置 \(x^* = \frac{\pi}{2(N+1)}\) ，並對 \(N = 31\) 算出 \(x^*\) 的數值。
  2. 證明 \(\int_0^{\pi}D_N(t)\,dt = 2\pi\) （ \(D_N\) 為 Dirichlet 核），並用此解釋「耳廓的質量來自 Dirichlet 核的副瓣」。
  3. 計算 \(\mathrm{Si}(\pi) = \int_0^\pi \frac{\sin t}{t}dt\) 的數值（用 Simpson 法， \(n = 8\) 個子區間），驗證約 1.8519，從而超調約 8.95%。
  4. 方波從 \(0\) 跳到 \(10\) （跳幅 10）的訊號，其 \(N\) 項傅立葉近似的峰值約多少？若示波器規格要求過衝小於 2%，截斷法可行嗎？
  5. 思考題：把方波先過「平滑器」（如移動平均）再取傅立葉級數，超調會消失。用「光滑度 \(\leftrightarrow\) 係數衰減」的價碼表（2.1 節）解釋原因，並指出代價。

進階研讀：超調的深化場景#

場景一：一階函數族的普遍超調。Gibbs 超調不是方波的特產：對任何在 \(x_0\) 有跳躍 \(\Delta = f(x_0^+) - f(x_0^-)\) 的有界變差函數，部分和在 \(x_0\) 側的極限峰高為 \(f(x_0^{\pm}) \pm 0.0895\,|\Delta|\) （Bôcher 定理）。證明思路：把 \(f\) 寫成「跳躍部分（方波型）+ 連續部分」，連續部分的貢獻在峰位附近趨於零，方波型的 8.95% 原樣保留——分解 + 主導項，又是模板式證明。

場景二：Lanczos \(\sigma\) -平滑。Cornelius Lanczos（1893–1974） 在 1950 年代提出：把 \(S_N(x)\) 的變數 \(x\) 乘上一個因子做局部平均（ \(\sigma\) 因子），等價於用 Fejér 型核重新加權係數——超調歸零，解析度損失一級。Lanczos 平滑至今在數值微分、譜方法求解 PDE 中是標準防振鈴工具（9.1 節的譜方法程式碼裡常見）。

場景三：Fejér 核的對比解剖。 \(F_N(t) = \frac{1}{N}\left(\frac{\sin(Nt/2)}{\sin(t/2)}\right)^2 \ge 0\) （非負！）， \(\int F_N = 2\pi\) ， \(\|F_N\|_{L^1} = 2\pi\) 不爆炸。非負核無超調：因為部分和變成「以 \(F_N\) 為權的局部加權平均」，平均值落在資料的最小最大值之間。Dirichlet 核有正有負（副瓣振盪），所以會過衝——超調的元兇是核的負部。這一視角統一了本節與 3.3 節。

場景四：Gibbs 現象在多維與小波世界。二維影像的截頻超調沿著邊緣呈「環狀紋」（MRI 的 Gibbs 偽影、JPEG 的振鈴）。小波方法（8.3 節）部分逃脫：因為小波基有局部支撐，跳躍的表示誤差只在跳躍附近的小區間累積，不會以全域振盪形式溢出——這是小波在邊緣表示上勝過傅立葉的核心理由之一。

場景五：測不準的下限感。把「峰位縮小 \(\frac{1}{N}\) 、峰高不變」重新敘述：部分和曲線在跳躍點附近的總振鈴能量 \(= O\!\left(\frac{1}{N}\right)\) （寬度）× 固定高度 \(= O\!\left(\frac1N\right) \to 0\) ——所以 \(L^2\) 收斂（3.1 節）依然成立！Gibbs 現象與 \(L^2\) 收斂不矛盾：超調是逐點/圖形意義的頑疾，能量意義上它逐漸無關緊要。「哪個度量在乎它」再次成為判斷標準。

補充：Michelson 諧波分析機的工程史#

Albert Michelson 與 Stratton 於 1898 年建造了一台能計算 75 個諧波的機械傅立葉分析/合成儀（齒輪 + 彈簧 + 拋光鏡面，精度驚人）。他用它做潮汐諧波分析、也嘗試合成方波，發現端點突刺無法消除，於是寫信給《Nature》質疑級數理論。這台機器是「傅立葉分析硬體化」的先驅：今日 GPU 上的 FFT、光學頻譜儀、手機的 DSP 晶片，都是同一條技術路線的後裔（6.3 節）。一場「機器 vs. 數學」的爭議，最後以「兩者都對，但極限順序要說清楚」落幕——工程與理論互相校準的典範案例。

名詞速查表（本節出現的核心術語）#

  • Gibbs 現象：傅立葉部分和在跳躍點附近固定約 8.95% 的超調與下衝。
  • 超調 / 下衝（overshoot / undershoot）：跳幅的 \(\pm 8.949\% = \dfrac{2}{\pi}\mathrm{Si}(\pi) - 1\) 。
  • 正弦積分 \(\mathrm{Si}(x) = \int_0^x \frac{\sin t}{t}dt\) ：超調值的來源。
  • 振鈴（ringing）：工程界對 Gibbs 型振盪的稱呼。
  • **峰位 \(x^* = \frac{\pi}{2(N+1)}\) **：超調峰值的位置，按 \(\frac1N\) 縮小。
  • 非一致收斂：Gibbs 現象的理論根源（極限不可交換）。
  • 非負核：Fejér 核型核函數，平均意義下無超調。
  • 加窗：以平滑窗函數壓制係數、抑制振鈴的工程手段。
  • ** \(\sigma\) -平滑（Lanczos）**：以重加權消除超調的數值技術。
  • 帶限訊號：頻譜受限於有限區間的訊號；不可能有真跳躍。

年表速查（本節涉及的關鍵年份）#

  • 1848：Wilbraham 首次精確計算超調，無人注意。
  • 1898：Michelson–Stratton 諧波分析機完成；Michelson 向《Nature》質疑。
  • 1898：Gibbs 第一封信（解釋錯誤）。
  • 1899：Gibbs 第二封信（正確的極限順序）。
  • 1906：Bôcher 系統化並命名「Gibbs 現象」，追認 Wilbraham。
  • 1900 / 1904：Fejér 定理——非負核與 Cesàro 收斂（3.3 節）。
  • 1950s：Lanczos \(\sigma\) -平滑，譜方法防振鈴標準化。
  • 1970s–：JPEG、MRI 等應用中 Gibbs 偽影成為行業常識。

深化問答三則#

問：既然加項不能消除超調，傅立葉級數是不是「壞」的表示法？ 答：不是壞，是有性格：它把跳躍表示成「無窮頻率的疊加」，因此在任何有限截斷處都要付 9% 的稅。若你的應用在乎波形形狀（顯示、控制），用 Fejér/加窗；若在乎能量（通訊、壓縮），Gibbs 超調無關緊要。工具沒有好壞，只有匹配。

問：超調 8.95% 這個數字裡有 \(\pi\) ，感覺很神秘——它從哪來？ 答：從 Dirichlet 核的形狀來。超調比 \(\frac{2}{\pi}\mathrm{Si}(\pi) - 1\) 中， \(\mathrm{Si}(\pi)\) 測量的是 \(\frac{\sin t}{t}\) 在一個半週期內的「積累盈餘」， \(\frac{1}{\pi}\) 來自核的歸一化質量。它不是巧合而是核函數形狀的積分泛函——換一個核（如 Dirichlet 核換成 Fejér 核），這個泛函就換一個值（Fejér 是 0）。神秘感的來源是不認識積分算子；認識了，9% 只是一個常數。

問：實務上如何向非數學同事解釋「為什麼不能靠加解析度消滅振鈴」？ 答：比喻：用有限支的筷子（有限頻率）夾一顆球（跳躍），筷子越多越貼合，但每支筷子都不可能貼出一個直角——貼出直角需要無窮支。加解析度（多加筷子）只會讓不合縫的區域更小、更靠近直角，但「縫的存在」不改。要「看不出縫」，靠的是距離夠遠（寬度 \(O(1/N)\) 縮小），不是縫消失。

（本節完）