[논문 리뷰] Exactification of Stirling's Approximation for the Logarithm of the Gamma Function

Victor Kowalenko|arXiv (Cornell University)|2014. 04. 10.

Iterative Methods for Nonlinear Equations참고 문헌 22인용 수 3

한 줄 요약

이 논문은 발산하는 점근급수를 Borel 및 Mellin-Barnes (MB) 정규화를 통해 스타링어의 감마 함수 로그에 대한 정확한 근사법을 엄밀하게 정의한다. 두 정규화 기법이 잘라낸 점근급수 전개와 로그 항을 결합할 경우 동일한 매우 정밀한 값을 도출하며, 이는 수학 소프트웨어의 정밀도에 도달한다. MB 정규화는 유한한 합의 무한급수를 피하는 컴act한 적분 형태를 통해 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)와 같은 다가값 함수의 다가값 성질을 유지할 수 있음을 보여주어 주 분 支를 초월한 정의역을 확장하는 것.
• Poincaré 점근해석의 전통적 교리에 도전하여, 완전한 점근급수 전개의 미분이 유효하며 정확한 결과를 도출할 수 있음을 보여주는 것.

제안 방법

• 논문은 Borel 합산을 활용해 점근급수의 발산 잔여항을 정규화하며, 이는 수렴하는 지수 적분의 무한합과 스토크스 현상에서 기인하는 불연속적인 로그 항으로 표현된다.
• Mellin-Barnes (MB) 정규화를 도입하여 잔여항을 단일 MB 적분으로 표현함으로써 수렴 영역이 겹치는 것을 허용하고 무한합을 피한다.
• Cauchy의 잔여정리에 기반하여 발산 급수를 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 합산된 잔여항은 MB 정규화와 비교해 유사한 정밀도를 확보하기 위해 매우 큰 값(예: 10^5)에서 잘라내야 하므로 훨씬 느리고 효율성이 떨어진다.
• 정규화 과정은 무한대를 제거하고 유한부를 통해 발산 급수를 다스림으로써 해드암 유한부와 유사한 방식으로 정의되며, 이는 점근 해석에서 이러한 기법의 유효성을 입증한다.
• 이 연구는 완전한 점근급수 전개의 미분이 정규화가 적용된 상태에서 허용되며 정확한 결과를 도출할 수 있음을 보여주며, Poincaré 점근해석의 오랜 가정에 도전한다.