贝塞尔函数

贡献者： addis

预备知识　二阶齐次变系数线性微分方程的幂级数解法

1. 贝塞尔函数

　　 贝塞尔微分方程（Bessel's differential equation）也被称为柱谐函数、圆柱函数或圆柱谐波，通常在使用分离变量法求解柱坐标中的拉普拉斯方程时得到：

\begin{equation} x \frac{\mathrm{d}}{\mathrm{d}{x}} \left(x \frac{\mathrm{d}{y}}{\mathrm{d}{x}} \right) + (x^2 - l ^2)y = 0~. \end{equation}

其中 $l$ 叫做阶数（order），一般来说可以是任意实数，但半整数和整数较为常见。两个线性无关的解分别是第一类贝塞尔函数（Bessel function of the first kind） $J_l(x)$ 和第二类贝塞尔函数（Bessel function of the second kind） $Y_l(x)$。这里只讨论 $x > 0$ 且 $l$ 为整数或半整数的情况¹。

图 1：第一类和第二类贝塞尔函数（非负整数阶）

　　 $J_l(x)$ 的级数形式为

\begin{equation} J_l(x) = \sum_{m=0}^\infty \frac{(-1)^m}{m!\Gamma(m+l+1)} \left(\frac x2 \right) ^{2m+l}~, \end{equation}

其中使用了 $\Gamma$ 函数。$Y_l(x)$ 可以通过 $J_l(x)$ 来定义

\begin{equation} Y_l(x) = \frac{J_l(x) \cos\left(l\pi\right) - J_{-l}(x)}{ \sin\left(l\pi\right) }~. \end{equation}

容易证明

\begin{equation} J_{1/2}(x) = \sqrt{\frac{2}{\pi x}} \sin x = -Y_{-1/2}(x)~, \end{equation}

\begin{equation} J_{-1/2}(x) = \sqrt{\frac{2}{\pi x}} \cos x = -Y_{1/2}(x)~. \end{equation}

变量替换

　　令 $u(x) = \sqrt{x} y(x)$，对式 1 替换变量可得更简洁的方程，使一阶导数项消失

\begin{equation} x^2 \frac{\mathrm{d}^{2}{u}}{\mathrm{d}{x}^{2}} + \left(x^2 - l^2 + \frac{1}{4} \right) u = 0~, \end{equation}

显然 $u(x)$ 的两组解就是 $\sqrt{x}J_l(x)$ 和 $\sqrt{x}Y_l(x)$。

汉克尔函数

　　贝塞尔方程的两个线性无关解也可以用第一类和第二类汉克尔函数（Hankel function）来表示²

\begin{equation} H_l ^{(1)}(x) = J_l(x) + \mathrm{i} Y_l(x)~, \qquad H_l ^{(2)}(x) = J_l(x) - \mathrm{i} Y_l(x)~. \end{equation}

2. 常用性质

　　令 $Z$ 为 $J, Y, H^{(1)}, H^{(2)}$ 的任意一种，则

\begin{equation} Z_{-l}(x) = (-1)^l Z_l(x)~. \end{equation}

递推关系

\begin{equation} \frac{2l}{x} Z_l(x) = Z_{l -1}(x) + Z_{l+1}(x)~. \end{equation}

一阶导数

\begin{equation} \frac{\mathrm{d}{Z_l}}{\mathrm{d}{x}} = \frac12 [Z_{l - 1}(x) - Z_{l +1}(x)]~. \end{equation}

渐进形式（$x \gg 1$）

\begin{equation} J_l(x) \to \sqrt{\frac{2}{\pi x}} \cos \left(x - \frac{l\pi}{2} - \frac{\pi}{4} \right) ~, \end{equation}

\begin{equation} Y_l(x) \to \sqrt{\frac{2}{\pi x}} \sin \left(x - \frac{l\pi}{2} - \frac{\pi}{4} \right) ~. \end{equation}

正交关系

\begin{equation} \int_0^1 J_l (u_{l ,m} x) J_l (u_{l ,n} x) x \,\mathrm{d}{x} = \frac{\delta_{m,n}}{2}[J_{l + 1} (u_{l ,m})]^2~. \end{equation}

其中 $u_{l, m}$ 是 $J_l(x)$ 的第 $m$ 个根。

　　连续正交关系（容易从渐进形式得到）

\begin{equation} \int_0^\infty \sqrt{k_1} J_l (k_1 x) \sqrt{k_2}J_l (k_2 x) x \,\mathrm{d}{x} = \delta(k_2 - k_1)~. \end{equation}

3. 修正贝塞尔函数

　　修正贝塞尔方程为

\begin{equation} x^2 \frac{\mathrm{d}^{2}{y}}{\mathrm{d}{x}^{2}} + x \frac{\mathrm{d}{y}}{\mathrm{d}{x}} - (x^2 + l ^2)y = 0~. \end{equation}

其解为两个修正贝塞尔函数（Modified Bessel Function），第一类为 $I_l(x)$，第二类为 $K_l(x)$，与贝塞尔函数的关系为

\begin{equation} I_l(x) = \mathrm{i} ^{-l} J_l( \mathrm{i} x)~, \qquad K_l(x) = \frac{\pi}{2} \mathrm{i} ^{l + 1} H_l ^{(1)}( \mathrm{i} x)~. \end{equation}

1. ^ 更详细的讨论请参见二阶齐次变系数线性微分方程的幂级数解法文章中的例 4 到例 5 部分。
2. ^ 可类比欧拉公式 $ \exp\left(\pm \mathrm{i} x\right) = \cos x \pm \mathrm{i} \sin x$（式 3 ）。

致读者：小时百科一直以来坚持所有内容免费无广告，这导致我们处于严重的亏损状态。长此以往很可能会最终导致我们不得不选择大量广告以及内容付费等。因此，我们请求广大读者热心打赏 ，使网站得以健康发展。如果看到这条信息的每位读者能慷慨打赏 20 元，我们一周就能脱离亏损，并在接下来的一年里向所有读者继续免费提供优质内容。但遗憾的是只有不到 1% 的读者愿意捐款，他们的付出帮助了 99% 的读者免费获取知识，我们在此表示感谢。