scipy.special.i1#

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

一阶修正贝塞尔函数。

定义为:

\[I_1(x) = \frac{1}{2}x \sum_{k=0}^\infty \frac{(x^2/4)^k}{k! (k + 1)!} = -\imath J_1(\imath x),\]

其中 \(J_1\) 是一阶第一类贝塞尔函数。

参数:
xarray_like

参数 (浮点数)

outndarray, optional

可选的输出数组,用于存储函数值

返回:
I标量或 ndarray

x 处的一阶修正贝塞尔函数的值。

另请参阅

iv

第一类修正贝塞尔函数

i1e

指数缩放的一阶修正贝塞尔函数

附注

范围被划分为两个区间 [0, 8] 和 (8, 无穷大)。在每个区间中使用切比雪夫多项式展开。

此函数是 Cephes [1] 例程 i1 的包装器。

数组 API 标准支持

i1 对 Python Array API 标准兼容的后端具有实验性支持,除了 NumPy 之外。请考虑通过设置环境变量 SCIPY_ARRAY_API=1 并提供 CuPy、PyTorch、JAX 或 Dask 数组作为数组参数来测试这些功能。支持以下后端和设备(或其他功能)的组合。

CPU

GPU

NumPy

不适用

CuPy

不适用

PyTorch

JAX

Dask

不适用

有关更多信息,请参阅 对数组 API 标准的支持

参考文献

[1]

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

示例

在一点处计算函数

>>> from scipy.special import i1
>>> i1(1.)
0.5651591039924851

在多个点计算函数

>>> import numpy as np
>>> i1(np.array([-2., 0., 6.]))
array([-1.59063685,  0.        , 61.34193678])

绘制函数在 -10 到 10 之间的图像。

>>> import matplotlib.pyplot as plt
>>> fig, ax = plt.subplots()
>>> x = np.linspace(-10., 10., 1000)
>>> y = i1(x)
>>> ax.plot(x, y)
>>> plt.show()
../../_images/scipy-special-i1-1.png