アフィリエイト広告を利用しています

2015年10月23日

Excel でスターリングの近似式を計算してみます

 ガンマ関数の表示形式は、このブログで扱っている積分形の他にもいくつか存在しています。ガンマ関数はある種の積分を機械的に実行できるという利点があるのですが、ガンマ関数それ自身の性質を調べる場合、どの表示形式を選んでもなかなか扱いにくい面があります。 Excel には予め GAMMA() という形でガンマ関数が用意されているのですが、 Excel など存在しなかった時代にはプログラムを組んで積分や無限乗積などをコンピューターに計算させなければなりませんでした。私の学生時代にはすでに Excel はありましたが、当時の私は全くその潜在能力に注目していませんでした(おそらく世間でも今ほど広く認知されていませんでした)。パソコンの隅で眠らせていただけです。実に悔やまれます。

Excel でスターリングの近似式を計算してみます

 ともかく、そういう時代には「多少は精度を犠牲にしてもいいから、扱いやすい近似式がほしいよ」ということになります。そこで登場するのが スターリングの近似式 で、次のような形で表されます:

スターリングの近似式

 この近似式と本物のΓ(x + 1) を同時にグラフで描いてみると・・・・・・

 ガンマ関数とスターリング近似式.gif

 赤で描かれた曲線が正確なΓ(x + 1) , 青がスターリングの近似式です。 x が大きくなるほど精度が良くなっていく(グラフが重なる)ことがわかりますね。スターリングの近似式は漸近展開と呼ばれる手法で導かれた式なのですが、一般に漸近展開はそれほど良い近似を与えません。その中でスターリングの近似式は優等生といっていいぐらいの精度があります。
 しかし現代の数値計算では、この見事な近似式も時代の遺物といってもいいかもしれません。極めて正確なガンマ関数がソフトに組み込まれている以上、もはやスターリングの公式を使う理由はないはずです。

 実用上はともかくとして、特殊な関数として眺めてみると、かなり面白い式です。 √x x x の部分で急激な x の増加を実現させ、ガンマ関数と比較して行き過ぎた部分を e−x で少し抑えておくという形になっています。変形できそうな部分も色々あり、さらに面白そうな関数を作れる予感があります。あれこれと試行錯誤して色々なグラフを描いてみたのですが、 x と e それぞれの肩に乗っている x を cosx と sinx に置き換えた関数が一番印象に残りました。

 スターリングの式変形.gif

 一定周期で鋭いピークが並んでいますね。x の増加と共にピークの値は急速に伸びていきます。しかしどれほど値を伸ばしても、x の肩に乗る三角関数が負の値を取った瞬間、一気に 0 付近まで落ち込むことになります。
   
posted by Blog Cat at 14:21 | Comment(0) | TrackBack(0) | 特殊関数
検索
Excel VBA 数学教室
数学問題集(解答付き)
下剋上算数
ベクトル解析
サッカーマティクス
Excelで学ぶ統計解析
和算的思考力
学び直し
整数論の理論と演習
大人が手こずる算数
東大生の知恵袋
フーリエ変換
インド式秒算術
Excelで学ぶ微分積分
Excel 数学シミュレーション
オイラーの贈物


ファン
最新記事
カテゴリーアーカイブ