Python SciPy special.kn用法及代码示例

本文简要介绍 python 语言中 scipy.special.kn 的用法。

用法:

scipy.special.kn(n, x, out=None) = <ufunc 'kn'>#

第二类整数阶n的修正贝塞尔函数

返回整数阶的第二类修正贝塞尔函数n在真实z.

这些有时也称为第三类函数、巴塞特函数或麦克唐纳函数。

参数 ：：

n： int的数组

贝塞尔函数的顺序(浮点数将截断并发出警告)

x： 类似浮点数的数组

评估贝塞尔函数的参数

out： ndarray，可选

函数结果的可选输出数组。

返回 ：：

标量或 ndarray

第二类修正贝塞尔函数的值，\(K_n(x)\)。

注意：

AMOS 的包装器[1]常规兹别斯克.有关所使用算法的讨论，请参见[2]以及其中的引用。

参考：

[1]

Donald E. Amos，“AMOS，用于复杂参数和非负阶贝塞尔函数的便携式软件包”，http://netlib.org/amos/

[2]

Donald E. Amos，“算法 644：复杂参数和非负阶贝塞尔函数的便携包”，ACM TOMS Vol。 12 第 3 期，1986 年 9 月，第 3 页。 265

例子：

绘制实际输入的几个订单的函数：

>>> import numpy as np
>>> from scipy.special import kn
>>> import matplotlib.pyplot as plt
>>> x = np.linspace(0, 5, 1000)
>>> for N in range(6):
...     plt.plot(x, kn(N, x), label='$K_{}(x)$'.format(N))
>>> plt.ylim(0, 10)
>>> plt.legend()
>>> plt.title(r'Modified Bessel function of the second kind $K_n(x)$')
>>> plt.show()

计算多个订单的单个值：

>>> kn([4, 5, 6], 1)
array([   44.23241585,   360.9605896 ,  3653.83831186])