高效评估超几何函数

时间:2009-01-25 21:16:22

标签: algorithm math numerical

有没有人有使用算法评估超几何函数的经验?我会对一般参考文献感兴趣,但我会描述我的特殊问题,以防有人处理它。

我的具体问题是评估形式3F2(a,b,1; c,d; 1)的函数,其中a,b,c和d都是正实数,c + d> 1。 A + B + 1。有许多特殊情况具有封闭形式的公式,但据我所知,一般没有这样的公式。以零为中心的幂级数收敛于1,但非常缓慢;连续系数的比率在限制中变为1。也许Aitken加速会有所帮助吗?

2 个答案:

答案 0 :(得分:6)

我测试了Aitken加速度,它似乎没有帮助解决这个问题(Richardson也没有推断)。这可能意味着Pade近似也不起作用。我可能做错了,所以一定要亲自试试。

我可以想到两种方法。

一种是在某些点评估序列,例如z = 0.5,其中收敛快速得到初始值,然后通过将hypergeometric differential equation插入ODE求解器前进到z = 1。我不知道这在实践中有多好用;它可能不会,因为z = 1是一个奇点(如果我没记错的话)。

第二种是根据Meijer G-function使用3F2的定义。定义Meijer G函数的轮廓积分可以通过将高斯或双指数正交应用于轮廓的分段来数值地评估。这不是非常有效,但它应该有效,它应该扩展到相对较高的精度。

答案 1 :(得分:1)

你想要一个你知道连续术语比例的系列并且它是一个合理的函数,这是正确的吗?

我认为Gosper's algorithm以及用于证明hypergeometric identities(并找到它们)的其他工具正是这样做的,对吗? (见Wilf和Zielberger的A=B book online.