3

一个简单的解析函数 f(s,Г) 怎么能近似逆不完全伽马函数Г(s,x)?这意味着写 x = f(s,Г) = 12*log(123.45*Г) + Г + 123.4^s 之类的东西。

(我至少需要想法或参考资料。)

4

3 回答 3

5

您可以查看 Boost 中的代码:http: //www.boost.org/doc/libs/1_35_0/libs/math/doc/sf_and_dist/html/math_toolkit/special/sf_gamma/igamma.html看看它们是什么使用。

编辑:他们也有逆: http: //www.boost.org/doc/libs/1_35_0/libs/math/doc/sf_and_dist/html/math_toolkit/special/sf_gamma/igamma_inv.html

于 2011-07-15T13:20:25.573 回答
1

我发现给定 s 的 x = f(s,Г) 可以很好地近似为 x = p0*(1-Г)^p1*ln(Г*p2)。至少它在 0.001 < Г < 0.999 区域内 s <= 15 对我有用。

这里 p0,p1,p2 - 是常数,在您选择 s 后通过 f(s,Г) 的近似值来选择。

于 2011-07-20T07:20:55.027 回答
0

Cephes有一个很好的实现。还有一个D 翻译,我认为它修复了 Cephes 版本中的一些错误。

于 2011-07-15T13:52:25.430 回答