免费贝塞尔函数计算器

贝塞尔函数是贝塞尔微分方程的解，这是一个出现在物理和工程许多领域的二阶线性常微分方程，从波传播到热传导。贝塞尔函数计算器是一个在线工具，用于计算实数阶 $\nu$ 和自变量 $x$（可为实数或复数）的 $J_\nu(x)$（第一类）和 $Y_\nu(x)$（第二类）。当 $x$ 为实数时，该工具还会生成贝塞尔函数图，直观展示函数的性质。

贝塞尔微分方程

贝塞尔方程写作

$$x^2 \frac{d^2 y}{dx^2} + x \frac{dy}{dx} + (x^2 - \nu^2) y = 0,$$

其中 $\nu$（阶数）是一个常数。这是一个二阶常微分方程，它的两个线性无关解是第一类和第二类贝塞尔函数，记为 $J_\nu(x)$ 和 $Y_\nu(x)$；它们统称为柱函数。

第一类贝塞尔函数

第一类函数 $J_\nu(x)$ 由幂级数定义

$$J_\nu(x) = \sum_{k=0}^{\infty} \frac{(-1)^k}{k! \, \Gamma(\nu + k + 1)} \left( \frac{x}{2} \right)^{\nu + 2k},$$

其中 $\Gamma(z)$ 是 Gamma 函数（阶乘在非整数上的推广）。对于非整数 $\nu$，函数 $J_\nu(x)$ 和 $J_{-\nu}(x)$ 是独立的；对于整数阶 $\nu = n$，有

$$J_{-n}(x) = (-1)^n J_n(x).$$

在实际中，该级数在适度数量的项后收敛，使得数值计算变得直截了当。

第二类贝塞尔函数

第二类函数 $Y_\nu(x)$（也记作 $N_\nu(x)$）是分段定义的，因为整数阶和非整数阶的公式不同。

非整数 $\nu$

对于 $\nu \notin \mathbb{Z}$，

$$Y_\nu(x) = \frac{J_\nu(x) \cos(\nu\pi) - J_{-\nu}(x)}{\sin(\nu\pi)}.$$

这个组合提供了与 $J_\nu(x)$ 线性无关的第二个独立解。

整数 $\nu = n$

当阶数为整数时，上述表达式变得不确定。正确形式通过取极限 $\nu \to n$ 得到：

$$\begin{aligned} Y_n(x) = &\; \frac{2}{\pi} J_n(x) \ln\frac{x}{2} \\ &- \frac{1}{\pi} \sum_{k=0}^{n-1} \frac{(n-k-1)!}{k!} \left(\frac{x}{2}\right)^{2k-n} \\ &- \frac{1}{\pi} \sum_{k=0}^{\infty} \frac{(-1)^k}{k! (n+k)!} \left(\frac{x}{2}\right)^{n+2k} \bigl[ \psi(k+1) + \psi(n+k+1) \bigr], \end{aligned}$$

其中 $\psi(z)$ 是 digamma 函数（$\Gamma(z)$ 的对数导数）。对于非负整数

$$\psi(n) = H_{n-1} - \gamma,$$

其中 $H_{n-1}$ 是第 $n-1$ 个调和数，$\gamma \approx 0.57721$ 是欧拉-马斯刻若尼常数。对于负整数，可利用对称性 $Y_{-n}(x) = (-1)^n Y_n(x)$。

汉克尔函数

第三类贝塞尔函数，即汉克尔函数，是 $J_\nu$ 和 $Y_\nu$ 的线性组合：

$$H_\nu^{(1)}(x) = J_\nu(x) + i Y_\nu(x), \qquad H_\nu^{(2)}(x) = J_\nu(x) - i Y_\nu(x),$$

其中 $i = \sqrt{-1}$。它们特别适用于描述向外和向内传播的波。

递推关系

所有柱函数（$C_\nu$ 表示 $J_\nu$ 或 $Y_\nu$）满足递推关系

$$C_{\nu-1}(x) + C_{\nu+1}(x) = \frac{2\nu}{x} C_\nu(x),$$ $$C_{\nu-1}(x) - C_{\nu+1}(x) = 2 C'_\nu(x),$$

其中撇号表示对 $x$ 的导数。这些恒等式允许在不重新计算级数的情况下快速计算导数或相邻阶数的值。

使用此计算器

要获取贝塞尔函数值，请输入阶数 $\nu$（实数）和点 $x$（实数或复数）。计算器同时返回 $J_\nu(x)$、$Y_\nu(x)$ 以及两个汉克尔函数 $H_\nu^{(1)}(x)$、$H_\nu^{(2)}(x)$。如果 $x$ 为实数，该工具还会显示第一类和第二类的贝塞尔函数图。

为了快速精确地计算：

• 阶数 $\nu$ 必须在 $[-99, 99]$ 内。
• $x$ 的实部必须在 $[-20, 20]$ 内（虚部不限）。
• 超出这些范围，计算可能会变慢或精度降低；对于这种情况，建议使用专门的算法。

此贝塞尔函数计算器提供了一种便捷方式，可在线计算 $J_\nu(x)$ 和 $Y_\nu(x)$ 值，并附带图表以便直观理解。