我正在估计两个伽玛函数的比率。两者的估计gamma(x) and gamma(y)非常大（> 10^300），但两者的比例应该相当小。
from scipy.special import gamma
gamma(x) / gamma(y)
不幸的是，有一个点gamma(x)太大，scipy 返回一个inf价值。有没有办法阻止这种情况，增加阈值，或者计算这个比率？
Thanks
假设 x 和 y 是实数且为正数，您可以使用以下恒等式：
a / b = e^(ln(a) - ln(b))
我可以建议使用gmpy2用于任意精度计算。它有gmpy2.lgamma，它返回 gamma 函数的对数，您可以使用gmpy2.exp从对数形式转换为您想要的比率。
本文内容由网友自发贡献，版权归原作者所有，本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容，请联系:hwhale#tublm.com(使用前将#替换为@)
