修正球贝塞尔函数 i_v(x)、k_v(x) 及其导数计算器

输入计算

阶 v（整数，无量纲）

x (argument, x > 0) For x ≤ 0 the real-valued results may be complex or undefined.

数学公式

Show calculation steps (2)

1. Upward recurrence to order v

   

   Applied for v = Order >= 1 to reach i_v and k_v; n runs from 1 up to v.

2. Derivatives

   

   First derivatives of the modified spherical Bessel functions at order v.

结果

Modified Spherical Bessel (first kind) i_v(x)

1.8134302039

无量纲

Second kind kv(x)	0.1062920829
Derivative i'v(x)	-0.9743827436
Derivative k'v(x)	-0.1594381243

这个计算器的用途

本工具用于计算第一类修正球贝塞尔函数 i_v(x) 与第二类 k_v(x)，以及它们的一阶导数 i'_v(x) 和 k'_v(x)，适用于非负整数阶 v 和正实数自变量 x。它们是纯粹的数学特殊函数——在任何地方的取值都完全一致，不涉及任何地区或单位上的约定。

原理与公式

这些函数是修正球贝塞尔方程 x²w″ + 2xw′ − (x² + v(v+1))w = 0 的解。它们通过半整数阶的平移与柱坐标下的修正贝塞尔函数相联系：i_v(x) = √(π/2x)·I_v+1/2(x)，k_v(x) = √(2/πx)·K_v+1/2(x)。由于这个 +1/2 的平移使整数 v 对应的阶变为半整数，因此函数可化简为含 sinh、cosh 与 exp 的初等表达式。我们以 i₀=sinh(x)/x、i₁=cosh(x)/x−sinh(x)/x²、k₀=(π/2x)e^−x、k₁=(π/2x)e^−x(1+1/x) 作为初始值，再逐阶向上递推到所需的阶数。导数则由 f'_v = −f_v+1 + (v/x)f_v 给出。

随着 x 增大，第一类函数 i_v(x) 增长，而第二类 k_v(x) 衰减。

使用方法

输入整数阶 v（0、1、2、…）以及满足 x > 0 的自变量 x，即可读取四个结果。请注意，这里采用的约定为 k_v(x)=√(2/πx)K_v+1/2(x)，它在 k₀ 中引入了 π/2 这一因子；部分文献会省略该因子。

计算示例（v = 0，x = 2）

i₀(2)=sinh(2)/2=3.6268604/2=1.8134302。i₁(2)=cosh(2)/2−sinh(2)/4=1.8810978−0.9067151=0.9743827，因此 i'₀(2)=−i₁(2)=−0.9743827。k₀(2)=(π/4)e⁻²=0.1062930，k₁(2)=k₀·1.5=0.1594394，因此 k'₀(2)=−k₁(2)=−0.1594394。

常见问题

可以使用非整数阶 v 吗？ 这段实值计算仅支持非负整数阶，此时函数为初等形式。非整数阶则需要完整地计算 I/K 贝塞尔函数。

为什么 x 必须为正？ 当 x→0 时 k_v(x) 会发散，而当 x<0 时结果可能变为复数，因此实值版本要求 x > 0。

i_v 与 k_v 有什么区别？ i_v 随 x 指数增长，并在原点处正则；k_v 随 x 指数衰减，并在原点处奇异。