n次元球の体積・表面積をExcelで計算してみる(GAMMALN関数)
数理科学のある号を読んでいると、学生さんにn次元球の体積にでてくる係数のガンマ関数の物理的意味は?と質問されたという話が書いてあった。でn次元球ってどうやって計算するかを思い出した。
n次元球の体積をVnとすると、Vn=Cn*r^nとかけるはず。表面積はこれを微分したもので、
Sn = n*Cn*r^(n-1).
ここで唐突だが
∫exp(-(x1^2 + x2^2+...+xn^2)) dx1dx2...dxn
を2通りの方法で計算する。(積分範囲は-∞~∞)
一つは、∫exp(-x1^2)dx1 * ∫exp(-x2^2)dx2 *...*∫exp(-xn^2)dxn =(√π)^n
もう一つは、r^2 = x1^2 + x2^2+...+xn^2 と置き、表面積にdrをかけた球殻の体積で積分するような形で計算すると、
∫exp(-r^2) *Sn dr =∫exp(-r^2) * n*Cn*r^(n-1)*drだが、r^2=s(ds=2rdr)と置くとn*Cn*∫s^(n/2-1) exp(-s)ds/2だが、Γ関数の定義から、これはCn = (n/2)*Γ(n/2) = Γ(n/2+1)となる。
なので、Vn = π^(n/2) * r^n / Γ(n/2+1)
Sn = n*π^(n/2) * r^(n-1) /Γ(n/2+1) = 2*π^(n/2) * r^(n-1) /Γ(n/2)
となる。で、これを計算するのに、Excelを使おう。GAMMALNという関数が最初から組み込まれていて、ガンマ関数の対数が簡単に計算できる。
それぞれの係数をプロットしたのがこちら。
これで行くと体積はn=5で最大、表面積はn=7で最大になる。nを整数に限らなければ(フラクタルが出てきてから非整数次元に抵抗がない、、、)
体積は n= 5.256946282, 表面積はn=7.256946337くらいで最大に。
複素数にしたらどうなるんだろうと思ったがそれは別の機会に。
« 大昔、ミュージックフェアで高中正義さんがTropic birdsを演奏していたその映像がYoutubeに! | トップページ | キンタロー。さんとももクロのかなこちゃんが並んだ写真が遠近法の概念を覆す。 »
「学問・資格」カテゴリの記事
- 「だまし絵と線形代数」を買った。あの錯視で超有名な杉原厚吉さんの著書。(2018.04.23)
- #仮面ライダービルド 第32話の話数を表す数式は、結晶族点群が32種類あるということ!(2018.04.22)
- 信頼性(寿命)評価に使われる”ワイブル分布のパラメータの算出”をカシオの高精度計算サイトに自作式として作った。(2018.04.20)
- ポアソン分布の例で馬に蹴られて死んだ兵士の例がよく出るが、原典見て実際に計算してみた。(2018.04.19)
- #仮面ライダービルド 第31話の話数を表す数式はメルセンヌ素数!M5=2^5-1=31(2018.04.15)
トラックバック
この記事のトラックバックURL:
http://app.cocolog-nifty.com/t/trackback/512682/57827858
この記事へのトラックバック一覧です: n次元球の体積・表面積をExcelで計算してみる(GAMMALN関数):
« 大昔、ミュージックフェアで高中正義さんがTropic birdsを演奏していたその映像がYoutubeに! | トップページ | キンタロー。さんとももクロのかなこちゃんが並んだ写真が遠近法の概念を覆す。 »
コメント