[论文解读] Rapidly Convergent Summation Formulas involving Stirling Series

Raphael Schumacher|arXiv (Cornell University)|Jan 31, 2016

Advanced Mathematical Identities参考文献 6被引用 3

一句话总结

本文提出了一种新颖的级数加速技术，利用第一类斯特林数将发散的欧拉-麦克拉劳林渐近展开式转化为有限解析函数求和的快速收敛级数。该方法生成了大量新颖且计算效率极高的公式，适用于包含平方根、调和级数、对数、黎曼ζ函数及交错级数的求和，数值计算中实现了数量级的性能提升。

ABSTRACT

This paper presents a family of rapidly convergent summation formulas for various finite sums of analytic functions. These summation formulas are obtained by applying a series acceleration transformation involving Stirling numbers of the first kind to the asymptotic, but divergent, expressions for the corresponding sums coming from the Euler-Maclaurin summation formula. While it is well-known that the expressions obtained from the Euler-Maclaurin summation formula diverge, our summation formulas are all very rapidly convergent and thus computationally efficient.

研究动机与目标

• 开发适用于有限解析函数求和的快速收敛级数展开式，以克服标准欧拉-麦克拉劳林渐近展开式的发散性。
• 通过第一类斯特林数，将反幂级数与斯特林级数关联的温尼格变换应用于收敛加速。
• 推导出∑√k、∑1/k²、∑log k及交错级数等求和的新颖、计算高效的公式，其中许多为首次提出。
• 展示该方法在高精度特殊函数（如digamma函数）计算中相比现有方法的显著计算优势。

提出的方法

• 利用欧拉-麦克拉劳林求和公式生成有限求和的渐近展开式，其本身具有发散性但渐近准确。
• 应用温尼格变换，通过伯努利数与第一类斯特林数的系数，将这些发散的反幂级数转换为收敛的斯特林级数（反阶乘级数）。
• 采用Pochhammer符号（上升阶乘）构造所得级数的分母，确保阶乘衰减以实现快速收敛。
• 通过第一类斯特林数Sₖ⁽¹⁾(l)与伯努利数Bₗ₊₁的嵌套求和表达系数，实现系统的计算。
• 通过与已知公式（如斯特林、丰塔纳及格雷戈里公式）比较验证结果，并识别出此前未知的展开式作为新贡献。
• 通过基准测试展示实际应用价值：新方法在3秒内完成ψ(10¹⁰)的400位高精度计算，而Mathematica 10需超过10分钟。

实验结果

研究问题

• RQ1 能否系统性地利用斯特林级数将有限求和的发散欧拉-麦克拉劳林展开式转化为快速收敛级数？
• RQ2 第一类斯特林数与温尼格变换在加速渐近级数收敛性方面发挥何种作用？
• RQ3 所推导的求和公式中有多少是新的？其计算优势相较于现有方法如何？
• RQ4 该方法能否推广至包括对数、反幂及交错求和在内的广泛解析函数类？
• RQ5 该技术在提升高精度特殊函数（如digamma函数）计算效率方面能达到何种程度？

主要发现

• 该方法为∑ₖ₌₀ⁿ√k生成了快速收敛级数，推导出两个不同公式——一个以√n为首项，另一个以1/√n为首项，均显著快于渐近展开式。
• 对于∑ₖ₌₁ⁿ1/k²，论文重新获得了1730年已知的斯特林型公式，证实了该方法与历史结果的一致性。
• 推导出∑ₖ₌₁ⁿlog k的新公式，与斯特林公式一致，但以包含伯努利数与斯特林数的收敛斯特林级数形式表达。
• 该方法生成了交错调和级数∑ₖ₌₁ⁿ(−1)ᵏ⁺¹/k的新收敛公式，其系数由伯努利数与带符号的斯特林数导出。
• 该技术使ψ(10¹⁰)的400位精度计算在3秒内完成，相较Mathematica 10原生实现提速10倍。
• 论文在16个求和公式中识别出13个新公式，仅有3个为已有，表明该方法具有高度新颖性与广泛适用性。