图

贝赛尔函数

   贝赛尔方程为

\begin{equation} x^2 \dvTwo[2]{y}{x} + x\dvTwo{y}{x} + (x^2 - \alpha ^2)y = 0 \end{equation}
其中 $\alpha$ 叫做阶数. 两个线性无关的解分别是第一类贝赛尔函数 $J_\alpha(x)$ 和第二类贝赛尔函数 $Y_\alpha(x)$. 这里只讨论 $x > 0$ 且 $\alpha$ 为整数或半整数的情况.

图
图1:第一类和第二类贝赛尔函数

   $J_\alpha(x)$ 的级数形式为

\begin{equation} J_\alpha(x) = \sum_{m=0}^\infty \frac{(-1)^m}{m!\Gamma(m+\alpha+1)} \qtyRound{\frac x2}^{2m+\alpha} \end{equation}
其中使用了 $\Gamma$ 函数. $Y_\alpha(x)$ 可以通过 $J_\alpha(x)$ 来定义
\begin{equation} Y_\alpha(x) = \frac{J_\alpha(x)\cos\alpha\pi - J_{-\alpha}(x)}{\sin\alpha\pi} \end{equation}

汉克尔函数

   贝赛尔方程的两个线性无关解也可以用第一类和第二类汉克尔函数来表示1

\begin{equation} H_\alpha ^{(1)}(x) = J_\alpha(x) + \I Y_\alpha(x) \qquad H_\alpha ^{(2)}(x) = J_\alpha(x) - \I Y_\alpha(x) \end{equation}

常用性质

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

\begin{equation} Z_{-\alpha}(x) = (-1)^\alpha Z_\alpha(x) \end{equation}
递推关系
\begin{equation} \frac{2\alpha}{x} Z_\alpha(x) = Z_{\alpha -1}(x) + Z_{\alpha+1}(x) \end{equation}
一阶导数
\begin{equation} \dvTwo{Z_\alpha}{x} = \frac12 [Z_{\alpha - 1}(x) - Z_{\alpha +1}(x)] \end{equation}
正交关系
\begin{equation} \int_0^1 x J_\alpha (u_{\alpha ,m} x) J_\alpha (u_{\alpha ,n} x) \dd{x} = \frac{\delta_{m,n}}{2}[J_{\alpha + 1} (u_{\alpha ,m})]^2 \end{equation}
其中 $u_{\alpha, m}$ 是 $J_\alpha(x)$ 的第 $m$ 个根.

修正贝赛尔函数

   修正贝赛尔方程为

\begin{equation} x^2 \dvTwo[2]{y}{x} + x \dvTwo{y}{x} - (x^2 + \alpha ^2)y = 0 \end{equation}
其解为两个修正贝赛尔函数(Modified Bessel Function),第一类为 $I_\alpha(x)$, 第二类为 $K_\alpha(x)$, 与贝赛尔函数的关系为
\begin{equation} I_\alpha(x) = \I^{-\alpha} J_\alpha(\I x) \qquad K_\alpha(x) = \frac{\pi}{2} \I^{\alpha + 1} H_\alpha ^{(1)}(\I x) \end{equation}


1. 可类比欧拉公式 $\expRound{\pm\I x} = \cos x \pm \I\sin x$ (式 3 ).

致读者: 小时物理百科一直以来坚持所有内容免费且不做广告,这导致我们处于日渐严重的亏损状态。长此以往很可能会最终导致我们不得不选择商业化,例如大量广告,内容付费,会员制,甚至被收购。因此,我们鼓起勇气在此请求广大读者热心捐款,使网站得以健康发展。如果看到这条信息的每位读者能慷慨捐助 10 元,我们几天内就能脱离亏损状态,并保证网站能在接下来的一整年里向所有读者继续免费提供优质内容。感谢您的支持。
—— 小时(项目创始人)

编辑词条 返回目录 返回主页 捐助项目 © 小时物理百科 保留一切权利