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

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

用法:

scipy.special.fresnel(z, out=None) = <ufunc 'fresnel'>#

菲涅耳积分。

菲涅耳积分定义为

\[\begin{split}S(z) &= \int_0^z \sin(\pi t^2 /2) dt \\ C(z) &= \int_0^z \cos(\pi t^2 /2) dt.\end{split}\]

有关详细信息，请参阅 [dlmf]。

参数 ：：

z： array_like

实值或复值参数

out： ndarray 的 2 元组，可选

函数结果的可选输出数组

返回 ：：

S, C： 标量或 ndarray 的 2 元组

菲涅耳积分的值

参考：

[dlmf]

NIST 数学函数数字 Library https://dlmf.nist.gov/7.2#iii

例子：

>>> import numpy as np
>>> import scipy.special as sc

随着 z 沿实轴趋于无穷大，S 和 C 收敛到 0.5。

>>> S, C = sc.fresnel([0.1, 1, 10, 100, np.inf])
>>> S
array([0.00052359, 0.43825915, 0.46816998, 0.4968169 , 0.5       ])
>>> C
array([0.09999753, 0.7798934 , 0.49989869, 0.4999999 , 0.5       ])

它们与错误函数 erf 相关。

>>> z = np.array([1, 2, 3, 4])
>>> zeta = 0.5 * np.sqrt(np.pi) * (1 - 1j) * z
>>> S, C = sc.fresnel(z)
>>> C + 1j*S
array([0.7798934 +0.43825915j, 0.48825341+0.34341568j,
       0.60572079+0.496313j  , 0.49842603+0.42051575j])
>>> 0.5 * (1 + 1j) * sc.erf(zeta)
array([0.7798934 +0.43825915j, 0.48825341+0.34341568j,
       0.60572079+0.496313j  , 0.49842603+0.42051575j])