超球の体積

初等幾何学における球体は決められた点から決められた距離以内にある点の全体が空間において占める領域であった。同様のことを テンプレート:Mvar-次元ユークリッド空間で行って テンプレート:Mvar-次元超球体が定義される。

本項ではテンプレート:Mvar-次元超球体の体積、および超球体の表面に相当する超球面の面積について述べる。

超球体の体積

テンプレート:Math theorem ここで$\left(\frac{n}{2}\right)!$はオイラーのガンマ函数テンプレート:Math（階乗函数の非整数引数への一般化）により定義される

$\left(\frac{n}{2}\right)!$$=\mathrm{\Gamma }(\frac{n}{2}+1)$$=\{\begin{array}{cc}(k+\frac{1}{2})((k-1)+\frac{1}{2})\cdots (1+\frac{1}{2})\sqrt{\pi }& \text{if }n=2k+1\\ k(k-1)\cdots 2\cdot 1& \text{if }n=2k\end{array}$

であり、「!!」は二重階乗

$n!!=\{\begin{array}{cc}n(n-2)(n-4)\cdots 3\cdot 1& \text{if }n\text{ is odd}\\ n(n-2)(n-4)\cdots 4\cdot 2& \text{if }n\text{ is even}\end{array}$

である。

具体的な値はテンプレート:Mvarの偶奇によって異なる。

最初のいくつかの次元

次元	半径 テンプレート:Mvar の球の体積	体積 テンプレート:Mvar の球の半径
0	$1$	全ての球の体積は1
1	$2R$	$V/2$
2	$\pi R^2$	$\frac{V^{1/2}}{\sqrt{\pi }}$
3	$\frac{4}{3}\pi R^3$	${\left(\frac{3V}{4\pi }\right)}^{1/3}$
4	$\frac{{\pi }^2}{2}{R}^4$	$\frac{(2V{)}^{1/4}}{\sqrt{\pi }}$
5	$\frac{8{\pi }^2}{15}{R}^5$	${\left(\frac{15V}{8{\pi }^2}\right)}^{1/5}$
6	$\frac{{\pi }^3}{6}{R}^6$	$\frac{(6V{)}^{1/6}}{\sqrt{\pi }}$
7	$\frac{16{\pi }^3}{105}{R}^7$	${\left(\frac{105V}{16{\pi }^3}\right)}^{1/7}$
8	$\frac{{\pi }^4}{24}{R}^8$	$\frac{(24V{)}^{1/8}}{\sqrt{\pi }}$
9	$\frac{32{\pi }^4}{945}{R}^9$	${\left(\frac{945V}{32{\pi }^4}\right)}^{1/9}$
10	$\frac{{\pi }^5}{120}{R}^{10}$	$\frac{(120V{)}^{1/10}}{\sqrt{\pi }}$
テンプレート:Mvar	テンプレート:Math	テンプレート:Math

奇数の場合

テンプレート:Mathの場合は、

$V_{2k+1}(R)$$=\frac{{\pi }^k\sqrt{\pi }}{(k+\frac{1}{2})!}{R}^{2k+1}$$=\frac{2(2\pi {)}^k}{(2k+1)!!}{R}^{2k+1}$$=\frac{2(2\pi {)}^k}{1\cdot 3\cdot \cdots \cdot (2k+1)}{R}^{2k+1}$

テンプレート:Mathである事を用いると、

$V_{2k+1}(R)$$=\frac{2(k!)(4\pi {)}^k}{(2k+1)!}{R}^{2k+1}$

とも表記できる。

偶数の場合

テンプレート:Mathの場合、

$V_{2k}(R)$$=\frac{{\pi }^k}{k!}{R}^{2k}$$=\frac{(2\pi {)}^k}{(2k)!!}{R}^{2k}$$=\frac{{\pi }^k}{1\cdot 2\cdot \cdots (k-1)\cdot k}{R}^{2k}$

漸近評価

前述のテンプレート:Mathの具体的な記述にスターリングの公式$$n!\sim \sqrt{2\pi n}{\left(\frac{n}{e}\right)}^n.$$とテンプレート:仮リンク

$x^{1-s}<\frac{\mathrm{\Gamma }(x+1)}{\mathrm{\Gamma }(x+s)}<(x+1{)}^{1-s}.$

を適用することで、以下が成立する事がわかる： テンプレート:Math theorem これは一辺の長さがテンプレート:Mvarの超立方体の体積テンプレート:Mvar対する超球体テンプレート:Mathがテンプレート:Mathのとき指数関数的に小さくなる事を意味する。

漸化式

超球体の体積を漸化式で表すと、下記のとおりである：

$V_n(R)=\{\begin{array}{cc}1& \text{if }n=0,\\ 2R& \text{if }n=1,\\ \frac{2\pi }{n}{R}^2\times V_{n-2}(R)& \text{otherwise}.\end{array}$

上述の漸化式はテンプレート:Mathをテンプレート:Mathと関係づけるが、テンプレート:Mathをテンプレート:Mathを関係づけると下記のようになる：

$V_n(R)=R\sqrt{\pi }\frac{(\frac{n-1}{2})!}{(\frac{n}{2})!}{V}_{n-1}(R)$$=a_{n}R\frac{(n-1)!!}{n!!}{V}_{n-1}(R)$

ここで、

$a_n=\{\begin{array}{cc}2& \text{if }n=2k+1\\ \pi & \text{if }n=2k\end{array}$

具体的に書けば、

$V_n(R)$$=\{\begin{array}{cc}2R\frac{(2k)(2k-2)\cdots 6\cdot 4\cdot 2}{(2k+1)(2k-1)\cdots 5\cdot 3\cdot 1}{V}_{2k}(R)& \text{if }n=2k+1\\ \pi R\frac{(2k-1)(2k-3)\cdots 5\cdot 3\cdot 1}{(2k)(2k-2)\cdots 6\cdot 4\cdot 2}{V}_{2k-1}(R)& \text{if }n=2k\end{array}$

体積から半径を求める公式

テンプレート:Mvar-次元球体の体積 テンプレート:Mvar をその半径 テンプレート:Mvar で表す代わりに、上記の公式を逆に解いて半径 テンプレート:Mvar を体積 テンプレート:Mvar の函数として表すこともできる：

$R_n(V)=\sqrt[n]{\frac{\left(\frac{n}{2}\right)!}{{\pi }^{n/2}}V}$$=\sqrt[n]{\frac{n!!}{2^{\lceil n/2\rceil }{\pi }^{\lfloor n/2\rfloor }V}}$$=\{\begin{array}{cc}\sqrt[2k+1]{\frac{(2k+1)!!}{2(2\pi {)}^k}V}& \text{if }n=2k+1\\ \sqrt[2k+1]{\frac{k!}{{\pi }^k}V}& \text{if }n=2k\end{array}$

超球面の面積

ユークリッド空間${ℝ}^{n+1}$における半径 テンプレート:Mvar の テンプレート:Mvar-次元超球面

$S^n(R)=\{(x_1,\dots ,x_n)\in {ℝ}^{n+1}|x_1{}^2+\cdots +x_{n+1}{}^2=R\}$

の面積テンプレート:Mathと${ℝ}^{n+1}$における半径 テンプレート:Mvar の テンプレート:Mvar-次元超球体

$B^{n+1}(R)=\{(x_1,\dots ,x_n)\in {ℝ}^n|x_1{}^2+\cdots +x_{n+1}{}^2\le R\}$

の体積テンプレート:Mathには

$A_n(R)=\frac{d}{dR}{V}_{n+1}(R)$$=(n+1)V_{n+1}(1)R^n$

という関係がある。これはテンプレート:Mathがテンプレート:Mathに対するテンプレート:Mathの和集合である事から従う。テンプレート:Mathの具体的表記から、以下が従う： テンプレート:Math theorem よってテンプレート:Mathの場合は、

$A_{2k+1}(R)$$=\frac{2{\pi }^k}{1\cdot 2\cdot \cdots (k-1)\cdot k}$$=\frac{2{\pi }^k}{k!}{R}^{2k+1}$

テンプレート:Mathの場合は、

$A_{2k}(R)$$=\frac{2(2\pi {)}^k}{1\cdot 3\cdot \cdots \cdot (2k-1)}{R}^{2k}$$=\frac{2(2\pi {)}^k}{(2k-1)!!}{R}^{2k}$$=\frac{2(k!)(4\pi {)}^k}{(2k-1)!}{R}^{2k}$

とも表記できる。

証明

上記の公式に関して多くの証明が存在する。

体積は半径の テンプレート:Mvar 乗に比例する

テンプレート:Mvar-次元球の体積についてのいくつかの証明においての重要なステップであり、それ以外にも有用性のある一般的な事実は、半径 テンプレート:Mvar の テンプレート:Mvar-次元球の体積は テンプレート:Mvar に比例すること、つまり $$V_n(R)\propto R^n$$ である。このときの比例定数は単位球の体積に等しい。

上記の関係は帰納法による簡単な証明がある。基底段階は テンプレート:Math であり、比例することは自明である。帰納段階は、次元 テンプレート:Math で比例することが真であると仮定する。テンプレート:Mvar-次元球体と一つの超平面との交わりは テンプレート:Math-次元球体であることに注意する。テンプレート:Mvar-次元球体の体積を テンプレート:Math-次元球体の体積の積分 $$V_n(R)=\underset{-R}{\overset{R}{\int }}{V}_{n-1}(\sqrt{R^2-x^2})dx$$ として書く時、帰納法の仮定により テンプレート:Math-次元球体の半径から テンプレート:Mvar-倍の因子を括りだして $$V_n(R)=R^{n-1}\underset{-R}{\overset{R}{\int }}{V}_{n-1}(\sqrt{1-(x/R{)}^2})dx$$ と書くことができる。変数変換 テンプレート:Math を施して導かれる $$V_n(R)=R^n\underset{-1}{\overset{1}{\int }}{V}_{n-1}(\sqrt{1-t^2})dt=R^n{V}_n(1)$$ は次元 テンプレート:Mvar における比例関係を示すものになっている。帰納法によって、全ての次元で比例関係は真である。

2次元漸化式

テンプレート:Mvar-次元球体と テンプレート:Math-次元球体の間の体積の漸化式の証明は、上記の比例式と円筒座標系における積分を用いて与えられる。球の中心を通る平面を固定する。テンプレート:Mvar を球面の中心と平面上の点との距離とし、テンプレート:Mvar を方位角とする。テンプレート:Mvar-次元球体と、半径と方位角を固定して定まる テンプレート:Math-次元平面とを交わらせれば、半径 テンプレート:Math の テンプレート:Math-次元球体が与えられる。球の体積は、従って テンプレート:Math-次元球体の体積の、取りうる半径および方位角に亘る逐次積分 $$V_n(R)=\underset{0}{\overset{2\pi }{\int }}\underset{0}{\overset{R}{\int }}{V}_{n-2}(\sqrt{R^2-r^2})rdrd\theta$$ として書くことができ、この方位角座標に関する積分は直ちに計算できる。比例関係を適用することで、この体積が $$V_n(R)=(2\pi )V_{n-2}(R)\underset{0}{\overset{R}{\int }}(1-(r/R{)}^2{)}^{(n-2)/2}rdr$$ に等しいことが示される。テンプレート:Math と置換することによって積分を評価することができ、$$\begin{array}{cc}{V}_n(R)& =(2\pi )V_{n-2}(R)(-\frac{R^2}{n}(1-(r/R{)}^2{)}^{n/2}){|}_{r=0}^{r=R}\\ & =\frac{2\pi R^2}{n}{V}_{n-2}(R)\end{array}$$ を得る。これが2次元漸化式である。

体積公式の帰納法による証明に同じ手法を用いることができる。帰納法の基底段階は テンプレート:Math-次元球体と テンプレート:Math-次元球体であり、ここで テンプレート:Math と テンプレート:Math という事実を用いて簡単に直接確認できる。再帰段階は上記と同様であるが、テンプレート:Math-次元球体の体積に比例関係を適用する代わりに、帰納法の仮定が適用される。

1次元漸化式

比例関係は テンプレート:Mvar-次元球体と テンプレート:Math-次元球体の体積の関係に関する漸化式の証明にも使われる。比例式の証明の際に見たように、テンプレート:Mvar-次元球の体積は テンプレート:Math-次元球体の体積の積分として書くことができる。置換の代わりに、比例関係を被積分関数に現れる テンプレート:Math-次元球体の体積に適用し、$$V_n(R)=V_{n-1}(R)\underset{-R}{\overset{R}{\int }}(1-(x/R{)}^2{)}^{(n-1)/2}dx$$ を得る。被積分関数は偶関数であるため、対称性によって積分区間を テンプレート:Closed-closed に制限することができる。区間 テンプレート:Closed-closed 上で テンプレート:Math なる置換を適用することができるから、式は $$V_{n-1}(R)R\underset{0}{\overset{1}{\int }}(1-u{)}^{(n-1)/2}{u}^{-1/2}du$$ と書き換えられる。この積分はベータ関数 テンプレート:Math と呼ばれるよく知られた特殊関数のある値に等しく、求める体積はベータ関数を用いて $$V_n(R)=V_{n-1}(R)R\mathrm{B}(\frac{n+1}{2},\frac{1}{2})$$ となる。階乗と二項係数との関係とほぼ同じ意味で、ベータ関数はガンマ関数を用いて表されるから、その関係式を適用して $$V_n(R)=V_{n-1}(R)R\frac{\mathrm{\Gamma }(\frac{n+1}{2})\mathrm{\Gamma }(\frac{1}{2})}{\mathrm{\Gamma }(\frac{n}{2}+1)}$$ が得られる。値 テンプレート:Math を用いて1次元漸化式 $$V_n(R)=R\sqrt{\pi }\frac{\mathrm{\Gamma }(\frac{n+1}{2})}{\mathrm{\Gamma }(\frac{n}{2}+1)}{V}_{n-1}(R)$$ が得られる。

2次元漸化式と同様に、体積公式の帰納法による証明を得るために同じ手法を使用することができる。

球座標における直接積分

体積を球座標における体積要素の積分によって計算することができる。球面座標系は動径座標 テンプレート:Mvar と偏角座標 テンプレート:Math を持つ。ここで テンプレート:Math を除く各 テンプレート:Mvar の変域は テンプレート:Closed-open であり、テンプレート:Math の変域は テンプレート:Closed-open である。球体積要素は $$dV=r^{n-1}{\sin }^{n-2}({\varphi }_1){\sin }^{n-3}({\varphi }_2)\cdots \sin ({\varphi }_{n-2})drd{\varphi }_{1}d{\varphi }_2\cdots d{\varphi }_{n-1}$$ で与えられる。そして求める体積は、テンプレート:Mvar は テンプレート:Math から テンプレート:Mvar までと、角は取りうるすべての値に亘って取った積分 $$V_n(R)=\underset{0}{\overset{R}{\int }}\underset{0}{\overset{\pi }{\int }}\cdots \underset{0}{\overset{2\pi }{\int }}{r}^{n-1}{\sin }^{n-2}({\varphi }_1)\cdots \sin ({\varphi }_{n-2})d{\varphi }_{n-1}\cdots d{\varphi }_{1}dr$$ に等しい。被積分関数の各因子は一変数のみに依存するため、従ってこの逐次積分は積分の積 $$V_n(R)=\left(\underset{0}{\overset{R}{\int }}{r}^{n-1}dr\right)(\underset{0}{\overset{\pi }{\int }}{\sin }^{n-2}({\varphi }_1)d{\varphi }_1)\cdots \left(\underset{0}{\overset{2\pi }{\int }}d{\varphi }_{n-1}\right)$$ として書くことができる。動径成分の積分は テンプレート:Mvar に等しく、また偏角成分の積分区間を対称性により テンプレート:Closed-closed と書き換えれば $$V_n(R)=\frac{R^n}{n}(2\underset{0}{\overset{\pi /2}{\int }}{\sin }^{n-2}({\varphi }_1)d{\varphi }_1)\cdots \left(4\underset{0}{\overset{\pi /2}{\int }}d{\varphi }_{n-1}\right)$$ を得る。残った各々の積分はいまやベータ関数の特定の値で、$$V_n(R)=\frac{R^n}{n}\mathrm{B}(\frac{n-1}{2},\frac{1}{2})\mathrm{B}(\frac{n-2}{2},\frac{1}{2})\cdots \mathrm{B}(\frac{2}{2},\frac{1}{2})\cdot 2\mathrm{B}(\frac{1}{2},\frac{1}{2})$$ となる。ベータ関数はガンマ関数に書き換えることができ、$$V_n(R)=\frac{R^n}{n}\frac{\mathrm{\Gamma }(\frac{n-1}{2})\mathrm{\Gamma }(\frac{1}{2})}{\mathrm{\Gamma }(\frac{n}{2})}\frac{\mathrm{\Gamma }(\frac{n-2}{2})\mathrm{\Gamma }(\frac{1}{2})}{\mathrm{\Gamma }(\frac{n-1}{2})}\cdots \frac{\mathrm{\Gamma }(\frac{2}{2})\mathrm{\Gamma }(\frac{1}{2})}{\mathrm{\Gamma }(\frac{3}{2})}\cdot 2\frac{\mathrm{\Gamma }(\frac{1}{2})\mathrm{\Gamma }(\frac{1}{2})}{\mathrm{\Gamma }(\frac{2}{2})}$$ を得るが、この積は連鎖的に約分して畳み込める。値 テンプレート:Math 関数等式 テンプレート:Math を組み合わせて $$V_n(R)=\frac{2{\pi }^{n/2}{R}^n}{n\mathrm{\Gamma }(\frac{n}{2})}=\frac{{\pi }^{n/2}{R}^n}{\mathrm{\Gamma }(\frac{n}{2}+1)}$$ が導かれる。

ガウス積分

体積公式はガウス積分を用いることにより直接証明することができる。関数 $f(x_1,\dots ,x_n)=\exp (-\frac{1}{2}{\sum }_{i=1}^n{x}_i^2)$ を考えると、この関数は回転不変かつ各々一変数の函数の積になっている。これが積に書けるという事実とガウス積分の公式を適用して $$\underset{{𝐑}^n}{\int }fdV=\prod _{i=1}^n(\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}\exp (-x_i^2/2)d{x}_i)=(2\pi {)}^{n/2}$$ が得られる。ここで テンプレート:Mvar は テンプレート:Mvar-次元体積要素である。回転不変性を用いれば、同じ積分を球座標に関して $$\underset{{𝐑}^n}{\int }fdV=\underset{0}{\overset{\mathrm{\infty }}{\int }}\underset{S^{n-1}(r)}{\int }\exp (-r^2/2)dAdr,$$ と計算できる。ここで テンプレート:Math は半径 テンプレート:Mvar の テンプレート:Math-次元球面であり、テンプレート:Mvar は表面積要素（すなわち テンプレート:Math-次元体積要素）である。球面の表面積は、球体の体積に関するのと同様の比例関係を満足する。すなわち テンプレート:Math を半径 テンプレート:Mvar のテンプレート:Math-次元球面の表面積とすれば $$A_{n-1}(r)=r^{n-1}{A}_{n-1}(1)$$ が成り立つ。上記の積分にこれを適用すると $$A_{n-1}(1)\underset{0}{\overset{\mathrm{\infty }}{\int }}\exp (-r^2/2)r^{n-1}dr$$ なる式を得る。置換 テンプレート:Math を適用すれば、この式は $$A_{n-1}(1)2^{n/2-1}\underset{0}{\overset{\mathrm{\infty }}{\int }}{e}^{-t}{t}^{n/2-1}dt$$ と変形でき、これはガンマ関数の テンプレート:Math における値である。

二つの積分を併せれば $$A_{n-1}(1)=\frac{2{\pi }^{n/2}}{\mathrm{\Gamma }(\frac{n}{2})}$$ が示される。この式から半径 テンプレート:Mvar の テンプレート:Mvar-次元球体の体積を導出するには、半径 テンプレート:Mvar (テンプレート:Math) の球面の表面積を積分し、関数等式 テンプレート:Math を適用すればよい。そうして $$V_n(R)=\underset{0}{\overset{R}{\int }}\frac{2{\pi }^{n/2}}{\mathrm{\Gamma }(\frac{n}{2})}{r}^{n-1}dr=\frac{2{\pi }^{n/2}}{n\mathrm{\Gamma }(\frac{n}{2})}{R}^n=\frac{{\pi }^{n/2}}{\mathrm{\Gamma }(\frac{n}{2}+1)}{R}^n$$ が得られる。

テンプレート:Mvar-ノルムに関する球体

次が知られていることがディリクレの時代から知られているテンプレート:Rpテンプレート:Rp：

よって特に次が成立する： テンプレート:Math theorem

テンプレート:Math の場合は通常のユークリッドノルムであり、それ以外の テンプレート:Mvar は情報理論、符号理論、テンプレート:仮リンクなどの様々な文脈において現れる。

特に テンプレート:Math および テンプレート:Math の場合の体積はそれぞれ

$V_1^n(R)=\frac{2^n}{n!}{R}^n,$

$V_{\mathrm{\infty }}^n(R)=(2R{)}^n$

で与えられる。これらは正軸体および超立方体の体積に一致する。

$V_p^n(R)$は以下の漸化式 を満たす：

$V_p^n(R)=\frac{2\left(\frac{1}{p}\right)!R\left(\frac{n-1}{p}\right)!}{\left(\frac{n}{p}\right)!}{V}_p^{n-1}(R)$

一部の例外的なテンプレート:Mvarを除き、テンプレート:Mvar-球面（テンプレート:Mvar-球体の境界）の表面積は、テンプレート:Mvar-球体の半径に関する微分として計算することはできないが、テンプレート:仮リンクを用いて、体積を表面積上の積分として表すことができる。

余面積公式には、点から点へ テンプレート:Mvar-ノルムがどのくらい変化するかを考慮した相関係数が含まれる。テンプレート:Math および テンプレート:Math に対してこの因子は テンプレート:Mvar だが、テンプレート:Math ならば相関因子は テンプレート:Math である（半径 テンプレート:Mvar の テンプレート:Math-次元 テンプレート:Math-球面の表面積は テンプレート:Math 掛ける テンプレート:Math-球体の体積の テンプレート:Mvar における微分係数）。ほとんどの テンプレート:Mvar の値に対してこの定数は複雑な積分になる。

関連項目

• 超球面
• 球充填
• ハミング限界

外部リンク

• http://www.brouty.fr/Maths/sphere.html (derivation in hyperspherical coordinates.)
• テンプレート:MathWorld
• http://www.mathreference.com/ca-int,hsp.html