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

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

用法:

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

伽马函数绝对值的对数。

定义为

\[\ln(\lvert\Gamma(x)\rvert)\]

其中\(\Gamma\) 是伽马函数。有关 gamma 函数的更多详细信息，请参阅 [dlmf]。

参数 ：：

x： array_like

真正的参数

out： ndarray，可选

函数结果的可选输出数组

返回 ：：

标量或 ndarray

伽马绝对值的对数值

注意：

它与 Python 标准库函数 math.lgamma 的函数相同。

当与 gammasgn 结合使用时，此函数对于在实轴上的日志空间中工作很有用，而无需通过关系 exp(gammaln(x)) = gammasgn(x) * gamma(x) 处理复数。

对于 complex-valued log-gamma，请使用 loggamma 而不是 gammaln 。

参考：

[dlmf]

NIST 数学函数数字 Library https://dlmf.nist.gov/5

例子：

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

它有两个正零。

>>> sc.gammaln([1, 2])
array([0., 0.])

它在非正整数处有极点。

>>> sc.gammaln([0, -1, -2, -3, -4])
array([inf, inf, inf, inf, inf])

它渐近地逼近x * log(x)(斯特林公式)。

>>> x = np.array([1e10, 1e20, 1e40, 1e80])
>>> sc.gammaln(x)
array([2.20258509e+11, 4.50517019e+21, 9.11034037e+41, 1.83206807e+82])
>>> x * np.log(x)
array([2.30258509e+11, 4.60517019e+21, 9.21034037e+41, 1.84206807e+82])