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

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

用法:

scipy.special.ellipj(u, m, out=None) = <ufunc 'ellipj'>#

雅可比椭圆函数

计算介于 0 和 1 之间的参数 m 和实参数 u 的雅可比椭圆函数。

参数 ：：

m： array_like

范围。

u： array_like

争论。

out： ndarray 的元组，可选

函数值的可选输出数组

返回 ：：

sn, cn, dn, ph： 标量或 ndarray 的 4 元组

返回的函数：

sn(u|m), cn(u|m), dn(u|m)

值酸碱度是这样的，如果u = 椭圆 (ph, m)， 然后sn(u|m) = sin(ph)和cn(u|m) = cos(ph).

注意：

Cephes 的包装器[1]常规埃尔普.

这些函数是周期性的，实轴上的quarter-period 等于完整的椭圆积分 ellipk(m)。

与不完全椭圆积分的关系：如果 u = ellipkinc(phi,m)，则 sn(u|m) = sin(phi)，且 cn(u|m) = cos(phi)。 phi 称为 u 的幅值。

计算是通过arithmetic-geometric 均值算法进行的，除非 m 在 0 或 1 的 1e-9 内。在后一种情况下 m 接近 1，近似值仅适用于 phi < pi/2。

参考：

[1]

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