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

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

用法:

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

第一类 1 阶贝塞尔函数。

参数 ：：

x： array_like

参数(浮点数)。

out： ndarray，可选

函数值的可选输出数组

返回 ：：

J： 标量或 ndarray

第一类 1 阶贝塞尔函数在 x 处的值。

注意：

该域分为区间 [0, 8] 和 (8, 无穷大)。在第一个区间中，使用 24 项切比雪夫展开式。在第二个例子中，使用两个 5/5 次有理函数来采用渐近三角表示。

该函数是 Cephes [1] 例程 j1 的包装器。不应将其与球面贝塞尔函数混淆(请参阅 spherical_jn )。

参考：

[1]

Cephes 数学函数库，http://www.netlib.org/cephes/

例子：

计算某一点的函数：

>>> from scipy.special import j1
>>> j1(1.)
0.44005058574493355

在几个点计算函数：

>>> import numpy as np
>>> j1(np.array([-2., 0., 4.]))
array([-0.57672481,  0.        , -0.06604333])

绘制 -20 到 20 范围内的函数。

>>> import matplotlib.pyplot as plt
>>> fig, ax = plt.subplots()
>>> x = np.linspace(-20., 20., 1000)
>>> y = j1(x)
>>> ax.plot(x, y)
>>> plt.show()