贝塞尔函数

贝塞尔函数与柱坐标波动

物理问题引入：圆形鼓膜为什么能敲出丰富的泛音？

一根一维琴弦振动时，两端固定边界条件简单地给出了正弦驻波 $\sin (n π x / L)$ sin(npi x/L)，泛音频率是基频的整倍数（ $1 : 2 : 3 : 4$ 1:2:3:4）．然而，当你敲击一个 圆形鼓膜 时，鼓声听起来却完全不像弦乐那样具有和谐的音高感．物理测量表明：圆形鼓的泛音频率比值约为：

1 : 1.593 : 2.135 : 2.295 : 2.653 : \dots

1 : 1.593 : 2.135 : 2.295 : 2.653 : dots

这些比值全是无理数！为什么圆形几何会产生这种非整倍数的分立振动模态？当我们把拉普拉斯算符改写为柱坐标系 $(r, θ, z)$ (r, theta, z) 时，径向微分方程自然引出了圆对称系统的标准数学工具——贝塞尔函数（Bessel Functions）．

1. 贝塞尔方程与两类解

在柱坐标系下，波动方程或亥姆霍兹方程分离变量后，径向函数 $R (r)$ R(r) 满足 $ν$ nu 阶贝塞尔方程（令无量纲宗量 $x = k r$ x = kr）：

x^{2} \frac{d^{2} y}{d x^{2}} + x \frac{d y}{d x} + (x^{2} - ν^{2}) y = 0

x^2 d^2 ydx^2 + x dydx + (x^2 - nu^2) y = 0

这是一个带有正则奇点 $x = 0$ x = 0 的二阶常微分方程，存在两个线性无关的特解：

1.1 第一类贝塞尔函数 $J_{ν} (x)$ J_nu(x)

由 Frobenius 广义幂级数解出，在原点 $x = 0$ x = 0 处 严格保持有限值：

J_{ν} (x) = \sum_{m = 0}^{\infty} \frac{(- 1)^{m}}{m! Γ (m + ν + 1)} {(\frac{x}{2})}^{2 m + ν}

J_nu(x) = sum_m=0^inf (-1)^mm! , Gamma(m + nu + 1) left(x2right)^2m + nu

当 $ν = 0$ nu = 0 时， $J_{0} (0) = 1$ J_0(0) = 1；当 $ν > 0$ nu > 0 时， $J_{ν} (0) = 0$ J_nu(0) = 0．

1.2 第二类贝塞尔函数（诺伊曼函数 $Y_{ν} (x)$ Y_nu(x)/ $N_{ν} (x)$ N_nu(x)）

第二个线性无关解在原点 $x \to 0$ x -> 0 处 发散趋向于负无穷大（含对数奇性）：

Y_{ν} (x) = \frac{J_{ν} (x) \cos (ν π) - J_{- ν} (x)}{\sin (ν π)} (x \to 0 时 Y_{0} (x) \propto \ln x)

Y_nu(x) = J_nu(x)cos(nupi) - J_-nu(x)sin(nupi) quad (x -> 0 时 Y_0(x) propto ln x)

如何为物理边界挑选特解？

实心内部问题（如实心圆形鼓膜、圆柱形激光光纤核心、实心导体圆棒）：因为原点包含在物理介质内部，物理量不能无端发散，强制要求第二类解系数为零，只保留有限解 $J_{ν} (x)$ J_nu(x)；
同轴/环状空心问题（如两同心导体圆筒之间的同轴电缆空腔、空心圆柱管道）：原点被排除在介质之外，通解必须同时保留 $J_{ν}$ J_nu 与 $Y_{ν}$ Y_nu 满足内壁与外壁的联合边界条件．

2. 关键微积分与渐近性质

2.1 递推关系

微商与不同阶贝塞尔函数之间满足精巧的互锁关系：

J_{ν - 1} (x) + J_{ν + 1} (x) = \frac{2 ν}{x} J_{ν} (x)

J_nu-1(x) + J_nu+1(x) = 2x J_nu(x)

\frac{d}{d x} [x^{ν} J_{ν} (x)] = x^{ν} J_{ν - 1} (x)

ddxleft[ x^nu J_nu(x) right] = x^nu J_nu-1(x)

2.2 远场大宗量渐近展开（退化为衰减平面波）

当远离轴心（ $x = k r ≫ 1$ x = kr gg 1）时，贝塞尔函数的行为退化为调幅正弦波：

J_{ν} (x) \approx \sqrt{\frac{2}{π x}} \cos (x - \frac{ν π}{2} - \frac{π}{4})

J_nu(x) ~ 2pi x cosleft( x - nu2 - 4 right)

物理图像：随着柱面波向外扩散，波前周长按 $2 π r$ 2pi r 线性扩展；根据能量守恒，振幅必须按 $1 / \sqrt{r}$ 1/r 衰减！渐近公式正是圆柱扩散能量守恒的直接体现．

3. 贝塞尔函数的带权正交性

在区间 $[0, a]$ [0, a] 上，设 $α_{ν, n}$ alpha_nu, n 为 $J_{ν} (x)$ J_nu(x) 的第 $n$ n 个正实数零点（即 $J_{ν} (α_{ν, n}) = 0$ J_nu(alpha_nu, n) = 0），令特征波矢 $k_{n} = α_{ν, n} / a$ k_n = alpha_nu, n/a．贝塞尔函数在带权函数 $w (r) = r$ w(r) = r 下构成 Sturm–Liouville 正交系：

\int_{0}^{a} J_{ν} (k_{n} r) J_{ν} (k_{m} r) r d r = \frac{a^{2}}{2} {[J_{ν + 1} (α_{ν, n})]}^{2} δ_{n m}

int_0^a J_nu(k_n r) J_nu(k_m r) , r , dr = a^22 left[ J_nu+1(alpha_nu, n) right]^2 delta_nm

这构成了圆盘区域内任何物理场（如鼓膜位移、波导截面场）展开为 傅里叶 - 贝塞尔级数 的数学基础．

4. 经典物理范例

4.1 圆形鼓膜的简正模态

两端固定圆形膜半径为 $a$ a，周边位移为零： $u (a, θ, t) = 0$ u(a, theta, t) = 0．分离变量解为 $u (r, θ, t) = J_{m} (k r) \cos (m θ) \cos (ω t)$ u(r, theta, t) = J_m(k r) cos(mtheta) cos(omega t)．边界条件强制要求： $J_{m} (k a) = 0 ⟹ k_{m n} = \frac{α_{m, n}}{a}$ J_m(k a) = 0 implies k_mn = alpha_m, na．固有角频率为 $ω_{m n} = v k_{m n} = \frac{v}{a} α_{m, n}$ omega_mn = v k_mn = va alpha_m, n．

基频（基模 $(0, 1)$ (0, 1)）： $α_{0, 1} \approx 2.4048$ alpha_0, 1 ~ 2.4048；
第一泛音（一节线模 $(1, 1)$ (1, 1)）： $α_{1, 1} \approx 3.8317$ alpha_1, 1 ~ 3.8317（比值 $3.8317 / 2.4048 \approx 1.593$ 3.8317 / 2.4048 ~ 1.593）；
第二泛音（一节圆模 $(0, 2)$ (0, 2)）： $α_{0, 2} \approx 5.5201$ alpha_0, 2 ~ 5.5201（比值 $5.5201 / 2.4048 \approx 2.295$ 5.5201 / 2.4048 ~ 2.295）．鼓面所有泛音比值，精准来自于贝塞尔函数的零点坐标！

4.2 夫琅禾费圆孔衍射与艾里斑（Airy Disk）

在物理光学中，波长为 $λ$ lambda 的平面波垂直照射直径为 $D = 2 a$ D = 2a 的圆孔时，远场衍射场强由二维傅里叶变换给出：

I (θ) = I_{0} {[\frac{2 J_{1} (k a \sin θ)}{k a \sin θ}]}^{2}

I(theta) = I_0 left[ 2 J_1(kasintheta)kasintheta right]^2

第一个暗环（强度降为零）对应 $J_{1} (x)$ J_1(x) 的第一个非零根： $x_{1} = α_{1, 1} \approx 3.8317$ x_1 = alpha_1, 1 ~ 3.8317．因为 $k a \sin θ = \frac{2 π a}{λ} \sin θ = \frac{π D}{λ} \sin θ = 3.8317$ kasintheta = 2pi a/lambdasintheta = pi D/lambdasintheta = 3.8317：

\sin θ_{1} = \frac{3.8317}{π} \frac{λ}{D} \approx 1.22 \frac{λ}{D}

sintheta_1 = 3.8317pi D ~ 1.22 D

这就是天文学与光学仪器中大名鼎鼎的 瑞利判据（Rayleigh Criterion）与艾里斑角分辨率常数 1.22 的纯数学来源！

5. 学习衔接

下一节：量子力学谐振子势能中的正交多项式，进入厄米多项式；
物理应用：在光学：衍射与干涉中掌握艾里斑衍射图样．

本页面最近更新：2026/10/1 10:49:01，更新历史
发现错误？想一起完善？在 GitHub 上编辑此页！
本页面贡献者：Leafuke, Physics Learning Wiki Team
本页面的全部内容在 CC BY-SA 4.0 和 SATA 协议之条款下提供，附加条款亦可能应用

贝塞尔函数