[論文レビュー] More accurate approximations for the Gamma function

Gergő Nemes|arXiv (Cornell University)|Mar 31, 2010

Mathematical functions and polynomials参考文献 11被引用数 24

ひとこと要約

本稿では、ゴスパーの公式にインspiredされた系列変換と中央二項係数の漸近展開を用いて、ガンマ関数の新しい漸近近似を提案する。この近似は、シフトされた変数の偶数次の項のみを含む展開形を採用し、特に偶数次近似において顕著な精度向上を達成しており、スターリング、ラ플ラス、ラマヌジャン、モルティチの公式を上回る数値的性能を示している。

ABSTRACT

A series transformation idea inspired by a formula of R. W. Gosper and some asymptotic expansions for the central binomial coefficients leads us to new accurate approximations for the Gamma function.

研究の動機と目的

• 大正の引数に対してより高精度なガンマ関数の漸近近似を開発すること。
• スターリングの公式やラ플ラスの公式といった既存の近似手法の限界、特に収束が遅いか精度が低いという問題を是正すること。
• 中央二項係数の展開構造を活用し、シフトされた変数の偶数次の項のみを含む新しい系列形を導出すること。
• ゴスパーの近似を変換して、係数を制御可能な高次の漸近系列に発展させ、数値的効率と精度を向上させること。
• 漸近的に高精度な近似をもたらす係数列を体系的に計算する手法を提供すること。

提案手法

• スターリングの公式における標準的な x を (x + 1/6) に置き換えることで初期精度を向上させるゴスパーの近似に基づく系列変換を用いる。
• 中央二項係数の 1/(n + 1/4) のべき級数としての漸近展開を応用し、シフトされた変数の偶数次の項のみを含む新しい系列形を考案する。
• Γ(x+1) ~ x^x e^{-x} √(2π(x + 1/6)) × Σ g_n / (x + v_n)^{2n} の形の新しい漸近展開を導出する。ここで g_n と v_n は再帰的に決定される。
• スターリングの系列から得られる係数 a_n と新しい係数 g_n、v_n の間の関係を、二項展開と漸近的マッチングを用いて確立する。式 (2.6) に示される再帰関係式を導出する。
• 漸近的級数展開の一意性を用いて、導出された形が x → ∞ のとき真のガンマ関数の挙動と一致することを保証する。
• 複数の x 値について、近似誤差の正確な小数桁数を数値的に計算し、既存の公式と比較することで、手法の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 シフトされた変数 (x + 1/6) の偶数次の項のみを含む新しいガンマ関数の漸近展開を構築可能か？ これにより収束性と精度が向上するか？
• RQ2 ゴスパーの近似を、係数を制御可能な高次の漸近系列に体系的に発展させることは可能か？ これにより数値的性能が向上するか？
• RQ3 大 x における正確な小数桁数の観点から、本稿の新しい展開はスターリング、ラプラス、ラマヌジャン、モルティチの公式を上回るか？
• RQ4 新しい系列における係数 g_n と v_n を支配する再帰関係式は何か？ これによりガンマ関数と漸近的に同等の挙動を示せるか？
• RQ5 特に偶数次近似において、本稿の新しい系列はラマヌジャンの公式を上回る高い精度を達成できるか？

主な発見

• 新しい漸近展開 (2.5) は、スターリング、ラプラス、ラマヌジャン、モルティチの公式を上回り、特に偶数次近似において顕著な精度向上を示している。
• x = 100 の場合、8次近似で19.2桁の正確な小数桁数を達成しており、ラマヌジャンの19.5桁、モルティチの19.4桁を上回っている。
• x = 1000 の場合、8次近似で27.2桁の正確な小数桁数を達成しており、ラマヌジャンの27.5桁、モルティチの27.2桁を上回っている。
• 再帰関係 (2.6) は、係数 g_n と v_n を効果的に生成できており、n=5 で v_n ≈ 0.249958 に収束しており、安定化していることが示された。
• 偶数次項のみを含む特別な系列 (2.5) は、x=1000 で10次近似で38.5桁の正確な小数桁数を達成しており、優れた収束性を示している。
• 数値的結果から、本稿の公式は奇数次および偶数次近似の両方において、特に高次項において、既存手法を一貫して上回ることが明らかになった。