n次元球面の体積

同義語:球面の体積volume of the n-spheresurface area of the n-sphere

概要

n次元球面の体積(volume of the n-sphere)とは、単位球面 $S^n\subset\mathbb R^{n+1}$ の Riemann 計量による $n$ 次元の体積 $\lvert S^n\rvert$ のことで、$\lvert S^n\rvert=2\pi^{(n+1)/2}/\Gamma(\frac{n+1}2)$ である。Stokes の定理から単位球の体積 $\omega_{n+1}$ との関係 $\lvert S^n\rvert=(n+1)\omega_{n+1}$ が出て、球対称な関数の積分公式を Gauss 積分に使うと値が決まる。単位球の体積は 5 次元で、球面の体積は $S^6$ で最大になり、次元を上げると $0$ に近づく。

$$\newcommand{C}[0]{\mathbb{C}} \newcommand{div}[0]{\mathbin{÷}} \newcommand{N}[0]{\mathbb{N}} \newcommand{Q}[0]{\mathbb{Q}} \newcommand{R}[0]{\mathbb{R}} \newcommand{Z}[0]{\mathbb{Z}} $$

前提知識: n次元球面, 変数変換公式(多変数積分), Γ関数, Stokesの定理

何を測る量か

単位円周 $S^1$ の長さは $2\pi$、単位球面 $S^2$ の面積は $4\pi$ である。一般の $n$ 次元球面 $S^n\subset\mathbb R^{n+1}$ についても、その $n$ 次元の体積(「表面積」)を Γ 関数で表す公式がある:
$$ \lvert S^n\rvert=\frac{2\pi^{(n+1)/2}}{\Gamma\bigl(\frac{n+1}2\bigr)},\qquad \lambda(B^{n+1})=\frac{\pi^{(n+1)/2}}{\Gamma\bigl(\frac{n+3}2\bigr)}=\frac{\lvert S^n\rvert}{n+1}. $$
右の式は、内部の球の体積との関係「体積を半径で微分すると表面積になる」を表す。
この記事では、球面の体積を Riemann 計量から定め、上の公式を 2 つの事実から導く。1 つは Stokes の定理による $\lvert S^n\rvert=(n+1)\lambda(B^{n+1})$ であり、もう 1 つは球対称な関数の積分を 1 変数の積分に直す公式である。後者を Gauss 積分 $\int_{\mathbb R^{n+1}}e^{-\lvert x\rvert^2}dx=\pi^{(n+1)/2}$ に使うと値が決まる。値を並べると、球の体積は 5 次元で、球面の体積は $S^6$ で最大になり、次元が上がると $0$ に近づく。

定義

$n\ge0$ とする。$\lvert x\rvert$ は Euclid ノルム、$\lambda$ は Lebesgue 測度で、
$$ S^n:=\{x\in\mathbb R^{n+1}\mid\lvert x\rvert=1\},\quad B^{n+1}:=\{x\mid\lvert x\rvert<1\},\quad D^{n+1}:=\{x\mid\lvert x\rvert\le1\} $$
とおく。$\omega_{n+1}:=\lambda(B^{n+1})$ を単位球の体積とする。$r>0$ について $\lambda(rB^{n+1})=r^{n+1}\omega_{n+1}$ である(変数変換公式(多変数積分) の命題「線形写像による体積の変化」を $x\mapsto rx$ に使う)。$D^{n+1}\setminus B^{n+1}=S^n$ は、すべての $\varepsilon\in(0,1)$ について $(1+\varepsilon)B^{n+1}\setminus(1-\varepsilon)B^{n+1}$ に含まれ、その体積 $\bigl((1+\varepsilon)^{n+1}-(1-\varepsilon)^{n+1}\bigr)\omega_{n+1}$ は $\varepsilon\to0$ で $0$ に近づくので、$\lambda(S^n)=0$、$\lambda(D^{n+1})=\omega_{n+1}$ である。

n次元球面の体積

$n\ge1$ のとき、$S^n$ は $\mathbb R^{n+1}$ のコンパクトな $n$ 次元部分多様体であり、$\mathbb R^{n+1}$ の内積から Riemann 計量を受け継ぐ。その Riemann 体積(Riemann多様体 の定義「リーマン体積」)による $S^n$ 全体の体積を $\lvert S^n\rvert$ と書き、$n$ 次元球面の体積 という。$n=0$ では $S^0=\{\pm1\}$ の点の個数をとり、$\lvert S^0\rvert:=2$ とする。

$n=1$ なら $\lvert S^1\rvert$ は単位円周の長さ、$n=2$ なら単位球面の面積である。半径 $r$ の球面 $S^n(r):=rS^n$ は、拡大 $x\mapsto rx$ によって計量を $r^2$ 倍した $S^n$ と同じ Riemann 多様体になるので、$\lvert S^n(r)\rvert=r^n\lvert S^n\rvert$ である(Riemann多様体 の命題「計量の一様な拡大と共形変形」)。

球面の体積と球の体積

$\mathbb R^{n+1}$ の座標を $x_1,\dots,x_{n+1}$ とし、$\Omega:=dx_1\wedge\cdots\wedge dx_{n+1}$、$E:=\sum_ix_i\,\partial/\partial x_i$ とおく。$\Omega$ に $E$ を代入した $n$ 形式を
$$ \sigma:=\iota_E\Omega=\sum_{i=1}^{n+1}(-1)^{i-1}x_i\,dx_1\wedge\cdots\wedge\widehat{dx_i}\wedge\cdots\wedge dx_{n+1} $$
とする($\widehat{\ \ }$ はその因子を除くことを表す)。$n=1$ なら $\sigma=x\,dy-y\,dx$ である。

球面の体積は球の体積のn+1倍

$n\ge0$ について $\lvert S^n\rvert=(n+1)\,\omega_{n+1}$ である。

$n=0$ では $\lvert S^0\rvert=2=\lambda((-1,1))$ である。$n\ge1$ とし、$S^n$ には閉球体 $D^{n+1}$ の境界としての向き(外向きの法線を先頭に置く向き。Stokesの定理 の定義「向き付けられた境界付き多様体と境界の向き」)を入れる。
$d\sigma=(n+1)\Omega$。 $d\bigl(x_i\,dx_1\wedge\cdots\widehat{dx_i}\cdots\wedge dx_{n+1}\bigr)=dx_i\wedge dx_1\wedge\cdots\widehat{dx_i}\cdots\wedge dx_{n+1}=(-1)^{i-1}\Omega$ なので、$d\sigma=\sum_i(-1)^{i-1}(-1)^{i-1}\Omega=(n+1)\Omega$ である。
$S^n$ 上で $\sigma$ は体積形式。 $p\in S^n$ で $E_p=p$ は長さ $1$ の外向き法線である。$T_pS^n$ の正の向きの正規直交基底 $(v_1,\dots,v_n)$ をとると、$(p,v_1,\dots,v_n)$ は $\mathbb R^{n+1}$ の正の向きの正規直交基底なので、$\sigma_p(v_1,\dots,v_n)=\Omega(p,v_1,\dots,v_n)=1$ である。一方、Riemann 体積形式 $\operatorname{vol}=\sqrt{\det(g_{ij})}\,dy^1\wedge\cdots\wedge dy^n$($y$ は正の向きの局所座標)も、$v_k=\sum_ia_{ik}\,\partial/\partial y^i$ と書くと、正規直交性 $a^{\mathsf T}(g_{ij})a=I$ から $(\det a)^2\det(g_{ij})=1$、$\det a>0$ なので $\operatorname{vol}(v_1,\dots,v_n)=\sqrt{\det(g_{ij})}\det a=1$ である。$n$ 次元空間の上の $n$ 次交代形式は 1 つの基底での値で決まるので、$S^n$ への引き戻しについて $\sigma|_{S^n}=\operatorname{vol}$ であり、$\int_{S^n}\sigma=\lvert S^n\rvert$ である。
Stokes の定理。 $D^{n+1}$ はコンパクトな境界付き多様体で、境界は $S^n$ である(Stokesの定理 の節「設定:境界付き多様体・向き・積分」)。Stokesの定理 の定理「Stokes の定理」により
$$ \lvert S^n\rvert=\int_{S^n}\sigma=\int_{D^{n+1}}d\sigma=(n+1)\int_{D^{n+1}}\Omega=(n+1)\,\lambda(D^{n+1})=(n+1)\,\omega_{n+1} $$
である。$\square$

この命題は「$\lambda(rB^{n+1})=r^{n+1}\omega_{n+1}$ を $r$ で微分すると $r^n\lvert S^n\rvert$」という関係の $r=1$ の場合である。

球対称な関数の積分

球対称な関数の積分

$n\ge0$ とし、$g\colon(0,\infty)\to[0,\infty]$ を Borel 可測関数とする。このとき
$$ \int_{\mathbb R^{n+1}}g(\lvert x\rvert)\,dx=\lvert S^n\rvert\int_0^\infty g(r)\,r^n\,dr $$
が成り立つ(両辺は $\infty$ でもよく、左辺の被積分関数は原点で $0$ とする)。$g$ が実数値で、$\int_0^\infty\lvert g(r)\rvert r^n\,dr<\infty$ のときも同じ等式が成り立つ。

$(0,\infty)$ の Borel 集合 $A$ について
$$ \nu_1(A):=\lambda\bigl(\{x\ne0\mid\lvert x\rvert\in A\}\bigr),\qquad\nu_2(A):=(n+1)\,\omega_{n+1}\int_Ar^n\,dr $$
とおく。$x\mapsto\lvert x\rvert$ は連続なので $\nu_1$ は測度であり、$\nu_2$ も測度である。$0\le a< b$ について $\{a<\lvert x\rvert\le b\}=bD^{n+1}\setminus aD^{n+1}$($a=0$ では原点を除く)なので
$$ \nu_1((a,b])=(b^{n+1}-a^{n+1})\,\omega_{n+1}=\nu_2((a,b]) $$
である。区間 $(a,b]$($0\le a< b$)と空集合の全体は π 系で、$(0,\infty)$ の Borel 集合全体を生成し、増加列 $(0,k]$ は $(0,\infty)$ を覆って $\nu_1((0,k])<\infty$ である。Dynkinのπ-λ定理 の定理「π系の上で一致する測度の一意性」の 2 により $\nu_1=\nu_2$ である。
$g=\chi_A$ なら主張は $\nu_1(A)=\nu_2(A)$ そのもの(prop-nsv-stokes で $(n+1)\omega_{n+1}=\lvert S^n\rvert$)である。両辺は $g$ について線形なので $0$ 以上の単関数でも成り立ち、単関数の増加列 $g_k\uparrow g$ に Lebesgue積分 の定理「単調収束定理」を両辺で使えば一般の $g\ge0$ で成り立つ。実数値の場合は、$\lvert g\rvert$ に使って両辺の可積分性を得たうえで、$g$ の正の部分と負の部分に分けて引けばよい。$\square$

主定理:Γ関数による値

n次元球面の体積

$n\ge0$ について
$$ \lvert S^n\rvert=\frac{2\pi^{(n+1)/2}}{\Gamma\bigl(\frac{n+1}2\bigr)},\qquad\omega_{n+1}=\frac{\pi^{(n+1)/2}}{\Gamma\bigl(\frac{n+3}2\bigr)} $$
である。半径 $r$ の球面と球については、それぞれ $r^n$ 倍、$r^{n+1}$ 倍である。

Gauss 積分を 2 通りに計算する

$I:=\int_{\mathbb R^{n+1}}e^{-\lvert x\rvert^2}\,dx$ を 2 通りに計算する。$e^{-\lvert x\rvert^2}=\prod_ie^{-x_i^2}$ なので、Tonelli の定理(Fubiniの定理 の定理「Tonelliの定理」)と 1 変数の値 $\int_{\mathbb R}e^{-t^2}dt=\sqrt\pi$(Gauss積分 の定理「Gauss積分の値」)から $I=\pi^{(n+1)/2}$ である。一方、prop-nsv-radial を $g(r)=e^{-r^2}$ に使い、$s=r^2$ と置換すると($r^n\,dr=\frac12s^{(n-1)/2}\,ds$)、
$$ I=\lvert S^n\rvert\int_0^\infty e^{-r^2}r^n\,dr=\frac{\lvert S^n\rvert}2\int_0^\infty s^{\frac{n+1}2-1}e^{-s}\,ds=\frac{\lvert S^n\rvert}2\,\Gamma\Bigl(\frac{n+1}2\Bigr) $$
である(Γ関数 の定義「$\Gamma$ 関数」)。両者を比べて $\lvert S^n\rvert$ の式を得る。prop-nsv-stokes と $\Gamma(s+1)=s\Gamma(s)$(Γ関数 の命題「関数等式と階乗」)を $s=\frac{n+1}2$ に使うと
$$ \omega_{n+1}=\frac{\lvert S^n\rvert}{n+1}=\frac{2\pi^{(n+1)/2}}{(n+1)\Gamma\bigl(\frac{n+1}2\bigr)}=\frac{\pi^{(n+1)/2}}{\Gamma\bigl(\frac{n+3}2\bigr)} $$
である。半径 $r$ の場合は定義の節で見た $r^n$、$r^{n+1}$ 倍である。$\square$

球の体積の式は Gauss積分 の定理「単位球の体積」にもあり、そこでは層ごとの積分 $\int f=\int_0^\infty\lambda(\{f>s\})\,ds$ で示されている。Axl26 5.44(p. 141)は次元を 2 つ下げる漸化式で、同じ値を Γ 関数を使わない形で示している。上の証明はどちらとも違い、球面の体積を先に求めて、球の体積を prop-nsv-stokes から得ている。
$\Gamma(\frac12)=\sqrt\pi$、$\Gamma(1)=1$ と関数等式から、次の値を得る。$\lvert S^n\rvert$ は $\pi^{\lceil n/2\rceil}$ の有理数倍である。

$n$$\lvert S^n\rvert$近似値$\omega_{n+1}=\lambda(B^{n+1})$近似値
$0$$2$$2$$2$$2$
$1$$2\pi$$6.2832$$\pi$$3.1416$
$2$$4\pi$$12.566$$\frac43\pi$$4.1888$
$3$$2\pi^2$$19.739$$\frac12\pi^2$$4.9348$
$4$$\frac83\pi^2$$26.319$$\frac8{15}\pi^2$$5.2638$
$5$$\pi^3$$31.006$$\frac16\pi^3$$5.1677$
$6$$\frac{16}{15}\pi^3$$33.073$$\frac{16}{105}\pi^3$$4.7248$
$7$$\frac13\pi^4$$32.470$$\frac1{24}\pi^4$$4.0587$

次元を上げたときのふるまい

$k\ge1$ について $\omega_k:=\lambda(B^k)=\pi^{k/2}/\Gamma(\frac k2+1)$ と書き、$\omega_0:=1$ とおく(同じ式で $k=0$ とした値)。

漸化式・最大値・極限
  1. $k\ge2$ について $\omega_k=\dfrac{2\pi}k\,\omega_{k-2}$ である。
  2. $n\ge0$ について $\lvert S^{n+2}\rvert=2\pi\,\omega_{n+1}$ である。
  3. $\omega_k$ は $k=5$ でただ 1 つの最大値 $\frac{8}{15}\pi^2$ をとり、$\lvert S^n\rvert$ は $n=6$ でただ 1 つの最大値 $\frac{16}{15}\pi^3$ をとる。
  4. $k\to\infty$ で $\omega_k\to0$、$n\to\infty$ で $\lvert S^n\rvert\to0$ である。

1:主定理の式と $\Gamma(\frac k2+1)=\frac k2\Gamma(\frac k2)$ から $\omega_k/\omega_{k-2}=\pi\,\Gamma(\frac k2)/\Gamma(\frac k2+1)=2\pi/k$ である($k=2$ では $\omega_2=\pi=\frac{2\pi}2\omega_0$)。
2:主定理から $\lvert S^{n+2}\rvert=2\pi^{(n+3)/2}/\Gamma(\frac{n+3}2)=2\pi\,\omega_{n+1}$ である。
3:1 により、$k$ と $k-2$ の比は $2\pi/k$ で、$k\le6$ なら $1$ より大きく、$k\ge7$ なら $1$ より小さい。よって偶数の $k$ では $\omega_0<\omega_2<\omega_4<\omega_6>\omega_8>\cdots$、奇数の $k$ では $\omega_1<\omega_3<\omega_5>\omega_7>\cdots$ である。$\omega_5=\frac8{15}\pi^2>\omega_6=\frac16\pi^3$($\pi<\frac{16}5$ による)なので、最大は $\omega_5$ だけである。2 により $\lvert S^n\rvert=2\pi\,\omega_{n-1}$($n\ge2$)であり、$\lvert S^0\rvert,\lvert S^1\rvert<2\pi\omega_5$ なので、$\lvert S^n\rvert$ の最大は $n-1=5$ のときだけである。
4:$k>4\pi$ なら $\omega_k<\frac12\omega_{k-2}$ なので、$\omega_k\to0$ である。2 により $\lvert S^n\rvert\to0$ である。$\square$

2 は Archimedes の「球面の面積は外接円柱の側面積に等しい」の高次元版にあたる($n=0$ では $\lvert S^2\rvert=2\pi\cdot2=4\pi$)。

高次元の球の体積は表面の近くに集まる

$0<\varepsilon<1$ とすると、半径 $1-\varepsilon$ の球と単位球の体積の比は $(1-\varepsilon)^k$ なので、殻 $B^k\setminus(1-\varepsilon)B^k$ が占める割合 $1-(1-\varepsilon)^k$ は $k\to\infty$ で $1$ に近づく。たとえば $\varepsilon=0.01$、$k=1000$ なら割合は $1-0.99^{1000}>0.9999$ である。また単位球を囲む立方体 $[-1,1]^k$ の体積は $2^k$ で、比 $\omega_k/2^k$ は 4 により $0$ に近づく。高次元では、球は立方体の中で無視できるほど小さい。

微分形式による見方

$\mathbb R^{n+1}\setminus\{0\}$ の上で $\tau:=\sigma/\lvert x\rvert^{n+1}$ とおく。$n=1$ なら $\tau=(x\,dy-y\,dx)/(x^2+y^2)$ は外微分 の例「原点を除いた平面の角度形式」である。

原点を除いた空間の閉形式

$n\ge1$ とする。$\tau$ は $\mathbb R^{n+1}\setminus\{0\}$ 上の閉形式で、$\int_{S^n}\tau=\lvert S^n\rvert>0$ である。したがって $\tau$ は完全形式でなく、$H^n_{dR}(\mathbb R^{n+1}\setminus\{0\})\ne0$ である。

1 形式 $\alpha$ について $\alpha\wedge\Omega=0$($n+2$ 次)なので、$0=\iota_E(\alpha\wedge\Omega)=\alpha(E)\Omega-\alpha\wedge\sigma$、すなわち $\alpha\wedge\sigma=\alpha(E)\Omega$ である。$\alpha=\sum_jx_j\,dx_j$ とすると $\alpha(E)=\lvert x\rvert^2$ である。$d(\lvert x\rvert^{-(n+1)})=-(n+1)\lvert x\rvert^{-(n+3)}\alpha$ なので、prf-nsv-stokes の $d\sigma=(n+1)\Omega$ とあわせて
$$ d\tau=-(n+1)\lvert x\rvert^{-(n+3)}\lvert x\rvert^2\,\Omega+\lvert x\rvert^{-(n+1)}(n+1)\Omega=0 $$
である。$S^n$ の上では $\lvert x\rvert=1$ なので $\tau=\sigma$ であり、prf-nsv-stokes により $\int_{S^n}\tau=\lvert S^n\rvert>0$ である。$\tau=d\eta$ なら、包含 $\iota\colon S^n\to\mathbb R^{n+1}\setminus\{0\}$ による引き戻しは $\iota^*\tau=d(\iota^*\eta)$ で、Stokesの定理 の系「閉多様体上の積分と完全形式」によりその積分は $0$ になり、矛盾する。$\square$

同じ計算から、$r:=\lvert x\rvert$ について $dr\wedge(\sigma/r)=\Omega$ も成り立つ($dr=\alpha/r$)。$\sigma/r$ を半径 $r$ の球面に制限したものはその体積形式であり、この等式は「体積要素 = 半径方向の長さ × 球面の体積要素」を形式の言葉で述べている。prop-nsv-radial はその測度論での表現である。

反例:条件を外すと崩れること

外す条件反例成り立たなくなること
prop-nsv-radial の「被積分関数が $\lvert x\rvert$ だけで決まる」$f(x)=\chi_{\{x_1>0\}}(x)\,e^{-\lvert x\rvert^2}$$\int f=\lvert S^n\rvert\int_0^\infty f(re_1)r^n\,dr$(左辺 $\frac12\pi^{(n+1)/2}$、右辺 $\pi^{(n+1)/2}$)
prop-nsv-radial の「$g\ge0$ または $\int_0^\infty\lvert g\rvert r^n\,dr<\infty$」$g(r)=r^{-(n+1)}\sin r$左辺が Lebesgue 積分として定まること(右辺の広義積分は $\lvert S^n\rvert\frac\pi2$ に収束する)
ノルムが Euclid ノルム$\lvert x\rvert_\infty=\max_i\lvert x_i\rvert$ の単位球 $(-1,1)^{n+1}$$\omega_{n+1}=\pi^{(n+1)/2}/\Gamma(\frac{n+3}2)$(体積は $2^{n+1}$)
反例の確認を開く

1 行目:半空間 $\{x_1>0\}$ と $\{x_1<0\}$ は鏡映 $x_1\mapsto-x_1$ で移りあい、境界の超平面は零集合なので、左辺は $I=\pi^{(n+1)/2}$ の半分である。右辺は $f(re_1)=e^{-r^2}$ を使った主定理の証明の計算そのもので $I$ に等しい。

2 行目:$g(r)r^n=\frac{\sin r}r$ で、$\int_0^\infty\frac{\sin r}r\,dr=\frac\pi2$ は広義積分として収束する。しかし $\int_0^\infty\frac{(\sin r)^\pm}{r}\,dr=\infty$($(\sin r)^+$ は $k$ が偶数の区間 $[k\pi,(k+1)\pi]$ で、$(\sin r)^-$ は $k$ が奇数の区間で、それぞれ寄与が $\frac{2}{(k+1)\pi}$ 以上であり、その和は発散する)なので、prop-nsv-radial を $g^\pm$ に使うと $\int g(\lvert x\rvert)^\pm\,dx=\infty$ であり、$g(\lvert x\rvert)$ の Lebesgue 積分は $\infty-\infty$ になって定まらない。

3 行目:$\{\lvert x\rvert_\infty<1\}=(-1,1)^{n+1}$ の体積は $2^{n+1}$ で、$n\ge1$ では $\omega_{n+1}$ と異なる($n=1$ で $4\ne\pi$)。$n=0$ では両方 $2$ である。

関連項目

参考文献

Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する