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