5.3 Bessel 方程#

本章導言#

Legendre 統治球面，Bessel 函數統治圓柱：圓膜振動、熱流散熱、光學衍射、天線輻射、貝塞爾濾波器——凡有圓的地方就有它。方程以其系統研究者 Friedrich Wilhelm Bessel（貝塞爾，1784–1846） 命名（1824 年行星擾動論文），但首度登場在 1732 年 Daniel Bernoulli 研究吊鏈振動、以及 1764 年 Euler 分析圓膜振動時——函數比名字早出生九十年。本節以 Frobenius 法（5.1 節）完整解出此方程，展示正則奇點處的三種劇本，並介紹 \(J_\nu\) 的振盪、衰減與正交性。

從圓膜振動出發#

圓形薄膜（半徑 \(R\) ）小振動滿足波方程（9.2 節） \(u_{tt} = c^2 \nabla^2 u\) 。分離變數 \(u = T(t)\, \Phi(\phi)\, Z(r)\) ：角度部分給 \(e^{im\phi}\) ，時間部分給 \(T'' + \omega^2 T = 0\) ，徑向部分化為

令 \(x = kr\) ，即得 Bessel 方程（ \(\nu = m\) ）。「柱面對稱 + 振盪 = Bessel」——正如「球對稱 + 振盪 = Legendre」（5.2 節）。

Frobenius 法完整演出#

\(x = 0\) 是正則奇點（3.3 節： \(x p = x \cdot \frac1x = 1\) 、 \(x^2 q = x^2 - \nu^2\) 皆解析）。設

代入方程（各項對齊 \(x^{r + n}\) ），最低次（ \(n = 0\) ）給指標方程：

其後遞推： \(a_n = -\frac{a_{n-2}}{(n + r)^2 - \nu^2}\) （ \(n \ge 2\) ），奇係數全零。取 \(r = \nu\) （ \(\nu \ge 0\) ），標準化 \(a_0 = \frac{1}{2^\nu \Gamma(\nu + 1)}\) ：

** \(\Gamma\) 函數（Euler 1729 年的階乘推廣： \(\Gamma(n+1) = n!\) ，全域定義）使 \(\nu\) 可以是任何實數——Bessel 函數天生就是連續階**，不似 Legendre 依賴整數截斷。

第二解的三種劇本（5.1 節場景一）：

Neumann 函數 \(Y_\nu\) （以 Karl Neumann 1867 命名，與 Heaviside/Weber 的先驅工作並存）統一第二解： \(x \to 0^+\) 時 \(Y_\nu \to -\infty\) （對數發散）。物理剔除：圓心處振幅有限 ⇒ 拋棄 \(Y_\nu\) （如同 5.2 節拋棄 \(Q_n\) ）——除非圓心有源（波導的芯線、渦流的核），那時 \(Y_\nu\) 重返舞台。

\(J_\nu\) 的長相與漸近行為#

• **小 \(x\) **： \(J_\nu(x) \sim \frac{1}{\Gamma(\nu+1)} \big(\frac{x}{2}\big)^\nu\) ——按 \(\nu\) 次冪貼地起飛（ \(J_0(0) = 1\) 、 \(J_1(x) \approx x/2\) ）。
• **大 \(x\) **：

振盪但衰減——與 \(\sin\) （不振盪衰減）和 \(e^{-x}\) （不振盪）都不同：柱面波的能量隨 \(1/r\) 稀釋，故振幅 \(\sim r^{-1/2}\) （能量 \(\sim\) 振幅 \(^2 \sim 1/r\) ）。漸近式是「柱面波」的身份證（漸近級數的精確形式見 5.1 節場景三；Hankel 函數 \(H_\nu^{(1,2)} = J_\nu \pm i Y_\nu\) 是「出射/入射柱面波」）。

零點： \(J_0\) 的零點 \(0, 2.405, 5.520, \ldots\) 決定圓膜振動模式：鼓面的同心圓節線。鼓手敲擊不同位置激發不同 \(J_m(k r) \cos m\phi\) 模式——鼓的音色是一張 Bessel 零點表（Rayleigh《聲學理論》1877）。

正交性與柱面調和分析#

正交關係： \(\displaystyle\int_0^R r\, J_\nu\big(\tfrac{\alpha_m r}{R}\big) J_\nu\big(\tfrac{\alpha_n r}{R}\big)\, dr = 0 \quad (m \neq n)\) ，其中 \(\alpha_m\) 是 \(J_\nu\) 的零點。

（證明與 5.2 節同型：兩個不同 \(k\) 的方程相減、分部積分、邊界由零點條件消去；額外的 \(r\) 權重來自極座標面積元 \(r\, dr\, d\phi\) 。）

由此建立Fourier–Bessel 級數： \(f(r) = \sum c_m J_\nu(\alpha_m r / R)\) ——圓盤上的「傅立葉分析」，是散熱片、圓板振動、光纖模式計算的標準工具。

積分表達式（Bessel 1824、Schläfli 1868）：

—— \(J_n\) 是「頻率調製」的信號成分：調頻廣播的邊帶功率分配就是 \(J_n\) 值的問題（FM 理論，Carson 1922）。一個十九世紀的擾動積分，成為二十世紀通訊工程的核心表。

一個具體例題#

例題 5.3.1：驗證 \(y = J_0(x)\) 滿足 \(\nu = 0\) 的 Bessel 方程（用級數前幾項直接驗證），並求 \(J_0\) 的第一個正零點的近似位置。

解： \(J_0(x) = 1 - \frac{x^2}{4} + \frac{x^4}{64} - \frac{x^6}{2304} + \cdots\) 。

計算導數： \(J_0' = -\frac{x}{2} + \frac{x^3}{16} - \frac{6x^5}{2304}\cdot\frac{?}{}\) （精確為 \(-\frac{x}{2} + \frac{x^3}{16} - \frac{x^5}{384} + \cdots\) ）， \(J_0'' = -\frac12 + \frac{3x^2}{16} - \frac{5x^4}{384} + \cdots\) 。

代入 \(x^2 J_0'' + x J_0' + x^2 J_0\) ，按冪次收集：

• \(x^0\) ： \(0 \cdot J_0'' + 0 + 0 = 0\) ✓（ \(x^2 J''\) 無常數項）
• \(x^2\) ： \(-\frac12 x^2 - \frac{x^2}{2} + x^2 = 0\) ✓
• \(x^4\) ： \(\frac{3}{16} x^4 + \frac{x^4}{16} - \frac{x^4}{4} = \frac{4}{16} x^4 - \frac{4}{16} x^4 = 0\) ✓

（一般項可由遞推保證。）故 \(J_0\) 是解。∎

零點估計：用漸近式 \(\cos(x - \pi/4)\) 的零點： \(x - \frac{\pi}{4} = \frac{\pi}{2} \Rightarrow x \approx \frac{3\pi}{4} = 2.356\) 。真值 \(\alpha_1 = 2.4048\) ——漸近式在第一個零點已給出 \(2\%\) 內誤差，第二零點誤差降到 \(0.4\%\) 。漸近公式「越遠越準」（與 Taylor 相反）——這是漸近分析的鮮明性格（5.1 節場景三）。

應用快照： \(\alpha_1 / \pi \approx 0.7655\) ——半徑 \(R\) 的圓膜基音頻率 \(f = \frac{2.405}{2\pi R}\sqrt{T/\rho}\) 。同尺寸方膜的基頻係數為 \(\sqrt{2}\pi / (2\pi) = \frac{\sqrt 2}{2} \approx 0.7071\) 。圓鼓比方鼓「高亢」約 8%——Bessel 零點與三角函數零點的競賽，就是鼓面形狀與音高的關係。

歷史啟示#

• 函數先於理論：1732 年 Bernoulli 的吊鏈、1764 年 Euler 的圓膜、1824 年 Bessel 的行星擾動，三個不相干問題逼出同一組級數；「特殊函數」作為學科的誕生，是問題驅動而非理論驅動的。
• 漸近與級數的分工：Taylor 級數管局部（ \(x\) 小），漸近展開管遠方（ \(x\) 大）；Bessel 函數的完整理解需要兩者。沒有一種展開能獨裁。
• 物理剔除原則的普適： \(Y_\nu\) 在原點發散被剔除，與 \(Q_n\) （5.2）在兩極發散被剔除如出一轍——「域的幾何決定解的許可證」。

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

• 5.1 節：Frobenius 法的招牌案例；指標方程三劇本全在此上演。
• 5.2 節：球面 vs 柱面的姊妹篇；正交性證明的同型結構。
• 9.1、9.2 節：圓柱內的熱傳導與波動問題，分離變數後徑向方程就是 Bessel。
• 4.1 節：正交性與 Sturm–Liouville 理論（Bessel 方程乘 \(\frac1x\) 後為自伴形式）。
• 光學：圓孔衍射（Airy 斑）的強度分佈 \(\big(\frac{2J_1(x)}{x}\big)^2\) ——望遠鏡解析度（Rayleigh 判據 1879）由 \(J_1\) 第一零點決定；天文學的「看得多清」，是一個 Bessel 零點。
• 訊號處理：FM 調變、Bessel 濾波器（最大平坦群延遲）。
• 數論與機率： \(e^{\frac{x}{2}(t - 1/t)}\) 生成函數連接整數分拆；隨機行走位移分佈的漸近是 \(J_0\) 。

練習#

1. 用遞推寫出 \(J_1(x)\) 的級數前三項，並驗證其滿足 \(\nu = 1\) 的方程（至 \(x^3\) 項）。
2. 證明 \(J_0'(x) = -J_1(x)\) （對級數逐項求導、指標平移比對）。
3. 由遞推關係 \(J_{\nu+1} = \frac{2\nu}{x} J_\nu - J_{\nu-1}\) 與 \(J_0' = -J_1\) ，導出 \(J_0'' + \frac{1}{x} J_0' + J_0 = 0\) 。
4. 對 \(\nu = 0\) 的重根情形，說明為何第二解必含 \(\ln x\) （用 4.1 節 Abel 恆等式： \(W = \frac{C}{x}\) ，若第二解也是純級數則 \(W \sim \frac{1}{x^k}\) 型，比對失敗）。
5. 思考題：圓膜的徑向模式 \(J_0(\alpha_m r / R)\) 有 \(m - 1\) 個同心圓節線。請解釋：為什麼節線半徑不是等距的（用 \(J_0\) 零點序列 \(2.405, 5.520, 8.654, \ldots\) 的間隔討論），並連結漸近式說明為什麼高階節線趨於等距。

進階研讀：Bessel 理論的五個深化場景#

場景一：修正 Bessel 函數。把方程中 \(x^2\) 換成 \(-x^2\) （ \(x^2 y'' + x y' - (x^2 + \nu^2) y = 0\) ），振盪變成指數型：解為 \(I_\nu(x) = i^{-\nu} J_\nu(ix)\) 、 \(K_\nu(x)\) 。應用：圓柱內的穩態熱傳、波的隧道效應（evanescent 波）、金融的 bond 定價積分。把振盪參數變號，函數族立即改行——特殊函數的「家族遺傳」現象。

場景二：WKB 與 Airy 連接。 \(J_\nu\) 在「轉折點」 \(x \approx \nu\) 兩側行為迥異（內側指數衰減、外側振盪衰減），衔接兩域的正是 Airy 函數（5.1 節）的均匀漸近（Langer 1934 的修正 \(x \to \sqrt{x^2 - \nu^2}\) ）。三個函數族在一個轉折點上交接班——這是數學物理中「uniform asymptotics」大業的樣板（Olver 1974 標準化）。

場景三：Bessel 的天文學現場。F. W. Bessel 是史上第一個測得恆星視差的人（1838 年，天鵝座 61，0.3 角秒）——證明地球繞日、確立宇宙距離尺度的人，同時把擾動函數展成 \(\frac{1}{\pi}\int \cos(n\tau - x\sin\tau)\, d\tau\) 。他的橢圓軌道擾動級數（今稱 Bessel 函數）計算的正是「開普勒方程的傅立葉係數」。觀測大師與函數發明者是同一人——十九世紀天文台是應用數學的引擎室。

場景四：隨機矩陣與黎曼 ζ。GUE 隨機矩陣特徵值的間距分佈、黎曼 ζ 零點的統計（Montgomery–Dyson 對話，1972），核心量是 Bessel 核 \(K(x) = \frac{\sin \pi(x - y)}{\pi (x - y)}\) 的變體與 \(J\) 函數的 Fredholm 行列式。Bessel 函數站在數學最深猜想的門口（黎曼假設的統計面向）。

場景五：現代計算中的 \(J_\nu\) 。快速 Bessel 變換（fast spherical/Bessel transforms）使天文映射（CMB 多極分析）、重力場建模在億萬模式規模下可行；不穩定遞推（向前算 \(J_n\) 會爆誤差—— \(J\) 遞推向下穩定、向上不穩定，10.1 節場景）要求向後遞推（Miller 演算法，1935）——計算特殊函數本身就是一堂數值穩定性課。

補充： \(\Gamma\) 函數——階乘的補票#

Euler 1729 年致 Goldbach 的信中定義 \(\Gamma(s) = \int_0^\infty t^{s-1} e^{-t}\, dt\) ，使階乘延拓到非整數： \(\Gamma(1/2) = \sqrt\pi\) 、 \(\Gamma(s + 1) = s \Gamma(s)\) 。Bessel 級數因 \(\Gamma\) 而擁有連續階 \(\nu\) ；量子場論的發散重整化、機率的 Beta 分佈都靠它。一個「補票」函數，成為特殊函數論的通用護照。

補充：為什麼柱面波衰減 \(\frac{1}{\sqrt r}\) 而球面波是 \(\frac{1}{r}\)#

能量守恆：波攜帶的能量通過圓柱面（面積 \(\propto r\) ）與球面（ \(\propto r^2\) ）。振幅平方 = 能量密度，故柱面振幅 \(\propto r^{-1/2}\) 、球面 \(\propto r^{-1}\) 。數學印證：漸近式 \(J_\nu \sim \sqrt{\frac{2}{\pi x}} \cos(\cdots)\) ；球面對應物（球 Bessel 函數） \(j_0(x) = \frac{\sin x}{x}\) 恰是 \(\frac{1}{x}\) 乘振盪。幾何的維數，直接寫在漸近式的指數上——這是漸近分析最壯觀的一次「物理直覺 = 數學結構」。

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

• Bessel 方程： \(x^2 y'' + x y' + (x^2 - \nu^2) y = 0\) ； \(x = 0\) 為正則奇點。
• ** \(J_\nu\) **：第一類 Bessel 函數；在原點有限。
• ** \(Y_\nu\) （Neumann 函數）**：第二類解；對數發散於原點。
• 指標方程： \(r^2 - \nu^2 = 0\) ；三種劇本（非整數差/重根/整數差）。
• ** \(\Gamma\) 函數**：階乘的全域推廣（Euler 1729）。
• 漸近式： \(J_\nu \sim \sqrt{\frac{2}{\pi x}} \cos(x - \frac{\nu\pi}{2} - \frac{\pi}{4})\) 。
• Fourier–Bessel 級數：以 \(J_\nu(\alpha_m r / R)\) 展開；圓盤的調和分析。
• Hankel 函數： \(H_\nu^{(1,2)} = J_\nu \pm i Y_\nu\) ；出射/入射柱面波。
• 修正 Bessel 函數： \(I_\nu, K_\nu\) ；虛宗量、指數型行為。
• Miller 演算法：向後遞推計算 \(J_n\) 的穩定方法。

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

• 1732：Daniel Bernoulli 吊鏈振動（ \(J\) 級數首現）。
• 1764：Euler 圓膜振動分析。
• 1729：Euler \(\Gamma\) 函數。
• 1824：Bessel 行星擾動論文（系統研究與命名起源）。
• 1838：Bessel 測得恆星視差。
• 1867：Neumann 函數。
• 1877：Rayleigh《聲學理論》，圓膜模式。
• 1879：Rayleigh 判據（衍射解析度， \(J_1\) 零點）。
• 1922：FM 理論的 Bessel 邊帶（Carson）。
• 1934：Langer 均匀漸近。
• 1972：Montgomery–Dyson（Bessel 核與 ζ 零點統計）。

深化問答三則#

問： \(J_0(0) = 1\) 而 \(J_1(0) = 0\) ，這對物理意味著什麼？ 答： \(\nu = 0\) 的模式在圓心處振幅非零（徑向無節點過圓心，模式「呼吸」）， \(\nu \ge 1\) 的模式在圓心為零（有 \(\phi\) 依賴 \(\cos m\phi\) ，圓心是 \(m\) 條節線交匯處，振幅必為零）。解在奇點處的行為（ \(J_\nu \sim x^\nu\) ）直接編碼了模式的幾何——Frobenius 指標 \(r = \nu\) 不是技術細節，是節線的出生證明。

問：為什麼 Bessel 方程的 \(x y'\) 項不可省略？ 答： \(x y'\) 項來自極座標的 Laplacian（ \(\frac1r \partial_r (r \partial_r)\) 的展開）——它是幾何的簽名。省略它得到的 \(y'' + y = 0\) 是平面波方程，振盪不衰減；加上它，能量守恆強制的 \(\frac{1}{\sqrt r}\) 衰減才出現（見補充）。每一個係數都是幾何的翻譯：這也是為什麼「變係數」不該被視為「常係數的麻煩版」——係數的變化就是空間的形狀。

問：向前遞推 \(J_n\) 為什麼會「爆」？ 答：Bessel 方程有兩個解 \(J_n\) （小）與 \(Y_n\) （大， \(\sim x^{-1}\) 相對而言其實是「大」的一支在遞推意義下），遞推 \(J_{n+1} = \frac{2n}{x} J_n - J_{n-1}\) 中捨入誤差會混入 \(Y_n\) 成分並按 \(n\) 放大——高階 \(n\) 時誤差壓倒真值。Miller 演算法反其道：從 \(n = N\) （誤差貢獻微小處）向下遞推，誤差按真解比例衰減。遞推方向決定生死——這是 3.3 節「剛性」思想的離散親戚，也是每個科學計算工程師的必修一課。

（本節完）