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

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

用法:

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

第二類1階貝塞爾函數。

參數 ：：

x： array_like

參數(浮點數)。

out： ndarray，可選

函數結果的可選輸出數組

返回 ：：

Y： 標量或 ndarray

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

注意：

該域分為區間 [0, 8] 和 (8, 無窮大)。在第一個區間中，使用 25 項切比雪夫展開式，並且需要計算\(J_1\)(第一類貝塞爾函數)。在第二個例子中，使用兩個 5/5 次有理函數來采用漸近三角表示。

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

參考：

[1]

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

例子：

計算某一點的函數：

>>> from scipy.special import y1
>>> y1(1.)
-0.7812128213002888

分幾個點計算：

>>> import numpy as np
>>> y1(np.array([0.5, 2., 3.]))
array([-1.47147239, -0.10703243,  0.32467442])

繪製從 0 到 10 的函數。

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