スターリングの近似

thumb|right|250px|log n! と n log n − n は n → ∞ のとき漸近するスターリングの近似（）またはスターリングの公式（）は、階乗、あるいはその拡張の一つであるガンマ関数の漸近近似である。名称は数学者にちなむ。

概要スターリングの近似は精度に応じていくつかの形がある。応用上よく使われる形の公式は、ランダウの記号を用いて、

である。における次の項はである。故に、次によい近似のは

である伏見康治「確率論及統計論」第I章　数学的補助手段 3節漸近展開 p.24 ISBN 9784874720127 http://ebsa.ism.ac.jp/ebooks/ebook/204。（ここで記号は両辺の比が（n → ∞ のとき） 1 に収束することを意味する。）

の漸近近似よりもむしろ上下からの評価が必要なことがある。そのような評価として、任意の正の整数に対して

が成り立ち、従って任意のに対して比はとの間にある。

スターリングの近似は階乗の複素引数への拡張の一つであるガンマ関数（正の整数に対しが成り立つ；ボーア・モレルップの定理も参照）に拡張することができ、

が成り立つ（ただし）。のときは収束が遅くなるため、応用上は相補公式などを用いて程度に制限することが多い。

導出初等的な導出スターリングの公式の厳密な証明にはオイラーの和公式、あるいは鞍点法といった複素解析の技法などを用いられることが多いが、初等的に導くことも可能である。まず階乗の対数を積分で近似する。logが凹関数であることから k-1<x<k (k=2,3,...) に対して

これを k-1 から k まで積分して

k=m+1,m+2,...,n に対して足し合わせると

ここでと定めると

m,n→∞のとき最左辺は1に収束するから、特に n=2m のとき

これとウォリスの公式の系と比較すると、

を得る。

精度の改善精度を改善するために an を評価する。

従って

オイラーの和公式による導出オイラーの乗積表示によるガンマ関数の定義の対数をとり

にオイラーの和公式を適用すれば

となる。右辺の定数を集めて

とすれば

となり、この主要部をガンマ関数の相補公式に代入してとすれば

となるが

であるから

を得る。剰余項については

として

である。故に

を得る。最初の数項を書き下せば

とやり、指数関数のテイラー展開により

となる。

鞍点法による導出スターリングの公式はの好適例とされることが多いが、実際に複素平面全体（負の実数を除く）で漸近近似が成立することを鞍点法によって示すのは困難であるから、ここではを正の実数に限定する。ガンマ関数はの置換により

となるが、が十分に大きければの付近が支配的であるから

という近似が許され、ガウス積分により

を得る。として、近似の誤差は

であり

であるから

を得る。これは

を示すに十分である。ただし、実際の誤差はであるが、それを鞍点法で示すのは困難である。

収束の速度と誤差見積もりより正確に記すと、次のようになるを見よ。も参照。。

ここで

スターリングの公式は以下の級数（スターリング級数）の近似（初項で打ち切ったもの）である。

としたとき、省かれた級数はその最初の項とそれ以降が相殺するように漸近していく。これは漸近展開の一例である。

以下のような階乗の対数の漸近展開も「スターリング級数」と呼ぶ。

この場合、誤差は打ち切った級数の初項と同じ符号で同程度の大きさであることが知られている。

ガンマ関数に対するスターリングの公式すべての正の整数に対して、

が成り立つ。ここで Γ はガンマ関数を表す。

しかしながら、は、階乗とは異なり、より広く、正でない整数を除いてすべての複素数に対して定義される。それにもかかわらず、スターリングの公式をなお適用することができる。Re z > 0 であれば

が成り立つ。部分積分を繰り返すことで次が得られる

ここで B は n 番目のベルヌーイ数である。（無限和は収束しないので、この公式は漸近展開にすぎないことに注意する。）公式はεを正数として |arg z| < − ε であるときに絶対値の十分大きい z に対して成り立つ。公式の右辺に現れる級数はスターリング級数と呼ばれる。最初の m 項が使われるとき誤差項はである。対応する近似は

のように書ける。この漸近展開のより進んだ応用は Re z が定数の複素変数 z に対してである。例えば直線 + it 上での Im z において適用されたスターリングの公式を見よ。

ビネーの公式スターリングの公式は収束しない級数を伴うので解析的に扱いづらいが、収束しない級数を収束する積分に換えたものとしてビネーの(第二)公式がある。

ビネーの公式は、スターリングの級数を形式的に(収束条件を無視して)操作することによっても導かれるが、厳密には対数ガンマ関数の導関数にアーベル・プラナの和公式を適用して得られる。

ならは右半平面において正則であるからプラナの和公式により

積分して

ならは有界であるから

である。スターリングの公式と比較して積分定数を求め

真数に直して

を得る。なお、ビネーの公式を元にして部分積分を繰り返すとスターリングの級数が得られる。

収束級数形式のスターリングの公式トーマス・ベイズの John Canton への書簡が1763年に王立協会により公表されている。それによると、スターリングの公式は収束級数ではないとされていたA letter from the late reverend Mr. Thomas Bayes, F. R. S. to John Canton, M. A. and F. R. S.。

スターリングの公式の収束級数形式を得るには以下を評価する。

一つの方法として、階乗冪の逆数の収束級数を使う方法がある。としたとき、次のようになる。

ここで

である。以上から次のようなスターリング級数が得られる。

これは、のとき収束する。

計算機向けの変形ガンマ関数の（関数電卓などの）計算機向けの近似として次の式がある。

これは、次と同等である。

これらはスターリングの公式を組み替えて、その結果生じる冪級数と双曲線正弦関数のテイラー展開の間の合致を観察することで得られる。この近似は z の実数部が 8 以上のとき、小数点以下 8 桁を超える精度を持つ。2002年、Robert H. Windschitl がリソースの制限された計算機（電卓など）でのそれなりの正確性を持った近似としてこれを示したToth, V. T. Programmable Calculators: Calculators and the Gamma Function. http://www.rskey.org/gamma.htm, modified 2006。

Gergő Nemes は 2007年にほぼ同程度の結果を与える近似式を提案した。こちらはより単純である。

これは、次と同等である。

歴史この公式は最初に次の形でアブラーム・ド・モアブルにより1730年に発見された。

スターリングの貢献は定数がであることを示したことである。より正確な形式はジャック・ビネが見出した。

スターリングの近似の「一次」バージョンは、マックス・プランクが1901年の黒体放射の論文で使用した。これは多量の光子や振動子についての黒体放射エネルギーの方程式にリンクしている。この近似は量子論でよく使われ、例えばピーター・デバイとルイ・ド・ブロイも使っている。アルベルト・アインシュタインとサティエンドラ・ボースは違う方式を採用した。非常に大きな n について確率分布をグラフに描画してみると、両者はほぼ平行になる。

脚注

参考文献

Paris, R. B., and Kaminsky, D., Asymptotics and the Mellin-Barnes Integrals, Cambridge University Press, 2001

スターリングの近似

Źródło i licencja