[论文解读] Exactification of Stirling's Approximation for the Logarithm of the Gamma Function

Victor Kowalenko|arXiv (Cornell University)|Apr 10, 2014

Iterative Methods for Nonlinear Equations参考文献 22被引用 3

一句话总结

本文提出了一种利用Borel和Mellin-Barnes（MB）正则化方法对伽马函数对数的斯特林近似进行严格精确化的技术。该方法表明，当与截断的渐近展开式及对数项结合时，两种正则化技术均能产生相同且高度精确的值，其精度与数学软件相当；其中MB正则化在计算效率上优于Borel求和，原因在于其紧凑的积分形式并避免了对指数积分无限级数的计算。

ABSTRACT

Exactification is the process of obtaining exact values of a function from its complete asymptotic expansion. Here Stirling's approximation for the logarithm of the gamma function or $\ln Γ(z)$ is derived completely whereby it is composed of the standard leading terms and an asymptotic series that is generally truncated. Nevertheless, to obtain values of $\ln Γ(z)$, the remainder must undergo regularization. Two regularization techniques are then applied: Borel summation and Mellin-Barnes (MB) regularization. The Borel-summed remainder possesses an infinite convergent sum of exponential integrals and discontinuous logarithmic terms across Stokes sectors and lines, while the MB-regularized remainders possess one MB integral, which is better to compute, and similar logarithmic terms. The MB integrals are valid over overlapping domains of convergence. Hence, two MB-regularized asymptotic forms can be used to evaluate $\ln Γ(z)$. Despite having to truncate the Borel-summed remainder, it is found that all the remainders combined with (1) the truncated asymptotic series, (2) the leading terms of Stirling's approximation and (3) their logarithmic terms yield identical values of $\ln Γ(z)$. In the few cases where the accuracy falls away, it is mostly due to a very high value for the truncation parameter, which results in the cancellation of redundant decimal places. Where possible, all the asymptotic forms yield the same values as the LogGamma routine in Mathematica.

研究动机与目标

• 开发一种从包含所有阶次以下主导项和超越级数的完整渐近展开式中获得伽马函数对数精确值的方法。
• 通过引入正则化技术，解决渐近展开式中发散余项的挑战，使其变为可求和且具有意义的表达式。
• 比较Borel求和与Mellin-Barnes正则化在计算伽马函数对数时的计算效率与数值稳定性。
• 证明正则化过程能够保留lnΓ(z)等多值函数的结构特征，特别是跨越不同黎曼面时的性质，从而扩展其适用范围。
• 挑战庞加莱渐近分析中的传统教条，表明对完整渐近展开式进行微分是有效的，并可得到精确结果。

提出的方法

• 本文采用Borel求和方法对渐近展开式的发散余项进行正则化，将其表示为指数积分的无限收敛级数，以及由斯托克斯现象引起的不连续对数项。
• 提出Mellin-Barnes（MB）正则化方法，通过单个MB积分表示余项，从而实现收敛域的重叠，并避免无限级数的出现。
• 该方法利用柯西留数定理，将发散级数表示为MB积分，舍弃发散的弧线围道积分，仅保留有限且正则化的值。
• 对Borel求和后的余项在极大值处（如10^5）进行截断，以实现高精度；而MB正则化形式则用黎曼ζ函数替代无限级数，提升计算紧凑性。
• 显式包含由斯托克斯不连续性产生的对数项——SD_M^{SL}(z)、SD_M^{SS}(z) 和 SD_MB(M,z)，以捕捉lnΓ(z)的多值特性，并实现对主支以外区域的扩展计算。
• 该方法将截断的渐近级数、主导斯特林项与正则化余项相结合，生成与高精度数学软件输出一致的结果。

实验结果

研究问题

• RQ1 斯特林近似展开式中lnΓ(z)的发散余项能否通过正则化方法转化为精确值？
• RQ2 在计算伽马函数对数时，Borel求和与Mellin-Barnes正则化在精度与计算效率方面如何比较？
• RQ3 正则化过程是否能保持lnΓ(z)的多值结构，特别是在不同黎曼面上的表现？
• RQ4 当正确正则化后，包含所有阶次以上项的完整渐近展开式是否可实现精确化？
• RQ5 即使渐近展开式本身发散，若施加正则化，是否仍可合法地对完整渐近展开式进行微分？

主要发现

• Borel求和与MB正则化余项在与截断渐近级数及主导斯特林项结合后，产生完全一致的数值结果，其精度与高精度数学软件包输出完全匹配。
• MB正则化形式在计算上更具优势，它用单个MB积分和黎曼ζ函数替代了指数积分的无限级数，实现了更快、更稳定的计算。
• 由斯托克斯不连续性产生的对数项对于捕捉lnΓ(z)的多值特性至关重要，使得在主支 −π < arg z ≤ π 之外的区域也能实现有效计算。
• Borel求和余项需在极大值（如10^5）处截断才能达到与MB正则化相当的精度，因此计算速度显著更慢，效率更低。
• 正则化过程通过有限部分（类似哈达玛有限部分）有效消除了发散级数中的无穷大，验证了此类方法在渐近分析中的有效性。
• 本研究证明，对完整渐近展开式进行微分是允许的，并能获得精确结果，从而挑战了庞加莱渐近分析中长期存在的假设。