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

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

用法:

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

0阶修正贝塞尔函数。

定义为，

\[I_0(x) = \sum_{k=0}^\infty \frac{(x^2/4)^k}{(k!)^2} = J_0(\imath x),\]

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

参数 ：：

x： array_like

参数(浮点数)

out： ndarray，可选

函数值的可选输出数组

返回 ：：

I： 标量或 ndarray

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

注意：

该范围分为两个区间 [0, 8] 和 (8, 无穷大)。每个区间均采用切比雪夫多项式展开式。

该函数是 Cephes [1] 例程 i0 的包装器。

参考：

[1]

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

例子：

计算某一点的函数：

>>> from scipy.special import i0
>>> i0(1.)
1.2660658777520082

分几个点计算：

>>> import numpy as np
>>> i0(np.array([-2., 0., 3.5]))
array([2.2795853 , 1.        , 7.37820343])

绘制 -10 到 10 之间的函数。

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