Python SciPy special.y0用法及代碼示例

本文簡要介紹 python 語言中 scipy.special.y0 的用法。

用法:

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

第二類0階貝塞爾函數。

參數 ：：

x： array_like

參數(浮點數)。

out： ndarray，可選

函數結果的可選輸出數組

返回 ：：

Y： 標量或 ndarray

第二類 0 階貝塞爾函數在 x 處的值。

注意：

該域分為區間 [0, 5] 和 (5, 無窮大)。在第一個區間中，采用有理近似\(R(x)\) 來計算，

\[Y_0(x) = R(x) + \frac{2 \log(x) J_0(x)}{\pi},\]

其中\(J_0\)是第一類0階貝塞爾函數。

在第二個區間中，漢克爾漸近展開式采用兩個 6/6 次和 7/7 次有理函數。

該函數是 Cephes [1] 例程 y0 的包裝器。

參考：

[1]

Cephes 數學函數庫，http://www.netlib.org/cephes/

例子：

計算某一點的函數：

>>> from scipy.special import y0
>>> y0(1.)
0.08825696421567697

分幾個點計算：

>>> import numpy as np
>>> y0(np.array([0.5, 2., 3.]))
array([-0.44451873,  0.51037567,  0.37685001])

繪製從 0 到 10 的函數。

>>> import matplotlib.pyplot as plt
>>> fig, ax = plt.subplots()
>>> x = np.linspace(0., 10., 1000)
>>> y = y0(x)
>>> ax.plot(x, y)
>>> plt.show()