這個方程從何而來
上一篇指南解了勒讓德方程，那是一切具有球對稱性之物的主方程。現在把球換成圓柱：一面振動的鼓膜、一根光纖的橫截面、熱量沿圓管向內滲透。在圓盤上寫下波動方程或熱方程，用分離變量法推下去，角向部分很容易（繞圓一圈的正弦與餘弦而已）——但徑向部分，那個依賴於到圓心距離 r 的部分，卻一再生出同一個頑固的方程。那個方程就是貝塞爾方程。勒讓德之於球，正如貝塞爾方程之於圓柱。
寫全了，n 階貝塞爾方程是 x^2 y'' + x y' + (x^2 - n^2) y = 0，其中 x 本質上是縮放後的半徑，參數 n（即階）是角向分離塞給你的那個數——它不必是整數。注意領頭係數 x^2：它在 x = 0 處為零，正是鼓的圓心。除過去把方程化成標準形式 y'' + P(x) y' + Q(x) y = 0，便得 P(x) = 1/x 與 Q(x) = 1 - n^2/x^2，二者都在原點爆掉。所以 x = 0 是一個奇點——而它恰好坐落在物理發生的地方，圓柱的軸線上。一個樸素的泰勒級數在那裡起步是靠不住的。
用弗羅貝尼烏斯方法造出第一個解
既然 x = 0 是正則奇點，樸素的冪級數就是錯的擬設——而弗羅貝尼烏斯方法恰恰對路。設一個帶未知領頭冪次的級數：y = x^r 乘 (a_0 + a_1 x + a_2 x^2 + ...)，其中 a_0 非零、指數 r 待定。把它代入 x^2 y'' + x y' + (x^2 - n^2) y = 0，盯住 x 的最低次冪。因為 a_0 非零，那個最低次冪的係數必須自行為零，於是給出指標方程 r^2 - n^2 = 0。它的兩個根是 r = +n 與 r = -n——解從奇異圓心冒出時兩種可能的領頭指數。
先取較大的根 r = +n；它總是給出一個乾淨的解。把它代回，要求 x 的每一個更高次冪也都為零，便產生一個把各係數串起來的遞推關係：a_2m = -a_{2m-2} / (4 m (m + n))。只有偶數下標的係數存活（奇數的全被逼為零），且每一個都是前一個除以一個齊整的乘積。從 a_0 搖起遞推，滾落出來的圖樣在分母裡含有 m! 與乘積 (n+1)(n+2)...(n+m)——這個乘積簡直求著要用伽馬函數來寫，因為 Gamma(m + n + 1) 恰好把這種類階乘的增長打包起來，即便 n 不是整數也照樣成立。
Bessel: x^2 y'' + x y' + (x^2 - n^2) y = 0, singular point x = 0
Frobenius ansatz: y = x^r * SUM_{m>=0} a_m x^m, a_0 != 0
lowest power of x -> indicial equation r^2 - n^2 = 0 -> r = +n, -n
take r = +n, match each power -> recurrence a_{2m} = - a_{2m-2} / (4 m (m+n))
(odd a's all zero)
choose a_0 = 1 / (2^n Gamma(n+1)) -> the standard first-kind solution
J_n(x) = SUM_{m>=0} (-1)^m / ( m! Gamma(m+n+1) ) * (x/2)^{2m+n}
leading term ~ (x/2)^n / Gamma(n+1) -> J_n(0)=0 for n>0, J_0(0)=1 (finite!)把自由常數 a_0 定到慣用值 1/(2^n Gamma(n+1))，這個級數就有了名字：J_n(x)，n 階第一類貝塞爾函數。它是有圓心世界裡餘弦與正弦的徑向對應物。要緊的是，每一項前面都拖著因子 x^n（即我們選的 +n 指數），所以在原點附近 J_n(x) 表現得像 (x/2)^n / Gamma(n+1)：當 n > 0 時它從零起步、平滑上升，而 J_0(0) = 1。無論哪種情形它都在圓心處有限——真實鼓膜在其軸線上必須具備的正是這種行為。
第二個解，以及它為何不老實
二階方程需要兩個獨立的解，那另一個在哪兒？指標方程交給我們第二個根 r = -n。對許多 n 值，你可以用這個指數跑同樣的遞推，得到一個真正獨立的級數 J_{-n}(x)。但指標根本身裡就埋著一個機關。兩根是 +n 與 -n，它們的差是 2n。當 n 是整數時（或半整數使 2n 成整數），這個差是一個整數——而指標方程的三分律警告：每當兩根相差一個整數，第二個弗羅貝尼烏斯級數就可能崩壞，一個對數項可能被逼進解裡。對整數階，J_{-n} 結果只是 J_n 的一個常數倍，根本不獨立——於是樸素的第二級數塌掉了。
補救之法是親手拼出一個真正獨立的第二解。定義 Y_n(x)，即第二類貝塞爾函數（有時稱諾伊曼函數），它是 J_n 與 J_{-n} 的一個特定組合，專門設計成能在整數階的極限下存活。這番設計的代價是一個對數：在原點附近 Y_n(x) 含有一項正比於 ln(x)（且當 n > 0 時，上面還疊加一塊像 1/x^n 那樣爆破的部分）。所以 Y_n(x) 在圓心處無界——當 x 趨於零時它向負無窮俯衝。那個對數式的爆破，正是兩根相等或相差整數這一情形可見的指紋，恰是弗羅貝尼烏斯理論早已警告過我們的那種情形。
它們究竟如何表現：衰減的餘弦
級數把 x = 0 附近的一切都告訴你，但要描繪遠處的函數它就無能為力了——你得用上幾百項。遠場行為來自一次單獨的漸近分析（與漸近那一級裡同樣的大宗量機械），而結果出奇地具體。當 x 很大時，J_n(x) 近似等於 sqrt(2/(pi x)) 乘 cos(x - n pi/2 - pi/4)，而 Y_n(x) 是同樣的振幅乘 sin(x - n pi/2 - pi/4)。把它念出來：一道餘弦（或正弦）波，錯開一個相位，其振幅像 1/sqrt(x) 那樣消退。它們是衰減的餘弦——恰是一圈圓形漣漪的畫面：它一邊向外行進一邊振盪，高度卻不斷下降，因為同樣的能量被攤到越來越大的圓周上。
對那個近似是什麼，要誠實。大 x 的餘弦公式是一個漸近展開，而非收斂級數：你不斷添上修正項，對任何固定的 x 它最終都會發散。然而在恰當之處截斷，它精確得驚人，而且 x 越大，單憑領頭項就越好。這是高等應用微積分反覆出現的教訓——一個發散的級數，用得明智，能比收斂的級數預測得更準。J_n 的冪級數是原點附近的精確工具；漸近形式是遠方的精確工具；二者都不是「那個公式」，真實的計算把它們縫在一起。
因為兩個函數都永遠振盪，每一個都無窮多次穿過零點，而那些零點正是物理被釘死的地方。把 J_n 的第 m 個正零點記作 j_{n,m}：J_n(j_{n,m}) = 0。在一面半徑為 R、邊緣被夾緊的鼓上，邊界條件是位移在那裡為零，這迫使縮放後的半徑落在某個零點上——於是挑出一組離散的允許振動頻率，即鼓的泛音。與弦不同，這些泛音不是某個基頻的整數倍（J_0 的零點間距並不均勻），這恰恰是為什麼鼓聽起來像鼓、而不像被撥動的弦。
一步步解一面真實的鼓
讓我們把各塊拼回到最初的問題上：一面半徑為 R、邊緣夾緊的圓鼓膜，以完全圓對稱（與角度無關）的模式振動。豎直位移 u 通過波動方程依賴於半徑 r 與時間 t，整樁任務就是找出哪些形狀與頻率被允許。下面是完整的弧線，從一個偏微分方程一路落到一個你能測量的數。
- 分離變量。寫 u(r, t) = F(r) G(t)，代入極坐標形式的波動方程。時間部分給出 G(t) = cos(omega t)（在頻率 omega 上的簡單振盪），而徑向部分被迫滿足 F'' + (1/r) F' + k^2 F = 0，其中 k = omega/c 把空間尺度與頻率綁在一起。
- 認出貝塞爾。代換 x = k r 來重新縮放，徑向方程就變成 x^2 F'' + x F' + x^2 F = 0——零階貝塞爾方程（n = 0，因為該模式沒有角向變化）。它的兩個解是 J_0(x) 與 Y_0(x)。
- 要求圓心處有限。區域包含 r = 0，鼓的正中，那裡位移必須有限。但 Y_0 在那裡對數式地爆破，所以它的係數必須為零。只有 F(r) = J_0(k r) 存活——幾何已經丟掉了一半的解空間。
- 套上夾緊的邊緣。邊緣 r = R 被固定不動，所以 F(R) = 0，即 J_0(k R) = 0。這迫使 k R 取 J_0 的某個零點：k R = j_{0,m}，m = 1, 2, 3, ...。允許的頻率是 omega_m = c j_{0,m} / R——一個離散譜，直接從貝塞爾函數的零點上讀出。
- 構造一般運動。每個 m 給出一個模態 J_0(j_{0,m} r/R) cos(omega_m t)。完整的振動是它們全體的疊加，而各權重由初始形狀經一個傅立葉-貝塞爾級數定下——它是普通傅立葉級數在圓柱情形下的表親，把起始輪廓在那些 J_0 模態上展開。
退一步，欣賞一下剛剛發生的事。一個關於實物的、亂糟糟的偏微分方程，在五個誠實的動作裡，溶解成了去查一個特殊函數的零點。鼓的音高之量子化——存在一個最低音、其上是一架離散的階梯——不是被強加的；它出自唯一的要求：解須在圓心處保持有限、在邊緣處為零。這正是整級最深的回報：物理那些著名的函數並不奇異，它們不過是當一個級數解必須服從它所棲居的幾何時，存活下來的那些自然形狀。