什么是贝塞尔 Y 函数数值表计算器？
这个工具用于列出第二类贝塞尔函数的数值表，该函数又称韦伯函数（Weber）或诺依曼函数（Neumann），记作 \(Y_{\nu}(x)\)。它是贝塞尔微分方程的第二个线性无关解。对于给定的实数阶数 \(\nu\)，计算器会在由起始值、步长和点数确定的一系列 \(x\) 取值上求出 \(Y_{\nu}(x)\)，从而生成一张完整的数值表。
使用方法
依次输入阶数 \(\nu\)（可以是非整数或负数）、\(x\) 的起始值、相邻点之间的步长（增量），以及迭代次数（即行数）。计算器按 $$x_i = \text{startX} + i \cdot \text{stepX}$$ （\(i\) 从 0 到 \(\text{pointCount} - 1\)）生成各个 \(x\) 值，并逐一列出对应的 \(Y_{\nu}(x)\)。需要注意的是，\(Y_{\nu}(x)\) 在 \(x = 0\) 处发散至负无穷，且只有当 \(x > 0\) 时才取实数值，因此任何 \(x \le 0\) 的行都会被标记为未定义。
计算公式
对于非整数阶： $$Y_{\nu}(x) = \frac{J_{\nu}(x)\cos(\nu\pi) - J_{-\nu}(x)}{\sin(\nu\pi)}$$ 对于整数阶 \(n\)，取极限可得到一个闭式表达式，其中包含一个含 \(J_n(x)\cdot\ln(x/2)\) 的对数项、一个有限的幂级数修正项，以及一个双伽马（digamma）级数。第一类贝塞尔函数 \(J_{\nu}(x)\) 由其幂级数求和得到，其中的伽马函数采用兰乔斯（Lanczos）近似计算。
实例演示
取 \(\nu = 0\)、\(\text{startX} = 0\)、\(\text{stepX} = 0.2\)、\(\text{pointCount} = 51\)，则各行的 \(x\) 从 \(0.0\) 一直到 \(10.0\)。其中 \(Y_0(0)\) 未定义（\(-\infty\)），\(Y_0(0.2) \approx -1.0811\)，\(Y_0(1.0) \approx 0.0883\)，\(Y_0(2.0) \approx 0.5104\)，\(Y_0(10.0) \approx 0.0557\)。结果中作为标题显示的"第一个有限值"即为 \(-1.0811\)。
定义与术语表
- 阶数 \(\nu\)
-
参数（
order字段）用于索引贝塞尔函数族。它可以是任意实数。整数阶（0、1、2、...）在具有圆柱对称性的物理问题中最常见；半整数阶给出球贝塞尔函数。 - 第二类贝塞尔函数 \(Y_\nu(x)\)
- 也称为韦伯函数或诺依曼函数（有时写作 \(N_\nu\)）。它是贝塞尔方程的解，在原点处无界（奇异）。对于非整数 \(\nu\) 定义为 \(Y_\nu(x) = \dfrac{J_\nu(x)\cos(\nu\pi) - J_{-\nu}(x)}{\sin(\nu\pi)}\)，整数情况通过极限得到。
- \(J_\nu\) 与 \(Y_\nu\) 的对比
- \(J_\nu(x)\)（第一类）在 \(x=0\) 处有限；\(Y_\nu(x)\)（第二类）当 \(x\to 0^+\) 时趋向于 \(-\infty\)。它们一起构成贝塞尔方程的完整独立解对。
- 贝塞尔微分方程
- 线性常微分方程 \(x^2 y'' + x y' + (x^2 - \nu^2) y = 0\)。其通解为 \(y = c_1 J_\nu(x) + c_2 Y_\nu(x)\)。
- 伽玛函数 \(\Gamma(z)\)
- 阶乘的连续扩张，\(\Gamma(n+1) = n!\)，出现在 \(J_\nu\) 和 \(Y_\nu\) 的级数系数中。
- digamma 函数 \(\psi(z)\)
- 对数导数 \(\psi(z) = \Gamma'(z)/\Gamma(z)\)。它明确出现在整数阶 \(Y_n(x)\) 的级数中，该级数包含对数项 \(\tfrac{2}{\pi}\ln(x/2)J_n(x)\) 加上 digamma 加权系数。
- Lanczos 逼近
- 用于计算伽玛函数 \(\Gamma(z)\) 对于复数或实数参数的高精度数值方法，通常在贝塞尔函数例程内使用以计算级数系数。
- 线性无关解
- 不能表示为第一个解的常数倍的第二个解。因为 \(J_\nu\) 本身无法表示在原点处奇异的解，\(Y_\nu\) 提供了通解所需的独立伴随函数。
常见问题
为什么第一行是未定义的？\(Y_{\nu}(x)\) 在 \(x = 0\) 处存在奇点，会发散到 \(-\infty\)，因此在该点没有有限值。
阶数可以取负数吗？可以。对于负整数阶，适用对称关系 \(Y_{-n}(x) = (-1)^n Y_n(x)\)；对于负的非整数阶，则直接使用上面的通用公式。
计算精度如何？级数会一直累加，直到各项小于机器精度为止，因此在中等大小的 \(x\) 范围内大约可达到 6—7 位有效数字。
