这个计算器能做什么
本工具用于求解贝塞尔微分方程 \(x^2y'' + xy' + (x^2 - v^2)y = 0\) 的两个线性无关解：第一类贝塞尔函数 \(J_{v}(x)\)、第二类贝塞尔函数 \(Y_{v}(x)\)（又称诺依曼函数 Neumann function），以及它们的一阶导数 \(J'_{v}(x)\) 和 \(Y'_{v}(x)\)。它支持任意实数阶 \(v\)（整数、分数或负数）和任意实数自变量 \(x\)。贝塞尔函数在物理与工程领域随处可见——圆形薄膜振动、圆柱体内的热传导、波导中的电磁波传播，以及信号处理等场景都会用到它。
使用方法
填入阶 \(v\)（例如 0、1 或 0.5）、自变量 \(x\)，并选择需要显示的小数位数。点击「计算」即可得到全部四个结果。当阶 \(v\) 为非整数时，请保证 \(x \geq 0\)，因为此时 \((x/2)^{v}\) 在 \(x\) 为负值时会变成复数。在 \(x = 0\) 处，第二类函数存在奇点，会被标记为「未定义」。
公式说明
\(J_{v}(x)\) 由其幂级数借助伽马函数（gamma function）逐项求和得到，采用稳定的递推算法，直到各项小于容差为止。 $$J_{\nu}(x) = \sum_{k=0}^{\infty} \frac{(-1)^{k}}{k!\,\Gamma(\nu+k+1)} \left(\frac{x}{2}\right)^{2k+\nu}, \qquad \nu = \text{Order }\nu,\; x = \text{Argument }x$$ 对于非整数阶，\(Y_{v}(x)\) 通过公式 $$Y_{\nu}(x) = \frac{J_{\nu}(x)\cos(\nu\pi) - J_{-\nu}(x)}{\sin(\nu\pi)}, \qquad \nu = \text{Order }\nu,\; x = \text{Argument }x$$ 计算；对于整数阶，由于 \(\sin(v\pi)=0\) 会导致除零，计算器会对 \(v\) 施加一个极小的扰动（约 \(1\mathrm{e}{-7}\)）来规避。导数则使用递推关系 $$C'_{\nu}(x) = \tfrac{1}{2}\bigl(C_{\nu-1}(x) - C_{\nu+1}(x)\bigr), \quad C \in \{J,\,Y\}, \quad \nu = \text{Order }\nu,\; x = \text{Argument }x$$ 其中特例为 \(C'_{0}(x) = -C_{1}(x)\)。
实例演算
当 \(v = 0\)、\(x = 1\) 时：\(J_{0}(1)\) 的级数展开为 $$1 - 0.25 + 0.015625 - 0.000434 + \ldots \approx 0.7651977$$ 其余已知值为 \(Y_{0}(1) \approx 0.0882570\)，\(J'_{0}(1) = -J_{1}(1) \approx -0.4400506\)，以及 \(Y'_{0}(1) = -Y_{1}(1) \approx 0.7812128\)。
关键术语和变量
- 阶 \(\nu\)
- 参数 \(\nu\) 在 \(J_\nu(x)\) 和 \(Y_\nu(x)\) 中，它确定微分方程 \(x^2 y'' + x y' + (x^2-\nu^2)y = 0\) 的形式。它可以是任意实数；整数阶（\(\nu = 0,1,2,\dots\)）源自圆柱坐标中的角度分离，而半整数阶给出可用初等函数表示的球形贝塞尔函数。
- 自变量 \(x\)
- 函数被求值的独立变量，通常是缩放径向距离 \(x = kr\)。对于实数 \(x\)，当 \(\nu\) 为整数时 \(J_\nu\) 为实值，而 \(Y_\nu\) 仅在 \(x>0\) 时定义。
- 第一类贝塞尔函数 \(J_\nu(x)\)
- 在原点处有限的解（对于 \(\nu\ge 0\)），由级数 \(J_\nu(x)=\sum_{k=0}^{\infty}\frac{(-1)^k}{k!\,\Gamma(\nu+k+1)}\left(\tfrac{x}{2}\right)^{2k+\nu}\) 定义。当 \(x\) 增大时，它以缓慢衰减的幅度振荡。
- 第二类贝塞尔函数 \(Y_\nu(x)\)
- 也称为诺依曼函数或韦伯函数，这是第二个线性独立解。它通过 \(Y_\nu(x)=\dfrac{J_\nu(x)\cos(\nu\pi)-J_{-\nu}(x)}{\sin(\nu\pi)}\) 定义（对整数 \(\nu\) 有极限形式），且在原点处以对数形式或 \(x\) 的幂次发散。
- 导数 \(J'_\nu(x)\)、\(Y'_\nu(x)\)
- 关于 \(x\) 的导数。它们满足递推关系 \(C'_\nu(x)=\tfrac{1}{2}\bigl(C_{\nu-1}(x)-C_{\nu+1}(x)\bigr)\) 和 \(C'_\nu(x)=C_{\nu-1}(x)-\tfrac{\nu}{x}C_\nu(x)\)，其中 \(C\) 代表 \(J\) 或 \(Y\)。特别地，\(J'_0(x)=-J_1(x)\)。
- 伽马函数 \(\Gamma(z)\)
- 阶乘的连续扩展，对于非负整数有 \(\Gamma(n+1)=n!\)，在 \(J_\nu\) 级数的分母中出现以允许非整数阶。请参阅伽马函数计算器获取个别值。
- 零点（根）
- 使得 \(J_\nu(x)=0\) 或 \(Y_\nu(x)=0\) 的值 \(j_{\nu,m}\) 和 \(y_{\nu,m}\)。这些在边值问题中用作特征值；例如，固定边界圆形膜的频率与零点 \(j_{\nu,m}\) 成正比。
常见问题
阶 \(v\) 可以是负数或分数吗？可以。级数与伽马函数能够处理任意实数 \(v\)。只需在 \(v\) 为非整数时保证 \(x \geq 0\) 即可。
为什么 \(x = 0\) 时 Y 未定义？所有第二类贝塞尔函数在 \(x \to 0\) 时都会发散至负无穷，因此在该点不存在有限值。
计算精度如何？计算采用双精度浮点数（约 15 位有效数字）。「显示位数」选项仅影响结果的呈现格式，不会改变底层运算精度。
