scipy.special.kv#
- scipy.special.kv(v, z, out=None) = <ufunc 'kv'>#
实数阶 v 的第二类修正贝塞尔函数
返回复数 z 处实数阶 v 的第二类修正贝塞尔函数。
这些有时也被称为第三类函数、巴塞特函数或麦克唐纳函数。它们被定义为修正贝塞尔方程的解,对于这些解,
\[K_v(x) \sim \sqrt{\pi/(2x)} \exp(-x)\]当 \(x \to \infty\) [3] 时。
- 参数:
- vfloat 类型的类数组
贝塞尔函数的阶数
- z复数类型的类数组
计算贝塞尔函数的参数
- outndarray,可选
函数结果的可选输出数组
- 返回:
- 标量或 ndarray
结果。请注意,输入必须是复数类型才能获得复数输出,例如,
kv(3, -2+0j)
而不是kv(3, -2)
。
注释
AMOS [1] 例程 zbesk 的包装器。有关所用算法的讨论,请参阅 [2] 及其中的参考文献。
参考文献
[1]Donald E. Amos,“AMOS,用于复数参数和非负阶的贝塞尔函数的便携式软件包”,http://netlib.org/amos/
[2]Donald E. Amos,“算法 644:用于复数参数和非负阶的贝塞尔函数的便携式软件包”,ACM TOMS 第 12 卷第 3 期,1986 年 9 月,第 265 页
[3]NIST 数字数学函数库,公式 10.25.E3。https://dlmf.nist.gov/10.25.E3
示例
绘制实数输入的多个阶数的函数
>>> import numpy as np >>> from scipy.special import kv >>> import matplotlib.pyplot as plt >>> x = np.linspace(0, 5, 1000) >>> for N in np.linspace(0, 6, 5): ... plt.plot(x, kv(N, x), label='$K_{{{}}}(x)$'.format(N)) >>> plt.ylim(0, 10) >>> plt.legend() >>> plt.title(r'Modified Bessel function of the second kind $K_\nu(x)$') >>> plt.show()
计算多个阶数的单个值
>>> kv([4, 4.5, 5], 1+2j) array([ 0.1992+2.3892j, 2.3493+3.6j , 7.2827+3.8104j])