估计两个伽马函数的比率,当各个函数 return inf 值
Estimate ratio of two gamma functions when individual functions return inf value
我正在估算两个伽马函数的比率。 gamma(x)
和 gamma(y)
的两个估计值都非常大 ( > 10^300) 但两者的比率应该相当小。
from scipy.special import gamma
gamma(x) / gamma(y)
不幸的是,有一点 gamma(x)
太大,scipy returns 是 inf
值。有没有办法阻止这种情况、提高阈值或计算这个比率?
谢谢
假设 x 和 y 是实数且为正,您可以使用以下恒等式:
a / b = e^(ln(a) - ln(b))
我建议使用 gmpy2
进行任意精度计算。它有 gmpy2.lgamma
,其中 returns 伽玛函数的对数,您可以使用 gmpy2.exp
将对数形式转换为您想要的比率。
我正在估算两个伽马函数的比率。 gamma(x)
和 gamma(y)
的两个估计值都非常大 ( > 10^300) 但两者的比率应该相当小。
from scipy.special import gamma
gamma(x) / gamma(y)
不幸的是,有一点 gamma(x)
太大,scipy returns 是 inf
值。有没有办法阻止这种情况、提高阈值或计算这个比率?
谢谢
假设 x 和 y 是实数且为正,您可以使用以下恒等式:
a / b = e^(ln(a) - ln(b))
我建议使用 gmpy2
进行任意精度计算。它有 gmpy2.lgamma
,其中 returns 伽玛函数的对数,您可以使用 gmpy2.exp
将对数形式转换为您想要的比率。