这个方程从何而来
上一篇指南解了勒让德方程，那是一切具有球对称性之物的主方程。现在把球换成圆柱：一面振动的鼓膜、一根光纤的横截面、热量沿圆管向内渗透。在圆盘上写下波动方程或热方程，用分离变量法推下去，角向部分很容易（绕圆一圈的正弦与余弦而已）——但径向部分，那个依赖于到圆心距离 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 模态上展开。
退一步，欣赏一下刚刚发生的事。一个关于实物的、乱糟糟的偏微分方程，在五个诚实的动作里，溶解成了去查一个特殊函数的零点。鼓的音高之量子化——存在一个最低音、其上是一架离散的阶梯——不是被强加的；它出自唯一的要求：解须在圆心处保持有限、在边缘处为零。这正是整级最深的回报：物理那些著名的函数并不奇异，它们不过是当一个级数解必须服从它所栖居的几何时，存活下来的那些自然形状。