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])