コンテンツにスキップ

ℝ^d の変数変換公式と極座標:|det A| から球面測度まで

対象読者: 本テーマの直積測度と Fubini の定理、Radon 測度の章、および逆関数定理と行列式を学んだ人。微分形式や多様体は仮定しません。

前提となる記事:直積測度と Fubini–Tonelli の定理:積分の順序交換と無限直積、Radon 測度と Riesz–Markov–角谷の表現定理:正の線形汎関数を測度で表す、逆関数定理・陰関数定理と条件付き極値:縮小写像で「方程式を解く」、行列式とその性質:符号付き体積としての det

生 Markdown

平面の積分は、座標を取り替えると計算しやすくなることがあります。たとえば円板の上の積分は、縦横の座標よりも、原点からの距離と角度で測るほうが楽です。このとき面積や体積が「何倍に伸び縮みするか」を正確に述べるのが変数変換公式です。この記事では、それを一般の次元で厳密に証明し、球の体積や球面の面積を求められるようにします。

  • Rd\mathbb{R}^d のルベーグ測度 mm は、正則な線形写像 TT によって一律に ∣det⁡T∣|\det T| 倍されます:m(T(E))=∣det⁡T∣ m(E)m(T(E))=|\det T|\,m(E)。証明の骨格は「並進不変でコンパクト集合で有限なボレル測度はルベーグ測度の定数倍しかない」という一意性で、定数は基本行列ごとに Fubini の定理で計算します。
  • C1C^1 微分同相 Φ:U→V\Phi:U\to V について、ルベーグ積分の意味で ∫Vf dm=∫Uf(Φ(x)) ∣det⁡DΦ(x)∣ dx\int_Vf\,dm=\int_Uf(\Phi(x))\,|\det D\Phi(x)|\,dx が成り立ちます。立方体を細かく分けて Φ\Phi を各小立方体で DΦ(xj)D\Phi(x_j) に置き換え、誤差を平均値不等式で抑え、外正則性(ボレル集合の測度を、それを含む開集合の測度で上からいくらでも近く近似できる性質)で一般のボレル集合(開集合から補集合と可算和をくり返して作れる集合)に広げます。
  • C1C^1 写像は測度 00 の集合を測度 00 の集合に写します。これを使うと、測度 00 の例外集合の上で単射性やヤコビアンの非退化性が崩れていても公式が使えます。極座標・球座標はこの形で使います。
  • 単位球面 Sd−1S^{d-1} 上の球面測度(球面の面積にあたるもの)を σ(A)=d⋅m({rω:0<r≤1, ω∈A})\sigma(A)=d\cdot m(\{r\omega:0<r\le1,\ \omega\in A\}) で定義すると、∫Rdf dm=∫0∞ ⁣∫Sd−1f(rω) dσ(ω) rd−1dr\int_{\mathbb{R}^d}f\,dm=\int_0^\infty\!\int_{S^{d-1}}f(r\omega)\,d\sigma(\omega)\,r^{d-1}dr が成り立ちます。σ\sigma は回転不変で、d=3d=3 では dσ=sin⁡θ dθ dφd\sigma=\sin\theta\,d\theta\,d\varphi です。
  • ガウス積分を二通りに計算して、σ(Sd−1)=2πd/2/Γ(d/2)\sigma(S^{d-1})=2\pi^{d/2}/\Gamma(d/2)、単位球の体積 πd/2/Γ(d/2+1)\pi^{d/2}/\Gamma(d/2+1) を得ます。

1. 動機 — 2 次元のリーマン積分から dd 次元のルベーグ積分へ

Section titled “1. 動機 — 2 次元のリーマン積分から ddd 次元のルベーグ積分へ”

重積分とフビニの定理 では、平面の有界閉集合上のリーマン積分について変数変換公式を証明しました。ところが解析学の多くの場面で使う変数変換は、この枠に収まりません。

  • 偏微分方程式では、ラプラス方程式の基本解 ∣x∣2−d|x|^{2-d} が原点の近くで可積分であることを、Rd\mathbb{R}^d の極座標で ∫01r2−drd−1 dr<∞\int_0^1r^{2-d}r^{d-1}\,dr<\infty と確かめます(ラプラス方程式)。被積分関数は非有界で、次元は一般の dd です。
  • フーリエ解析では、⟨x⟩−s=(1+∣x∣2)−s/2\langle x\rangle^{-s}=(1+|x|^2)^{-s/2} が Rd\mathbb{R}^d 全体で可積分になる条件 s>ds>d を、やはり極座標で判定します。積分範囲は有界ではありません。
  • 確率論では、正規分布 N(μ,Σ)N(\mu,\Sigma) の密度を、標準正規分布を線形写像 x↦Σ1/2x+μx\mapsto\Sigma^{1/2}x+\mu で写して求めます。必要なのは「線形写像はルベーグ測度を ∣det⁡∣|\det| 倍する」という事実で、関数は一般の可測関数です。
  • 物理では、3 次元の球面上の積分を sin⁡θ dθ dφ\sin\theta\,d\theta\,d\varphi で計算します。球面上の「面積測度」とは何か、なぜ回転不変なのかを、一度は確かめておく必要があります。

この記事は、これらをすべてルベーグ積分の枠で一度に正当化します。リーマン積分版との違いは、(i) 次元が一般の dd、(ii) 被積分関数が一般のボレル関数(非有界・非有界領域でよい)、(iii) 線形写像の場合も含めて証明が測度論の道具(一意性定理・外正則性・単調収束定理)で閉じている、の三点です。

Rd\mathbb{R}^d の点 x=(x1,…,xd)x=(x_1,\ldots,x_d) に、ユークリッドノルム ∣x∣=(∑ixi2)1/2|x|=(\sum_ix_i^2)^{1/2} と最大値ノルム(各成分の絶対値の最大)∥x∥∞=max⁡i∣xi∣\|x\|_\infty=\max_i|x_i| を使い分けます。行列 A=(aij)∈Md(R)A=(a_{ij})\in M_d(\mathbb{R}) には

∥A∥∞=max⁡1≤i≤d∑j=1d∣aij∣\|A\|_\infty=\max_{1\le i\le d}\sum_{j=1}^d|a_{ij}|

を使います。∣(Ax)i∣≤∑j∣aij∣∣xj∣≤∥A∥∞∥x∥∞|(Ax)_i|\le\sum_j|a_{ij}||x_j|\le\|A\|_\infty\|x\|_\infty なので ∥Ax∥∞≤∥A∥∞∥x∥∞\|Ax\|_\infty\le\|A\|_\infty\|x\|_\infty、また ∥AB∥∞≤∥A∥∞∥B∥∞\|AB\|_\infty\le\|A\|_\infty\|B\|_\infty です。

a∈Rda\in\mathbb{R}^d、r>0r>0 に対し、閉立方体を Q(a,r)={x:∥x−a∥∞≤r}Q(a,r)=\{x:\|x-a\|_\infty\le r\} と書きます(辺の長さ 2r2r、中心 aa)。箱とは区間の直積 I1×⋯×IdI_1\times\cdots\times I_d のことです。

2.2. Rd\mathbb{R}^d のルベーグ測度

Section titled “2.2. Rd\mathbb{R}^dRd のルベーグ測度”

定義 2.1(Rdℝ^d のルベーグ測度)↗

B(Rd)\mathcal{B}(\mathbb{R}^d) を Rd\mathbb{R}^d のボレル集合族、m1m_1 を R\mathbb{R} のルベーグ測度のボレル集合への制限とする。d≥2d\ge2 について帰納的に

md=m1×md−1(B(R)⊗B(Rd−1)=B(Rd) 上)m_d=m_1\times m_{d-1}\qquad\bigl(\mathcal{B}(\mathbb{R})\otimes\mathcal{B}(\mathbb{R}^{d-1})=\mathcal{B}(\mathbb{R}^d)\text{ 上}\bigr)

と定め、Rd\mathbb{R}^d のルベーグ測度と呼ぶ。次元が明らかなときは mm と書く。さらに、ボレル集合 BB と、ある m(Z)=0m(Z)=0 のボレル集合 ZZ の部分集合 NN によって E=B∪NE=B\cup N と書ける集合 EE をルベーグ可測といい、その全体を Ld\mathcal{L}_d、m(E):=m(B)m(E):=m(B) と定める。

σ-加法族の等式は 命題 2.3[直積測度と Fubini–Tonelli の定理] によります。直積測度は 定理 3.2[直積測度と Fubini–Tonelli の定理] により存在して一意です。d=1d=1 のとき L1\mathcal{L}_1 はカラテオドリの条件による可測集合族と一致します(命題 6.6[可測集合とルベーグ測度])。Radon 測度と Riesz–Markov の表現定理 では mdm_d を m1×⋯×m1m_1\times\cdots\times m_1 と書きましたが、括弧の付け方によらず箱の上で同じ値をとるので、次の補題と測度の一意性(定理 4.5[測度空間と Carathéodory の拡張定理])により同じ測度です。

補題 2.2(箱の測度と超平面)✓

  1. 区間 I1,…,IdI_1,\ldots,I_d(端点を含むかどうかは任意)について m(I1×⋯×Id)=ℓ(I1)⋯ℓ(Id)m(I_1\times\cdots\times I_d)=\ell(I_1)\cdots\ell(I_d)。ここで ℓ\ell は区間の長さ。
  2. 座標超平面 H={x:xi=c}H=\{x:x_i=c\} は m(H)=0m(H)=0 を満たす。とくに閉立方体の境界は測度 00 である。
証明(補題 2.2)
  1. dd についての帰納法です。d=1d=1 は 命題 3.5[可測集合とルベーグ測度] です。d≥2d\ge2 では、I2×⋯×IdI_2\times\cdots\times I_d がボレル集合であることに注意して 定理 3.2[直積測度と Fubini–Tonelli の定理] の長方形の値の式を使うと、md(I1×(I2×⋯×Id))=ℓ(I1) md−1(I2×⋯×Id)m_d(I_1\times(I_2\times\cdots\times I_d))=\ell(I_1)\,m_{d-1}(I_2\times\cdots\times I_d) で、帰納法の仮定から主張が出ます。

  2. H=⋃kHkH=\bigcup_{k}H_k、Hk={x∈H:∥x∥∞≤k}H_k=\{x\in H:\|x\|_\infty\le k\} と書けます。HkH_k は第 ii 因子が 1 点 [c,c][c,c](長さ 00)の箱なので、1 から m(Hk)=0m(H_k)=0 です。可算劣加法性から m(H)=0m(H)=0。閉立方体の境界は 2d2d 枚の座標超平面の和集合に含まれます。

k≥0k\ge0 について、∏i=1d[ji2−k,(ji+1)2−k)\prod_{i=1}^d[j_i2^{-k},(j_i+1)2^{-k})(ji∈Zj_i\in\mathbb{Z})の形の集合を第 kk 世代の 2 進立方体と呼びます。同じ世代の 2 進立方体は互いに素で Rd\mathbb{R}^d を覆い、異なる世代の 2 進立方体は「一方が他方を含む」か「互いに素」のどちらかです(各座標で 2 進区間について同じことが成り立つためです)。

次の補題のねらいは、形の複雑な開集合を、扱いやすい 2 進立方体の可算個の貼り合わせに分解することです。立方体について示したことを、開集合へ、さらに外正則性を通じて一般のボレル集合へ広げるために、この後で何度も使います。証明は、開集合の中に入る 2 進立方体のうち、これ以上大きい 2 進立方体に含まれないものを集める、という方針です。

補題 2.3(開集合の 2 進立方体分解)✓

W⊆RdW\subseteq\mathbb{R}^d を開集合とする。互いに素な 2 進立方体の可算族 {Dn}\{D_n\} で、各閉包 Dn‾\overline{D_n} が WW に含まれ、W=⋃nDn=⋃nDn‾W=\bigcup_nD_n=\bigcup_n\overline{D_n} となるものが存在する。とくに m(W)=∑nm(Dn)=∑nm(Dn‾)m(W)=\sum_nm(D_n)=\sum_nm(\overline{D_n}) であり、WW は内部が互いに素な閉立方体の可算和である。

証明(補題 2.3)

C\mathcal{C} を、世代 k≥0k\ge0 の 2 進立方体 DD で D‾⊆W\overline D\subseteq W となるもの全体とします。x∈Wx\in W なら、WW が開なので ∥y−x∥∞≤δ\|y-x\|_\infty\le\delta ならば y∈Wy\in W となる δ>0\delta>0 があり、2−k≤δ2^{-k}\le\delta となる世代 kk で xx を含む 2 進立方体 DD を取ると、D‾\overline D の点 yy は ∥y−x∥∞≤2−k≤δ\|y-x\|_\infty\le2^{-k}\le\delta を満たすので D∈CD\in\mathcal{C} です。よって W=⋃CW=\bigcup\mathcal{C}。D∈CD\in\mathcal{C} を含む世代 ≥0\ge0 の 2 進立方体は有限個(各世代に一つ)なので、C\mathcal{C} の中で包含に関して極大な元が DD を含む形で存在します。極大元どうしは、一方が他方を含めば極大性から一致するので、異なる極大元は互いに素です。2 進立方体全体は可算なので極大元も可算個で、これを DnD_n とします。W⊆⋃Dn⊆⋃Dn‾⊆WW\subseteq\bigcup D_n\subseteq\bigcup\overline{D_n}\subseteq W です。測度の式は可算加法性と、Dn‾∖Dn\overline{D_n}\setminus D_n が立方体の境界に含まれ測度 00(補題 2.2)であることから従います。

外正則性も使います。命題 5.4[Radon 測度と Riesz–Markov–角谷の表現定理] により、任意のボレル集合 EE について m(E)=inf⁡{m(W):W⊇E は開}m(E)=\inf\{m(W):W\supseteq E\text{ は開}\} です。

3. 線形写像によるルベーグ測度の変換

Section titled “3. 線形写像によるルベーグ測度の変換”

この節で使う論法を先に述べておきます。二つの測度が等しいことを示すには、すべてのボレル集合を調べる代わりに、共通部分で閉じた集合族(π-系)の上で値が一致し、有限測度の集合で全体を覆えることを確かめれば十分です(測度の一意性定理)。箱や 2 進立方体の測度は計算できるので、そこから全体へ一致が広がります。

補題 3.1(並進不変性)✓

a∈Rda\in\mathbb{R}^d とボレル集合 EE について、E+a={x+a:x∈E}E+a=\{x+a:x\in E\} はボレル集合で m(E+a)=m(E)m(E+a)=m(E) である。

証明(補題 3.1)

τ(x)=x−a\tau(x)=x-a は連続なので、E+a=τ−1(E)E+a=\tau^{-1}(E) はボレル集合です(連続写像によるボレル集合の逆像はボレル集合:逆像がボレル集合になる集合の全体は σ-加法族で、開集合を含むためです)。ν(E)=m(E+a)\nu(E)=m(E+a) は B(Rd)\mathcal{B}(\mathbb{R}^d) 上の測度です。半開箱 P=∏i[ai,bi)P=\prod_i[a_i,b_i) の全体に ∅\emptyset を加えた族 P\mathcal{P} は共通部分で閉じた π-系で、補題 2.3 により開集合はすべて σ(P)\sigma(\mathcal{P}) に入るので σ(P)=B(Rd)\sigma(\mathcal{P})=\mathcal{B}(\mathbb{R}^d) です。P+aP+a は同じ辺の長さの半開箱なので、補題 2.2 から ν(P)=m(P)\nu(P)=m(P)。[−n,n)d∈P[-n,n)^d\in\mathcal{P} は増大して Rd\mathbb{R}^d を覆い測度は有限です。定理 4.5[測度空間と Carathéodory の拡張定理] により ν=m\nu=m です。

補題 3.2(並進不変測度の一意性)✓

μ\mu を B(Rd)\mathcal{B}(\mathbb{R}^d) 上の測度とし、すべてのボレル集合 EE と a∈Rda\in\mathbb{R}^d について μ(E+a)=μ(E)\mu(E+a)=\mu(E)、かつ c:=μ([0,1)d)<∞c:=\mu([0,1)^d)<\infty とする。このとき μ=c m\mu=c\,m である。

証明(補題 3.2)

第 kk 世代の 2 進立方体はどれも [0,2−k)d[0,2^{-k})^d の平行移動で、[0,1)d[0,1)^d はそのうち 2kd2^{kd} 個の互いに素な和です。並進不変性と加法性から c=2kdμ([0,2−k)d)c=2^{kd}\mu([0,2^{-k})^d)、つまり第 kk 世代の立方体 DD について μ(D)=c 2−kd=c m(D)\mu(D)=c\,2^{-kd}=c\,m(D) です。

P2\mathcal{P}_2 を、端点がすべて 2 進有理数(j2−kj2^{-k} の形)の半開箱の全体に ∅\emptyset を加えた族とします。共通部分で閉じているので π-系です。P∈P2P\in\mathcal{P}_2 は、端点の分母を 2k2^k にそろえれば第 kk 世代の 2 進立方体の有限個の互いに素な和なので、μ(P)=c m(P)\mu(P)=c\,m(P) です。補題 2.3 により σ(P2)=B(Rd)\sigma(\mathcal{P}_2)=\mathcal{B}(\mathbb{R}^d) で、[−n,n)d∈P2[-n,n)^d\in\mathcal{P}_2 の測度は有限です。定理 4.5[測度空間と Carathéodory の拡張定理] により μ=c m\mu=c\,m です(c=0c=0 でも同じ議論が通ります)。

この一意性が、次の定理の証明の中心です。「並進不変」と「単位立方体で有限」だけで測度が決まってしまうので、線形写像で写した測度も mm の定数倍になり、残るのは定数の計算だけになります。

補題 3.3(正則行列は基本行列の積)✓

正則行列 T∈Md(R)T\in M_d(\mathbb{R}) は、次の 3 種類の基本行列の有限個の積である:(I) 2 つの座標を入れ替える置換行列、(II) 1 つの座標を λ≠0\lambda\ne0 倍する対角行列、(III) I+cEpqI+cE_{pq}(p≠qp\ne q、EpqE_{pq} は (p,q)(p,q) 成分だけが 11 の行列)。

証明(補題 3.3)

命題 8.7[行列と連立一次方程式] により TT と行同値な簡約階段形は II で、補題 7.2[行列と連立一次方程式] の 1 により基本行列 E1,…,EkE_1,\ldots,E_k を使って Ek⋯E1T=IE_k\cdots E_1T=I と書けます。よって T=E1−1⋯Ek−1T=E_1^{-1}\cdots E_k^{-1} で、同補題の 2 により各 Ei−1E_i^{-1} も基本行列です。行変形の 3 つの型に対応する基本行列は、それぞれ (I)(II)(III) の形です。

定理 3.4(線形写像によるルベーグ測度の変換)✓

T∈Md(R)T\in M_d(\mathbb{R}) とする。

  1. det⁡T≠0\det T\ne0 ならば、任意のボレル集合 EE について T(E)T(E) はボレル集合で m(T(E))=∣det⁡T∣ m(E).m(T(E))=|\det T|\,m(E).
  2. det⁡T=0\det T=0 ならば、T(Rd)T(\mathbb{R}^d) は測度 00 の閉集合である。したがって任意の E⊆RdE\subseteq\mathbb{R}^d について T(E)T(E) は測度 00 のボレル集合に含まれ、上の等式は(両辺 00 として)成り立つ。
  3. 1 の等式はルベーグ可測集合 E∈LdE\in\mathcal{L}_d についても成り立ち、T(E)∈LdT(E)\in\mathcal{L}_d である(det⁡T≠0\det T\ne0 のとき)。
証明(定理 3.4)

1. TT は連続な逆写像 T−1T^{-1} をもつので、T(E)=(T−1)−1(E)T(E)=(T^{-1})^{-1}(E) はボレル集合です。μT(E)=m(T(E))\mu_T(E)=m(T(E)) とおくと、TT は単射なので互いに素な集合を互いに素な集合に写し、μT\mu_T は B(Rd)\mathcal{B}(\mathbb{R}^d) 上の測度です。T(E+a)=T(E)+TaT(E+a)=T(E)+Ta と 補題 3.1 から μT\mu_T は並進不変です。T([0,1)d)T([0,1)^d) は有界(∥Tx∥∞≤∥T∥∞\|Tx\|_\infty\le\|T\|_\infty)なので、ある立方体に含まれ μT([0,1)d)<∞\mu_T([0,1)^d)<\infty です。補題 3.2 により、定数 c(T)≥0c(T)\ge0 が存在して

m(T(E))=c(T) m(E)(E∈B(Rd)).m(T(E))=c(T)\,m(E)\qquad(E\in\mathcal{B}(\mathbb{R}^d)).

正則な S,TS,T について m(ST(E))=c(S)m(T(E))=c(S)c(T)m(E)m(ST(E))=c(S)m(T(E))=c(S)c(T)m(E) なので c(ST)=c(S)c(T)c(ST)=c(S)c(T) です。一方 ∣det⁡(ST)∣=∣det⁡S∣∣det⁡T∣|\det(ST)|=|\det S||\det T|(定理 6.1[行列式とその性質])です。補題 3.3 により、基本行列 TT について c(T)=∣det⁡T∣c(T)=|\det T| を示せば十分です。Q0=[0,1)dQ_0=[0,1)^d とします。

  • (I) 座標の入れ替え:T(Q0)=Q0T(Q_0)=Q_0 なので c(T)=1c(T)=1。行列式は −1-1 です(命題 4.3[行列式とその性質] の 1)。
  • (II) 第 pp 座標の λ\lambda 倍:T(Q0)T(Q_0) は第 pp 因子が長さ ∣λ∣|\lambda| の区間([0,λ)[0,\lambda) または (λ,0](\lambda,0])で他が [0,1)[0,1) の箱なので、補題 2.2 から c(T)=∣λ∣=∣det⁡T∣c(T)=|\lambda|=|\det T|。
  • (III) T=I+cEpqT=I+cE_{pq}、つまり (Tx)p=xp+cxq(Tx)_p=x_p+cx_q で他の座標は不変。置換行列 PP を、Pep=e1Pe_p=e_1、Peq=e2Pe_q=e_2 となるように選ぶと PEpqP−1=E12PE_{pq}P^{-1}=E_{12}、つまり P T P−1=I+cE12P\,T\,P^{-1}=I+cE_{12} です。PP と P−1P^{-1} は座標の入れ替え(互換)の積なので、(I) と乗法性から c(P)=c(P−1)=1c(P)=c(P^{-1})=1 で、 c(T)=c(I+cE12)c(T)=c(I+cE_{12}) です。T′=I+cE12T'=I+cE_{12} について T′(Q0)={y:y′∈[0,1)d−1, y1−cy2∈[0,1)}T'(Q_0)=\{y:y'\in[0,1)^{d-1},\ y_1-cy_2\in[0,1)\}(y′=(y2,…,yd)y'=(y_2,\ldots,y_d))です。md=m1×md−1m_d=m_1\times m_{d-1} に 定理 3.2[直積測度と Fubini–Tonelli の定理] の切断の式を使うと、y′∈[0,1)d−1y'\in[0,1)^{d-1} での切断は区間 [cy2,cy2+1)[cy_2,cy_2+1)(測度 11)、それ以外では空集合なので m(T′(Q0))=md−1([0,1)d−1)=1m(T'(Q_0))=m_{d-1}([0,1)^{d-1})=1。行列式は 11 です(命題 4.3[行列式とその性質] の 3)。

2. W=T(Rd)W=T(\mathbb{R}^d) は次元 k<dk<d の部分空間です。WW の基底 w1,…,wkw_1,\ldots,w_k を Rd\mathbb{R}^d の基底 w1,…,wdw_1,\ldots,w_d に延長し、Swi=eiSw_i=e_i で正則な SS を定めると、S(W)=span⁡(e1,…,ek)⊆H:={x:xd=0}S(W)=\operatorname{span}(e_1,\ldots,e_k)\subseteq H:=\{x:x_d=0\} です。HH は閉集合で m(H)=0m(H)=0(補題 2.2)、W=S−1(S(W))W=S^{-1}(S(W)) は閉集合で、1 から m(W)≤m(S−1(H))=∣det⁡S−1∣ m(H)=0m(W)\le m(S^{-1}(H))=|\det S^{-1}|\,m(H)=0 です。

3. E=B∪NE=B\cup N、N⊆ZN\subseteq Z(B,ZB,Z はボレル集合、m(Z)=0m(Z)=0)とすると、T(E)=T(B)∪T(N)T(E)=T(B)\cup T(N)、T(N)⊆T(Z)T(N)\subseteq T(Z)、m(T(Z))=∣det⁡T∣⋅0=0m(T(Z))=|\det T|\cdot0=0 なので T(E)∈LdT(E)\in\mathcal{L}_d で m(T(E))=m(T(B))=∣det⁡T∣m(B)=∣det⁡T∣m(E)m(T(E))=m(T(B))=|\det T|m(B)=|\det T|m(E) です。

系 3.5(回転不変性と相似)✓

  1. 直交行列 OO(OTO=IO^{\mathsf{T}}O=I)について m(O(E))=m(E)m(O(E))=m(E)。
  2. t>0t>0 について m(tE)=tdm(E)m(tE)=t^dm(E)。
  3. TT が正則、f≥0f\ge0 がボレル関数ならば ∫Rdf(Tx) dx=∣det⁡T∣−1∫Rdf(y) dy\int_{\mathbb{R}^d}f(Tx)\,dx=|\det T|^{-1}\int_{\mathbb{R}^d}f(y)\,dy。f∈L1f\in L^1 のときも同じ等式が成り立つ。
証明(系 3.5)
  1. 1=det⁡(OTO)=det⁡(OT)det⁡O=(det⁡O)21=\det(O^{\mathsf{T}}O)=\det(O^{\mathsf{T}})\det O=(\det O)^2(定理 6.1[行列式とその性質] と 命題 3.3[行列式とその性質] の転置不変性)なので ∣det⁡O∣=1|\det O|=1 です。2. det⁡(tI)=td\det(tI)=t^d です。3. f=1Af=\mathbf{1}_A なら f(Tx)=1T−1(A)(x)f(Tx)=\mathbf{1}_{T^{-1}(A)}(x) で、左辺は m(T−1(A))=∣det⁡T∣−1m(A)m(T^{-1}(A))=|\det T|^{-1}m(A) です。両辺とも ff について非負係数の線形性をもつので非負単関数で成り立ち、定理 3.7[ルベーグ積分の定義と収束定理] と単調収束定理(定理 5.2[ルベーグ積分の定義と収束定理])で一般の f≥0f\ge0 に移ります。f∈L1f\in L^1 なら f=f+−f−f=f^+-f^- に分けて適用します。

例 3.6(平行六面体の体積)

ベクトル v1,…,vdv_1,\ldots,v_d の張る平行六面体 {∑itivi:0≤ti≤1}\{\sum_it_iv_i:0\le t_i\le1\} は、viv_i を列とする行列 TT による [0,1]d[0,1]^d の像なので、体積は ∣det⁡(v1 ⋯ vd)∣|\det(v_1\ \cdots\ v_d)| です。d=2d=2 では 行列式 の 命題 8.1[行列式とその性質] の面積公式を、ルベーグ測度の言葉で一般次元に広げたことになります。viv_i が一次従属なら体積は 00 で、これは 2 に対応します。

例 3.7(正定値 2 次形式のガウス積分)

AA を正定値な実対称行列とします。系 4.3[スペクトル定理] により A=OΛOTA=O\Lambda O^{\mathsf{T}}(OO は直交行列、Λ=diag⁡(λ1,…,λd)\Lambda=\operatorname{diag}(\lambda_1,\ldots,\lambda_d)、λi>0\lambda_i>0)と書け、S=Λ1/2OTS=\Lambda^{1/2}O^{\mathsf{T}} とおくと ⟨Ax,x⟩=∣Sx∣2\langle Ax,x\rangle=|Sx|^2、∣det⁡S∣=(λ1⋯λd)1/2=(det⁡A)1/2|\det S|=(\lambda_1\cdots\lambda_d)^{1/2}=(\det A)^{1/2} です。系 3.5 の 3 と、e−∣y∣2=∏ie−yi2e^{-|y|^2}=\prod_ie^{-y_i^2} に Tonelli の定理(定理 3.3[直積測度と Fubini–Tonelli の定理])を d−1d-1 回使い、例 4.2[直積測度と Fubini–Tonelli の定理] の ∫Re−t2dt=π\int_{\mathbb{R}}e^{-t^2}dt=\sqrt\pi を使うと

∫Rde−⟨Ax,x⟩dx=1∣det⁡S∣∫Rde−∣y∣2dy=πd/2det⁡A.\int_{\mathbb{R}^d}e^{-\langle Ax,x\rangle}dx=\frac{1}{|\det S|}\int_{\mathbb{R}^d}e^{-|y|^2}dy=\frac{\pi^{d/2}}{\sqrt{\det A}}.

多変量正規分布の密度の正規化定数 (2π)−d/2(det⁡Σ)−1/2(2\pi)^{-d/2}(\det\Sigma)^{-1/2} は、A=12Σ−1A=\frac12\Sigma^{-1} としたこの式です。

4. C1C^1 微分同相に対する変数変換公式

Section titled “4. C1C^1C1 微分同相に対する変数変換公式”

定義 4.1(C¹ 微分同相)↗

U,V⊆RdU,V\subseteq\mathbb{R}^d を開集合とする。全単射 Φ:U→V\Phi:U\to V が C1C^1 微分同相であるとは、Φ\Phi と逆写像 Φ−1:V→U\Phi^{-1}:V\to U がともに C1C^1 級であることをいう。DΦ(x)D\Phi(x) は Φ\Phi のヤコビ行列、det⁡DΦ(x)\det D\Phi(x) をヤコビアンと呼ぶ。

Φ−1∘Φ=idU\Phi^{-1}\circ\Phi=\mathrm{id}_U に連鎖律(定理 9.4[多変数関数の微分と偏微分])を使うと D(Φ−1)(Φ(x)) DΦ(x)=ID(\Phi^{-1})(\Phi(x))\,D\Phi(x)=I なので、DΦ(x)D\Phi(x) はつねに正則です。逆に、次の補題により「単射でヤコビアンが消えない」ことを確かめれば微分同相が得られます。

補題 4.2(単射な局所微分同相は微分同相)✓

U⊆RdU\subseteq\mathbb{R}^d を開集合、Φ:U→Rd\Phi:U\to\mathbb{R}^d を単射な C1C^1 写像で、すべての x∈Ux\in U で det⁡DΦ(x)≠0\det D\Phi(x)\ne0 とする。このとき V=Φ(U)V=\Phi(U) は開集合で、Φ:U→V\Phi:U\to V は C1C^1 微分同相である。

証明(補題 4.2)

x∈Ux\in U に逆関数定理(定理 4.1[逆関数定理・陰関数定理と条件付き極値])を使うと、xx の開近傍 Vx⊆UV_x\subseteq U と Φ(x)\Phi(x) の開近傍 WxW_x があって、Φ∣Vx:Vx→Wx\Phi|_{V_x}:V_x\to W_x は全単射で逆写像 GxG_x は C1C^1 級です。Wx⊆Φ(U)W_x\subseteq\Phi(U) なので VV は開集合です。Ψ:V→U\Psi:V\to U を Φ\Phi の逆写像とすると、y∈Wxy\in W_x について Φ(Gx(y))=y=Φ(Ψ(y))\Phi(G_x(y))=y=\Phi(\Psi(y)) で、Φ\Phi の単射性から Ψ(y)=Gx(y)\Psi(y)=G_x(y) です。よって Ψ\Psi は各点の近傍で C1C^1 級です。

立方体の像を立方体で押さえるために、最大値ノルムでの平均値不等式を用意します。ユークリッドノルムでの形は 定理 9.7[多変数関数の微分と偏微分] ですが、そのまま使うとノルムの換算で d\sqrt d の因子が出て、下の証明の「(1+ε)d(1+\varepsilon)^d」が 11 に近づかなくなります。

補題 4.3(最大値ノルムでの平均値不等式)✓

U⊆RdU\subseteq\mathbb{R}^d を開集合、Φ:U→Rd\Phi:U\to\mathbb{R}^d を C1C^1 写像とし、x,y∈Ux,y\in U を結ぶ線分 [x,y][x,y] が UU に含まれるとする。このとき

∥Φ(x)−Φ(y)∥∞≤sup⁡z∈[x,y]∥DΦ(z)∥∞  ∥x−y∥∞.\|\Phi(x)-\Phi(y)\|_\infty\le\sup_{z\in[x,y]}\|D\Phi(z)\|_\infty\;\|x-y\|_\infty .
証明(補題 4.3)

成分 ii を固定し、g(t)=Φi(y+t(x−y))g(t)=\Phi_i(y+t(x-y))(0≤t≤10\le t\le1)とおきます。連鎖律により gg は微分可能で g′(t)=∑j∂jΦi(zt)(xj−yj)g'(t)=\sum_j\partial_j\Phi_i(z_t)(x_j-y_j)(zt=y+t(x−y)z_t=y+t(x-y))です。ラグランジュの平均値の定理(定理 3.3[平均値の定理とテイラーの定理])から、ある θ∈(0,1)\theta\in(0,1) で

∣Φi(x)−Φi(y)∣=∣g′(θ)∣≤∑j∣∂jΦi(zθ)∣ ∣xj−yj∣≤∥DΦ(zθ)∥∞∥x−y∥∞.|\Phi_i(x)-\Phi_i(y)|=|g'(\theta)|\le\sum_j|\partial_j\Phi_i(z_\theta)|\,|x_j-y_j|\le\|D\Phi(z_\theta)\|_\infty\|x-y\|_\infty .

ii について最大をとれば主張です。

補題 4.4(立方体の像の測度)✓

Φ:U→Rd\Phi:U\to\mathbb{R}^d を C1C^1 写像、Q=Q(a,r)⊆UQ=Q(a,r)\subseteq U を閉立方体、TT を正則行列とする。このとき Φ(Q)\Phi(Q) はコンパクトで

m(Φ(Q))≤∣det⁡T∣(sup⁡z∈Q∥T−1DΦ(z)∥∞)dm(Q).m(\Phi(Q))\le|\det T|\Bigl(\sup_{z\in Q}\|T^{-1}D\Phi(z)\|_\infty\Bigr)^d m(Q).
証明(補題 4.4)

Φ(Q)\Phi(Q) はコンパクト集合の連続像なのでコンパクトで、とくに閉集合なのでボレル集合です。G=T−1∘ΦG=T^{-1}\circ\Phi は C1C^1 級で、連鎖律から DG=T−1DΦDG=T^{-1}D\Phi です。右辺の上限を MM とします。QQ は凸なので、x∈Qx\in Q について [x,a]⊆Q[x,a]\subseteq Q で、補題 4.3 から ∥G(x)−G(a)∥∞≤M∥x−a∥∞≤Mr\|G(x)-G(a)\|_\infty\le M\|x-a\|_\infty\le Mr です。よって G(Q)⊆Q(G(a),Mr)G(Q)\subseteq Q(G(a),Mr) で、補題 2.2 から m(G(Q))≤(2Mr)d=Mdm(Q)m(G(Q))\le(2Mr)^d=M^dm(Q)。Φ(Q)=T(G(Q))\Phi(Q)=T(G(Q)) なので、定理 3.4 から m(Φ(Q))=∣det⁡T∣ m(G(Q))m(\Phi(Q))=|\det T|\,m(G(Q)) です。

この補題で T=DΦ(xj)T=D\Phi(x_j) と選ぶのが、証明の核心です。xjx_j の近くでは T−1DΦ(z)≈IT^{-1}D\Phi(z)\approx I なので、上限はほぼ 11 になり、「小さな立方体の像の体積 ≈∣det⁡DΦ(xj)∣×\approx|\det D\Phi(x_j)|\times 立方体の体積」が誤差 (1+ε)d(1+\varepsilon)^d の範囲で正確に言えます。

x_jQ を N^d 個の Q_j に分割ΦΦ(Q_j) ≈ Φ(x_j) + DΦ(x_j)(Q_j − x_j)Φ(Q)
変数変換公式の証明の第 1 段。立方体 Q を小立方体 Q_j に分けると、Φ は各 Q_j の上で線形写像 DΦ(x_j) でよく近似され、Φ(Q_j) はほぼ平行六面体になる。その体積は |det DΦ(x_j)| m(Q_j) の (1+ε)^d 倍以下に収まる。

定理 4.5(変数変換公式(ルベーグ積分版))✓

U,V⊆RdU,V\subseteq\mathbb{R}^d を開集合、Φ:U→V\Phi:U\to V を C1C^1 微分同相とする。

  1. 任意のボレル集合 E⊆UE\subseteq U について Φ(E)\Phi(E) はボレル集合で、m(Φ(E))=∫E∣det⁡DΦ(x)∣ dxm(\Phi(E))=\int_E|\det D\Phi(x)|\,dx。
  2. 任意のボレル関数 f:V→[0,∞]f:V\to[0,\infty] について ∫Vf(y) dy=∫Uf(Φ(x)) ∣det⁡DΦ(x)∣ dx.\int_Vf(y)\,dy=\int_Uf(\Phi(x))\,|\det D\Phi(x)|\,dx .
  3. ボレル関数 f:V→Rf:V\to\mathbb{R}(または C\mathbb{C})について、f∈L1(V)f\in L^1(V) と (f∘Φ) ∣det⁡DΦ∣∈L1(U)(f\circ\Phi)\,|\det D\Phi|\in L^1(U) は同値で、そのとき 2 の等式が成り立つ。

証明は 5 段で進みます。第 1〜3 段では、立方体、開集合、一般のボレル集合の順に、J=∣det⁡DΦ∣J=|\det D\Phi| として m(Φ(E))≤∫EJm(\Phi(E))\le\int_EJ を示します。第 1 段では、小さな立方体の上で Φ\Phi を線形写像 DΦ(xj)D\Phi(x_j) に置き換え、その誤差が ε\varepsilon で一様に抑えられるよう、立方体の分け方 NN を選びます。第 3 段の UkU_k は、UU の内側にあって有界で、縁から離れた開集合の増加列です。その上では JJ が最大値 CkC_k をもつので、誤差を評価できます。第 4 段で関数の不等式にし、第 5 段で逆写像 Φ−1\Phi^{-1} に同じ不等式を使って逆向きを得て、等号にします。

証明(定理 4.5)

J(x)=∣det⁡DΦ(x)∣J(x)=|\det D\Phi(x)| とおきます。JJ は連続(行列式は成分の多項式)です。E⊆UE\subseteq U がボレル集合なら、Φ(E)\Phi(E) は連続写像 Φ−1:V→U\Phi^{-1}:V\to U による EE の逆像なので VV の相対ボレル集合で、VV が開集合なので Rd\mathbb{R}^d のボレル集合です。

第 1 段:閉立方体 Q⊆UQ\subseteq U について m(Φ(Q))≤∫QJm(\Phi(Q))\le\int_QJ。 Q=Q(a,r)Q=Q(a,r) とし、0<ε<10<\varepsilon<1 を固定します。z↦DΦ(z)−1=(det⁡DΦ(z))−1adj⁡DΦ(z)z\mapsto D\Phi(z)^{-1}=(\det D\Phi(z))^{-1}\operatorname{adj}D\Phi(z)(命題 5.4[行列式とその性質])は連続なので、コンパクト集合 QQ 上で K=max⁡z∈Q∥DΦ(z)−1∥∞<∞K=\max_{z\in Q}\|D\Phi(z)^{-1}\|_\infty<\infty です(系 6.2[コンパクト性])。DΦD\Phi の各成分と JJ はコンパクト集合 QQ 上で一様連続(補題 2.2[Arzelà–Ascoli と Stone–Weierstrass])なので、δ>0\delta>0 を、z,z′∈Qz,z'\in Q、∥z−z′∥∞<δ\|z-z'\|_\infty<\delta ならば ∥DΦ(z)−DΦ(z′)∥∞≤ε/K\|D\Phi(z)-D\Phi(z')\|_\infty\le\varepsilon/K かつ ∣J(z)−J(z′)∣≤ε|J(z)-J(z')|\le\varepsilon となるように取れます。

r/N<δr/N<\delta となる自然数 NN を取り、QQ を辺の長さ 2r/N2r/N の NdN^d 個の閉立方体 Qj=Q(xj,r/N)Q_j=Q(x_j,r/N) に分けます。z∈Qjz\in Q_j なら ∥z−xj∥∞≤r/N<δ\|z-x_j\|_\infty\le r/N<\delta なので

∥DΦ(xj)−1DΦ(z)∥∞≤∥I∥∞+∥DΦ(xj)−1∥∞∥DΦ(z)−DΦ(xj)∥∞≤1+K⋅εK=1+ε.\|D\Phi(x_j)^{-1}D\Phi(z)\|_\infty\le\|I\|_\infty+\|D\Phi(x_j)^{-1}\|_\infty\|D\Phi(z)-D\Phi(x_j)\|_\infty\le1+K\cdot\frac{\varepsilon}{K}=1+\varepsilon .

補題 4.4 を T=DΦ(xj)T=D\Phi(x_j) で使い、J(xj)≤J(z)+εJ(x_j)\le J(z)+\varepsilon(z∈Qjz\in Q_j)を QjQ_j 上で積分すると

m(Φ(Qj))≤(1+ε)dJ(xj) m(Qj)≤(1+ε)d∫Qj(J+ε) dm.m(\Phi(Q_j))\le(1+\varepsilon)^dJ(x_j)\,m(Q_j)\le(1+\varepsilon)^d\int_{Q_j}(J+\varepsilon)\,dm .

QjQ_j どうしは境界でしか重ならず、境界は測度 00(補題 2.2)なので、∑j1Qj=1Q\sum_j\mathbf{1}_{Q_j}=\mathbf{1}_Q がほとんど至るところ成り立ち、∑j∫Qjg=∫Qg\sum_j\int_{Q_j}g=\int_Qg(g≥0g\ge0)です。Φ(Q)=⋃jΦ(Qj)\Phi(Q)=\bigcup_j\Phi(Q_j) と劣加法性から

m(Φ(Q))≤(1+ε)d(∫QJ dm+ε m(Q)).m(\Phi(Q))\le(1+\varepsilon)^d\Bigl(\int_QJ\,dm+\varepsilon\,m(Q)\Bigr).

ε→0\varepsilon\to0 として第 1 段の不等式を得ます。

第 2 段:開集合 W⊆UW\subseteq U について m(Φ(W))≤∫WJm(\Phi(W))\le\int_WJ。 補題 2.3 により W=⋃nDn‾W=\bigcup_n\overline{D_n}(DnD_n は互いに素な 2 進立方体、Dn‾⊆W\overline{D_n}\subseteq W は閉立方体)と書けます。第 1 段と、Dn‾∖Dn\overline{D_n}\setminus D_n が測度 00 であることから

m(Φ(W))≤∑nm(Φ(Dn‾))≤∑n∫Dn‾J=∑n∫DnJ=∫WJ.m(\Phi(W))\le\sum_nm(\Phi(\overline{D_n}))\le\sum_n\int_{\overline{D_n}}J=\sum_n\int_{D_n}J=\int_WJ .

第 3 段:ボレル集合 E⊆UE\subseteq U について m(Φ(E))≤∫EJm(\Phi(E))\le\int_EJ。 dist⁡∞(x,Uc)=inf⁡{∥x−y∥∞:y∉U}\operatorname{dist}_\infty(x,U^{c})=\inf\{\|x-y\|_\infty:y\notin U\}(U=RdU=\mathbb{R}^d なら +∞+\infty)とし、

Uk={x∈U:∥x∥∞<k, dist⁡∞(x,Uc)>1/k}U_k=\{x\in U:\|x\|_\infty<k,\ \operatorname{dist}_\infty(x,U^{c})>1/k\}

とおきます。距離関数は連続なので UkU_k は開集合で、その閉包は有界閉集合 {∥x∥∞≤k, dist⁡∞(x,Uc)≥1/k}⊆U\{\|x\|_\infty\le k,\ \operatorname{dist}_\infty(x,U^{c})\ge1/k\}\subseteq U に含まれます。そこで JJ の最大値を CkC_k とします。

まず E⊆UkE\subseteq U_k とします。m(E)<∞m(E)<\infty で、外正則性から開集合 W′⊇EW'\supseteq E で m(W′)<m(E)+ηm(W')<m(E)+\eta となるものがあります。W=W′∩UkW=W'\cap U_k は E⊆W⊆UkE\subseteq W\subseteq U_k を満たす開集合で m(W∖E)<ηm(W\setminus E)<\eta です。第 2 段から

m(Φ(E))≤m(Φ(W))≤∫WJ=∫EJ+∫W∖EJ≤∫EJ+Ckη,m(\Phi(E))\le m(\Phi(W))\le\int_WJ=\int_EJ+\int_{W\setminus E}J\le\int_EJ+C_k\eta ,

η→0\eta\to0 で m(Φ(E))≤∫EJm(\Phi(E))\le\int_EJ。一般の EE では E∩Uk↑EE\cap U_k\uparrow E なので、測度の下からの連続性(命題 2.1[ルベーグ積分の定義と収束定理])と単調収束定理で k→∞k\to\infty とします。

第 4 段:関数の不等式。 ボレル関数 f≥0f\ge0 について ∫Vf≤∫U(f∘Φ)J\int_Vf\le\int_U(f\circ\Phi)J を示します。f=1Af=\mathbf{1}_A(A⊆VA\subseteq V ボレル)なら、E=Φ−1(A)E=\Phi^{-1}(A) に第 3 段を使って ∫V1A=m(Φ(E))≤∫EJ=∫U(1A∘Φ)J\int_V\mathbf{1}_A=m(\Phi(E))\le\int_EJ=\int_U(\mathbf{1}_A\circ\Phi)J。非負単関数には線形性で、一般の ff には単関数の増加列(定理 3.7[ルベーグ積分の定義と収束定理])と単調収束定理(定理 5.2[ルベーグ積分の定義と収束定理])で移ります。f∘Φf\circ\Phi は連続写像との合成なのでボレル関数です。

第 5 段:逆向きの不等式と 1, 2。 第 4 段は任意の C1C^1 微分同相について成り立つので、Ψ=Φ−1:V→U\Psi=\Phi^{-1}:V\to U とボレル関数 g=(f∘Φ)J≥0g=(f\circ\Phi)J\ge0 に使えます。連鎖律から DΨ(y)=DΦ(Ψ(y))−1D\Psi(y)=D\Phi(\Psi(y))^{-1}、行列式の積の定理から ∣det⁡DΨ(y)∣=J(Ψ(y))−1|\det D\Psi(y)|=J(\Psi(y))^{-1} なので

∫U(f∘Φ)J≤∫Vg(Ψ(y)) ∣det⁡DΨ(y)∣ dy=∫Vf(y) J(Ψ(y)) J(Ψ(y))−1 dy=∫Vf.\int_U(f\circ\Phi)J\le\int_Vg(\Psi(y))\,|\det D\Psi(y)|\,dy=\int_Vf(y)\,J(\Psi(y))\,J(\Psi(y))^{-1}\,dy=\int_Vf .

第 4 段と合わせて 2 の等式です。1 は f=1Φ(E)f=\mathbf{1}_{\Phi(E)} の場合です。

3. 2 を ∣f∣|f| に使うと ∫V∣f∣=∫U∣f∘Φ∣J\int_V|f|=\int_U|f\circ\Phi|J なので可積分性は同値です。実数値なら f=f+−f−f=f^+-f^-、複素数値なら実部と虚部に分けて 2 を使えば等式が出ます。

系 4.6(ルベーグ可測集合と零集合の保存)✓

Φ:U→V\Phi:U\to V を C1C^1 微分同相とする。E∈LdE\in\mathcal{L}_d、E⊆UE\subseteq U ならば Φ(E)∈Ld\Phi(E)\in\mathcal{L}_d であり、m(E)=0m(E)=0 ならば m(Φ(E))=0m(\Phi(E))=0 である。定理 4.5 の 2, 3 は Ld\mathcal{L}_d-可測な関数 ff についても成り立つ。

証明(系 4.6)

E=B∪NE=B\cup N、N⊆ZN\subseteq Z(B,ZB,Z はボレル集合、m(Z)=0m(Z)=0)とし、Z⊆UZ\subseteq U としてよいとします(Z∩UZ\cap U に取り替えます)。定理 4.5 の 1 から m(Φ(Z))=∫ZJ=0m(\Phi(Z))=\int_ZJ=0 なので、Φ(E)=Φ(B)∪Φ(N)\Phi(E)=\Phi(B)\cup\Phi(N)、Φ(N)⊆Φ(Z)\Phi(N)\subseteq\Phi(Z) により Φ(E)∈Ld\Phi(E)\in\mathcal{L}_d で、m(Φ(E))=m(Φ(B))=∫BJ=∫EJm(\Phi(E))=m(\Phi(B))=\int_BJ=\int_EJ です。Φ−1\Phi^{-1} にも同じことが言えるので、Ld\mathcal{L}_d-可測な ff について {f∘Φ>t}=Φ−1({f>t})∈Ld\{f\circ\Phi>t\}=\Phi^{-1}(\{f>t\})\in\mathcal{L}_d です。あとは第 4 段と同じく、指示関数、単関数、単調収束定理の順に等式を広げます。

命題 4.7(1 次元の置換積分(二乗の場合))✓

ボレル関数 g:(0,∞)→[0,∞]g:(0,\infty)\to[0,\infty] について ∫0∞g(t) dt=∫0∞g(r2) 2r dr\displaystyle\int_0^\infty g(t)\,dt=\int_0^\infty g(r^2)\,2r\,dr が成り立つ。gg が非負連続なら、両辺はそれぞれ広義リーマン積分とも一致する。

証明(命題 4.7)

Φ(r)=r2\Phi(r)=r^2 は (0,∞)(0,\infty) から (0,∞)(0,\infty) への C1C^1 微分同相(逆写像 t\sqrt t は t>0t>0 で C1C^1 級)で、Φ′(r)=2r>0\Phi'(r)=2r>0 です。定理 4.5 から等式が従います。gg が非負連続なら、[1/n,n][1/n,n] 上ではリーマン積分とルベーグ積分が一致し(定理 8.1[ルベーグ積分の定義と収束定理])、n→∞n\to\infty で単調収束定理を使えば、ルベーグ積分は広義リーマン積分に等しくなります。右辺の g(r2) 2rg(r^2)\,2r も非負連続なので同様です。

例 4.8(1 次元の置換積分)

d=1d=1 で Φ(r)=r2\Phi(r)=r^2 は (0,∞)(0,\infty) から (0,∞)(0,\infty) への C1C^1 微分同相(逆写像 t\sqrt t は t>0t>0 で C1C^1 級)で、Φ′(r)=2r\Phi'(r)=2r です。定理 4.5 から、ボレル関数 g≥0g\ge0 について

∫0∞g(t) dt=∫0∞g(r2) 2r dr.\int_0^\infty g(t)\,dt=\int_0^\infty g(r^2)\,2r\,dr .

被積分関数が非負で連続なら、左右の積分は広義リーマン積分と一致します([1/n,n][1/n,n] 上ではリーマン積分とルベーグ積分が一致し(定理 8.1[ルベーグ積分の定義と収束定理])、n→∞n\to\infty で単調収束定理を使います)。

一方 Φ(t)=t3\Phi(t)=t^3 は (−1,1)(-1,1) から (−1,1)(-1,1) への C1C^1 級の全単射ですが、逆写像 t1/3t^{1/3} は 00 で微分可能でなく、微分同相ではありません。ヤコビアン 3t23t^2 が t=0t=0 で消えるためです。この場合も公式自体は成り立ちますが、それを保証するのは次節の例外集合を許す形です。

極座標 (r,θ)↦(rcos⁡θ,rsin⁡θ)(r,\theta)\mapsto(r\cos\theta,r\sin\theta) は、r=0r=0 ではヤコビアン rr が消え、θ=0\theta=0 と θ=2π\theta=2\pi で同じ点に行きます。実際の計算では、こうした「悪い集合」が測度 00 であることを使って無視します。その根拠が次の補題と系です。

補題 5.1(C¹ 写像は零集合を零集合に写す)✓

U⊆RdU\subseteq\mathbb{R}^d を開集合、Φ:U→Rd\Phi:U\to\mathbb{R}^d を C1C^1 写像、Z⊆UZ\subseteq U をボレル集合で m(Z)=0m(Z)=0 とする。このとき Φ(Z)\Phi(Z) は測度 00 のボレル集合に含まれる。ZZ が UU の相対閉集合なら Φ(Z)\Phi(Z) 自身がボレル集合で、m(Φ(Z))=0m(\Phi(Z))=0 である。

証明(補題 5.1)

補題 2.3 により U=⋃kQkU=\bigcup_kQ_k(Qk⊆UQ_k\subseteq U は閉立方体)と書けるので、各 kk について Φ(Z∩Qk)\Phi(Z\cap Q_k) が測度 00 のボレル集合に含まれることを示せば十分です。M=max⁡Qk∥DΦ∥∞M=\max_{Q_k}\|D\Phi\|_\infty とし、η>0\eta>0 とします。外正則性から開集合 W⊇Z∩QkW\supseteq Z\cap Q_k で m(W)<ηm(W)<\eta となるものを取り、補題 2.3 で W=⋃iCiW=\bigcup_iC_i(CiC_i は内部が互いに素な閉立方体、辺の長さ sis_i、∑im(Ci)=m(W)\sum_im(C_i)=m(W))と書きます。Ci∩Qk≠∅C_i\cap Q_k\ne\emptyset のとき、Ci∩QkC_i\cap Q_k は箱なので凸で Qk⊆UQ_k\subseteq U に含まれ、その点 pip_i を一つ選ぶと、補題 4.3 から x∈Ci∩Qkx\in C_i\cap Q_k について ∥Φ(x)−Φ(pi)∥∞≤Msi\|\Phi(x)-\Phi(p_i)\|_\infty\le Ms_i です。したがって

Φ(Z∩Qk)⊆⋃iQ(Φ(pi),Msi)=:Bη,m(Bη)≤∑i(2Msi)d=(2M)d m(W)<(2M)dη.\Phi(Z\cap Q_k)\subseteq\bigcup_iQ(\Phi(p_i),Ms_i)=:B_\eta,\qquad m(B_\eta)\le\sum_i(2Ms_i)^d=(2M)^d\,m(W)<(2M)^d\eta .

B=⋂nB1/nB=\bigcap_nB_{1/n} は Φ(Z∩Qk)\Phi(Z\cap Q_k) を含む測度 00 のボレル集合です。

ZZ が UU の相対閉集合なら、Z∩QkZ\cap Q_k はコンパクトで、Φ(Z)=⋃kΦ(Z∩Qk)\Phi(Z)=\bigcup_k\Phi(Z\cap Q_k) はコンパクト集合の可算和なのでボレル集合です。

系 5.2(例外集合を許す変数変換公式)✓

U⊆RdU\subseteq\mathbb{R}^d を開集合、Φ:U→Rd\Phi:U\to\mathbb{R}^d を C1C^1 写像とする。UU の相対閉集合 ZZ で m(Z)=0m(Z)=0 となるものがあり、U∖ZU\setminus Z の上で Φ\Phi は単射かつ det⁡DΦ≠0\det D\Phi\ne0 であるとする。このとき Φ(U)\Phi(U) はボレル集合で、任意のボレル関数 f:Rd→[0,∞]f:\mathbb{R}^d\to[0,\infty] について

∫Φ(U)f(y) dy=∫Uf(Φ(x)) ∣det⁡DΦ(x)∣ dx.\int_{\Phi(U)}f(y)\,dy=\int_Uf(\Phi(x))\,|\det D\Phi(x)|\,dx .

ff が実数値(複素数値)なら、片方の積分が絶対収束すれば他方も絶対収束して等号が成り立つ。

証明(系 5.2)

U0=U∖ZU_0=U\setminus Z は開集合です。補題 4.2 により V0=Φ(U0)V_0=\Phi(U_0) は開集合で、Φ:U0→V0\Phi:U_0\to V_0 は C1C^1 微分同相です。定理 4.5 から ∫V0f=∫U0(f∘Φ)∣det⁡DΦ∣\int_{V_0}f=\int_{U_0}(f\circ\Phi)|\det D\Phi| で、m(Z)=0m(Z)=0 なので右辺は ∫U(f∘Φ)∣det⁡DΦ∣\int_U(f\circ\Phi)|\det D\Phi| に等しくなります。一方 Φ(U)=V0∪Φ(Z)\Phi(U)=V_0\cup\Phi(Z) で、補題 5.1 により Φ(Z)\Phi(Z) は測度 00 のボレル集合です。よって Φ(U)\Phi(U) はボレル集合で ∫Φ(U)f=∫V0f\int_{\Phi(U)}f=\int_{V_0}f です。後半は ∣f∣|f| に適用してから実部・虚部の正負の部分に分けます。

例 5.3(平面の極座標)

開集合 U=(0,∞)×(0,2π)U=(0,\infty)\times(0,2\pi) と Φ(r,θ)=(rcos⁡θ,rsin⁡θ)\Phi(r,\theta)=(r\cos\theta,r\sin\theta) を考えます(θ\theta の範囲に端点を含めると UU は開集合でなくなり、補題 4.2 が使えません)。det⁡DΦ=r>0\det D\Phi=r>0 で、Φ\Phi は UU 上で単射です(r=∣Φ∣r=|\Phi| が決まり、θ∈(0,2π)\theta\in(0,2\pi) は (cos⁡θ,sin⁡θ)(\cos\theta,\sin\theta) から一つに決まります)。像は R2∖{(x,0):x≥0}\mathbb{R}^2\setminus\{(x,0):x\ge0\} で、除かれた半直線は直線 {y=0}\{y=0\} に含まれるので測度 00 です(補題 2.2)。系 5.2(ここでは Z=∅Z=\emptyset)から、ボレル関数 f≥0f\ge0 について

∫R2f dm=∫0∞ ⁣ ⁣∫02πf(rcos⁡θ,rsin⁡θ) r dθ dr.\int_{\mathbb{R}^2}f\,dm=\int_0^\infty\!\!\int_0^{2\pi}f(r\cos\theta,r\sin\theta)\,r\,d\theta\,dr .

右辺の累次積分への書き換えは Tonelli の定理(定理 3.3[直積測度と Fubini–Tonelli の定理])です。積分範囲を閉じた [0,∞)×[0,2π][0,\infty)\times[0,2\pi] と書いても、付け加わった境界は測度 00 なので値は変わりません。このように、悪い点を定義域から除いて像の側の測度 00 の集合に押し付けるのが、例外集合の典型的な扱い方です。

6. Rd\mathbb{R}^d の極座標と球面測度

Section titled “6. Rd\mathbb{R}^dRd の極座標と球面測度”

Sd−1={x∈Rd:∣x∣=1}S^{d-1}=\{x\in\mathbb{R}^d:|x|=1\} とし、そのボレル集合族(Rd\mathbb{R}^d からの相対位相に関するもの)を B(Sd−1)\mathcal{B}(S^{d-1}) と書きます。Sd−1S^{d-1} は閉集合なので、B(Sd−1)\mathcal{B}(S^{d-1}) の元は Rd\mathbb{R}^d のボレル集合です({A⊆Sd−1:A∈B(Rd)}\{A\subseteq S^{d-1}:A\in\mathcal{B}(\mathbb{R}^d)\} は相対開集合を含む σ-加法族だからです)。写像

Ψ:Rd∖{0}→(0,∞)×Sd−1,Ψ(x)=(∣x∣,x∣x∣)\Psi:\mathbb{R}^d\setminus\{0\}\to(0,\infty)\times S^{d-1},\qquad\Psi(x)=\Bigl(|x|,\frac{x}{|x|}\Bigr)

は連続な逆写像 (r,ω)↦rω(r,\omega)\mapsto r\omega をもつ同相写像です。A∈B(Sd−1)A\in\mathcal{B}(S^{d-1}) と t>0t>0 について、「AA の上の錐」を

At={rω:0<r≤t, ω∈A}=Ψ−1((0,t]×A)A_t=\{r\omega:0<r\le t,\ \omega\in A\}=\Psi^{-1}\bigl((0,t]\times A\bigr)

とおきます。これはボレル集合で、At=t A1A_t=t\,A_1 です。

定義 6.1(球面測度)↗

A∈B(Sd−1)A\in\mathcal{B}(S^{d-1}) について

σ(A)=d⋅m(A1)=d⋅m({rω:0<r≤1, ω∈A})\sigma(A)=d\cdot m(A_1)=d\cdot m\bigl(\{r\omega:0<r\le1,\ \omega\in A\}\bigr)

と定める。σ\sigma を Sd−1S^{d-1} 上の球面測度(面積測度)という。

A↦A1A\mapsto A_1 は互いに素な集合を互いに素な集合に写し、可算和を保つので、σ\sigma は B(Sd−1)\mathcal{B}(S^{d-1}) 上の測度です。σ(Sd−1)=d m({0<∣x∣≤1})≤d 2d<∞\sigma(S^{d-1})=d\,m(\{0<|x|\le1\})\le d\,2^d<\infty です。係数 dd は、半径 11 の球の体積を「面積 ×\times 半径 ÷d\div d」と見る(錐の体積の公式)ための正規化で、次の定理で rd−1drr^{d-1}dr と組み合わせたときにちょうど辻褄が合います。

たとえば d=2d=2、A=S1A=S^1 なら、A1A_1 は原点を除いた単位円板で m(A1)=πm(A_1)=\pi なので、σ(S1)=2π\sigma(S^1)=2\pi となり、円周の長さと一致します。

0A ⊂ S^(d−1)a半径 1濃い部分は半径 a から 1 の殻
球面測度の定義。球面上の集合 A の上の錐 A_1 の体積に d を掛けたものを σ(A) とする。極座標の定理の証明では、殻状の集合 {rω : a < r ≤ b, ω ∈ A}(A_b から A_a を除いた集合)の体積が (b^d − a^d) m(A_1) であることを使う。

補題 6.2(極座標側のボレル集合族)✓

B((0,∞)×Sd−1)=B((0,∞))⊗B(Sd−1)\mathcal{B}\bigl((0,\infty)\times S^{d-1}\bigr)=\mathcal{B}((0,\infty))\otimes\mathcal{B}(S^{d-1}) であり、この σ-加法族は {(a,b]×A:0<a<b, A∈B(Sd−1)}∪{∅}\{(a,b]\times A:0<a<b,\ A\in\mathcal{B}(S^{d-1})\}\cup\{\emptyset\} という π-系で生成される。

証明(補題 6.2)

(⊇\supseteq) 射影は連続なので、命題 2.3[直積測度と Fubini–Tonelli の定理] の証明と同じく、I×Sd−1I\times S^{d-1} と (0,∞)×A(0,\infty)\times A はボレル集合で、その共通部分 I×AI\times A もボレル集合です。(⊆\subseteq) (0,∞)(0,\infty) は有理数端点の開区間、Sd−1S^{d-1} は有理数端点の開箱との共通部分という可算な開基をもちます。直積の開集合 GG の各点は、開基の元の直積で GG に含まれるものの中にあるので、GG はそのような直積の可算和で、B((0,∞))⊗B(Sd−1)\mathcal{B}((0,\infty))\otimes\mathcal{B}(S^{d-1}) に入ります。生成については、(a,b]×A(a,b]\times A の形の集合を含む σ-加法族 S\mathcal{S} を考えると、{I⊆(0,∞):I×A∈S}\{I\subseteq(0,\infty):I\times A\in\mathcal{S}\} は σ-加法族で、開区間 (a,b)=⋃n(a,b−1/n](a,b)=\bigcup_n(a,b-1/n] を含むので (0,∞)(0,\infty) のボレル集合をすべて含みます。よって S\mathcal{S} は可測長方形をすべて含みます。

次の定理のねらいは、Rd\mathbb{R}^d の点を(原点からの距離,方向)の組に対応させる写像 Ψ\Psi で、体積を測る測度 mm を移すと、距離の方向の測度 ρ\rho と球面の方向の測度 σ\sigma の直積測度になる、と言うことです。直積測度になるので、積分は距離と方向で別々に計算できます。

定理 6.3(Rdℝ^d の極座標分解)✓

(0,∞)(0,\infty) 上の測度 ρ\rho を ρ(I)=∫Ird−1 dr\rho(I)=\int_Ir^{d-1}\,dr で定める。B(Rd∖{0})\mathcal{B}(\mathbb{R}^d\setminus\{0\}) 上の mm を Ψ\Psi で写した測度 m∘Ψ−1m\circ\Psi^{-1} は、直積測度 ρ×σ\rho\times\sigma に等しい。したがって、ボレル関数 f:Rd→[0,∞]f:\mathbb{R}^d\to[0,\infty] について、r↦∫Sd−1f(rω) dσ(ω)r\mapsto\int_{S^{d-1}}f(r\omega)\,d\sigma(\omega) はボレル関数で

∫Rdf(x) dx=∫0∞(∫Sd−1f(rω) dσ(ω))rd−1 dr.\int_{\mathbb{R}^d}f(x)\,dx=\int_0^\infty\Bigl(\int_{S^{d-1}}f(r\omega)\,d\sigma(\omega)\Bigr)r^{d-1}\,dr .

f∈L1(Rd)f\in L^1(\mathbb{R}^d) なら、ほとんどすべての r>0r>0 で f(r ⋅)∈L1(σ)f(r\,\cdot)\in L^1(\sigma) で、同じ等式が成り立つ。

証明(定理 6.3)

測度の一致。 0<a<b0<a<b、A∈B(Sd−1)A\in\mathcal{B}(S^{d-1}) について、Ψ−1((a,b]×A)=Ab∖Aa\Psi^{-1}((a,b]\times A)=A_b\setminus A_a で Aa⊆AbA_a\subseteq A_b です。At=tA1A_t=tA_1 と 系 3.5 の 2 から m(At)=tdm(A1)m(A_t)=t^dm(A_1) なので

m(Ψ−1((a,b]×A))=(bd−ad) m(A1)=∫abrd−1dr⋅d m(A1)=ρ((a,b]) σ(A).m\bigl(\Psi^{-1}((a,b]\times A)\bigr)=(b^d-a^d)\,m(A_1)=\int_a^br^{d-1}dr\cdot d\,m(A_1)=\rho((a,b])\,\sigma(A).

右辺は 定理 3.2[直積測度と Fubini–Tonelli の定理] により (ρ×σ)((a,b]×A)(\rho\times\sigma)((a,b]\times A) です(ρ\rho も σ\sigma も σ-有限)。補題 6.2 の π-系の上で二つの測度は一致し、(1/n,n]×Sd−1(1/n,n]\times S^{d-1} は増大して全体を覆い測度有限です。定理 4.5[測度空間と Carathéodory の拡張定理] により m∘Ψ−1=ρ×σm\circ\Psi^{-1}=\rho\times\sigma です。

積分の式。 m({0})=0m(\{0\})=0 なので ∫Rdf=∫Rd∖{0}f\int_{\mathbb{R}^d}f=\int_{\mathbb{R}^d\setminus\{0\}}f です。F(r,ω)=f(rω)F(r,\omega)=f(r\omega) は連続写像 Ψ−1\Psi^{-1} とボレル関数の合成なので、補題 6.2 により B((0,∞))⊗B(Sd−1)\mathcal{B}((0,\infty))\otimes\mathcal{B}(S^{d-1})-可測です。指示関数 f=1Bf=\mathbf{1}_B については ∫f dm=m(Ψ−1(Ψ(B∖{0})))=(ρ×σ)(Ψ(B∖{0}))=∫F d(ρ×σ)\int f\,dm=m(\Psi^{-1}(\Psi(B\setminus\{0\})))=(\rho\times\sigma)(\Psi(B\setminus\{0\}))=\int F\,d(\rho\times\sigma) で、単関数と単調収束定理により一般の f≥0f\ge0 で ∫Rdf dm=∫F d(ρ×σ)\int_{\mathbb{R}^d}f\,dm=\int F\,d(\rho\times\sigma) です。Tonelli の定理(定理 3.3[直積測度と Fubini–Tonelli の定理])で右辺を ∫0∞(∫F(r,ω) dσ)dρ(r)\int_0^\infty\bigl(\int F(r,\omega)\,d\sigma\bigr)d\rho(r) と書き、最後に ∫g dρ=∫0∞g(r)rd−1dr\int g\,d\rho=\int_0^\infty g(r)r^{d-1}dr(指示関数では ρ\rho の定義そのもの。単関数と単調収束定理で一般の g≥0g\ge0 へ)を使います。L1L^1 の場合は Fubini の定理(定理 3.4[直積測度と Fubini–Tonelli の定理])を使います。

この定理は「Rd∖{0}\mathbb{R}^d\setminus\{0\} は、測度まで込めて、半直線 (0,∞)(0,\infty)(密度 rd−1r^{d-1})と球面 Sd−1S^{d-1}(測度 σ\sigma)の直積である」と言っています。角度座標を一切使っていないので、どの次元でも同じ形で成り立ちます。

定理はボレル関数について述べましたが、ルベーグ可測な ff にも使えます。f≥0f\ge0 がルベーグ可測なら、ボレル関数 hh と測度 00 のボレル集合 ZZ で f=hf=h(ZZ の外)となるものがあります。Ψ(Z∖{0})\Psi(Z\setminus\{0\}) は ρ×σ\rho\times\sigma-零集合なので(上の測度の一致)、Tonelli の定理により、ほとんどすべての rr で {ω:rω∈Z}\{\omega:r\omega\in Z\} は σ\sigma-零集合です。そのような rr では f(r ⋅)=h(r ⋅)f(r\,\cdot)=h(r\,\cdot) が σ\sigma-ほとんど至るところ成り立つので、∫f dm=∫h dm\int f\,dm=\int h\,dm と内側の積分の値は、ff と hh で一致します。L1L^1 の場合も同様です。

命題 6.4(球面測度の回転不変性)✓

OO を dd 次の直交行列とする。A∈B(Sd−1)A\in\mathcal{B}(S^{d-1}) について O(A)∈B(Sd−1)O(A)\in\mathcal{B}(S^{d-1}) かつ σ(O(A))=σ(A)\sigma(O(A))=\sigma(A) であり、ボレル関数 g:Sd−1→[0,∞]g:S^{d-1}\to[0,\infty](または g∈L1(σ)g\in L^1(\sigma))について ∫Sd−1g(Oω) dσ(ω)=∫Sd−1g dσ\int_{S^{d-1}}g(O\omega)\,d\sigma(\omega)=\int_{S^{d-1}}g\,d\sigma である。

証明(命題 6.4)

∣Ox∣=∣x∣|Ox|=|x| なので OO は Sd−1S^{d-1} をそれ自身に同相に写し、O(A)O(A) はボレル集合です。O(rω)=rOωO(r\omega)=rO\omega から (O(A))1=O(A1)(O(A))_1=O(A_1) で、系 3.5 の 1 により σ(O(A))=d m(O(A1))=d m(A1)=σ(A)\sigma(O(A))=d\,m(O(A_1))=d\,m(A_1)=\sigma(A) です。関数については、g=1Ag=\mathbf{1}_A なら g∘O=1O−1(A)g\circ O=\mathbf{1}_{O^{-1}(A)} で、O−1O^{-1} も直交行列なので両辺は等しく、単関数・単調収束定理・正負の部分への分解で一般の gg に移ります。

系 6.5(動径関数の積分)✓

ボレル関数 g:(0,∞)→[0,∞]g:(0,\infty)\to[0,\infty] について

∫Rdg(∣x∣) dx=σ(Sd−1)∫0∞g(r) rd−1 dr.\int_{\mathbb{R}^d}g(|x|)\,dx=\sigma(S^{d-1})\int_0^\infty g(r)\,r^{d-1}\,dr .
証明(系 6.5)

定理 6.3 を f(x)=g(∣x∣)f(x)=g(|x|)(x≠0x\ne0)に使うと、内側の積分は ∫Sd−1g(r) dσ=g(r) σ(Sd−1)\int_{S^{d-1}}g(r)\,d\sigma=g(r)\,\sigma(S^{d-1}) です。

命題 6.6(べき関数の可積分性)✓

a∈Ra\in\mathbb{R} とする。

∫∣x∣≤1∣x∣−adx<∞  ⟺  a<d,∫∣x∣>1∣x∣−adx<∞  ⟺  a>d.\int_{|x|\le1}|x|^{-a}dx<\infty\iff a<d,\qquad\int_{|x|>1}|x|^{-a}dx<\infty\iff a>d .
証明(命題 6.6)

系 6.5 から ∫∣x∣≤1∣x∣−adx=σ(Sd−1)∫01rd−1−adr\int_{|x|\le1}|x|^{-a}dx=\sigma(S^{d-1})\int_0^1r^{d-1-a}dr、∫∣x∣>1∣x∣−adx=σ(Sd−1)∫1∞rd−1−adr\int_{|x|>1}|x|^{-a}dx=\sigma(S^{d-1})\int_1^\infty r^{d-1-a}dr です。σ(Sd−1)∈(0,∞)\sigma(S^{d-1})\in(0,\infty) なので、右辺の積分の有限性を見ればよく、非負連続な被積分関数のルベーグ積分は、区間 [1/n,1][1/n,1] や [1,n][1,n] でのリーマン積分(定理 8.1[ルベーグ積分の定義と収束定理])の n→∞n\to\infty での極限(単調収束定理)、すなわち広義リーマン積分に等しいので、命題 2.4[広義積分・パラメータ積分とガンマ関数] を p=a+1−dp=a+1-d に使えます。∫01\int_0^1 が有限なのは p<1p<1、すなわち a<da<d のときに限り、∫1∞\int_1^\infty が有限なのは p>1p>1、すなわち a>da>d のときに限ります。

例 6.7(べき関数の可積分性)

a∈Ra\in\mathbb{R} とします。系 6.5 から

∫∣x∣≤1∣x∣−adx=σ(Sd−1)∫01rd−1−adr<∞  ⟺  a<d,∫∣x∣>1∣x∣−adx<∞  ⟺  a>d.\int_{|x|\le1}|x|^{-a}dx=\sigma(S^{d-1})\int_0^1r^{d-1-a}dr<\infty\iff a<d,\qquad\int_{|x|>1}|x|^{-a}dx<\infty\iff a>d .

1 次元の積分 ∫01r−pdr\int_0^1r^{-p}dr、∫1∞r−pdr\int_1^\infty r^{-p}dr の収束判定(例 2.5[広義積分・パラメータ積分とガンマ関数])に帰着しました。原点での特異性は「次元 dd より弱ければ可積分」、無限遠での減衰は「次元 dd より速ければ可積分」です。

  • d≥3d\ge3 のラプラス方程式の基本解 ∣x∣2−d|x|^{2-d} は a=d−2<da=d-2<d なので局所可積分です。
  • ⟨x⟩−s=(1+∣x∣2)−s/2\langle x\rangle^{-s}=(1+|x|^2)^{-s/2} は、∣x∣≤1|x|\le1 で有界、∣x∣≥1|x|\ge1 で 2−s/2∣x∣−s≤⟨x⟩−s≤∣x∣−s2^{-s/2}|x|^{-s}\le\langle x\rangle^{-s}\le|x|^{-s}(s≥0s\ge0)なので、Rd\mathbb{R}^d 上可積分であることと s>ds>d は同値です。

定理 6.8(球の体積と球面の面積)✓

d≥1d\ge1 について

σ(Sd−1)=2πd/2Γ(d/2),m({∣x∣<R})=m({∣x∣≤R})=πd/2Γ(d/2+1)Rd.\sigma(S^{d-1})=\frac{2\pi^{d/2}}{\Gamma(d/2)},\qquad m(\{|x|<R\})=m(\{|x|\le R\})=\frac{\pi^{d/2}}{\Gamma(d/2+1)}R^d .

ここで Γ\Gamma はガンマ関数(定義 5.1[広義積分・パラメータ積分とガンマ関数])である。

証明(定理 6.8)

ガウス積分を二通りに計算する。 e−∣x∣2=∏i=1de−xi2e^{-|x|^2}=\prod_{i=1}^de^{-x_i^2} に Tonelli の定理を d−1d-1 回使い、∫Re−t2dt=π\int_{\mathbb{R}}e^{-t^2}dt=\sqrt\pi(定理 3.3[広義積分・パラメータ積分とガンマ関数]。非負連続関数なので、定理 8.1[ルベーグ積分の定義と収束定理] と単調収束定理により広義リーマン積分とルベーグ積分は一致します)から ∫Rde−∣x∣2dx=πd/2\int_{\mathbb{R}^d}e^{-|x|^2}dx=\pi^{d/2} です。他方 系 6.5 から

∫Rde−∣x∣2dx=σ(Sd−1)∫0∞e−r2rd−1dr.\int_{\mathbb{R}^d}e^{-|x|^2}dx=\sigma(S^{d-1})\int_0^\infty e^{-r^2}r^{d-1}dr .

命題 4.7 を g(t)=12e−ttd/2−1g(t)=\frac12e^{-t}t^{d/2-1} に使うと(g(r2) 2r=e−r2rd−1g(r^2)\,2r=e^{-r^2}r^{d-1})、右辺の積分は 12∫0∞e−ttd/2−1dt=12Γ(d/2)\frac12\int_0^\infty e^{-t}t^{d/2-1}dt=\frac12\Gamma(d/2) です(被積分関数は非負連続なので、広義リーマン積分で定義された Γ\Gamma と一致します)。よって πd/2=σ(Sd−1)Γ(d/2)/2\pi^{d/2}=\sigma(S^{d-1})\Gamma(d/2)/2。

体積。 系 6.5 を g=1(0,R)g=\mathbf{1}_{(0,R)} と g=1(0,R]g=\mathbf{1}_{(0,R]} に使うと、どちらも m=σ(Sd−1)Rd/dm=\sigma(S^{d-1})R^d/d です(m({0})=0m(\{0\})=0)。関数等式 d2Γ(d2)=Γ(d2+1)\frac d2\Gamma(\frac d2)=\Gamma(\frac d2+1)(定理 5.3[広義積分・パラメータ積分とガンマ関数])から主張の形になります。

Γ(1/2)=π\Gamma(1/2)=\sqrt\pi、Γ(1)=1\Gamma(1)=1、Γ(3/2)=π/2\Gamma(3/2)=\sqrt\pi/2、Γ(2)=1\Gamma(2)=1、Γ(5/2)=3π/4\Gamma(5/2)=3\sqrt\pi/4、Γ(3)=2\Gamma(3)=2 を代入すると次の表になります(体積は Γ(d/2+1)\Gamma(d/2+1) なので d=3d=3 で Γ(5/2)\Gamma(5/2)、d=4d=4 で Γ(3)\Gamma(3) が現れます)。

ddσ(Sd−1)\sigma(S^{d-1})単位球の体積
122(2 点 ±1\pm1 の数え上げ)22(区間 [−1,1][-1,1] の長さ)
22π2\pi(円周の長さ)π\pi
34π4\pi4π/34\pi/3
42π22\pi^2π2/2\pi^2/2

d=1d=1 では S0={±1}S^0=\{\pm1\} で、σ({1})=1⋅m((0,1])=1\sigma(\{1\})=1\cdot m((0,1])=1 なので σ\sigma は数え上げ測度です。極座標の式は ∫Rf=∫0∞(f(r)+f(−r)) dr\int_{\mathbb{R}}f=\int_0^\infty(f(r)+f(-r))\,dr となり、当たり前の式に戻ります。

注意 6.9(高次元の球はやせる)

単位球の体積 πd/2/Γ(d/2+1)\pi^{d/2}/\Gamma(d/2+1) は d=5d=5 で最大(約 5.265.26)になり、その後 d→∞d\to\infty で 00 に収束します。分母の Γ(d/2+1)\Gamma(d/2+1) がスターリングの公式(定理 7.1[広義積分・パラメータ積分とガンマ関数])により (d/2e)d/2(d/2e)^{d/2} 程度で増え、分子の πd/2\pi^{d/2} を圧倒するからです。単位立方体 [−1,1]d[-1,1]^d の体積 2d2^d との比は急速に 00 に近づき、「高次元の立方体の体積はほとんど角に集まっている」ことを示しています。

0246d=5(約 5.26)12345678910次元 d
次元 d と単位球の体積 π^(d/2)/Γ(d/2+1)。d = 5 で最大(約 5.26)になり、その後は 0 に向かって減っていく。

6.5. 3 次元の球座標と dσ=sin⁡θ dθ dφd\sigma=\sin\theta\,d\theta\,d\varphi

Section titled “6.5. 3 次元の球座標と dσ=sin⁡θ dθ dφd\sigma=\sin\theta\,d\theta\,d\varphidσ=sinθdθdφ”

命題 6.10(S² の球面測度の球座標表示)✓

ω(θ,φ)=(sin⁡θcos⁡φ, sin⁡θsin⁡φ, cos⁡θ)\omega(\theta,\varphi)=(\sin\theta\cos\varphi,\ \sin\theta\sin\varphi,\ \cos\theta) とする。ボレル関数 g:S2→[0,∞]g:S^2\to[0,\infty] について

∫S2g dσ=∫02π ⁣ ⁣∫0πg(ω(θ,φ)) sin⁡θ dθ dφ.\int_{S^2}g\,d\sigma=\int_0^{2\pi}\!\!\int_0^\pi g(\omega(\theta,\varphi))\,\sin\theta\,d\theta\,d\varphi .

とくに σ(S2)=4π\sigma(S^2)=4\pi で、g∈L1(σ)g\in L^1(\sigma) についても同じ等式が成り立つ。

証明(命題 6.10)

球座標の変数変換。 U=(0,∞)×(0,π)×(0,2π)U=(0,\infty)\times(0,\pi)\times(0,2\pi) 上で Φ(r,θ,φ)=r ω(θ,φ)\Phi(r,\theta,\varphi)=r\,\omega(\theta,\varphi) とします。ヤコビ行列の第 3 行 (cos⁡θ, −rsin⁡θ, 0)(\cos\theta,\,-r\sin\theta,\,0) に沿って余因子展開(定理 5.2[行列式とその性質])すると、(3,1)(3,1) 余因子は r2sin⁡θcos⁡θr^2\sin\theta\cos\theta、(3,2)(3,2) 余因子は −rsin⁡2θ-r\sin^2\theta なので

det⁡DΦ=cos⁡θ⋅r2sin⁡θcos⁡θ+(−rsin⁡θ)(−rsin⁡2θ)=r2sin⁡θ>0.\det D\Phi=\cos\theta\cdot r^2\sin\theta\cos\theta+(-r\sin\theta)(-r\sin^2\theta)=r^2\sin\theta>0 .

Φ\Phi は UU 上で単射です:r=∣Φ∣r=|\Phi|、cos⁡θ=Φ3/r\cos\theta=\Phi_3/r で θ∈(0,π)\theta\in(0,\pi) が決まり、sin⁡θ>0\sin\theta>0 なので (cos⁡φ,sin⁡φ)(\cos\varphi,\sin\varphi) が決まり φ∈(0,2π)\varphi\in(0,2\pi) が決まります。同じ計算から、Φ(U)\Phi(U) は R3\mathbb{R}^3 から半平面 H={(x,0,z):x≥0}H=\{(x,0,z):x\ge0\} を除いた集合で、H⊆{y=0}H\subseteq\{y=0\} は測度 00 です(補題 2.2)。系 5.2(Z=∅Z=\emptyset)から、ボレル関数 f≥0f\ge0 について

∫R3f dm=∫Uf(r ω(θ,φ)) r2sin⁡θ  d(r,θ,φ).\int_{\mathbb{R}^3}f\,dm=\int_Uf(r\,\omega(\theta,\varphi))\,r^2\sin\theta\;d(r,\theta,\varphi).

二つの計算の比較。 f(x)=g(x/∣x∣) 1(0,1](∣x∣)f(x)=g(x/|x|)\,\mathbf{1}_{(0,1]}(|x|)(x≠0x\ne0)とします。定理 6.3 から ∫f dm=∫01r2dr∫S2g dσ=13∫S2g dσ\int f\,dm=\int_0^1r^2dr\int_{S^2}g\,d\sigma=\frac13\int_{S^2}g\,d\sigma。上の式の右辺は、f(rω(θ,φ))=g(ω(θ,φ))1(0,1](r)f(r\omega(\theta,\varphi))=g(\omega(\theta,\varphi))\mathbf{1}_{(0,1]}(r) と Tonelli の定理から 13∫02π∫0πg(ω(θ,φ))sin⁡θ dθ dφ\frac13\int_0^{2\pi}\int_0^\pi g(\omega(\theta,\varphi))\sin\theta\,d\theta\,d\varphi です。両者を比べて主張を得ます。g=1g=1 なら ∫02π∫0πsin⁡θ dθ dφ=4π\int_0^{2\pi}\int_0^\pi\sin\theta\,d\theta\,d\varphi=4\pi で、定理 6.8 とも一致します。

注意 6.11(一般次元の角度座標)

同じ方法で、Rd\mathbb{R}^d の超球座標 x=r(cos⁡θ1, sin⁡θ1cos⁡θ2, …)x=r(\cos\theta_1,\ \sin\theta_1\cos\theta_2,\ \ldots) のヤコビアン rd−1sin⁡d−2θ1sin⁡d−3θ2⋯sin⁡θd−2r^{d-1}\sin^{d-2}\theta_1\sin^{d-3}\theta_2\cdots\sin\theta_{d-2} を計算すれば、dσ=sin⁡d−2θ1⋯sin⁡θd−2 dθ1⋯dθd−1d\sigma=\sin^{d-2}\theta_1\cdots\sin\theta_{d-2}\,d\theta_1\cdots d\theta_{d-1} が得られます。ただし実際の計算では、定理 6.3・命題 6.4・系 6.5 のように角度座標を使わずに済ませるほうが、見通しがよいことが多いです。また σ\sigma は、Sd−1S^{d-1} を Rd\mathbb{R}^d の部分多様体と見たときの誘導計量の体積測度とも一致します(d=3d=3 で確かめたのが上の命題です)。

演習 7.1易

a1,…,ad>0a_1,\ldots,a_d>0 とする。楕円体 E={x:∑ixi2/ai2≤1}\mathcal{E}=\{x:\sum_ix_i^2/a_i^2\le1\} の体積を求めよ。

解答

T=diag⁡(a1,…,ad)T=\operatorname{diag}(a_1,\ldots,a_d) とすると E=T({∣x∣≤1})\mathcal{E}=T(\{|x|\le1\}) です(y=Txy=Tx なら ∑iyi2/ai2=∣x∣2\sum_iy_i^2/a_i^2=|x|^2)。定理 3.4 と 定理 6.8 から

m(E)=a1⋯ad⋅πd/2Γ(d/2+1).m(\mathcal{E})=a_1\cdots a_d\cdot\frac{\pi^{d/2}}{\Gamma(d/2+1)} .

d=3d=3 なら 43πa1a2a3\frac43\pi a_1a_2a_3 です。

演習 7.2標準

d=2d=2 のとき、ボレル関数 g:S1→[0,∞]g:S^1\to[0,\infty] について ∫S1g dσ=∫02πg(cos⁡φ,sin⁡φ) dφ\int_{S^1}g\,d\sigma=\int_0^{2\pi}g(\cos\varphi,\sin\varphi)\,d\varphi を示せ。

解答

f(x)=g(x/∣x∣)1(0,1](∣x∣)f(x)=g(x/|x|)\mathbf{1}_{(0,1]}(|x|) とします。定理 6.3 から ∫R2f=∫01r dr∫S1g dσ=12∫S1g dσ\int_{\mathbb{R}^2}f=\int_0^1r\,dr\int_{S^1}g\,d\sigma=\frac12\int_{S^1}g\,d\sigma。一方 例 5.3 の式と Tonelli の定理から ∫R2f=∫01r dr∫02πg(cos⁡φ,sin⁡φ) dφ=12∫02πg(cos⁡φ,sin⁡φ) dφ\int_{\mathbb{R}^2}f=\int_0^1r\,dr\int_0^{2\pi}g(\cos\varphi,\sin\varphi)\,d\varphi=\frac12\int_0^{2\pi}g(\cos\varphi,\sin\varphi)\,d\varphi。両者を比べればよいです。g=1g=1 で σ(S1)=2π\sigma(S^1)=2\pi となり、定理 6.8 と一致します。

演習 7.3標準

∫Sd−1ωiωj dσ(ω)=σ(Sd−1)dδij\int_{S^{d-1}}\omega_i\omega_j\,d\sigma(\omega)=\dfrac{\sigma(S^{d-1})}{d}\delta_{ij} を示せ。また、これを使って単位球 B={∣x∣≤1}B=\{|x|\le1\} 上の ∫Bxixj dx\int_Bx_ix_j\,dx を求めよ。

解答

i≠ji\ne j のとき、第 ii 座標の符号だけを変える直交行列 OO について (Oω)i(Oω)j=−ωiωj(O\omega)_i(O\omega)_j=-\omega_i\omega_j なので、命題 6.4 から積分 IijI_{ij} は Iij=−IijI_{ij}=-I_{ij}、すなわち 00 です(∣ωiωj∣≤1|\omega_i\omega_j|\le1 で σ\sigma は有限なので可積分です)。i=ji=j のとき、座標 11 と ii を入れ替える置換行列も直交行列なので Iii=I11I_{ii}=I_{11} で、∑iIii=∫∣ω∣2dσ=σ(Sd−1)\sum_iI_{ii}=\int|\omega|^2d\sigma=\sigma(S^{d-1}) から Iii=σ(Sd−1)/dI_{ii}=\sigma(S^{d-1})/d です。

定理 6.3 により

∫Bxixj dx=∫01r2 rd−1dr⋅Iij=σ(Sd−1)d(d+2)δij=m(B)d+2δij.\int_Bx_ix_j\,dx=\int_0^1r^2\,r^{d-1}dr\cdot I_{ij}=\frac{\sigma(S^{d-1})}{d(d+2)}\delta_{ij}=\frac{m(B)}{d+2}\delta_{ij}.

d=3d=3 なら ∫Bx12 dx=4π15\int_Bx_1^2\,dx=\frac{4\pi}{15} です。

演習 7.4標準

∫Rde−∣x∣ dx\int_{\mathbb{R}^d}e^{-|x|}\,dx を求めよ。d=3d=3 ではいくつか。

解答

系 6.5 と 定理 6.8 から

∫Rde−∣x∣dx=σ(Sd−1)∫0∞e−rrd−1dr=2πd/2Γ(d/2) Γ(d)=2πd/2(d−1)!Γ(d/2).\int_{\mathbb{R}^d}e^{-|x|}dx=\sigma(S^{d-1})\int_0^\infty e^{-r}r^{d-1}dr=\frac{2\pi^{d/2}}{\Gamma(d/2)}\,\Gamma(d)=\frac{2\pi^{d/2}(d-1)!}{\Gamma(d/2)} .

rr の積分は Γ(d)\Gamma(d) の定義(定義 5.1[広義積分・パラメータ積分とガンマ関数])そのもので、Γ(d)=(d−1)!\Gamma(d)=(d-1)! は 定理 5.3[広義積分・パラメータ積分とガンマ関数] によります。d=3d=3 では Γ(3/2)=π/2\Gamma(3/2)=\sqrt\pi/2 から 2π3/2⋅2π/2=8π\frac{2\pi^{3/2}\cdot2}{\sqrt\pi/2}=8\pi です。水素原子の基底状態の波動関数 e−r/ae^{-r/a} を正規化するときに現れるのが、この種の積分です。

  • G. B. Folland, Real Analysis: Modern Techniques and Their Applications, 2nd ed., Wiley, 1999 — §2.6(Rn\mathbb{R}^n のルベーグ測度、線形写像と C1C^1 微分同相による変換)、§2.7(極座標での積分)。この記事の証明の骨格はここに沿っています。
  • W. Rudin, Real and Complex Analysis, 3rd ed., McGraw-Hill, 1987 — Chapter 2(Rk\mathbb{R}^k のルベーグ測度と線形変換)、Chapter 7(微分と変数変換定理)。
  • E. M. Stein and R. Shakarchi, Real Analysis: Measure Theory, Integration, and Hilbert Spaces, Princeton University Press, 2005 — Chapter 2(線形変換とルベーグ測度)、Chapter 6(極座標の積分公式)。
  • 伊藤清三『ルベーグ積分入門』裳華房、1963 — ユークリッド空間のルベーグ測度と変数変換を扱う章。

このテーマの 100 本ノック(100 問)

執筆: Claude Code(Opus 5.5)、Claude Code(Sonnet 5.5)・監査 3 回(Claude Code)・最終 2026-10-12

この記事の誤りを報告する ・運営: 夢現技研合同会社 ・料金 ・利用条件 ・特定商取引法に基づく表記

© 2026 夢現技研合同会社 ・コード例は MIT ライセンスです。