丸みのある屋根の下の空間の広さを知りたいとき、床を細かい升目に区切り、それぞれの升の上に細い柱を立てて、柱の体積を足し合わせるのが素朴な方法です。この記事は、その足し算を「積分」として厳密に定義し、実際の計算は「1 変数の積分を 2 回行う」だけで済む理由を説明します。さらに、座標を取り替えて円板などの積分を楽にする方法も扱います。読み終えると、三角形や円板の上の積分を自分で式に直して計算できるようになります。
2 重積分は、長方形を細かく切り、その上に立つ「柱」の体積を足し上げたリーマン和の極限として定義します。1 変数の定積分が符号付き面積だったのに対し、2 重積分は符号付き体積です。
有界閉長方形の上の連続関数は可積分です。鍵になるのは一様連続性 (関数の値の差を、2 点の距離だけで場所によらず小さく抑えられる性質)で、1 変数のときとまったく同じ筋道です。
フビニの定理は「重積分(体積)=累次積分(1 変数の積分を 2 回)」を保証します。これがないと、重積分は定義できても計算できません。
一般の領域では、領域を縦線集合として書くか横線集合として書くかが、そのまま積分順序の選択になります。順序の取り替えで値が変わるのは、重積分そのものが存在しないとき(この記事の例では被積分関数が非有界なとき)です。
変数変換公式 d x d y = ∣ det J Φ ∣ d u d v dx\,dy = \lvert\det J_\Phi\rvert\,du\,dv d x d y = ∣ det J Φ ∣ d u d v に現れるヤコビアン det J Φ \det J_\Phi det J Φ は、写像 Φ \Phi Φ が微小な面積を何倍に引き伸ばすかという「面積の拡大率」です。極座標では r r r になります。
応用として、ガウス積分 ∫ − ∞ ∞ e − x 2 d x = π \int_{-\infty}^{\infty}e^{-x^2}\,dx=\sqrt{\pi} ∫ − ∞ ∞ e − x 2 d x = π を極座標変換で計算します。
この記事で使う事前知識
上積分・下積分 : 区間を細かく切ったときの「高さの最大値で作った和」の下限と「最小値で作った和」の上限のこと。両者が一致するとき 1 変数の関数は可積分です。→ 上積分・下積分とリーマン可積分性
一様連続性 : 値の差を、2 点の距離だけで(場所によらず)小さく抑えられるという性質。→ 一様連続性
偏微分 : 他の変数を固定して、1 つの変数だけで微分したもの。→ 偏微分係数と偏導関数
行列式 : 正方行列に決まる数で、2 次では a d − b c ad-bc a d − b c のこと。→ 行列式
広義積分 : 積分区間が無限であったり、被積分関数が非有界であったりするときに、極限として定める積分。→ 広義積分
1 変数の定積分 ∫ a b f ( x ) d x \int_a^b f(x)\,dx ∫ a b f ( x ) d x は、グラフの下の符号付き面積でした。定義は「区間を細かく切り、細長い長方形の面積を足し上げ、細かくする極限を取る」というものです。詳しくは 積分の基本定理と定積分 の リーマン可積分性の定義(定義 3.1)[積分の基本定理と定積分] を参照してください。
2 変数関数 z = f ( x , y ) z = f(x,y) z = f ( x , y ) のグラフは空間内の曲面です。その下の立体の体積を測りたい、というのが 2 重積分の出発点です。定義の作り方は 1 変数と同じで、平面領域を細かい長方形に切り、その上に立つ直方体(柱)の体積 f ( 代表点 ) × ( 底面積 ) f(\text{代表点}) \times (\text{底面積}) f ( 代表点 ) × ( 底面積 ) を足し上げ、細かくする極限を取ります。
ここで素朴な疑問が出ます。「x x x で積分してから y y y で積分すればいいのではないか」。この操作、つまり 1 変数の積分を 2 回繰り返すことを累次積分 といいます。実際これが唯一の実用的な計算法です。しかし累次積分は「操作」であって、体積の定義ではありません。定義(体積)と計算法(2 回の積分)が一致する保証は、まったく自明ではないのです。
自明でないことは、次の 2 点からわかります。第一に、累次積分には d y d x dy\,dx d y d x の順と d x d y dx\,dy d x d y の順の 2 通りがあり、この 2 つの値が食い違う関数が実在します(例 5.8 )。もし累次積分がいつでも体積に等しいなら、そんなことは起こりえません。第二に、累次積分は「各断面の面積 F ( x ) = ∫ f ( x , y ) d y F(x)=\int f(x,y)\,dy F ( x ) = ∫ f ( x , y ) d y を求めてから、それを x x x 方向に積み上げる」という操作ですが、断面の面積が定まっても立体の体積が定まる保証はありません。この橋渡しをするのがフビニの定理 です。
歴史的にも、この橋渡しは繊細な問題として扱われてきました。リーマンが 1854 年の教授資格論文で積分を厳密に定義したあと、多変数への拡張と累次積分との関係が整理され、グイド・フビニが 1907 年に、ルベーグ積分の枠組みで一般的な定理を証明しました。翌々年にはトネリが、非負関数に対する使いやすい形(可積分性を仮定せずに済む形)を与えています。
もう一つの主役が変数変換 です。1 変数では置換積分 ∫ f ( g ( t ) ) g ′ ( t ) d t \int f(g(t))g'(t)\,dt ∫ f ( g ( t )) g ′ ( t ) d t が積分計算の要でした。円板や楕円板の上の積分では、直交座標のままでは領域の記述が − a 2 − x 2 ≤ y ≤ a 2 − x 2 -\sqrt{a^2-x^2} \le y \le \sqrt{a^2-x^2} − a 2 − x 2 ≤ y ≤ a 2 − x 2 のように醜くなります。極座標に取り替えれば領域は長方形になり、計算が一気に楽になります。そのとき d x d y dx\,dy d x d y が r d r d θ r\,dr\,d\theta r d r d θ に化ける理由が、ヤコビ(1841 年に関数行列式を組織的に論じました)の名を冠したヤコビアンです。
この記事では、定義(リーマン和)→ フビニの定理 → 一般の領域と積分順序 → 変数変換とヤコビアン、の順に進みます。全体の見取り図は次のとおりです。
2 重積分を計算するときの流れ
この記事では次の記号を使います。
R = [ a , b ] × [ c , d ] R = [a,b]\times[c,d] R = [ a , b ] × [ c , d ] (a < b a < b a < b , c < d c < d c < d )を R 2 \mathbb{R}^2 R 2 の閉長方形とし、その面積を ∣ R ∣ = ( b − a ) ( d − c ) \lvert R\rvert = (b-a)(d-c) ∣ R ∣ = ( b − a ) ( d − c ) と書きます。
点は P = ( x , y ) P=(x,y) P = ( x , y ) のように書き、∥ P − Q ∥ \lVert P - Q\rVert ∥ P − Q ∥ はユークリッド距離です。
f f f が集合 S S S 上有界 とは、ある K K K があって S S S 上つねに ∣ f ∣ ≤ K \lvert f\rvert \le K ∣ f ∣ ≤ K となることをいいます。
前提として次の 3 つを使います。
1 変数のリーマン積分(定義、区間加法性、単調性、線形性)。積分の基本定理と定積分 の 定積分の基本性質(命題 4.1)[積分の基本定理と定積分] を参照してください。
有界閉集合上の連続関数は一様連続であり、最大値・最小値を取ること(ハイネ・カントールの定理、最大値原理)。有界閉区間では ハイネ・カントールの定理(定理 7.1)[連続関数と一様連続性] 、閉長方形では 長方形上の連続関数の一様連続性(補題 3.1)[広義積分・パラメータ積分とガンマ関数] 、最大値・最小値は 最大値・最小値の定理(系 6.2)[コンパクト性] を参照してください。
偏微分とヤコビ行列。多変数関数の微分と偏微分 の 偏微分係数と偏導関数の定義(定義 3.1)[多変数関数の微分と偏微分] を参照してください。
まず、いちばん素直な領域である長方形の上で積分を定義します。一般の領域は第 5 節で扱います。
定義 3.1 (長方形の分割と分割の幅 ) ↗
R = [ a , b ] × [ c , d ] R=[a,b]\times[c,d] R = [ a , b ] × [ c , d ] とする。[ a , b ] [a,b] [ a , b ] の分点 a = x 0 < x 1 < ⋯ < x m = b a = x_0 < x_1 < \cdots < x_m = b a = x 0 < x 1 < ⋯ < x m = b と [ c , d ] [c,d] [ c , d ] の分点 c = y 0 < y 1 < ⋯ < y n = d c = y_0 < y_1 < \cdots < y_n = d c = y 0 < y 1 < ⋯ < y n = d の組を R R R の分割 Δ \Delta Δ と呼び、R R R を小長方形
R i j = [ x i − 1 , x i ] × [ y j − 1 , y j ] ( 1 ≤ i ≤ m , 1 ≤ j ≤ n ) R_{ij} = [x_{i-1}, x_i]\times[y_{j-1}, y_j] \qquad (1\le i\le m,\ 1\le j\le n) R ij = [ x i − 1 , x i ] × [ y j − 1 , y j ] ( 1 ≤ i ≤ m , 1 ≤ j ≤ n ) に分ける。Δ x i = x i − x i − 1 \Delta x_i = x_i - x_{i-1} Δ x i = x i − x i − 1 、Δ y j = y j − y j − 1 \Delta y_j = y_j - y_{j-1} Δ y j = y j − y j − 1 、∣ R i j ∣ = Δ x i Δ y j \lvert R_{ij}\rvert = \Delta x_i \Delta y_j ∣ R ij ∣ = Δ x i Δ y j と書く。また
∣ Δ ∣ = max i , j ( Δ x i ) 2 + ( Δ y j ) 2 \lvert \Delta\rvert = \max_{i,j}\sqrt{(\Delta x_i)^2 + (\Delta y_j)^2} ∣ Δ ∣ = i , j max ( Δ x i ) 2 + ( Δ y j ) 2 (小長方形の対角線の長さの最大値)を分割の幅 という。
幅を対角線で測っておくと、「同じ小長方形に属する 2 点の距離は ∣ Δ ∣ \lvert\Delta\rvert ∣ Δ ∣ 以下」がそのまま言えるので、あとの一様連続性の議論が楽になります。
定義 3.2 (長方形上の 2 重積分 ) ↗
f : R → R f\colon R\to\mathbb{R} f : R → R を有界関数とする。分割 Δ \Delta Δ と各小長方形から選んだ代表点 ξ i j ∈ R i j \xi_{ij}\in R_{ij} ξ ij ∈ R ij に対し
Σ ( f , Δ , ξ ) = ∑ i = 1 m ∑ j = 1 n f ( ξ i j ) Δ x i Δ y j \Sigma(f,\Delta,\xi) = \sum_{i=1}^{m}\sum_{j=1}^{n} f(\xi_{ij})\,\Delta x_i\,\Delta y_j Σ ( f , Δ , ξ ) = i = 1 ∑ m j = 1 ∑ n f ( ξ ij ) Δ x i Δ y j をリーマン和 という。ある実数 I I I が存在して、任意の ε > 0 \varepsilon > 0 ε > 0 に対しある δ > 0 \delta > 0 δ > 0 があり、∣ Δ ∣ < δ \lvert\Delta\rvert < \delta ∣ Δ ∣ < δ を満たすすべての 分割 Δ \Delta Δ とすべての 代表点の選び方 ξ \xi ξ について
∣ Σ ( f , Δ , ξ ) − I ∣ < ε \lvert \Sigma(f,\Delta,\xi) - I\rvert < \varepsilon ∣ Σ ( f , Δ , ξ ) − I ∣ < ε が成り立つとき、f f f は R R R 上**(リーマン)可積分**であるといい、I I I を f f f の R R R 上の 2 重積分 と呼んで
∬ R f ( x , y ) d x d y = I \iint_R f(x,y)\,dx\,dy = I ∬ R f ( x , y ) d x d y = I と書く。
f ≥ 0 f \ge 0 f ≥ 0 のとき、f ( ξ i j ) Δ x i Δ y j f(\xi_{ij})\Delta x_i\Delta y_j f ( ξ ij ) Δ x i Δ y j は底面 R i j R_{ij} R ij 、高さ f ( ξ i j ) f(\xi_{ij}) f ( ξ ij ) の直方体の体積です。リーマン和はこの柱をすべて足したもの、2 重積分はそれを細かくした極限で、これが「グラフの下の立体の体積」の定義になります。
定義はこれでよいのですが、このままでは「すべての分割とすべての代表点」を相手にしなければならず、可積分性の判定に使えません。1 変数のときと同じく、代表点を inf \inf inf と sup \sup sup に振り切った和を考えると扱いやすくなります。分割 Δ \Delta Δ に対し
m i j = inf R i j f , M i j = sup R i j f , s ( f , Δ ) = ∑ i , j m i j ∣ R i j ∣ , S ( f , Δ ) = ∑ i , j M i j ∣ R i j ∣ m_{ij} = \inf_{R_{ij}} f,\quad M_{ij} = \sup_{R_{ij}} f,\qquad
s(f,\Delta) = \sum_{i,j} m_{ij}\lvert R_{ij}\rvert,\quad
S(f,\Delta) = \sum_{i,j} M_{ij}\lvert R_{ij}\rvert m ij = R ij inf f , M ij = R ij sup f , s ( f , Δ ) = i , j ∑ m ij ∣ R ij ∣ , S ( f , Δ ) = i , j ∑ M ij ∣ R ij ∣
とおき、それぞれ不足和 (下ダルブー和)、過剰和 (上ダルブー和)と呼びます。inf \inf inf と sup \sup sup は、小長方形の上での f f f の最小値・最大値のことです(最小値・最大値が実際には取られない場合も含めて言えるように、下限・上限という言葉を使います)。不足和は各小長方形で「いちばん低い柱」を、過剰和は「いちばん高い柱」を立てた体積の合計ですから、本当の体積があるとすれば、必ず不足和以上、過剰和以下になります。
次の命題のねらいは、「すべての分割とすべての代表点を調べる」という重い条件を、「不足和と過剰和の差を小さくできる分割が 1 つあるか」という確かめやすい条件に言い換えることです。証明の前半は、リーマン和が不足和と過剰和に挟まれることを使い、後半は、分割を細かくすると不足和は増え過剰和は減ることを使います。
命題 3.3 (ダルブーの判定法 ) ✓
f : R → R f\colon R\to\mathbb{R} f : R → R を有界関数とする。f f f が R R R 上可積分であるための必要十分条件は、任意の ε > 0 \varepsilon > 0 ε > 0 に対して
S ( f , Δ ) − s ( f , Δ ) < ε S(f,\Delta) - s(f,\Delta) < \varepsilon S ( f , Δ ) − s ( f , Δ ) < ε を満たす R R R の分割 Δ \Delta Δ が存在することである。またこのとき
∬ R f ( x , y ) d x d y = sup Δ s ( f , Δ ) = inf Δ S ( f , Δ ) \iint_R f(x,y)\,dx\,dy = \sup_{\Delta} s(f,\Delta) = \inf_{\Delta} S(f,\Delta) ∬ R f ( x , y ) d x d y = Δ sup s ( f , Δ ) = Δ inf S ( f , Δ ) が成り立つ。
証明(命題 3.3) 準備。 任意の分割 Δ \Delta Δ と代表点 ξ \xi ξ に対し、ξ i j ∈ R i j \xi_{ij}\in R_{ij} ξ ij ∈ R ij だから m i j ≤ f ( ξ i j ) ≤ M i j m_{ij}\le f(\xi_{ij})\le M_{ij} m ij ≤ f ( ξ ij ) ≤ M ij であり、正の数 ∣ R i j ∣ \lvert R_{ij}\rvert ∣ R ij ∣ を掛けて足せば
s ( f , Δ ) ≤ Σ ( f , Δ , ξ ) ≤ S ( f , Δ ) s(f,\Delta) \le \Sigma(f,\Delta,\xi)\le S(f,\Delta) s ( f , Δ ) ≤ Σ ( f , Δ , ξ ) ≤ S ( f , Δ ) を得ます。さらに η > 0 \eta > 0 η > 0 を任意に取ると、inf \inf inf の定義から各 ( i , j ) (i,j) ( i , j ) について f ( ξ i j ) < m i j + η f(\xi_{ij}) < m_{ij}+\eta f ( ξ ij ) < m ij + η となる ξ i j ∈ R i j \xi_{ij}\in R_{ij} ξ ij ∈ R ij が選べ、そのとき Σ < s ( f , Δ ) + η ∣ R ∣ \Sigma < s(f,\Delta) + \eta\lvert R\rvert Σ < s ( f , Δ ) + η ∣ R ∣ です。η \eta η は任意なので、代表点をすべて動かしたときのリーマン和の下限はちょうど s ( f , Δ ) s(f,\Delta) s ( f , Δ ) です。同様に上限は S ( f , Δ ) S(f,\Delta) S ( f , Δ ) です。
(必要性) f f f が可積分で I = ∬ R f I = \iint_R f I = ∬ R f とします。ε > 0 \varepsilon > 0 ε > 0 に対し 定義 3.2 の δ \delta δ を ε / 3 \varepsilon/3 ε /3 に対して取り、∣ Δ ∣ < δ \lvert\Delta\rvert < \delta ∣ Δ ∣ < δ なる分割を 1 つ固定します。すべての代表点で I − ε / 3 < Σ < I + ε / 3 I - \varepsilon/3 < \Sigma < I + \varepsilon/3 I − ε /3 < Σ < I + ε /3 ですから、代表点について上限・下限を取れば、上の準備より S ( f , Δ ) ≤ I + ε / 3 S(f,\Delta)\le I+\varepsilon/3 S ( f , Δ ) ≤ I + ε /3 かつ s ( f , Δ ) ≥ I − ε / 3 s(f,\Delta)\ge I-\varepsilon/3 s ( f , Δ ) ≥ I − ε /3 。よって S − s ≤ 2 ε / 3 < ε S - s \le 2\varepsilon/3 < \varepsilon S − s ≤ 2 ε /3 < ε です。
(十分性) まず細分についての単調性を確かめます。分割 Δ ′ \Delta' Δ ′ が Δ \Delta Δ の細分 (分点を追加して得られる分割)であるとし、Δ \Delta Δ の小長方形 R i j R_{ij} R ij が Δ ′ \Delta' Δ ′ の小長方形 R 1 ′ , … , R k ′ R'_1,\dots,R'_k R 1 ′ , … , R k ′ に分かれたとします。R l ′ ⊂ R i j R'_l\subset R_{ij} R l ′ ⊂ R ij なので inf R l ′ f ≥ inf R i j f = m i j \inf_{R'_l} f \ge \inf_{R_{ij}} f = m_{ij} inf R l ′ f ≥ inf R ij f = m ij であり、∑ l ∣ R l ′ ∣ = ∣ R i j ∣ \sum_l \lvert R'_l\rvert = \lvert R_{ij}\rvert ∑ l ∣ R l ′ ∣ = ∣ R ij ∣ だから
∑ l ( inf R l ′ f ) ∣ R l ′ ∣ ≥ m i j ∣ R i j ∣ . \sum_{l}\Bigl(\inf_{R'_l} f\Bigr)\lvert R'_l\rvert \ \ge\ m_{ij}\lvert R_{ij}\rvert . l ∑ ( R l ′ inf f ) ∣ R l ′ ∣ ≥ m ij ∣ R ij ∣ . ( i , j ) (i,j) ( i , j ) について足せば s ( f , Δ ′ ) ≥ s ( f , Δ ) s(f,\Delta')\ge s(f,\Delta) s ( f , Δ ′ ) ≥ s ( f , Δ ) 、同様に S ( f , Δ ′ ) ≤ S ( f , Δ ) S(f,\Delta')\le S(f,\Delta) S ( f , Δ ′ ) ≤ S ( f , Δ ) です。ゆえに任意の 2 つの分割 Δ 1 , Δ 2 \Delta_1,\Delta_2 Δ 1 , Δ 2 に対し、両方の分点を合わせた共通細分 Δ 1 ∨ Δ 2 \Delta_1\vee\Delta_2 Δ 1 ∨ Δ 2 を挟んで
s ( f , Δ 1 ) ≤ s ( f , Δ 1 ∨ Δ 2 ) ≤ S ( f , Δ 1 ∨ Δ 2 ) ≤ S ( f , Δ 2 ) s(f,\Delta_1)\le s(f,\Delta_1\vee\Delta_2)\le S(f,\Delta_1\vee\Delta_2)\le S(f,\Delta_2) s ( f , Δ 1 ) ≤ s ( f , Δ 1 ∨ Δ 2 ) ≤ S ( f , Δ 1 ∨ Δ 2 ) ≤ S ( f , Δ 2 ) が成り立ち、I ‾ : = sup Δ s ( f , Δ ) ≤ inf Δ S ( f , Δ ) = : I ‾ \underline{I} := \sup_\Delta s(f,\Delta) \le \inf_\Delta S(f,\Delta) =: \overline{I} I := sup Δ s ( f , Δ ) ≤ inf Δ S ( f , Δ ) =: I を得ます。仮定よりどんな ε > 0 \varepsilon>0 ε > 0 に対しても I ‾ − I ‾ ≤ S ( f , Δ ) − s ( f , Δ ) < ε \overline{I}-\underline{I}\le S(f,\Delta)-s(f,\Delta)<\varepsilon I − I ≤ S ( f , Δ ) − s ( f , Δ ) < ε となる Δ \Delta Δ が取れるので、I ‾ = I ‾ \underline I = \overline I I = I 。これを J J J とおきます。
残るは、リーマン和が ∣ Δ ∣ → 0 \lvert\Delta\rvert\to 0 ∣ Δ ∣ → 0 で J J J に収束することです。K = sup R ∣ f ∣ K = \sup_R \lvert f\rvert K = sup R ∣ f ∣ とおき(f f f は有界)、ε > 0 \varepsilon>0 ε > 0 に対し S ( f , Δ 0 ) − s ( f , Δ 0 ) < ε / 2 S(f,\Delta_0)-s(f,\Delta_0)<\varepsilon/2 S ( f , Δ 0 ) − s ( f , Δ 0 ) < ε /2 となる Δ 0 \Delta_0 Δ 0 を取り、その内部の分点が x x x 方向に p p p 本、y y y 方向に q q q 本あるとします。C = p ( d − c ) + q ( b − a ) C = p(d-c)+q(b-a) C = p ( d − c ) + q ( b − a ) とおきます。任意の分割 Δ \Delta Δ に対し共通細分 Δ ∨ Δ 0 \Delta\vee\Delta_0 Δ ∨ Δ 0 を考えると、Δ \Delta Δ の小長方形のうち実際に切られるのは Δ 0 \Delta_0 Δ 0 の分点を通る直線と交わるものだけで、それらは幅 ∣ Δ ∣ \lvert\Delta\rvert ∣ Δ ∣ 以下の帯に含まれるので、面積の合計は C ∣ Δ ∣ C\lvert\Delta\rvert C ∣ Δ ∣ 以下です。切られた小長方形 1 つあたり過剰和の差は 2 K ∣ R i j ∣ 2K\lvert R_{ij}\rvert 2 K ∣ R ij ∣ 以下(sup \sup sup の値はどれも絶対値 K K K 以下)だから
S ( f , Δ ) − S ( f , Δ ∨ Δ 0 ) ≤ 2 K C ∣ Δ ∣ . S(f,\Delta) - S(f,\Delta\vee\Delta_0) \le 2KC\lvert\Delta\rvert . S ( f , Δ ) − S ( f , Δ ∨ Δ 0 ) ≤ 2 K C ∣ Δ ∣ . そこで δ = ε 8 K C + 1 \delta = \dfrac{\varepsilon}{8KC+1} δ = 8 K C + 1 ε と取れば、∣ Δ ∣ < δ \lvert\Delta\rvert<\delta ∣ Δ ∣ < δ のとき右辺は ε / 4 \varepsilon/4 ε /4 未満で、細分の単調性 S ( f , Δ ∨ Δ 0 ) ≤ S ( f , Δ 0 ) S(f,\Delta\vee\Delta_0)\le S(f,\Delta_0) S ( f , Δ ∨ Δ 0 ) ≤ S ( f , Δ 0 ) と合わせて
S ( f , Δ ) ≤ S ( f , Δ 0 ) + ε 4 < s ( f , Δ 0 ) + ε 2 + ε 4 ≤ J + 3 ε 4 S(f,\Delta) \le S(f,\Delta_0) + \frac{\varepsilon}{4} < s(f,\Delta_0)+\frac{\varepsilon}{2}+\frac{\varepsilon}{4} \le J + \frac{3\varepsilon}{4} S ( f , Δ ) ≤ S ( f , Δ 0 ) + 4 ε < s ( f , Δ 0 ) + 2 ε + 4 ε ≤ J + 4 3 ε (最後は s ( f , Δ 0 ) ≤ I ‾ = J s(f,\Delta_0)\le \underline I = J s ( f , Δ 0 ) ≤ I = J )。同じ議論で s ( f , Δ ) > J − 3 ε / 4 s(f,\Delta) > J - 3\varepsilon/4 s ( f , Δ ) > J − 3 ε /4 です。準備の不等式 s ≤ Σ ≤ S s\le\Sigma\le S s ≤ Σ ≤ S より、∣ Δ ∣ < δ \lvert\Delta\rvert<\delta ∣ Δ ∣ < δ なるすべての分割と代表点で ∣ Σ − J ∣ < ε \lvert \Sigma - J\rvert < \varepsilon ∣ Σ − J ∣ < ε 。よって f f f は可積分で積分値は J J J です。
∎
定理 3.4 (連続関数の可積分性 ) ✓
R = [ a , b ] × [ c , d ] R = [a,b]\times[c,d] R = [ a , b ] × [ c , d ] (a < b a<b a < b , c < d c<d c < d )上の連続関数 f f f は R R R 上可積分である。
証明(定理 3.4) R R R は有界閉集合なので、f f f は R R R 上一様連続です(第 2 節の前提 2)。ε > 0 \varepsilon>0 ε > 0 を与え、
∥ P − Q ∥ < δ ⟹ ∣ f ( P ) − f ( Q ) ∣ < ε 2 ∣ R ∣ \lVert P-Q\rVert < \delta \implies \lvert f(P)-f(Q)\rvert < \frac{\varepsilon}{2\lvert R\rvert} ∥ P − Q ∥ < δ ⟹ ∣ f ( P ) − f ( Q )∣ < 2 ∣ R ∣ ε となる δ > 0 \delta>0 δ > 0 を取ります。∣ Δ ∣ < δ \lvert\Delta\rvert<\delta ∣ Δ ∣ < δ なる分割を 1 つ取ると、定義 3.1 の幅の定義より、同じ小長方形 R i j R_{ij} R ij に属する 2 点の距離は ∣ Δ ∣ \lvert\Delta\rvert ∣ Δ ∣ 以下、したがって δ \delta δ 未満です。また R i j R_{ij} R ij は有界閉集合なので f f f はそこで最大値と最小値を取り、M i j = f ( Q i j ) M_{ij}=f(Q_{ij}) M ij = f ( Q ij ) 、m i j = f ( P i j ) m_{ij}=f(P_{ij}) m ij = f ( P ij ) となる点 P i j , Q i j ∈ R i j P_{ij},Q_{ij}\in R_{ij} P ij , Q ij ∈ R ij が存在します。ゆえに
M i j − m i j = f ( Q i j ) − f ( P i j ) < ε 2 ∣ R ∣ M_{ij}-m_{ij} = f(Q_{ij})-f(P_{ij}) < \frac{\varepsilon}{2\lvert R\rvert} M ij − m ij = f ( Q ij ) − f ( P ij ) < 2 ∣ R ∣ ε であり、∑ i , j ∣ R i j ∣ = ∣ R ∣ \sum_{i,j}\lvert R_{ij}\rvert = \lvert R\rvert ∑ i , j ∣ R ij ∣ = ∣ R ∣ を使って
S ( f , Δ ) − s ( f , Δ ) = ∑ i , j ( M i j − m i j ) ∣ R i j ∣ ≤ ε 2 ∣ R ∣ ⋅ ∣ R ∣ = ε 2 < ε . S(f,\Delta)-s(f,\Delta) = \sum_{i,j}(M_{ij}-m_{ij})\lvert R_{ij}\rvert \le \frac{\varepsilon}{2\lvert R\rvert}\cdot\lvert R\rvert = \frac{\varepsilon}{2} < \varepsilon . S ( f , Δ ) − s ( f , Δ ) = i , j ∑ ( M ij − m ij ) ∣ R ij ∣ ≤ 2 ∣ R ∣ ε ⋅ ∣ R ∣ = 2 ε < ε . 命題 3.3 より f f f は可積分です。
∎
例 3.5 (リーマン和から直接計算する )
f ( x , y ) = x y f(x,y)=xy f ( x , y ) = x y 、R = [ 0 , 1 ] × [ 0 , 1 ] R=[0,1]\times[0,1] R = [ 0 , 1 ] × [ 0 , 1 ] とします。f f f は連続なので 定理 3.4 より可積分です。可積分だとわかっていれば、幅が 0 0 0 に近づく都合のよい 分割・代表点の列を 1 つ選んで極限を計算すれば、それが積分値です。
x x x 方向も y y y 方向も n n n 等分し、代表点を小長方形の右上頂点 ξ i j = ( i / n , j / n ) \xi_{ij}=(i/n,\ j/n) ξ ij = ( i / n , j / n ) に取ります。Δ x i = Δ y j = 1 / n \Delta x_i = \Delta y_j = 1/n Δ x i = Δ y j = 1/ n なので
Σ = ∑ i = 1 n ∑ j = 1 n i n ⋅ j n ⋅ 1 n 2 = 1 n 4 ( ∑ i = 1 n i ) ( ∑ j = 1 n j ) = 1 n 4 ⋅ n ( n + 1 ) 2 ⋅ n ( n + 1 ) 2 = ( n + 1 ) 2 4 n 2 . \Sigma = \sum_{i=1}^{n}\sum_{j=1}^{n}\frac{i}{n}\cdot\frac{j}{n}\cdot\frac{1}{n^2}
= \frac{1}{n^4}\Bigl(\sum_{i=1}^{n} i\Bigr)\Bigl(\sum_{j=1}^{n} j\Bigr)
= \frac{1}{n^4}\cdot\frac{n(n+1)}{2}\cdot\frac{n(n+1)}{2}
= \frac{(n+1)^2}{4n^2}. Σ = i = 1 ∑ n j = 1 ∑ n n i ⋅ n j ⋅ n 2 1 = n 4 1 ( i = 1 ∑ n i ) ( j = 1 ∑ n j ) = n 4 1 ⋅ 2 n ( n + 1 ) ⋅ 2 n ( n + 1 ) = 4 n 2 ( n + 1 ) 2 . 分割の幅は 2 / n → 0 \sqrt{2}/n\to 0 2 / n → 0 で、( n + 1 ) 2 4 n 2 → 1 4 \dfrac{(n+1)^2}{4n^2}\to\dfrac14 4 n 2 ( n + 1 ) 2 → 4 1 です。したがって
∬ [ 0 , 1 ] 2 x y d x d y = 1 4 . \iint_{[0,1]^2} xy\,dx\,dy = \frac14 . ∬ [ 0 , 1 ] 2 x y d x d y = 4 1 . 可積分性を先に確かめた点が大事です。可積分性を知らずに 1 つの列だけを計算しても、他の分割で別の値になる可能性を排除できません。
命題 3.6 (重積分の基本性質 ) ✓
f , g f,g f , g が R R R 上可積分、α , β ∈ R \alpha,\beta\in\mathbb{R} α , β ∈ R とする。
(1)(線形性) α f + β g \alpha f+\beta g α f + β g は可積分で ∬ R ( α f + β g ) = α ∬ R f + β ∬ R g \iint_R(\alpha f+\beta g)=\alpha\iint_R f+\beta\iint_R g ∬ R ( α f + β g ) = α ∬ R f + β ∬ R g 。
(2)(単調性) R R R 上つねに f ≤ g f\le g f ≤ g なら ∬ R f ≤ ∬ R g \iint_R f\le\iint_R g ∬ R f ≤ ∬ R g 。とくに f ≥ 0 f\ge0 f ≥ 0 なら ∬ R f ≥ 0 \iint_R f\ge0 ∬ R f ≥ 0 。
(3)(領域についての加法性) Δ 0 \Delta_0 Δ 0 を R R R の分割とし、その小長方形を R 1 , … , R N R_1,\dots,R_N R 1 , … , R N とする。有界関数 f f f が R R R 上可積分であることと、各 R k R_k R k 上可積分であることは同値で、そのとき ∬ R f = ∑ k ∬ R k f \iint_R f=\sum_k\iint_{R_k}f ∬ R f = ∑ k ∬ R k f である。
証明(命題 3.6) (1) 同じ分割 Δ \Delta Δ と代表点 ξ \xi ξ に対し Σ ( α f + β g , Δ , ξ ) = α Σ ( f , Δ , ξ ) + β Σ ( g , Δ , ξ ) \Sigma(\alpha f+\beta g,\Delta,\xi)=\alpha\Sigma(f,\Delta,\xi)+\beta\Sigma(g,\Delta,\xi) Σ ( α f + β g , Δ , ξ ) = α Σ ( f , Δ , ξ ) + β Σ ( g , Δ , ξ ) です。I f , I g I_f,I_g I f , I g を f , g f,g f , g の積分とし、ε > 0 \varepsilon>0 ε > 0 に対し 定義 3.2 の δ \delta δ を f f f と g g g のそれぞれに ε / ( ∣ α ∣ + ∣ β ∣ + 1 ) \varepsilon/(|\alpha|+|\beta|+1) ε / ( ∣ α ∣ + ∣ β ∣ + 1 ) で取って小さいほうを採れば、∣ Δ ∣ < δ \lvert\Delta\rvert<\delta ∣ Δ ∣ < δ なるすべての Δ , ξ \Delta,\xi Δ , ξ で ∣ Σ ( α f + β g ) − α I f − β I g ∣ < ε \lvert\Sigma(\alpha f+\beta g)-\alpha I_f-\beta I_g\rvert<\varepsilon ∣ Σ ( α f + β g ) − α I f − β I g ∣ < ε です。
(2) f ≤ g f\le g f ≤ g なら Σ ( f , Δ , ξ ) ≤ Σ ( g , Δ , ξ ) \Sigma(f,\Delta,\xi)\le\Sigma(g,\Delta,\xi) Σ ( f , Δ , ξ ) ≤ Σ ( g , Δ , ξ ) です。ε > 0 \varepsilon>0 ε > 0 に対し ∣ Δ ∣ \lvert\Delta\rvert ∣ Δ ∣ が十分小さい分割を取ると I f − ε < Σ ( f ) ≤ Σ ( g ) < I g + ε I_f-\varepsilon<\Sigma(f)\le\Sigma(g)<I_g+\varepsilon I f − ε < Σ ( f ) ≤ Σ ( g ) < I g + ε なので I f < I g + 2 ε I_f<I_g+2\varepsilon I f < I g + 2 ε で、ε \varepsilon ε は任意です。
(3) 小長方形の inf \inf inf は細分で増え sup \sup sup は減るので、細分すると S S S は増えず s s s は減りません。したがって S − s S-s S − s も増えません。Δ 0 \Delta_0 Δ 0 の分点をすべて含む R R R の分割 Δ \Delta Δ は、各 R k R_k R k の分割 Δ k \Delta_k Δ k を定め、S ( f , Δ ) = ∑ k S ( f , Δ k ) S(f,\Delta)=\sum_kS(f,\Delta_k) S ( f , Δ ) = ∑ k S ( f , Δ k ) 、s ( f , Δ ) = ∑ k s ( f , Δ k ) s(f,\Delta)=\sum_ks(f,\Delta_k) s ( f , Δ ) = ∑ k s ( f , Δ k ) なので
S ( f , Δ ) − s ( f , Δ ) = ∑ k ( S ( f , Δ k ) − s ( f , Δ k ) ) , 各項 ≥ 0 S(f,\Delta)-s(f,\Delta)=\sum_k\bigl(S(f,\Delta_k)-s(f,\Delta_k)\bigr),\qquad\text{各項}\ge0 S ( f , Δ ) − s ( f , Δ ) = k ∑ ( S ( f , Δ k ) − s ( f , Δ k ) ) , 各項 ≥ 0 です。R R R 上可積分とし ε > 0 \varepsilon>0 ε > 0 とします。命題 3.3 で S − s < ε S-s<\varepsilon S − s < ε の分割を取り、Δ 0 \Delta_0 Δ 0 の分点を加えて細分すると、各 k k k で S ( f , Δ k ) − s ( f , Δ k ) < ε S(f,\Delta_k)-s(f,\Delta_k)<\varepsilon S ( f , Δ k ) − s ( f , Δ k ) < ε なので、命題 3.3 より f f f は R k R_k R k 上可積分です。逆に各 R k R_k R k 上可積分とし、各 k k k で S ( f , Δ k ′ ) − s ( f , Δ k ′ ) < ε / N S(f,\Delta_k')-s(f,\Delta_k')<\varepsilon/N S ( f , Δ k ′ ) − s ( f , Δ k ′ ) < ε / N となる R k R_k R k の分割 Δ k ′ \Delta_k' Δ k ′ を取ります。Δ 0 \Delta_0 Δ 0 とすべての Δ k ′ \Delta_k' Δ k ′ の分点を合わせた R R R の分割 Δ \Delta Δ が R k R_k R k に定める分割 Δ k \Delta_k Δ k は Δ k ′ \Delta_k' Δ k ′ の細分なので S ( f , Δ k ) − s ( f , Δ k ) < ε / N S(f,\Delta_k)-s(f,\Delta_k)<\varepsilon/N S ( f , Δ k ) − s ( f , Δ k ) < ε / N で、和を取って S ( f , Δ ) − s ( f , Δ ) < ε S(f,\Delta)-s(f,\Delta)<\varepsilon S ( f , Δ ) − s ( f , Δ ) < ε 、よって R R R 上可積分です。
積分値。上のような Δ \Delta Δ では s ( f , Δ k ) ≤ ∬ R k f ≤ S ( f , Δ k ) s(f,\Delta_k)\le\iint_{R_k}f\le S(f,\Delta_k) s ( f , Δ k ) ≤ ∬ R k f ≤ S ( f , Δ k ) (命題 3.3 の sup s = inf S = \sup s=\inf S= sup s = inf S = 積分)で、足すと s ( f , Δ ) ≤ ∑ k ∬ R k f ≤ S ( f , Δ ) s(f,\Delta)\le\sum_k\iint_{R_k}f\le S(f,\Delta) s ( f , Δ ) ≤ ∑ k ∬ R k f ≤ S ( f , Δ ) 、また s ( f , Δ ) ≤ ∬ R f ≤ S ( f , Δ ) s(f,\Delta)\le\iint_Rf\le S(f,\Delta) s ( f , Δ ) ≤ ∬ R f ≤ S ( f , Δ ) です。S − s < ε S-s<\varepsilon S − s < ε なので両者の差は ε \varepsilon ε 未満で、ε \varepsilon ε は任意だから ∬ R f = ∑ k ∬ R k f \iint_Rf=\sum_k\iint_{R_k}f ∬ R f = ∑ k ∬ R k f です。
∎
定義 4.1 (累次積分 ) ↗
f : R = [ a , b ] × [ c , d ] → R f\colon R=[a,b]\times[c,d]\to\mathbb{R} f : R = [ a , b ] × [ c , d ] → R とする。各 x ∈ [ a , b ] x\in[a,b] x ∈ [ a , b ] を固定したとき 1 変数関数 y ↦ f ( x , y ) y\mapsto f(x,y) y ↦ f ( x , y ) が [ c , d ] [c,d] [ c , d ] 上可積分で、そうして定まる関数
F ( x ) = ∫ c d f ( x , y ) d y F(x) = \int_c^d f(x,y)\,dy F ( x ) = ∫ c d f ( x , y ) d y が [ a , b ] [a,b] [ a , b ] 上可積分であるとき、
∫ a b ( ∫ c d f ( x , y ) d y ) d x = ∫ a b F ( x ) d x \int_a^b\!\left(\int_c^d f(x,y)\,dy\right)dx = \int_a^b F(x)\,dx ∫ a b ( ∫ c d f ( x , y ) d y ) d x = ∫ a b F ( x ) d x を f f f の累次積分 (反復積分)という。x x x と y y y の役割を入れ替えた ∫ c d ( ∫ a b f ( x , y ) d x ) d y \displaystyle\int_c^d\!\left(\int_a^b f(x,y)\,dx\right)dy ∫ c d ( ∫ a b f ( x , y ) d x ) d y も同様に定める。
2 つの積分の性格の違いを並べておきます。
重積分 ∬ R f d x d y \displaystyle\iint_R f\,dx\,dy ∬ R f d x d y 累次積分 ∫ a b ( ∫ c d f d y ) d x \displaystyle\int_a^b\!\left(\int_c^d f\,dy\right)dx ∫ a b ( ∫ c d f d y ) d x 定義 平面の分割によるリーマン和の極限 1 変数の定積分を 2 回行う操作 幾何的意味 グラフの下の符号付き体積 断面積 F ( x ) F(x) F ( x ) を x x x 方向に積み上げたもの 存在条件 有界性 + 命題 3.3 各断面の可積分性 + F F F の可積分性 計算 ふつう直接は計算できない 1 変数の技法がそのまま使える
この 2 つを結ぶのがフビニの定理です。仮定 (ii) は「切り口ごとの面積がそもそも定義できる」ことの要請で、これを落とすと右辺の意味が失われます。
定理 4.2 (フビニの定理(長方形上) ) ✓
R = [ a , b ] × [ c , d ] R=[a,b]\times[c,d] R = [ a , b ] × [ c , d ] とし、f : R → R f\colon R\to\mathbb{R} f : R → R が次の 2 条件を満たすとする。
f f f は R R R 上有界かつ可積分である。
各 x ∈ [ a , b ] x\in[a,b] x ∈ [ a , b ] に対し、1 変数関数 y ↦ f ( x , y ) y\mapsto f(x,y) y ↦ f ( x , y ) は [ c , d ] [c,d] [ c , d ] 上可積分である。
このとき F ( x ) = ∫ c d f ( x , y ) d y F(x)=\displaystyle\int_c^d f(x,y)\,dy F ( x ) = ∫ c d f ( x , y ) d y は [ a , b ] [a,b] [ a , b ] 上可積分で、
∫ a b ( ∫ c d f ( x , y ) d y ) d x = ∬ R f ( x , y ) d x d y \int_a^b\!\left(\int_c^d f(x,y)\,dy\right)dx = \iint_R f(x,y)\,dx\,dy ∫ a b ( ∫ c d f ( x , y ) d y ) d x = ∬ R f ( x , y ) d x d y が成り立つ。x x x と y y y の役割を入れ替えた主張(各 y y y について x ↦ f ( x , y ) x\mapsto f(x,y) x ↦ f ( x , y ) が可積分なら ∫ c d ( ∫ a b f d x ) d y = ∬ R f \int_c^d\!\left(\int_a^b f\,dx\right)dy = \iint_R f ∫ c d ( ∫ a b f d x ) d y = ∬ R f )も同様に成り立つ。
証明の方針は次のとおりです。各 x x x での断面積 F ( x ) F(x) F ( x ) は、小長方形ごとの「いちばん低い柱」と「いちばん高い柱」の高さの和で上下から挟めます。この挟み込みを x x x 方向にも行うと、F F F の積分が不足和と過剰和に挟まれます。ただし F F F が可積分かどうかはまだ分かっていないので、1 変数の上積分・下積分 (可積分でなくても定まる、積分の下からの値と上からの値)を使い、それらが一致することを示して可積分性も同時に得ます。
証明(定理 4.2) R R R の分割 Δ \Delta Δ を任意に取り、その x x x 方向の分点だけを取り出した [ a , b ] [a,b] [ a , b ] の分割を Δ x \Delta_x Δ x と書きます。
x ∈ [ x i − 1 , x i ] x\in[x_{i-1},x_i] x ∈ [ x i − 1 , x i ] を固定します。仮定 2 と 1 変数の積分の区間加法性から
F ( x ) = ∫ c d f ( x , y ) d y = ∑ j = 1 n ∫ y j − 1 y j f ( x , y ) d y . F(x) = \int_c^d f(x,y)\,dy = \sum_{j=1}^{n}\int_{y_{j-1}}^{y_j} f(x,y)\,dy . F ( x ) = ∫ c d f ( x , y ) d y = j = 1 ∑ n ∫ y j − 1 y j f ( x , y ) d y . y ∈ [ y j − 1 , y j ] y\in[y_{j-1},y_j] y ∈ [ y j − 1 , y j ] のとき ( x , y ) ∈ R i j (x,y)\in R_{ij} ( x , y ) ∈ R ij なので m i j ≤ f ( x , y ) ≤ M i j m_{ij}\le f(x,y)\le M_{ij} m ij ≤ f ( x , y ) ≤ M ij であり、1 変数の積分の単調性から
m i j Δ y j ≤ ∫ y j − 1 y j f ( x , y ) d y ≤ M i j Δ y j . m_{ij}\,\Delta y_j \le \int_{y_{j-1}}^{y_j} f(x,y)\,dy \le M_{ij}\,\Delta y_j . m ij Δ y j ≤ ∫ y j − 1 y j f ( x , y ) d y ≤ M ij Δ y j . j j j について足すと、[ x i − 1 , x i ] [x_{i-1},x_i] [ x i − 1 , x i ] のすべての x x x について
∑ j = 1 n m i j Δ y j ≤ F ( x ) ≤ ∑ j = 1 n M i j Δ y j \sum_{j=1}^{n} m_{ij}\Delta y_j \ \le\ F(x)\ \le\ \sum_{j=1}^{n} M_{ij}\Delta y_j j = 1 ∑ n m ij Δ y j ≤ F ( x ) ≤ j = 1 ∑ n M ij Δ y j が成り立ちます。左辺・右辺は x x x によらない定数なので、[ x i − 1 , x i ] [x_{i-1},x_i] [ x i − 1 , x i ] 上での F F F の下限・上限について
inf [ x i − 1 , x i ] F ≥ ∑ j m i j Δ y j , sup [ x i − 1 , x i ] F ≤ ∑ j M i j Δ y j \inf_{[x_{i-1},x_i]} F \ \ge\ \sum_{j} m_{ij}\Delta y_j,
\qquad
\sup_{[x_{i-1},x_i]} F \ \le\ \sum_{j} M_{ij}\Delta y_j [ x i − 1 , x i ] inf F ≥ j ∑ m ij Δ y j , [ x i − 1 , x i ] sup F ≤ j ∑ M ij Δ y j が従います。Δ x i > 0 \Delta x_i>0 Δ x i > 0 を掛けて i i i について足せば
s ( f , Δ ) = ∑ i , j m i j Δ y j Δ x i ≤ s ( F , Δ x ) ≤ S ( F , Δ x ) ≤ ∑ i , j M i j Δ y j Δ x i = S ( f , Δ ) s(f,\Delta) = \sum_{i,j} m_{ij}\Delta y_j\Delta x_i \ \le\ s(F,\Delta_x)
\ \le\ S(F,\Delta_x)\ \le\ \sum_{i,j} M_{ij}\Delta y_j \Delta x_i = S(f,\Delta) s ( f , Δ ) = i , j ∑ m ij Δ y j Δ x i ≤ s ( F , Δ x ) ≤ S ( F , Δ x ) ≤ i , j ∑ M ij Δ y j Δ x i = S ( f , Δ ) を得ます(真ん中の不等号は不足和 ≤ \le ≤ 過剰和)。ここで 1 変数の下積分・上積分は s ( F , Δ x ) ≤ ∫ a b ‾ F ≤ ∫ a b ‾ F ≤ S ( F , Δ x ) s(F,\Delta_x)\le \underline{\int_a^b} F \le \overline{\int_a^b} F\le S(F,\Delta_x) s ( F , Δ x ) ≤ ∫ a b F ≤ ∫ a b F ≤ S ( F , Δ x ) を満たすので、結局すべての分割 Δ \Delta Δ に対して
s ( f , Δ ) ≤ ∫ a b ‾ F ≤ ∫ a b ‾ F ≤ S ( f , Δ ) . s(f,\Delta)\ \le\ \underline{\int_a^b} F \ \le\ \overline{\int_a^b} F\ \le\ S(f,\Delta) . s ( f , Δ ) ≤ ∫ a b F ≤ ∫ a b F ≤ S ( f , Δ ) . 仮定 1 と 命題 3.3 より sup Δ s ( f , Δ ) = inf Δ S ( f , Δ ) = ∬ R f \sup_\Delta s(f,\Delta)=\inf_\Delta S(f,\Delta)=\iint_R f sup Δ s ( f , Δ ) = inf Δ S ( f , Δ ) = ∬ R f ですから、Δ \Delta Δ について上限・下限を取ると挟み撃ちで
∫ a b ‾ F = ∫ a b ‾ F = ∬ R f ( x , y ) d x d y . \underline{\int_a^b} F = \overline{\int_a^b} F = \iint_R f(x,y)\,dx\,dy . ∫ a b F = ∫ a b F = ∬ R f ( x , y ) d x d y . 上積分と下積分が一致するので F F F は [ a , b ] [a,b] [ a , b ] 上可積分であり(1 変数版の リーマンの可積分条件(定理 3.4)[積分の基本定理と定積分] )、その積分値は ∬ R f \iint_R f ∬ R f です。x x x と y y y を入れ替えても議論はまったく同じです。
∎
系 4.3 (連続関数に対するフビニの定理 ) ✓
f f f が R = [ a , b ] × [ c , d ] R=[a,b]\times[c,d] R = [ a , b ] × [ c , d ] 上連続ならば、2 つの累次積分はともに存在して等しく、重積分に一致する。すなわち
∫ a b ( ∫ c d f ( x , y ) d y ) d x = ∫ c d ( ∫ a b f ( x , y ) d x ) d y = ∬ R f ( x , y ) d x d y . \int_a^b\!\left(\int_c^d f(x,y)\,dy\right)dx
= \int_c^d\!\left(\int_a^b f(x,y)\,dx\right)dy
= \iint_R f(x,y)\,dx\,dy . ∫ a b ( ∫ c d f ( x , y ) d y ) d x = ∫ c d ( ∫ a b f ( x , y ) d x ) d y = ∬ R f ( x , y ) d x d y .
証明(系 4.3) f f f は有界閉集合 R R R 上の連続関数なので有界で、定理 3.4 より可積分です。これで 定理 4.2 の仮定 1 が満たされます。次に x ∈ [ a , b ] x\in[a,b] x ∈ [ a , b ] を固定すると、y ↦ f ( x , y ) y\mapsto f(x,y) y ↦ f ( x , y ) は [ c , d ] [c,d] [ c , d ] 上の連続関数なので 1 変数の意味で可積分です。これで仮定 2 も満たされ、定理 4.2 から第 1 の式と第 3 の式が等しいことがわかります。同じ議論を x x x と y y y を入れ替えて行えば、第 2 の式も第 3 の式に等しくなります。
∎
例 4.4 (順序の選び方で手間が変わる )
∬ [ 0 , 1 ] × [ 0 , π ] x cos ( x y ) d x d y \iint_{[0,1]\times[0,\pi]} x\cos(xy)\,dx\,dy ∬ [ 0 , 1 ] × [ 0 , π ] x cos ( x y ) d x d y を計算します。被積分関数は連続なので 系 4.3 よりどちらの順でもよく、値は一致します。
先に y y y で積分する場合。 x x x を固定すると、y y y の関数として cos ( x y ) \cos(xy) cos ( x y ) の原始関数は sin ( x y ) / x \sin(xy)/x sin ( x y ) / x (x ≠ 0 x\neq0 x = 0 )ですから
∫ 0 π x cos ( x y ) d y = x ⋅ [ sin ( x y ) x ] y = 0 y = π = sin ( π x ) . \int_0^{\pi} x\cos(xy)\,dy = x\cdot\left[\frac{\sin(xy)}{x}\right]_{y=0}^{y=\pi} = \sin(\pi x). ∫ 0 π x cos ( x y ) d y = x ⋅ [ x sin ( x y ) ] y = 0 y = π = sin ( π x ) . x = 0 x=0 x = 0 のときは被積分関数が 0 0 0 なので値は 0 0 0 で、これも sin ( π ⋅ 0 ) = 0 \sin(\pi\cdot 0)=0 sin ( π ⋅ 0 ) = 0 と一致します。よって
∫ 0 1 sin ( π x ) d x = [ − cos ( π x ) π ] 0 1 = 1 + 1 π = 2 π . \int_0^1 \sin(\pi x)\,dx = \left[-\frac{\cos(\pi x)}{\pi}\right]_0^1 = \frac{1+1}{\pi} = \frac{2}{\pi}. ∫ 0 1 sin ( π x ) d x = [ − π cos ( π x ) ] 0 1 = π 1 + 1 = π 2 . 先に x x x で積分する場合。 y > 0 y>0 y > 0 を固定して部分積分(u = x u=x u = x , d v = cos ( x y ) d x dv=\cos(xy)\,dx d v = cos ( x y ) d x , v = sin ( x y ) / y v=\sin(xy)/y v = sin ( x y ) / y )を行うと
∫ 0 1 x cos ( x y ) d x = [ x sin ( x y ) y ] 0 1 − ∫ 0 1 sin ( x y ) y d x = sin y y − 1 y [ − cos ( x y ) y ] 0 1 = sin y y + cos y − 1 y 2 . \int_0^1 x\cos(xy)\,dx = \left[\frac{x\sin(xy)}{y}\right]_0^1 - \int_0^1\frac{\sin(xy)}{y}\,dx
= \frac{\sin y}{y} - \frac{1}{y}\left[-\frac{\cos(xy)}{y}\right]_0^1
= \frac{\sin y}{y}+\frac{\cos y - 1}{y^2}. ∫ 0 1 x cos ( x y ) d x = [ y x sin ( x y ) ] 0 1 − ∫ 0 1 y sin ( x y ) d x = y sin y − y 1 [ − y cos ( x y ) ] 0 1 = y sin y + y 2 cos y − 1 . これを G ( y ) G(y) G ( y ) と書きます。y → 0 y\to 0 y → 0 では sin y / y → 1 \sin y/y\to1 sin y / y → 1 、( cos y − 1 ) / y 2 → − 1 / 2 (\cos y-1)/y^2\to-1/2 ( cos y − 1 ) / y 2 → − 1/2 なので G ( 0 + ) = 1 / 2 G(0^+)=1/2 G ( 0 + ) = 1/2 となり、y = 0 y=0 y = 0 での直接計算 ∫ 0 1 x d x = 1 / 2 \int_0^1 x\,dx = 1/2 ∫ 0 1 x d x = 1/2 と一致します。次に
d d y ( 1 − cos y y ) = y sin y − ( 1 − cos y ) y 2 = sin y y − 1 − cos y y 2 = G ( y ) \frac{d}{dy}\left(\frac{1-\cos y}{y}\right)
= \frac{y\sin y - (1-\cos y)}{y^2}
= \frac{\sin y}{y} - \frac{1-\cos y}{y^2} = G(y) d y d ( y 1 − cos y ) = y 2 y sin y − ( 1 − cos y ) = y sin y − y 2 1 − cos y = G ( y ) なので、( 1 − cos y ) / y (1-\cos y)/y ( 1 − cos y ) / y が G G G の原始関数です(y → 0 y\to0 y → 0 での値は 0 0 0 )。したがって
∫ 0 π G ( y ) d y = [ 1 − cos y y ] 0 π = 2 π − 0 = 2 π . \int_0^{\pi} G(y)\,dy = \left[\frac{1-\cos y}{y}\right]_{0}^{\pi} = \frac{2}{\pi} - 0 = \frac{2}{\pi}. ∫ 0 π G ( y ) d y = [ y 1 − cos y ] 0 π = π 2 − 0 = π 2 . 両者は確かに一致しました。ただし手間はまるで違います。順序を選べること自体がフビニの定理の実用的な価値です。
現実に出てくる領域は長方形ではなく、三角形・円板・2 曲線に挟まれた領域などです。これらを扱う標準的な手口は「長方形まで 0 0 0 で延長する」ことです。
定義 5.1 (一般の有界領域上の重積分と面積確定 ) ↗
D ⊂ R 2 D\subset\mathbb{R}^2 D ⊂ R 2 を有界集合、f : D → R f\colon D\to\mathbb{R} f : D → R を有界関数とする。D ⊂ R D\subset R D ⊂ R となる閉長方形 R R R を 1 つ取り、
f ~ ( x , y ) = { f ( x , y ) ( ( x , y ) ∈ D ) 0 ( ( x , y ) ∈ R ∖ D ) \tilde f(x,y) = \begin{cases} f(x,y) & ((x,y)\in D)\\ 0 & ((x,y)\in R\setminus D)\end{cases} f ~ ( x , y ) = { f ( x , y ) 0 (( x , y ) ∈ D ) (( x , y ) ∈ R ∖ D ) とおく。f ~ \tilde f f ~ が R R R 上可積分であるとき、f f f は D D D 上可積分であるといい
∬ D f ( x , y ) d x d y : = ∬ R f ~ ( x , y ) d x d y \iint_D f(x,y)\,dx\,dy := \iint_R \tilde f(x,y)\,dx\,dy ∬ D f ( x , y ) d x d y := ∬ R f ~ ( x , y ) d x d y と定める。とくに f ≡ 1 f\equiv 1 f ≡ 1 が D D D 上可積分であるとき D D D は面積確定 であるといい、∣ D ∣ = ∬ D 1 d x d y \lvert D\rvert = \iint_D 1\,dx\,dy ∣ D ∣ = ∬ D 1 d x d y を D D D の面積 という。
この定義が R R R の取り方によらないことは次のようにわかります。R ⊂ R ′ R\subset R' R ⊂ R ′ なる長方形を取り、R ′ R' R ′ の分割で R R R の 4 辺を分点に含むものを考えると、R ′ ∖ R R'\setminus R R ′ ∖ R の側の小長方形では f ~ ≡ 0 \tilde f\equiv 0 f ~ ≡ 0 なのでリーマン和に寄与せず、R R R 側の和とちょうど一致します。一般の 2 つの長方形は、両方を含む長方形を経由して比べればよいです。
ここから 定理 5.6 までの補題は、すべて 1 つの目的のためにあります。D D D の外を 0 0 0 で埋めた f ~ \tilde f f ~ は、D D D の境界の上で値が飛ぶので、そのままでは可積分かどうか分かりません。そこで、境界が「面積 0 0 0 」の集合であれば無視できることを示します。流れは、(1) 連続関数のグラフは面積 0 0 0 である、(2) 面積 0 0 0 とは「総面積がいくらでも小さい長方形で覆える」ことである、(3) 不連続点がそのような集合に入っていれば可積分である、の 3 段階です。なお (2) と (3) の証明に出てくる 1 Z \mathbf 1_Z 1 Z は、Z Z Z の上で 1 1 1 、外で 0 0 0 とする関数(指示関数)のことです。また、D D D を縦に切るように表せる領域を縦線集合 、横に切るように表せる領域を横線集合 といい、次の定理で正確に定義します。
補題 5.2 (連続関数のグラフは面積 0 ) ✓
φ : [ a , b ] → R \varphi\colon[a,b]\to\mathbb{R} φ : [ a , b ] → R を連続関数とすると、そのグラフ G = { ( t , φ ( t ) ) ∣ t ∈ [ a , b ] } G=\{(t,\varphi(t)) \mid t\in[a,b]\} G = {( t , φ ( t )) ∣ t ∈ [ a , b ]} は面積確定で ∣ G ∣ = 0 \lvert G\rvert = 0 ∣ G ∣ = 0 である。同じことは横向きのグラフ { ( ψ ( t ) , t ) } \{(\psi(t),t)\} {( ψ ( t ) , t )} についても成り立つ。
証明(補題 5.2) ε > 0 \varepsilon>0 ε > 0 を取ります。φ \varphi φ は有界閉区間上の連続関数なので一様連続で、∣ t − t ′ ∣ < δ ⇒ ∣ φ ( t ) − φ ( t ′ ) ∣ < ε \lvert t-t'\rvert<\delta \Rightarrow \lvert\varphi(t)-\varphi(t')\rvert<\varepsilon ∣ t − t ′ ∣ < δ ⇒ ∣ φ ( t ) − φ ( t ′ )∣ < ε となる δ > 0 \delta>0 δ > 0 があります。[ a , b ] [a,b] [ a , b ] を n n n 等分して ( b − a ) / n < δ (b-a)/n<\delta ( b − a ) / n < δ とし、小区間を I 1 , … , I n I_1,\dots,I_n I 1 , … , I n とします。φ \varphi φ は各 I k I_k I k 上で最大値 β k \beta_k β k と最小値 α k \alpha_k α k を取り、I k I_k I k 内の 2 点での値の差は ε \varepsilon ε 未満なので β k − α k ≤ ε \beta_k-\alpha_k\le\varepsilon β k − α k ≤ ε です。
G G G を含む閉長方形 R = [ a , b ] × [ c , d ] R=[a,b]\times[c,d] R = [ a , b ] × [ c , d ] を取り、G G G 上で 1 1 1 、R ∖ G R\setminus G R ∖ G で 0 0 0 となる関数を g g g とします。R R R の分割 Δ \Delta Δ として、x x x 方向は上の n n n 等分、y y y 方向は c , d c,d c , d と 2 n 2n 2 n 個の値 α 1 , β 1 , … , α n , β n \alpha_1,\beta_1,\dots,\alpha_n,\beta_n α 1 , β 1 , … , α n , β n を分点に持つものを取ります。I k I_k I k の上にある G G G の部分は帯 I k × [ α k , β k ] I_k\times[\alpha_k,\beta_k] I k × [ α k , β k ] に含まれ、α k , β k \alpha_k,\beta_k α k , β k は分点なので、この帯は第 k k k 列の小長方形いくつかの合併になっています。それ以外の小長方形では g ≡ 0 g\equiv0 g ≡ 0 、すなわち sup g = 0 \sup g=0 sup g = 0 ですから
S ( g , Δ ) ≤ ∑ k = 1 n b − a n ( β k − α k ) ≤ ∑ k = 1 n b − a n ε = ( b − a ) ε . S(g,\Delta) \le \sum_{k=1}^{n} \frac{b-a}{n}\,(\beta_k-\alpha_k) \le \sum_{k=1}^{n}\frac{b-a}{n}\,\varepsilon = (b-a)\varepsilon . S ( g , Δ ) ≤ k = 1 ∑ n n b − a ( β k − α k ) ≤ k = 1 ∑ n n b − a ε = ( b − a ) ε . 一方、どの小長方形も幅と高さが正であり、各 x x x 座標に対して G G G の点は高々 1 個しかないので、小長方形が G G G に含まれることはなく、必ず G G G の外の点を含みます。そこでは g = 0 g=0 g = 0 なので inf g = 0 \inf g=0 inf g = 0 となり、s ( g , Δ ) = 0 s(g,\Delta)=0 s ( g , Δ ) = 0 です。したがって S − s ≤ ( b − a ) ε S-s\le(b-a)\varepsilon S − s ≤ ( b − a ) ε が任意の ε > 0 \varepsilon>0 ε > 0 で成り立ち、命題 3.3 より g g g は可積分で、積分値は inf Δ S = 0 \inf_\Delta S=0 inf Δ S = 0 です。横向きのグラフも x x x と y y y を入れ替えれば同じです。
∎
補題 5.3 (面積 0 の集合の被覆による特徴づけ ) ✓
有界集合 Z ⊂ R 2 Z\subset\mathbb{R}^2 Z ⊂ R 2 について、次の 2 条件は同値である。
Z Z Z は面積確定で ∣ Z ∣ = 0 \lvert Z\rvert=0 ∣ Z ∣ = 0 (定義 5.1 の意味)。
任意の ε > 0 \varepsilon>0 ε > 0 に対し、総面積が ε \varepsilon ε 未満の有限個の閉長方形で Z Z Z を覆える。
したがって、2 を満たす集合の部分集合、およびそのような集合の有限個の合併は 2 を満たし、座標軸に平行な線分も 2 を満たす。
証明(補題 5.3) (1)⇒ \Rightarrow ⇒ (2)。Z Z Z を含む閉長方形 R 0 R_0 R 0 をとります。1 Z \mathbf 1_Z 1 Z は R 0 R_0 R 0 上可積分で積分値は 0 0 0 なので、命題 3.3 により S ( 1 Z , Δ ) < ε S(\mathbf 1_Z,\Delta)<\varepsilon S ( 1 Z , Δ ) < ε となる分割 Δ \Delta Δ があります。sup 1 Z = 1 \sup\mathbf 1_Z=1 sup 1 Z = 1 となる小長方形は Z Z Z と交わるものだけで、それらが Z Z Z を覆い、総面積は S ( 1 Z , Δ ) S(\mathbf 1_Z,\Delta) S ( 1 Z , Δ ) に等しいので ε \varepsilon ε 未満です。
(2)⇒ \Rightarrow ⇒ (1)。ε > 0 \varepsilon>0 ε > 0 とし、総面積 ε \varepsilon ε 未満の閉長方形 P 1 , … , P m P_1,\dots,P_m P 1 , … , P m が Z Z Z を覆うとします。各 P l P_l P l を、面積が 2 2 2 倍未満になる範囲で少し広げた閉長方形 P l ′ P_l' P l ′ で、P l P_l P l が P l ′ P_l' P l ′ の内部に入るものに置き換えます。P l ′ P_l' P l ′ をすべて含む閉長方形 R 0 ⊃ Z R_0\supset Z R 0 ⊃ Z を取り、その分割 Δ \Delta Δ として、すべての P l ′ P_l' P l ′ の辺(R 0 R_0 R 0 の内部を通るもの)を分点に含むものを取ります。Z Z Z と交わる小長方形は、どれかの P l P_l P l 、したがって P l ′ P_l' P l ′ の内部と交わりますが、小長方形の内部には分点が無いので、その小長方形は P l ′ P_l' P l ′ に含まれます。ゆえに S ( 1 Z , Δ ) ≤ ∑ l ∣ P l ′ ∣ < 2 ε S(\mathbf 1_Z,\Delta)\le\sum_l\lvert P_l'\rvert<2\varepsilon S ( 1 Z , Δ ) ≤ ∑ l ∣ P l ′ ∣ < 2 ε で、s ( 1 Z , Δ ) ≥ 0 s(\mathbf 1_Z,\Delta)\ge0 s ( 1 Z , Δ ) ≥ 0 です。ε \varepsilon ε は任意なので 命題 3.3 より 1 Z \mathbf 1_Z 1 Z は可積分で、積分値は inf Δ S = 0 \inf_\Delta S=0 inf Δ S = 0 です。
最後の主張。部分集合は同じ被覆でよく、2 つの集合の合併は ε / 2 \varepsilon/2 ε /2 ずつの被覆を併せればよいです。線分 [ p , q ] × { y 0 } [p,q]\times\{y_0\} [ p , q ] × { y 0 } は長方形 [ p , q ] × [ y 0 − τ , y 0 + τ ] [p,q]\times[y_0-\tau,y_0+\tau] [ p , q ] × [ y 0 − τ , y 0 + τ ] (面積 2 τ ( q − p ) 2\tau(q-p) 2 τ ( q − p ) )で覆え、τ \tau τ はいくらでも小さく取れます。縦の線分も同様です。
∎
命題 5.4 (リーマンの可積分条件(ジョルダン版) ) ✓
g : R → R g\colon R\to\mathbb{R} g : R → R を有界関数とし、その不連続点全体を Z Z Z とする。Z Z Z が 補題 5.3 の条件 2(総面積がいくらでも小さい有限個の閉長方形で覆える。Z Z Z が面積 0 0 0 の集合に含まれればこれを満たす)を満たすならば、g g g は R R R 上可積分である。
証明(命題 5.4) M = sup R ∣ g ∣ M=\sup_R\lvert g\rvert M = sup R ∣ g ∣ とし、ε > 0 \varepsilon>0 ε > 0 を与えます。δ = ε / ( 2 M + 1 ) \delta=\varepsilon/(2M+1) δ = ε / ( 2 M + 1 ) とおきます。
被覆。 仮定により、総面積 δ / 2 \delta/2 δ /2 未満の閉長方形 P 1 , … , P m P_1,\dots,P_m P 1 , … , P m が Z Z Z を覆います。補題 5.3 の証明と同じく、各 P l P_l P l を少し広げた閉長方形 P l ′ P_l' P l ′ (P l P_l P l は P l ′ P_l' P l ′ の内部に入る)で ∑ l ∣ P l ′ ∣ < δ \sum_l\lvert P_l'\rvert<\delta ∑ l ∣ P l ′ ∣ < δ となるものを取ります。R R R の分割 Δ 1 \Delta_1 Δ 1 として、R R R の内部を通る P l ′ P_l' P l ′ のすべての辺を分点に含むものを取ります。Z Z Z と交わる Δ 1 \Delta_1 Δ 1 の小長方形は、前の補題の証明と同じ理由でどれかの P l ′ P_l' P l ′ に含まれるので、それらの総面積は ∑ l ∣ P l ′ ∣ < δ \sum_l\lvert P_l'\rvert<\delta ∑ l ∣ P l ′ ∣ < δ です(小長方形の内部は互いに交わりません)。これらを「悪い小長方形」、残りを「良い小長方形」と呼びます。
良い小長方形。 良い小長方形 Q Q Q は Z Z Z と交わらないので、g g g は Q Q Q のすべての点で連続です(R R R 上の関数として連続なので、Q Q Q への制限も連続です)。Q Q Q は閉長方形なので、補題 3.1[広義積分・パラメータ積分とガンマ関数] より g ∣ Q g|_Q g ∣ Q は一様連続です。良い小長方形は有限個なので、それらすべてに共通の η > 0 \eta>0 η > 0 を、「p , q ∈ Q p,q\in Q p , q ∈ Q 、∥ p − q ∥ < η \lVert p-q\rVert<\eta ∥ p − q ∥ < η ならば ∣ g ( p ) − g ( q ) ∣ < ε \lvert g(p)-g(q)\rvert<\varepsilon ∣ g ( p ) − g ( q )∣ < ε 」となるように取れます。
細分。 Δ 1 \Delta_1 Δ 1 に等間隔の分点を加えて、幅 ∣ Δ 2 ∣ < η \lvert\Delta_2\rvert<\eta ∣ Δ 2 ∣ < η の分割 Δ 2 \Delta_2 Δ 2 を作ります。Δ 2 \Delta_2 Δ 2 の小長方形はどれも Δ 1 \Delta_1 Δ 1 のただ 1 つの小長方形に含まれます。悪い小長方形に含まれるものでは M i j − m i j ≤ 2 M M_{ij}-m_{ij}\le2M M ij − m ij ≤ 2 M で、その面積の合計は δ \delta δ 未満です。良い小長方形に含まれるものでは、小長方形内の 2 点の距離が η \eta η 未満なので M i j − m i j ≤ ε M_{ij}-m_{ij}\le\varepsilon M ij − m ij ≤ ε です。したがって
S ( g , Δ 2 ) − s ( g , Δ 2 ) ≤ 2 M δ + ε ∣ R ∣ ≤ ε ( 1 + ∣ R ∣ ) . S(g,\Delta_2)-s(g,\Delta_2)\le2M\delta+\varepsilon\lvert R\rvert\le\varepsilon\,(1+\lvert R\rvert). S ( g , Δ 2 ) − s ( g , Δ 2 ) ≤ 2 M δ + ε ∣ R ∣ ≤ ε ( 1 + ∣ R ∣) . ε \varepsilon ε は任意なので、命題 3.3 より g g g は可積分です。
∎
定理 5.6 (縦線集合上の累次積分 ) ✓
φ 1 , φ 2 : [ a , b ] → R \varphi_1,\varphi_2\colon[a,b]\to\mathbb{R} φ 1 , φ 2 : [ a , b ] → R を連続関数で φ 1 ≤ φ 2 \varphi_1\le\varphi_2 φ 1 ≤ φ 2 を満たすものとし、
D = { ( x , y ) ∈ R 2 ∣ a ≤ x ≤ b , φ 1 ( x ) ≤ y ≤ φ 2 ( x ) } D = \{(x,y)\in\mathbb{R}^2 \mid a\le x\le b,\ \varphi_1(x)\le y\le \varphi_2(x)\} D = {( x , y ) ∈ R 2 ∣ a ≤ x ≤ b , φ 1 ( x ) ≤ y ≤ φ 2 ( x )} とおく(このような D D D を縦線集合 という)。f : D → R f\colon D\to\mathbb{R} f : D → R が連続ならば、f f f は D D D 上可積分で
∬ D f ( x , y ) d x d y = ∫ a b ( ∫ φ 1 ( x ) φ 2 ( x ) f ( x , y ) d y ) d x \iint_D f(x,y)\,dx\,dy = \int_a^b\!\left(\int_{\varphi_1(x)}^{\varphi_2(x)} f(x,y)\,dy\right)dx ∬ D f ( x , y ) d x d y = ∫ a b ( ∫ φ 1 ( x ) φ 2 ( x ) f ( x , y ) d y ) d x が成り立つ。同様に、ψ 1 ≤ ψ 2 \psi_1\le\psi_2 ψ 1 ≤ ψ 2 を連続関数として E = { ( x , y ) ∣ c ≤ y ≤ d , ψ 1 ( y ) ≤ x ≤ ψ 2 ( y ) } E=\{(x,y)\mid c\le y\le d,\ \psi_1(y)\le x\le\psi_2(y)\} E = {( x , y ) ∣ c ≤ y ≤ d , ψ 1 ( y ) ≤ x ≤ ψ 2 ( y )} (横線集合 )と E E E 上連続な f f f に対しては
∬ E f ( x , y ) d x d y = ∫ c d ( ∫ ψ 1 ( y ) ψ 2 ( y ) f ( x , y ) d x ) d y . \iint_E f(x,y)\,dx\,dy = \int_c^d\!\left(\int_{\psi_1(y)}^{\psi_2(y)} f(x,y)\,dx\right)dy . ∬ E f ( x , y ) d x d y = ∫ c d ( ∫ ψ 1 ( y ) ψ 2 ( y ) f ( x , y ) d x ) d y .
証明(定理 5.6) φ 1 , φ 2 \varphi_1,\varphi_2 φ 1 , φ 2 は連続だから有界で、D D D は有界です。さらに D D D は有界閉集合なので、その上の連続関数 f f f は有界です。そこで D ⊂ R = [ a , b ] × [ c , d ] D\subset R=[a,b]\times[c,d] D ⊂ R = [ a , b ] × [ c , d ] となる閉長方形を取り、定義 5.1 の f ~ \tilde f f ~ を考えます。
(仮定 1 の確認) f ~ \tilde f f ~ の不連続点を調べます。D D D の内点では、その近傍で f ~ = f \tilde f=f f ~ = f であり f f f は連続です。D D D の外部(補集合の内点)では近傍で f ~ ≡ 0 \tilde f\equiv0 f ~ ≡ 0 なので連続です。したがって不連続点は境界 ∂ D \partial D ∂ D に含まれます。∂ D \partial D ∂ D は 2 つのグラフ y = φ 1 ( x ) y=\varphi_1(x) y = φ 1 ( x ) 、y = φ 2 ( x ) y=\varphi_2(x) y = φ 2 ( x ) と、2 本の縦線分 x = a x=a x = a 、x = b x=b x = b の部分集合の合併に含まれます。補題 5.2 により 2 つのグラフは面積 0 0 0 で、補題 5.3 によりそれらと 2 本の縦線分の合併、したがってその部分集合 ∂ D \partial D ∂ D は、総面積がいくらでも小さい有限個の閉長方形で覆えます。よって 命題 5.4 より f ~ \tilde f f ~ は R R R 上可積分です。これで 定理 4.2 の仮定 1 が満たされました。
(仮定 2 の確認) x ∈ [ a , b ] x\in[a,b] x ∈ [ a , b ] を固定します。y ↦ f ~ ( x , y ) y\mapsto\tilde f(x,y) y ↦ f ~ ( x , y ) は、[ c , φ 1 ( x ) ) [c,\varphi_1(x)) [ c , φ 1 ( x )) と ( φ 2 ( x ) , d ] (\varphi_2(x),d] ( φ 2 ( x ) , d ] では 0 0 0 、[ φ 1 ( x ) , φ 2 ( x ) ] [\varphi_1(x),\varphi_2(x)] [ φ 1 ( x ) , φ 2 ( x )] では y ↦ f ( x , y ) y\mapsto f(x,y) y ↦ f ( x , y ) に等しく、f f f が D D D 上連続なのでこの区間上連続です。つまり有界で、不連続点は高々 2 点 y = φ 1 ( x ) , φ 2 ( x ) y=\varphi_1(x),\varphi_2(x) y = φ 1 ( x ) , φ 2 ( x ) しかありません。1 変数の有界関数は不連続点が有限個なら可積分ですから、仮定 2 も満たされ、さらに区間加法性から
∫ c d f ~ ( x , y ) d y = ∫ φ 1 ( x ) φ 2 ( x ) f ( x , y ) d y \int_c^d \tilde f(x,y)\,dy = \int_{\varphi_1(x)}^{\varphi_2(x)} f(x,y)\,dy ∫ c d f ~ ( x , y ) d y = ∫ φ 1 ( x ) φ 2 ( x ) f ( x , y ) d y です(φ 1 ( x ) = φ 2 ( x ) \varphi_1(x)=\varphi_2(x) φ 1 ( x ) = φ 2 ( x ) のときは両辺 0 0 0 )。
以上より 定理 4.2 が適用でき、
∬ D f = ∬ R f ~ = ∫ a b ( ∫ c d f ~ ( x , y ) d y ) d x = ∫ a b ( ∫ φ 1 ( x ) φ 2 ( x ) f ( x , y ) d y ) d x \iint_D f = \iint_R\tilde f = \int_a^b\!\left(\int_c^d \tilde f(x,y)\,dy\right)dx = \int_a^b\!\left(\int_{\varphi_1(x)}^{\varphi_2(x)} f(x,y)\,dy\right)dx ∬ D f = ∬ R f ~ = ∫ a b ( ∫ c d f ~ ( x , y ) d y ) d x = ∫ a b ( ∫ φ 1 ( x ) φ 2 ( x ) f ( x , y ) d y ) d x を得ます。横線集合の場合は x x x と y y y の役割を入れ替えれば同じ議論です。
∎
次の例の三角形 0 ≤ y ≤ x ≤ 1 0\le y\le x\le1 0 ≤ y ≤ x ≤ 1 を、縦に切る場合と横に切る場合で並べると、内側の積分の範囲が違って見えます。
0 1 1 x y は 0 から x まで x を固定して縦に切る 0 1 1 y x は y から 1 まで y を固定して横に切る 同じ三角形 0 ≤ y ≤ x ≤ 1 を縦に切る場合(左)と横に切る場合(右)。太線が内側の積分の範囲
例 5.7 (順序を換えないと計算できない積分 )
∫ 0 1 ( ∫ y 1 e x 2 d x ) d y \int_0^1\!\left(\int_y^1 e^{x^2}\,dx\right)dy ∫ 0 1 ( ∫ y 1 e x 2 d x ) d y を計算します。内側の ∫ y 1 e x 2 d x \int_y^1 e^{x^2}dx ∫ y 1 e x 2 d x は、e x 2 e^{x^2} e x 2 の原始関数が初等関数で書けないため、この順序のままでは手が出ません。
積分領域を読み取ります。外側が 0 ≤ y ≤ 1 0\le y\le 1 0 ≤ y ≤ 1 、内側が y ≤ x ≤ 1 y\le x\le 1 y ≤ x ≤ 1 ですから、
D = { ( x , y ) ∣ 0 ≤ y ≤ 1 , y ≤ x ≤ 1 } D = \{(x,y)\mid 0\le y\le 1,\ y\le x\le 1\} D = {( x , y ) ∣ 0 ≤ y ≤ 1 , y ≤ x ≤ 1 } で、これは ( 0 , 0 ) , ( 1 , 0 ) , ( 1 , 1 ) (0,0),(1,0),(1,1) ( 0 , 0 ) , ( 1 , 0 ) , ( 1 , 1 ) を頂点とする三角形です(0 ≤ y ≤ x ≤ 1 0\le y\le x\le 1 0 ≤ y ≤ x ≤ 1 と書くと見やすいです)。同じ三角形を縦線集合として書き直すと、各 x ∈ [ 0 , 1 ] x\in[0,1] x ∈ [ 0 , 1 ] に対して y y y は 0 0 0 から x x x まで動くので
D = { ( x , y ) ∣ 0 ≤ x ≤ 1 , 0 ≤ y ≤ x } . D = \{(x,y)\mid 0\le x\le 1,\ 0\le y\le x\}. D = {( x , y ) ∣ 0 ≤ x ≤ 1 , 0 ≤ y ≤ x } . 被積分関数 e x 2 e^{x^2} e x 2 は D D D 上連続なので 定理 5.6 が両方の表示に使え、どちらの累次積分も ∬ D e x 2 d x d y \iint_D e^{x^2}dx\,dy ∬ D e x 2 d x d y に等しくなります。したがって
∫ 0 1 ( ∫ y 1 e x 2 d x ) d y = ∫ 0 1 ( ∫ 0 x e x 2 d y ) d x = ∫ 0 1 e x 2 [ y ] 0 x d x = ∫ 0 1 x e x 2 d x . \int_0^1\!\left(\int_y^1 e^{x^2}dx\right)dy
= \int_0^1\!\left(\int_0^x e^{x^2}\,dy\right)dx
= \int_0^1 e^{x^2}\bigl[y\bigr]_{0}^{x}\,dx
= \int_0^1 x e^{x^2}\,dx . ∫ 0 1 ( ∫ y 1 e x 2 d x ) d y = ∫ 0 1 ( ∫ 0 x e x 2 d y ) d x = ∫ 0 1 e x 2 [ y ] 0 x d x = ∫ 0 1 x e x 2 d x . 最後は t = x 2 t=x^2 t = x 2 と置換して(d t = 2 x d x dt=2x\,dx d t = 2 x d x )
∫ 0 1 x e x 2 d x = [ e x 2 2 ] 0 1 = e − 1 2 . \int_0^1 xe^{x^2}dx = \left[\frac{e^{x^2}}{2}\right]_0^1 = \frac{e-1}{2}. ∫ 0 1 x e x 2 d x = [ 2 e x 2 ] 0 1 = 2 e − 1 . 内側の積分に現れた x x x という因子は、三角形の縦の長さ(断面の長さ)です。これが x e x 2 xe^{x^2} x e x 2 の x x x の正体で、順序交換の御利益はここに集約されています。
例 5.8 (2 つの累次積分が食い違う例 )
f ( x , y ) = x 2 − y 2 ( x 2 + y 2 ) 2 f(x,y) = \frac{x^2-y^2}{(x^2+y^2)^2} f ( x , y ) = ( x 2 + y 2 ) 2 x 2 − y 2 を [ 0 , 1 ] 2 ∖ { ( 0 , 0 ) } [0,1]^2\setminus\{(0,0)\} [ 0 , 1 ] 2 ∖ {( 0 , 0 )} 上で考えます。この f f f は原点に近づけると非有界です。実際 y y y を固定して x → 0 + x\to 0^+ x → 0 + とすると f ( x , y ) → − 1 / y 2 f(x,y)\to -1/y^2 f ( x , y ) → − 1/ y 2 ですが、y = x y=x y = x 上では f = 0 f=0 f = 0 、y → 0 + y\to0^+ y → 0 + を先にすると f ( x , y ) → 1 / x 2 f(x,y)\to 1/x^2 f ( x , y ) → 1/ x 2 で、f ( x , x / 2 ) = x 2 − x 2 / 4 ( 5 x 2 / 4 ) 2 = 12 25 x 2 → + ∞ f(x,x/2)= \frac{x^2-x^2/4}{(5x^2/4)^2}=\frac{12}{25x^2}\to+\infty f ( x , x /2 ) = ( 5 x 2 /4 ) 2 x 2 − x 2 /4 = 25 x 2 12 → + ∞ です。したがって f f f は [ 0 , 1 ] 2 [0,1]^2 [ 0 , 1 ] 2 上有界でなく、定義 3.2 の意味ではそもそも重積分が定義されません 。
それでも 2 つの累次積分は計算できます。x > 0 x>0 x > 0 を固定すると y ↦ f ( x , y ) y\mapsto f(x,y) y ↦ f ( x , y ) は [ 0 , 1 ] [0,1] [ 0 , 1 ] 上連続です(分母は x 2 x^2 x 2 以上)。商の微分から
∂ ∂ y ( y x 2 + y 2 ) = ( x 2 + y 2 ) − y ⋅ 2 y ( x 2 + y 2 ) 2 = x 2 − y 2 ( x 2 + y 2 ) 2 = f ( x , y ) \frac{\partial}{\partial y}\left(\frac{y}{x^2+y^2}\right)
= \frac{(x^2+y^2)-y\cdot 2y}{(x^2+y^2)^2}
= \frac{x^2-y^2}{(x^2+y^2)^2} = f(x,y) ∂ y ∂ ( x 2 + y 2 y ) = ( x 2 + y 2 ) 2 ( x 2 + y 2 ) − y ⋅ 2 y = ( x 2 + y 2 ) 2 x 2 − y 2 = f ( x , y ) なので
∫ 0 1 f ( x , y ) d y = [ y x 2 + y 2 ] y = 0 y = 1 = 1 x 2 + 1 . \int_0^1 f(x,y)\,dy = \left[\frac{y}{x^2+y^2}\right]_{y=0}^{y=1} = \frac{1}{x^2+1}. ∫ 0 1 f ( x , y ) d y = [ x 2 + y 2 y ] y = 0 y = 1 = x 2 + 1 1 . これを x x x で積分して
∫ 0 1 ( ∫ 0 1 f ( x , y ) d y ) d x = ∫ 0 1 d x 1 + x 2 = [ arctan x ] 0 1 = π 4 . \int_0^1\!\left(\int_0^1 f(x,y)\,dy\right)dx = \int_0^1\frac{dx}{1+x^2} = \Bigl[\arctan x\Bigr]_0^1 = \frac{\pi}{4}. ∫ 0 1 ( ∫ 0 1 f ( x , y ) d y ) d x = ∫ 0 1 1 + x 2 d x = [ arctan x ] 0 1 = 4 π . 一方 f ( y , x ) = y 2 − x 2 ( x 2 + y 2 ) 2 = − f ( x , y ) f(y,x) = \dfrac{y^2-x^2}{(x^2+y^2)^2} = -f(x,y) f ( y , x ) = ( x 2 + y 2 ) 2 y 2 − x 2 = − f ( x , y ) という反対称性があるので、x x x と y y y の名前を入れ替えて同じ計算をすれば
∫ 0 1 ( ∫ 0 1 f ( x , y ) d x ) d y = − π 4 . \int_0^1\!\left(\int_0^1 f(x,y)\,dx\right)dy = -\frac{\pi}{4}. ∫ 0 1 ( ∫ 0 1 f ( x , y ) d x ) d y = − 4 π . (x = 0 x=0 x = 0 では内側の積分 ∫ 0 1 f ( 0 , y ) d y = ∫ 0 1 ( − y − 2 ) d y \int_0^1 f(0,y)\,dy=\int_0^1(-y^{-2})dy ∫ 0 1 f ( 0 , y ) d y = ∫ 0 1 ( − y − 2 ) d y が発散するので、外側の積分は x ∈ ( 0 , 1 ] x\in(0,1] x ∈ ( 0 , 1 ] 上の広義積分 として読みます。1 / ( 1 + x 2 ) 1/(1+x^2) 1/ ( 1 + x 2 ) は x = 0 x=0 x = 0 まで連続に延びるので値は π / 4 \pi/4 π /4 のままです。)
π / 4 ≠ − π / 4 \pi/4 \ne -\pi/4 π /4 = − π /4 ですから、2 つの累次積分は食い違います。これは 定理 4.2 と矛盾しません。仮定 1(R R R 上可積分)が破れているからです。フビニの定理は「重積分が存在すること」を前提に、それを累次積分に翻訳する定理であって、2 つの累次積分が勝手に一致することを主張するものではありません。
1 変数の 置換積分(定理 6.1)[積分の基本定理と定積分] は、g g g が C 1 C^1 C 1 級で単調なら
∫ g ( α ) g ( β ) f ( x ) d x = ∫ α β f ( g ( t ) ) g ′ ( t ) d t \int_{g(\alpha)}^{g(\beta)} f(x)\,dx = \int_{\alpha}^{\beta} f(g(t))\,g'(t)\,dt ∫ g ( α ) g ( β ) f ( x ) d x = ∫ α β f ( g ( t )) g ′ ( t ) d t
でした。向きを気にせず「長さ」だけで読むと、g g g は微小区間 [ t , t + d t ] [t,t+dt] [ t , t + d t ] を長さ ∣ g ′ ( t ) ∣ d t \lvert g'(t)\rvert\,dt ∣ g ′ ( t )∣ d t の区間に写すので、d x = ∣ g ′ ( t ) ∣ d t dx = \lvert g'(t)\rvert\,dt d x = ∣ g ′ ( t )∣ d t という置き換えが起きています。2 変数では d x d y dx\,dy d x d y が面積なので、「微小な面積を何倍にするか」を表す量が必要になります。それがヤコビアンです。
一次変換で面積が何倍になるかを、定義 5.1 の面積の定義と 定理 5.6 だけから導きます。「平行四辺形の面積 = 底辺 × 高さ」や、なす角の sin \sin sin は使いません(三角関数を級数から定義し、扇形の面積を重積分で計算する 初等関数を級数から定義する の議論がこの節に依存するので、ここで幾何の公式を使うと循環になります)。準備として、平行移動・拡大と面積の関係から始めます。
行列式 の絶対値が面積の拡大率になることは既に知っている方も多いと思いますが、以下の補題が続くのは、それを重積分の定義だけから示すためです。方針は、(1) 平行移動・拡大・線分の面積を押さえる、(2) 平行四辺形が縦線集合として書けることを示し、縦方向・横方向の「ずらし」で面積が変わらないことを確かめる、(3) ずらしを 2 回行って平行四辺形を長方形に変え、その面積を直接計算する、というものです。
補題 6.1 (平行移動と拡大による面積の変化 ) ✓
E ⊂ R 2 E\subset\mathbb{R}^2 E ⊂ R 2 を面積確定な有界集合、c ∈ R 2 \boldsymbol c\in\mathbb{R}^2 c ∈ R 2 、h > 0 h>0 h > 0 とし、Ψ ( z ) = c + h z \Psi(\boldsymbol z)=\boldsymbol c+h\boldsymbol z Ψ ( z ) = c + h z とおく。このとき Ψ ( E ) = { c + h z ∣ z ∈ E } \Psi(E)=\{\boldsymbol c+h\boldsymbol z\mid\boldsymbol z\in E\} Ψ ( E ) = { c + h z ∣ z ∈ E } は面積確定で、∣ Ψ ( E ) ∣ = h 2 ∣ E ∣ \lvert\Psi(E)\rvert=h^2\lvert E\rvert ∣ Ψ ( E )∣ = h 2 ∣ E ∣ である。とくに(h = 1 h=1 h = 1 )平行移動 c + E \boldsymbol c+E c + E は面積確定で ∣ c + E ∣ = ∣ E ∣ \lvert\boldsymbol c+E\rvert=\lvert E\rvert ∣ c + E ∣ = ∣ E ∣ 、(c = 0 \boldsymbol c=\boldsymbol 0 c = 0 )拡大 h E hE h E は面積確定で ∣ h E ∣ = h 2 ∣ E ∣ \lvert hE\rvert=h^2\lvert E\rvert ∣ h E ∣ = h 2 ∣ E ∣ である。
証明(補題 6.1) E ⊂ R = [ a , b ] × [ c ′ , d ′ ] E\subset R=[a,b]\times[c',d'] E ⊂ R = [ a , b ] × [ c ′ , d ′ ] となる閉長方形を取り、c = ( c 1 , c 2 ) \boldsymbol c=(c_1,c_2) c = ( c 1 , c 2 ) とします。Ψ \Psi Ψ は各座標を t ↦ c i + h t t\mapsto c_i+ht t ↦ c i + h t で写す(h > 0 h>0 h > 0 なので順序を保つ)全単射なので、Ψ ( R ) = [ c 1 + h a , c 1 + h b ] × [ c 2 + h c ′ , c 2 + h d ′ ] \Psi(R)=[c_1+ha,\,c_1+hb]\times[c_2+hc',\,c_2+hd'] Ψ ( R ) = [ c 1 + ha , c 1 + hb ] × [ c 2 + h c ′ , c 2 + h d ′ ] は Ψ ( E ) \Psi(E) Ψ ( E ) を含む閉長方形です。
R R R の分割 Δ \Delta Δ (分点 x i x_i x i 、y j y_j y j )に対し、分点 c 1 + h x i c_1+hx_i c 1 + h x i 、c 2 + h y j c_2+hy_j c 2 + h y j で決まる Ψ ( R ) \Psi(R) Ψ ( R ) の分割を Ψ Δ \Psi\Delta ΨΔ と書きます。逆に Ψ ( R ) \Psi(R) Ψ ( R ) の分割は分点を Ψ − 1 \Psi^{-1} Ψ − 1 で戻せば R R R の分割になるので、Ψ ( R ) \Psi(R) Ψ ( R ) の分割はすべて Ψ Δ \Psi\Delta ΨΔ の形です。Ψ Δ \Psi\Delta ΨΔ の小長方形は Ψ ( R i j ) \Psi(R_{ij}) Ψ ( R ij ) で、辺の長さが h Δ x i h\Delta x_i h Δ x i 、h Δ y j h\Delta y_j h Δ y j なので面積は h 2 ∣ R i j ∣ h^2\lvert R_{ij}\rvert h 2 ∣ R ij ∣ です。
g = 1 E g=\mathbf 1_E g = 1 E (R R R 上)、g ′ = 1 Ψ ( E ) g'=\mathbf 1_{\Psi(E)} g ′ = 1 Ψ ( E ) (Ψ ( R ) \Psi(R) Ψ ( R ) 上)とおきます。Ψ \Psi Ψ は全単射なので g ′ ( Ψ ( z ) ) = g ( z ) g'(\Psi(\boldsymbol z))=g(\boldsymbol z) g ′ ( Ψ ( z )) = g ( z ) で、g ′ g' g ′ の Ψ ( R i j ) \Psi(R_{ij}) Ψ ( R ij ) 上の sup \sup sup ・inf \inf inf は g g g の R i j R_{ij} R ij 上の sup \sup sup ・inf \inf inf に等しい。したがって
S ( g ′ , Ψ Δ ) = h 2 S ( g , Δ ) , s ( g ′ , Ψ Δ ) = h 2 s ( g , Δ ) S(g',\Psi\Delta)=h^2S(g,\Delta),\qquad s(g',\Psi\Delta)=h^2s(g,\Delta) S ( g ′ , ΨΔ ) = h 2 S ( g , Δ ) , s ( g ′ , ΨΔ ) = h 2 s ( g , Δ ) です。g g g は可積分なので、命題 3.3 より任意の ε > 0 \varepsilon>0 ε > 0 に対し S ( g , Δ ) − s ( g , Δ ) < ε / h 2 S(g,\Delta)-s(g,\Delta)<\varepsilon/h^2 S ( g , Δ ) − s ( g , Δ ) < ε / h 2 となる Δ \Delta Δ があり、そのとき S ( g ′ , Ψ Δ ) − s ( g ′ , Ψ Δ ) < ε S(g',\Psi\Delta)-s(g',\Psi\Delta)<\varepsilon S ( g ′ , ΨΔ ) − s ( g ′ , ΨΔ ) < ε です。再び 命題 3.3 より g ′ g' g ′ は可積分、すなわち Ψ ( E ) \Psi(E) Ψ ( E ) は面積確定で、
∣ Ψ ( E ) ∣ = inf Δ ′ S ( g ′ , Δ ′ ) = inf Δ h 2 S ( g , Δ ) = h 2 ∣ E ∣ \lvert\Psi(E)\rvert=\inf_{\Delta'}S(g',\Delta')=\inf_{\Delta}h^2S(g,\Delta)=h^2\lvert E\rvert ∣ Ψ ( E )∣ = Δ ′ inf S ( g ′ , Δ ′ ) = Δ inf h 2 S ( g , Δ ) = h 2 ∣ E ∣ です(Δ ′ \Delta' Δ ′ は Ψ ( R ) \Psi(R) Ψ ( R ) の分割全体を動き、それは Ψ Δ \Psi\Delta ΨΔ の全体と一致する)。
∎
補題 6.2 (線分は面積 0 ) ✓
p , q ∈ R 2 \boldsymbol p,\boldsymbol q\in\mathbb{R}^2 p , q ∈ R 2 を結ぶ線分 [ p , q ] = { p + t ( q − p ) ∣ 0 ≤ t ≤ 1 } [\boldsymbol p,\boldsymbol q]=\{\boldsymbol p+t(\boldsymbol q-\boldsymbol p)\mid 0\le t\le1\} [ p , q ] = { p + t ( q − p ) ∣ 0 ≤ t ≤ 1 } は面積確定で、その面積は 0 0 0 である。その部分集合も同様である。
証明(補題 6.2) p = ( p 1 , p 2 ) \boldsymbol p=(p_1,p_2) p = ( p 1 , p 2 ) 、q − p = ( d 1 , d 2 ) \boldsymbol q-\boldsymbol p=(d_1,d_2) q − p = ( d 1 , d 2 ) とします。d 1 ≠ 0 d_1\ne0 d 1 = 0 なら、x = p 1 + t d 1 x=p_1+td_1 x = p 1 + t d 1 と t t t が 1 対 1 に対応するので、線分は 1 次関数 φ ( x ) = p 2 + d 2 d 1 ( x − p 1 ) \varphi(x)=p_2+\dfrac{d_2}{d_1}(x-p_1) φ ( x ) = p 2 + d 1 d 2 ( x − p 1 ) (多項式なので連続。系 5.2[極限と連続性] )の、p 1 p_1 p 1 と q 1 q_1 q 1 を端点とする閉区間上のグラフです。補題 5.2 より面積確定で面積 0 0 0 です。d 1 = 0 d_1=0 d 1 = 0 なら、線分は縦の線分 { p 1 } × [ min ( p 2 , q 2 ) , max ( p 2 , q 2 ) ] \{p_1\}\times[\min(p_2,q_2),\max(p_2,q_2)] { p 1 } × [ min ( p 2 , q 2 ) , max ( p 2 , q 2 )] (d 2 = 0 d_2=0 d 2 = 0 なら 1 点)で、一辺 τ \tau τ の正方形を縦に ⌈ ∣ d 2 ∣ / τ ⌉ + 1 \lceil\lvert d_2\rvert/\tau\rceil+1 ⌈∣ d 2 ∣ / τ ⌉ + 1 個並べれば覆え、総面積 τ ( ∣ d 2 ∣ + 2 τ ) \tau(\lvert d_2\rvert+2\tau) τ (∣ d 2 ∣ + 2 τ ) は τ \tau τ を小さくすればいくらでも小さくなります。補題 5.3 の (2)⇒ \Rightarrow ⇒ (1) より面積確定で面積 0 0 0 です。部分集合は同じ被覆で覆えるので、補題 5.3 より同じ結論です。
∎
以下、2 次正方行列 B = ( b 11 b 12 b 21 b 22 ) B=\begin{pmatrix}b_{11}&b_{12}\\ b_{21}&b_{22}\end{pmatrix} B = ( b 11 b 21 b 12 b 22 ) の列ベクトル b 1 = ( b 11 , b 21 ) \boldsymbol b_1=(b_{11},b_{21}) b 1 = ( b 11 , b 21 ) 、b 2 = ( b 12 , b 22 ) \boldsymbol b_2=(b_{12},b_{22}) b 2 = ( b 12 , b 22 ) が張る平行四辺形を
P ( B ) = { u b 1 + v b 2 ∣ 0 ≤ u ≤ 1 , 0 ≤ v ≤ 1 } = B ( [ 0 , 1 ] 2 ) P(B)=\{u\boldsymbol b_1+v\boldsymbol b_2\mid 0\le u\le1,\ 0\le v\le1\}=B\bigl([0,1]^2\bigr) P ( B ) = { u b 1 + v b 2 ∣ 0 ≤ u ≤ 1 , 0 ≤ v ≤ 1 } = B ( [ 0 , 1 ] 2 )
と書きます。また、縦方向のずらし S k ( x , y ) = ( x , y + k x ) S_k(x,y)=(x,\,y+kx) S k ( x , y ) = ( x , y + k x ) と横方向のずらし H m ( x , y ) = ( x + m y , y ) H_m(x,y)=(x+my,\,y) H m ( x , y ) = ( x + m y , y ) (k , m ∈ R k,m\in\mathbb{R} k , m ∈ R )を使います。どちらも一次変換なので、S k ( P ( B ) ) = P ( S k B ) S_k(P(B))=P(S_kB) S k ( P ( B )) = P ( S k B ) 、H m ( P ( B ) ) = P ( H m B ) H_m(P(B))=P(H_mB) H m ( P ( B )) = P ( H m B ) です(S k B S_kB S k B は列ベクトルを S k S_k S k で写した行列)。
補題 6.3 (平行四辺形は縦線集合であり、ずらしで面積は変わらない ) ✓
det B ≠ 0 \det B\ne0 det B = 0 とする。
実数 a ≤ b a\le b a ≤ b と [ a , b ] [a,b] [ a , b ] 上の連続関数 φ 1 ≤ φ 2 \varphi_1\le\varphi_2 φ 1 ≤ φ 2 があって、P ( B ) = { ( x , y ) ∣ a ≤ x ≤ b , φ 1 ( x ) ≤ y ≤ φ 2 ( x ) } P(B)=\{(x,y)\mid a\le x\le b,\ \varphi_1(x)\le y\le\varphi_2(x)\} P ( B ) = {( x , y ) ∣ a ≤ x ≤ b , φ 1 ( x ) ≤ y ≤ φ 2 ( x )} (縦線集合)と書ける。同様に P ( B ) P(B) P ( B ) は横線集合としても書ける。
P ( B ) P(B) P ( B ) は面積確定である。さらに任意の k , m k,m k , m に対し S k ( P ( B ) ) S_k(P(B)) S k ( P ( B )) 、H m ( P ( B ) ) H_m(P(B)) H m ( P ( B )) も面積確定で、∣ S k ( P ( B ) ) ∣ = ∣ H m ( P ( B ) ) ∣ = ∣ P ( B ) ∣ \lvert S_k(P(B))\rvert=\lvert H_m(P(B))\rvert=\lvert P(B)\rvert ∣ S k ( P ( B ))∣ = ∣ H m ( P ( B ))∣ = ∣ P ( B )∣ である。
証明(補題 6.3) D = det B = b 11 b 22 − b 12 b 21 ≠ 0 D=\det B=b_{11}b_{22}-b_{12}b_{21}\ne0 D = det B = b 11 b 22 − b 12 b 21 = 0 とし、I D I_D I D を 0 0 0 と D D D を端点とする閉区間とします。
(1) ( x , y ) = u b 1 + v b 2 (x,y)=u\boldsymbol b_1+v\boldsymbol b_2 ( x , y ) = u b 1 + v b 2 すなわち x = u b 11 + v b 12 x=ub_{11}+vb_{12} x = u b 11 + v b 12 、y = u b 21 + v b 22 y=ub_{21}+vb_{22} y = u b 21 + v b 22 を満たす ( u , v ) (u,v) ( u , v ) は、第 1 式の b 22 b_{22} b 22 倍から第 2 式の b 12 b_{12} b 12 倍を引き、また第 2 式の b 11 b_{11} b 11 倍から第 1 式の b 21 b_{21} b 21 倍を引くと
u = b 22 x − b 12 y D , v = − b 21 x + b 11 y D u=\frac{b_{22}x-b_{12}y}{D},\qquad v=\frac{-b_{21}x+b_{11}y}{D} u = D b 22 x − b 12 y , v = D − b 21 x + b 11 y に限られ、逆にこの ( u , v ) (u,v) ( u , v ) は 2 式を満たします(代入して確かめられる)。したがって
( x , y ) ∈ P ( B ) ⟺ b 22 x − b 12 y ∈ I D かつ − b 21 x + b 11 y ∈ I D (x,y)\in P(B)\iff b_{22}x-b_{12}y\in I_D\ \text{かつ}\ -b_{21}x+b_{11}y\in I_D ( x , y ) ∈ P ( B ) ⟺ b 22 x − b 12 y ∈ I D かつ − b 21 x + b 11 y ∈ I D です。それぞれの条件を書き直します。
第 1 の条件。b 12 ≠ 0 b_{12}\ne0 b 12 = 0 なら、y y y が b 22 x b 12 \dfrac{b_{22}x}{b_{12}} b 12 b 22 x と b 22 x − D b 12 \dfrac{b_{22}x-D}{b_{12}} b 12 b 22 x − D を端点とする閉区間にあることと同値です。2 つの端点の差は定数 D / b 12 D/b_{12} D / b 12 なので、小さいほうは L 1 ( x ) = b 22 x b 12 − max ( 0 , D b 12 ) L_1(x)=\dfrac{b_{22}x}{b_{12}}-\max\!\left(0,\dfrac D{b_{12}}\right) L 1 ( x ) = b 12 b 22 x − max ( 0 , b 12 D ) 、大きいほうは U 1 ( x ) = L 1 ( x ) + ∣ D ∣ ∣ b 12 ∣ U_1(x)=L_1(x)+\dfrac{\lvert D\rvert}{\lvert b_{12}\rvert} U 1 ( x ) = L 1 ( x ) + ∣ b 12 ∣ ∣ D ∣ で、どちらも x x x の 1 次式です。b 12 = 0 b_{12}=0 b 12 = 0 なら D = b 11 b 22 D=b_{11}b_{22} D = b 11 b 22 で b 22 ≠ 0 b_{22}\ne0 b 22 = 0 なので、条件は x x x が 0 0 0 と b 11 b_{11} b 11 を端点とする閉区間 J 1 J_1 J 1 にあることと同値です。
第 2 の条件。b 11 ≠ 0 b_{11}\ne0 b 11 = 0 なら同様に、y y y が 1 次式 L 2 ( x ) L_2(x) L 2 ( x ) と U 2 ( x ) = L 2 ( x ) + ∣ D ∣ ∣ b 11 ∣ U_2(x)=L_2(x)+\dfrac{\lvert D\rvert}{\lvert b_{11}\rvert} U 2 ( x ) = L 2 ( x ) + ∣ b 11 ∣ ∣ D ∣ の間にあることと同値です(端点は b 21 x b 11 \dfrac{b_{21}x}{b_{11}} b 11 b 21 x と b 21 x + D b 11 \dfrac{b_{21}x+D}{b_{11}} b 11 b 21 x + D )。b 11 = 0 b_{11}=0 b 11 = 0 なら D = − b 12 b 21 D=-b_{12}b_{21} D = − b 12 b 21 で b 21 ≠ 0 b_{21}\ne0 b 21 = 0 なので、条件は x x x が 0 0 0 と b 12 b_{12} b 12 を端点とする閉区間 J 2 J_2 J 2 にあることと同値です。
b 11 = b 12 = 0 b_{11}=b_{12}=0 b 11 = b 12 = 0 なら D = 0 D=0 D = 0 なので、少なくとも一方の条件は y y y の区間の条件です。
b 12 = 0 b_{12}=0 b 12 = 0 のときは P ( B ) = { x ∈ J 1 , L 2 ( x ) ≤ y ≤ U 2 ( x ) } P(B)=\{x\in J_1,\ L_2(x)\le y\le U_2(x)\} P ( B ) = { x ∈ J 1 , L 2 ( x ) ≤ y ≤ U 2 ( x )} 、b 11 = 0 b_{11}=0 b 11 = 0 のときは P ( B ) = { x ∈ J 2 , L 1 ( x ) ≤ y ≤ U 1 ( x ) } P(B)=\{x\in J_2,\ L_1(x)\le y\le U_1(x)\} P ( B ) = { x ∈ J 2 , L 1 ( x ) ≤ y ≤ U 1 ( x )} で、主張の形です。b 11 , b 12 b_{11},b_{12} b 11 , b 12 がともに 0 0 0 でないときは、φ 1 = max ( L 1 , L 2 ) \varphi_1=\max(L_1,L_2) φ 1 = max ( L 1 , L 2 ) 、φ 2 = min ( U 1 , U 2 ) \varphi_2=\min(U_1,U_2) φ 2 = min ( U 1 , U 2 ) (連続関数の最大・最小なので連続)とおけば P ( B ) = { ( x , y ) ∣ x ∈ X , φ 1 ( x ) ≤ y ≤ φ 2 ( x ) } P(B)=\{(x,y)\mid x\in X,\ \varphi_1(x)\le y\le\varphi_2(x)\} P ( B ) = {( x , y ) ∣ x ∈ X , φ 1 ( x ) ≤ y ≤ φ 2 ( x )} 、X = { x ∣ φ 1 ( x ) ≤ φ 2 ( x ) } X=\{x\mid\varphi_1(x)\le\varphi_2(x)\} X = { x ∣ φ 1 ( x ) ≤ φ 2 ( x )} です。L i ≤ U i L_i\le U_i L i ≤ U i なので X = { x ∣ L 1 ( x ) ≤ U 2 ( x ) } ∩ { x ∣ L 2 ( x ) ≤ U 1 ( x ) } X=\{x\mid L_1(x)\le U_2(x)\}\cap\{x\mid L_2(x)\le U_1(x)\} X = { x ∣ L 1 ( x ) ≤ U 2 ( x )} ∩ { x ∣ L 2 ( x ) ≤ U 1 ( x )} で、各集合は 1 次式の不等式の解なので R \mathbb{R} R 、空集合、または ( − ∞ , t ] (-\infty,t] ( − ∞ , t ] ・[ t , ∞ ) [t,\infty) [ t , ∞ ) の形の閉半直線です。X X X は x = 0 x=0 x = 0 を含み(u = v = 0 u=v=0 u = v = 0 の点 0 ∈ P ( B ) \boldsymbol 0\in P(B) 0 ∈ P ( B ) )、x ∈ X x\in X x ∈ X なら ( x , φ 1 ( x ) ) ∈ P ( B ) (x,\varphi_1(x))\in P(B) ( x , φ 1 ( x )) ∈ P ( B ) なので ∣ x ∣ ≤ ∣ b 11 ∣ + ∣ b 12 ∣ \lvert x\rvert\le\lvert b_{11}\rvert+\lvert b_{12}\rvert ∣ x ∣ ≤ ∣ b 11 ∣ + ∣ b 12 ∣ と有界です。有界で空でない、2 つの閉半直線(または R \mathbb{R} R )の共通部分は閉区間 [ a , b ] [a,b] [ a , b ] です。
横線集合としての表示は、x x x と y y y の役割を入れ替えた行列 ( b 21 b 22 b 11 b 12 ) \begin{pmatrix}b_{21}&b_{22}\\ b_{11}&b_{12}\end{pmatrix} ( b 21 b 11 b 22 b 12 ) (行列式は − D ≠ 0 -D\ne0 − D = 0 )に同じ議論を使えば得られます。
(2) (1) の縦線集合の表示と 定理 5.6 (f ≡ 1 f\equiv1 f ≡ 1 )により、P ( B ) P(B) P ( B ) は面積確定で
∣ P ( B ) ∣ = ∫ a b ( ∫ φ 1 ( x ) φ 2 ( x ) 1 d y ) d x = ∫ a b ( φ 2 ( x ) − φ 1 ( x ) ) d x \lvert P(B)\rvert=\int_a^b\Bigl(\int_{\varphi_1(x)}^{\varphi_2(x)}1\,dy\Bigr)dx=\int_a^b\bigl(\varphi_2(x)-\varphi_1(x)\bigr)\,dx ∣ P ( B )∣ = ∫ a b ( ∫ φ 1 ( x ) φ 2 ( x ) 1 d y ) d x = ∫ a b ( φ 2 ( x ) − φ 1 ( x ) ) d x です(定数の積分 ∫ α β c d t = c ( β − α ) \int_\alpha^\beta c\,dt=c(\beta-\alpha) ∫ α β c d t = c ( β − α ) は、F ( t ) = c t F(t)=ct F ( t ) = c t に 定理 5.4[積分の基本定理と定積分] を使えば得られる)。( x , y ) ∈ P ( B ) (x,y)\in P(B) ( x , y ) ∈ P ( B ) と ( x , y + k x ) ∈ S k ( P ( B ) ) (x,y+kx)\in S_k(P(B)) ( x , y + k x ) ∈ S k ( P ( B )) は同値なので
S k ( P ( B ) ) = { ( x , y ) ∣ a ≤ x ≤ b , φ 1 ( x ) + k x ≤ y ≤ φ 2 ( x ) + k x } S_k(P(B))=\{(x,y)\mid a\le x\le b,\ \varphi_1(x)+kx\le y\le\varphi_2(x)+kx\} S k ( P ( B )) = {( x , y ) ∣ a ≤ x ≤ b , φ 1 ( x ) + k x ≤ y ≤ φ 2 ( x ) + k x } で、これも連続関数で挟まれた縦線集合です。定理 5.6 から面積確定で、面積は ∫ a b ( ( φ 2 ( x ) + k x ) − ( φ 1 ( x ) + k x ) ) d x = ∣ P ( B ) ∣ \int_a^b\bigl((\varphi_2(x)+kx)-(\varphi_1(x)+kx)\bigr)dx=\lvert P(B)\rvert ∫ a b ( ( φ 2 ( x ) + k x ) − ( φ 1 ( x ) + k x ) ) d x = ∣ P ( B )∣ です。H m H_m H m については、(1) の横線集合の表示を使って x x x と y y y の役割を入れ替えた同じ議論をします。
∎
補題 6.4 (平行四辺形の面積は行列式の絶対値 ) ✓
任意の 2 次正方行列 B B B に対し、P ( B ) P(B) P ( B ) は面積確定で ∣ P ( B ) ∣ = ∣ det B ∣ \lvert P(B)\rvert=\lvert\det B\rvert ∣ P ( B )∣ = ∣ det B ∣ である。
証明(補題 6.4) det B = 0 \det B=0 det B = 0 の場合。 b 1 = b 2 = 0 \boldsymbol b_1=\boldsymbol b_2=\boldsymbol 0 b 1 = b 2 = 0 なら P ( B ) = { 0 } P(B)=\{\boldsymbol 0\} P ( B ) = { 0 } で、1 点は線分の部分集合なので 補題 6.2 より面積 0 0 0 です。そうでなければ、b 1 , b 2 \boldsymbol b_1,\boldsymbol b_2 b 1 , b 2 のうち 0 \boldsymbol 0 0 でないほうを w = ( w 1 , w 2 ) \boldsymbol w=(w_1,w_2) w = ( w 1 , w 2 ) 、他方を w ′ = ( w 1 ′ , w 2 ′ ) \boldsymbol w'=(w_1',w_2') w ′ = ( w 1 ′ , w 2 ′ ) とします。行列式が 0 0 0 なので w 1 w 2 ′ = w 2 w 1 ′ w_1w_2'=w_2w_1' w 1 w 2 ′ = w 2 w 1 ′ です。w 1 ≠ 0 w_1\ne0 w 1 = 0 なら λ = w 1 ′ / w 1 \lambda=w_1'/w_1 λ = w 1 ′ / w 1 とおくと λ w 1 = w 1 ′ \lambda w_1=w_1' λ w 1 = w 1 ′ 、λ w 2 = w 2 w 1 ′ / w 1 = w 2 ′ \lambda w_2=w_2w_1'/w_1=w_2' λ w 2 = w 2 w 1 ′ / w 1 = w 2 ′ です。w 1 = 0 w_1=0 w 1 = 0 なら w 2 ≠ 0 w_2\ne0 w 2 = 0 で、w 2 w 1 ′ = 0 w_2w_1'=0 w 2 w 1 ′ = 0 から w 1 ′ = 0 w_1'=0 w 1 ′ = 0 、λ = w 2 ′ / w 2 \lambda=w_2'/w_2 λ = w 2 ′ / w 2 とおけば同様です。いずれにせよ w ′ = λ w \boldsymbol w'=\lambda\boldsymbol w w ′ = λ w なので、P ( B ) P(B) P ( B ) の点は t w t\boldsymbol w t w (∣ t ∣ ≤ 1 + ∣ λ ∣ = : M \lvert t\rvert\le1+\lvert\lambda\rvert=:M ∣ t ∣ ≤ 1 + ∣ λ ∣ =: M )の形で、P ( B ) P(B) P ( B ) は線分 [ − M w , M w ] [-M\boldsymbol w,M\boldsymbol w] [ − M w , M w ] に含まれます。補題 6.2 より P ( B ) P(B) P ( B ) は面積確定で、面積は 0 = ∣ det B ∣ 0=\lvert\det B\rvert 0 = ∣ det B ∣ です。
det B ≠ 0 \det B\ne0 det B = 0 の場合。 D = det B D=\det B D = det B とします。b 11 = 0 b_{11}=0 b 11 = 0 なら b 1 \boldsymbol b_1 b 1 と b 2 \boldsymbol b_2 b 2 を入れ替えます。P ( B ) P(B) P ( B ) は変わらず、行列式は符号だけが変わり、入れ替えた後の ( 1 , 1 ) (1,1) ( 1 , 1 ) 成分 b 12 b_{12} b 12 は 0 0 0 ではありません(b 11 = b 12 = 0 b_{11}=b_{12}=0 b 11 = b 12 = 0 なら D = 0 D=0 D = 0 )。そこで b 11 ≠ 0 b_{11}\ne0 b 11 = 0 と仮定してよい。
k = − b 21 / b 11 k=-b_{21}/b_{11} k = − b 21 / b 11 とすると S k b 1 = ( b 11 , 0 ) S_k\boldsymbol b_1=(b_{11},0) S k b 1 = ( b 11 , 0 ) 、S k b 2 = ( b 12 , b 22 + k b 12 ) = ( b 12 , D / b 11 ) S_k\boldsymbol b_2=(b_{12},\,b_{22}+kb_{12})=(b_{12},\,D/b_{11}) S k b 2 = ( b 12 , b 22 + k b 12 ) = ( b 12 , D / b 11 ) です。B 1 = S k B B_1=S_kB B 1 = S k B の行列式は b 11 ⋅ D / b 11 = D ≠ 0 b_{11}\cdot D/b_{11}=D\ne0 b 11 ⋅ D / b 11 = D = 0 で、補題 6.3 (2) より ∣ P ( B 1 ) ∣ = ∣ S k ( P ( B ) ) ∣ = ∣ P ( B ) ∣ \lvert P(B_1)\rvert=\lvert S_k(P(B))\rvert=\lvert P(B)\rvert ∣ P ( B 1 )∣ = ∣ S k ( P ( B ))∣ = ∣ P ( B )∣ です。
次に d = D / b 11 ≠ 0 d=D/b_{11}\ne0 d = D / b 11 = 0 、m = − b 12 / d m=-b_{12}/d m = − b 12 / d とすると H m ( b 11 , 0 ) = ( b 11 , 0 ) H_m(b_{11},0)=(b_{11},0) H m ( b 11 , 0 ) = ( b 11 , 0 ) 、H m ( b 12 , d ) = ( b 12 + m d , d ) = ( 0 , d ) H_m(b_{12},d)=(b_{12}+md,\,d)=(0,d) H m ( b 12 , d ) = ( b 12 + m d , d ) = ( 0 , d ) です。補題 6.3 (2) を B 1 B_1 B 1 に使うと ∣ P ( B 2 ) ∣ = ∣ P ( B 1 ) ∣ \lvert P(B_2)\rvert=\lvert P(B_1)\rvert ∣ P ( B 2 )∣ = ∣ P ( B 1 )∣ 、B 2 = ( b 11 0 0 d ) B_2=\begin{pmatrix}b_{11}&0\\0&d\end{pmatrix} B 2 = ( b 11 0 0 d ) です。
P ( B 2 ) = { ( u b 11 , v d ) } P(B_2)=\{(ub_{11},vd)\} P ( B 2 ) = {( u b 11 , v d )} は長方形 [ min ( 0 , b 11 ) , max ( 0 , b 11 ) ] × [ min ( 0 , d ) , max ( 0 , d ) ] [\min(0,b_{11}),\max(0,b_{11})]\times[\min(0,d),\max(0,d)] [ min ( 0 , b 11 ) , max ( 0 , b 11 )] × [ min ( 0 , d ) , max ( 0 , d )] で、φ 1 , φ 2 \varphi_1,\varphi_2 φ 1 , φ 2 を定数とする縦線集合なので、定理 5.6 (f ≡ 1 f\equiv1 f ≡ 1 )と定数の積分(定理 5.4[積分の基本定理と定積分] を F ( t ) = c t F(t)=ct F ( t ) = c t に使う)から ∣ P ( B 2 ) ∣ = ∣ b 11 ∣ ∣ d ∣ = ∣ D ∣ \lvert P(B_2)\rvert=\lvert b_{11}\rvert\,\lvert d\rvert=\lvert D\rvert ∣ P ( B 2 )∣ = ∣ b 11 ∣ ∣ d ∣ = ∣ D ∣ です。以上で ∣ P ( B ) ∣ = ∣ det B ∣ \lvert P(B)\rvert=\lvert\det B\rvert ∣ P ( B )∣ = ∣ det B ∣ です。
∎
命題 6.5 (一次変換による面積の拡大率 ) ✓
A = ( a 11 a 12 a 21 a 22 ) A=\begin{pmatrix} a_{11} & a_{12}\\ a_{21} & a_{22}\end{pmatrix} A = ( a 11 a 21 a 12 a 22 ) を実 2 次正方行列とし、T ( u ) = A u T(\boldsymbol u)=A\boldsymbol u T ( u ) = A u とする。
単位正方形 [ 0 , 1 ] × [ 0 , 1 ] [0,1]\times[0,1] [ 0 , 1 ] × [ 0 , 1 ] の像 T ( [ 0 , 1 ] 2 ) T([0,1]^2) T ([ 0 , 1 ] 2 ) は A A A の 2 つの列ベクトル a 1 = ( a 11 , a 21 ) \boldsymbol a_1=(a_{11},a_{21}) a 1 = ( a 11 , a 21 ) 、a 2 = ( a 12 , a 22 ) \boldsymbol a_2=(a_{12},a_{22}) a 2 = ( a 12 , a 22 ) の張る平行四辺形であり、その面積は ∣ det A ∣ \lvert\det A\rvert ∣ det A ∣ に等しい。
さらに det A ≠ 0 \det A\ne0 det A = 0 とすると、面積確定な有界集合 E ⊂ R 2 E\subset\mathbb{R}^2 E ⊂ R 2 に対し T ( E ) T(E) T ( E ) も面積確定で ∣ T ( E ) ∣ = ∣ det A ∣ ⋅ ∣ E ∣ \lvert T(E)\rvert = \lvert\det A\rvert\cdot\lvert E\rvert ∣ T ( E )∣ = ∣ det A ∣ ⋅ ∣ E ∣ が成り立つ。
証明(命題 6.5) (1) u = ( u , v ) ∈ [ 0 , 1 ] 2 \boldsymbol u=(u,v)\in[0,1]^2 u = ( u , v ) ∈ [ 0 , 1 ] 2 に対し T ( u ) = u a 1 + v a 2 T(\boldsymbol u)=u\boldsymbol a_1+v\boldsymbol a_2 T ( u ) = u a 1 + v a 2 ですから、像は a 1 , a 2 \boldsymbol a_1,\boldsymbol a_2 a 1 , a 2 が張る平行四辺形 P ( A ) P(A) P ( A ) そのものです。補題 6.4 より、これは面積確定で面積は ∣ det A ∣ \lvert\det A\rvert ∣ det A ∣ です。
(2) E ⊂ R = [ a , b ] × [ c , d ] E\subset R=[a,b]\times[c,d] E ⊂ R = [ a , b ] × [ c , d ] となる閉長方形を取り、ε > 0 \varepsilon>0 ε > 0 とします。1 E \mathbf 1_E 1 E は R R R 上可積分なので、命題 3.3 より S ( 1 E , Δ ) − s ( 1 E , Δ ) < ε S(\mathbf 1_E,\Delta)-s(\mathbf 1_E,\Delta)<\varepsilon S ( 1 E , Δ ) − s ( 1 E , Δ ) < ε となる分割 Δ \Delta Δ があり、s ( 1 E , Δ ) ≤ ∣ E ∣ ≤ S ( 1 E , Δ ) s(\mathbf 1_E,\Delta)\le\lvert E\rvert\le S(\mathbf 1_E,\Delta) s ( 1 E , Δ ) ≤ ∣ E ∣ ≤ S ( 1 E , Δ ) です。1 E \mathbf 1_E 1 E の小長方形 R i j R_{ij} R ij 上の sup \sup sup は R i j R_{ij} R ij が E E E と交わるとき 1 1 1 、そうでないとき 0 0 0 、inf \inf inf は R i j ⊂ E R_{ij}\subset E R ij ⊂ E のとき 1 1 1 、そうでないとき 0 0 0 なので
s ( 1 E , Δ ) = ∑ R i j ⊂ E ∣ R i j ∣ , S ( 1 E , Δ ) = ∑ R i j ∩ E ≠ ∅ ∣ R i j ∣ s(\mathbf 1_E,\Delta)=\sum_{R_{ij}\subset E}\lvert R_{ij}\rvert,\qquad S(\mathbf 1_E,\Delta)=\sum_{R_{ij}\cap E\ne\emptyset}\lvert R_{ij}\rvert s ( 1 E , Δ ) = R ij ⊂ E ∑ ∣ R ij ∣ , S ( 1 E , Δ ) = R ij ∩ E = ∅ ∑ ∣ R ij ∣ です。
小長方形の像。 R i j = p + [ 0 , w ] × [ 0 , k ] R_{ij}=\boldsymbol p+[0,w]\times[0,k] R ij = p + [ 0 , w ] × [ 0 , k ] (w = Δ x i w=\Delta x_i w = Δ x i 、k = Δ y j k=\Delta y_j k = Δ y j )と書くと、T T T の線形性から T ( p + ( u w , v k ) ) = T ( p ) + u w a 1 + v k a 2 T(\boldsymbol p+(uw,vk))=T(\boldsymbol p)+u\,w\boldsymbol a_1+v\,k\boldsymbol a_2 T ( p + ( u w , v k )) = T ( p ) + u w a 1 + v k a 2 なので、T ( R i j ) = T ( p ) + P ( A d i a g ( w , k ) ) T(R_{ij})=T(\boldsymbol p)+P\bigl(A\,\mathrm{diag}(w,k)\bigr) T ( R ij ) = T ( p ) + P ( A diag ( w , k ) ) です。補題 6.4 と 補題 6.1 (h = 1 h=1 h = 1 )より T ( R i j ) T(R_{ij}) T ( R ij ) は面積確定で、∣ T ( R i j ) ∣ = ∣ det A ∣ w k = ∣ det A ∣ ∣ R i j ∣ \lvert T(R_{ij})\rvert=\lvert\det A\rvert\,wk=\lvert\det A\rvert\,\lvert R_{ij}\rvert ∣ T ( R ij )∣ = ∣ det A ∣ w k = ∣ det A ∣ ∣ R ij ∣ です。また R i j R_{ij} R ij の周は 4 本の線分で、T T T は線分 [ q , q ′ ] [\boldsymbol q,\boldsymbol q'] [ q , q ′ ] を線分 [ T q , T q ′ ] [T\boldsymbol q,T\boldsymbol q'] [ T q , T q ′ ] に写すので、補題 6.2 と 補題 5.3 (有限個の合併)より T ( ∂ R i j ) T(\partial R_{ij}) T ( ∂ R ij ) は面積 0 0 0 です。det A ≠ 0 \det A\ne0 det A = 0 より T T T は単射(補題 6.3 の証明の (1) と同じ計算で、T u T\boldsymbol u T u から u \boldsymbol u u が決まる)なので T ( i n t R i j ) = T ( R i j ) ∖ T ( ∂ R i j ) T(\mathrm{int}\,R_{ij})=T(R_{ij})\setminus T(\partial R_{ij}) T ( int R ij ) = T ( R ij ) ∖ T ( ∂ R ij ) 、したがって 1 T ( i n t R i j ) = 1 T ( R i j ) − 1 T ( ∂ R i j ) \mathbf 1_{T(\mathrm{int}\,R_{ij})}=\mathbf 1_{T(R_{ij})}-\mathbf 1_{T(\partial R_{ij})} 1 T ( int R ij ) = 1 T ( R ij ) − 1 T ( ∂ R ij ) は可積分で、積分は ∣ det A ∣ ∣ R i j ∣ \lvert\det A\rvert\,\lvert R_{ij}\rvert ∣ det A ∣ ∣ R ij ∣ です(命題 3.6 (1))。
挟み撃ち。 T ( R ) T(R) T ( R ) は平行四辺形の平行移動なので有界です。T ( R ) T(R) T ( R ) を含む閉長方形 R ′ R' R ′ を取り、以下の関数はすべて R ′ R' R ′ 上で考えます(集合の外では 0 0 0 。定義 5.1 )。
g 1 = ∑ R i j ⊂ E 1 T ( i n t R i j ) , f = 1 T ( E ) , g 2 = ∑ R i j ∩ E ≠ ∅ 1 T ( R i j ) g_1=\sum_{R_{ij}\subset E}\mathbf 1_{T(\mathrm{int}\,R_{ij})},\qquad f=\mathbf 1_{T(E)},\qquad g_2=\sum_{R_{ij}\cap E\ne\emptyset}\mathbf 1_{T(R_{ij})} g 1 = R ij ⊂ E ∑ 1 T ( int R ij ) , f = 1 T ( E ) , g 2 = R ij ∩ E = ∅ ∑ 1 T ( R ij ) とおきます。異なる小長方形の内部は交わらず T T T は単射なので、g 1 g_1 g 1 の値は 0 0 0 か 1 1 1 で、g 1 = 1 g_1=1 g 1 = 1 となる点は E E E に含まれる R i j R_{ij} R ij の点の像、つまり T ( E ) T(E) T ( E ) の点です。よって g 1 ≤ f g_1\le f g 1 ≤ f です。T ( E ) T(E) T ( E ) の点は E E E の点 z \boldsymbol z z の像で、z \boldsymbol z z を含む R i j R_{ij} R ij は E E E と交わるので g 2 ≥ 1 g_2\ge1 g 2 ≥ 1 です。よって f ≤ g 2 f\le g_2 f ≤ g 2 です。命題 3.6 (1) と上の計算から、L = ∣ det A ∣ ∣ E ∣ L=\lvert\det A\rvert\,\lvert E\rvert L = ∣ det A ∣ ∣ E ∣ とおくと
∬ R ′ g 1 = ∣ det A ∣ s ( 1 E , Δ ) > L − ∣ det A ∣ ε , ∬ R ′ g 2 = ∣ det A ∣ S ( 1 E , Δ ) < L + ∣ det A ∣ ε \iint_{R'}g_1=\lvert\det A\rvert\,s(\mathbf 1_E,\Delta)>L-\lvert\det A\rvert\varepsilon,\qquad
\iint_{R'}g_2=\lvert\det A\rvert\,S(\mathbf 1_E,\Delta)<L+\lvert\det A\rvert\varepsilon ∬ R ′ g 1 = ∣ det A ∣ s ( 1 E , Δ ) > L − ∣ det A ∣ ε , ∬ R ′ g 2 = ∣ det A ∣ S ( 1 E , Δ ) < L + ∣ det A ∣ ε です(s > S − ε ≥ ∣ E ∣ − ε s>S-\varepsilon\ge\lvert E\rvert-\varepsilon s > S − ε ≥ ∣ E ∣ − ε 、S < s + ε ≤ ∣ E ∣ + ε S<s+\varepsilon\le\lvert E\rvert+\varepsilon S < s + ε ≤ ∣ E ∣ + ε )。定義 3.2 を g 1 , g 2 g_1,g_2 g 1 , g 2 に使うと、ある δ > 0 \delta>0 δ > 0 があって、R ′ R' R ′ の分割 Δ ′ \Delta' Δ ′ で ∣ Δ ′ ∣ < δ \lvert\Delta'\rvert<\delta ∣ Δ ′ ∣ < δ のものとすべての代表点 ξ \xi ξ について Σ ( g 1 , Δ ′ , ξ ) > ∬ g 1 − ε \Sigma(g_1,\Delta',\xi)>\iint g_1-\varepsilon Σ ( g 1 , Δ ′ , ξ ) > ∬ g 1 − ε 、Σ ( g 2 , Δ ′ , ξ ) < ∬ g 2 + ε \Sigma(g_2,\Delta',\xi)<\iint g_2+\varepsilon Σ ( g 2 , Δ ′ , ξ ) < ∬ g 2 + ε です。g 1 ≤ f ≤ g 2 g_1\le f\le g_2 g 1 ≤ f ≤ g 2 よりリーマン和も Σ ( g 1 ) ≤ Σ ( f ) ≤ Σ ( g 2 ) \Sigma(g_1)\le\Sigma(f)\le\Sigma(g_2) Σ ( g 1 ) ≤ Σ ( f ) ≤ Σ ( g 2 ) なので
∣ Σ ( f , Δ ′ , ξ ) − L ∣ < ( ∣ det A ∣ + 1 ) ε \bigl\lvert\Sigma(f,\Delta',\xi)-L\bigr\rvert<(\lvert\det A\rvert+1)\,\varepsilon Σ ( f , Δ ′ , ξ ) − L < (∣ det A ∣ + 1 ) ε です。ε > 0 \varepsilon>0 ε > 0 は任意なので、定義 3.2 により f f f は R ′ R' R ′ 上可積分で積分は L L L 、すなわち T ( E ) T(E) T ( E ) は面積確定で ∣ T ( E ) ∣ = ∣ det A ∣ ∣ E ∣ \lvert T(E)\rvert=\lvert\det A\rvert\,\lvert E\rvert ∣ T ( E )∣ = ∣ det A ∣ ∣ E ∣ です。
∎
定義 6.6 (ヤコビ行列とヤコビアン ) ↗
U ⊂ R 2 U\subset\mathbb{R}^2 U ⊂ R 2 を開集合、Φ : U → R 2 \Phi\colon U\to\mathbb{R}^2 Φ : U → R 2 、Φ ( u , v ) = ( x ( u , v ) , y ( u , v ) ) \Phi(u,v)=\bigl(x(u,v),\,y(u,v)\bigr) Φ ( u , v ) = ( x ( u , v ) , y ( u , v ) ) を C 1 C^1 C 1 級写像とする。
J Φ ( u , v ) = ( ∂ x ∂ u ∂ x ∂ v ∂ y ∂ u ∂ y ∂ v ) J_\Phi(u,v) = \begin{pmatrix}
\dfrac{\partial x}{\partial u} & \dfrac{\partial x}{\partial v}\\[10pt]
\dfrac{\partial y}{\partial u} & \dfrac{\partial y}{\partial v}
\end{pmatrix} J Φ ( u , v ) = ∂ u ∂ x ∂ u ∂ y ∂ v ∂ x ∂ v ∂ y を Φ \Phi Φ のヤコビ行列 、その行列式
det J Φ ( u , v ) = ∂ x ∂ u ∂ y ∂ v − ∂ x ∂ v ∂ y ∂ u ( = ∂ ( x , y ) ∂ ( u , v ) とも書く ) \det J_\Phi(u,v) = \frac{\partial x}{\partial u}\frac{\partial y}{\partial v}-\frac{\partial x}{\partial v}\frac{\partial y}{\partial u}
\ \ \left(=\frac{\partial(x,y)}{\partial(u,v)} \text{ とも書く}\right) det J Φ ( u , v ) = ∂ u ∂ x ∂ v ∂ y − ∂ v ∂ x ∂ u ∂ y ( = ∂ ( u , v ) ∂ ( x , y ) とも書く ) をヤコビアン (関数行列式)という。
ヤコビ行列は Φ \Phi Φ の全微分(1 次近似)を表す行列でした。多変数関数の微分と偏微分 の C¹ 級ならば全微分可能(定理 4.5)[多変数関数の微分と偏微分] で見たとおり、C 1 C^1 C 1 級写像は各点で
Φ ( u + Δ u , v + Δ v ) = Φ ( u , v ) + J Φ ( u , v ) ( Δ u Δ v ) + o ( ∥ ( Δ u , Δ v ) ∥ ) \Phi(u+\Delta u,\ v+\Delta v) = \Phi(u,v) + J_\Phi(u,v)\begin{pmatrix}\Delta u\\ \Delta v\end{pmatrix} + o\bigl(\lVert(\Delta u,\Delta v)\rVert\bigr) Φ ( u + Δ u , v + Δ v ) = Φ ( u , v ) + J Φ ( u , v ) ( Δ u Δ v ) + o ( ∥( Δ u , Δ v )∥ )
と 1 次近似できます。そこで u v uv uv 平面の微小な長方形 [ u , u + Δ u ] × [ v , v + Δ v ] [u,u+\Delta u]\times[v,v+\Delta v] [ u , u + Δ u ] × [ v , v + Δ v ] を Φ \Phi Φ で写すと、誤差項を無視すれば、Φ ( u , v ) \Phi(u,v) Φ ( u , v ) を頂点として 2 つのベクトル
∂ Φ ∂ u Δ u , ∂ Φ ∂ v Δ v \frac{\partial\Phi}{\partial u}\Delta u,\qquad \frac{\partial\Phi}{\partial v}\Delta v ∂ u ∂ Φ Δ u , ∂ v ∂ Φ Δ v
が張る平行四辺形になります。この 2 ベクトルを列とする行列は J Φ d i a g ( Δ u , Δ v ) J_\Phi\,\mathrm{diag}(\Delta u,\Delta v) J Φ diag ( Δ u , Δ v ) ですから、命題 6.5 よりその面積は
∣ det J Φ ( u , v ) ∣ Δ u Δ v \lvert\det J_\Phi(u,v)\rvert\,\Delta u\,\Delta v ∣ det J Φ ( u , v )∣ Δ u Δ v
です。「微小面積 Δ u Δ v \Delta u\Delta v Δ u Δ v がヤコビアンの絶対値倍に化ける」——これが変数変換公式の正体です。
uv 平面 xy 平面 u v x y Φ Δu Δv (u, v) Φu Δu Φv Δv 面積 = Δu Δv 面積 ≒ |det J| Δu Δv ヤコビアンは微小な面積の拡大率。uv 平面の小長方形は、1 次近似の下で偏微分ベクトルが張る平行四辺形に写る
定理 6.7 (重積分の変数変換公式 ) ✓
U , V ⊂ R 2 U,V\subset\mathbb{R}^2 U , V ⊂ R 2 を開集合、Φ : U → V \Phi\colon U\to V Φ : U → V を全単射な C 1 C^1 C 1 級写像で、すべての ( u , v ) ∈ U (u,v)\in U ( u , v ) ∈ U において det J Φ ( u , v ) ≠ 0 \det J_\Phi(u,v)\ne0 det J Φ ( u , v ) = 0 を満たすものとする。E E E を有界閉かつ面積確定な集合で E ⊂ U E\subset U E ⊂ U を満たすものとすると、Φ ( E ) \Phi(E) Φ ( E ) も有界閉かつ面積確定であり、Φ ( E ) \Phi(E) Φ ( E ) 上の任意の連続関数 f f f に対して
∬ Φ ( E ) f ( x , y ) d x d y = ∬ E f ( Φ ( u , v ) ) ∣ det J Φ ( u , v ) ∣ d u d v \iint_{\Phi(E)} f(x,y)\,dx\,dy = \iint_{E} f\bigl(\Phi(u,v)\bigr)\,\bigl\lvert\det J_\Phi(u,v)\bigr\rvert\,du\,dv ∬ Φ ( E ) f ( x , y ) d x d y = ∬ E f ( Φ ( u , v ) ) det J Φ ( u , v ) d u d v が成り立つ。
系 6.9 (面積 0 の例外集合を許す変数変換公式 ) ✓
W ⊂ R 2 W\subset\mathbb{R}^2 W ⊂ R 2 を開集合、Φ : W → R 2 \Phi\colon W\to\mathbb{R}^2 Φ : W → R 2 を C 1 C^1 C 1 級写像、E ⊂ W E\subset W E ⊂ W を有界閉かつ面積確定な集合、Z ⊂ E Z\subset E Z ⊂ E を面積 0 0 0 の有界閉集合とする。Φ \Phi Φ が i n t E ∖ Z \mathrm{int}\,E\setminus Z int E ∖ Z の上で単射で、そこで det J Φ ≠ 0 \det J_\Phi\ne0 det J Φ = 0 ならば、Φ ( E ) \Phi(E) Φ ( E ) は有界閉かつ面積確定であり、Φ ( E ) \Phi(E) Φ ( E ) 上の任意の連続関数 f f f に対して
∬ Φ ( E ) f ( x , y ) d x d y = ∬ E f ( Φ ( u , v ) ) ∣ det J Φ ( u , v ) ∣ d u d v \iint_{\Phi(E)} f(x,y)\,dx\,dy = \iint_{E} f\bigl(\Phi(u,v)\bigr)\,\bigl\lvert\det J_\Phi(u,v)\bigr\rvert\,du\,dv ∬ Φ ( E ) f ( x , y ) d x d y = ∬ E f ( Φ ( u , v ) ) det J Φ ( u , v ) d u d v が成り立つ。
証明は 定理 6.7 の証明のあと、Appendix に置きます。
例 6.11 (極座標変換とそのヤコビアン )
Φ ( r , θ ) = ( r cos θ , r sin θ ) \Phi(r,\theta) = (r\cos\theta,\ r\sin\theta) Φ ( r , θ ) = ( r cos θ , r sin θ ) とします。x = r cos θ x=r\cos\theta x = r cos θ 、y = r sin θ y=r\sin\theta y = r sin θ の偏微分を並べると
J Φ = ( cos θ − r sin θ sin θ r cos θ ) , det J Φ = r cos 2 θ + r sin 2 θ = r . J_\Phi = \begin{pmatrix}\cos\theta & -r\sin\theta\\ \sin\theta & r\cos\theta\end{pmatrix},
\qquad
\det J_\Phi = r\cos^2\theta + r\sin^2\theta = r . J Φ = ( cos θ sin θ − r sin θ r cos θ ) , det J Φ = r cos 2 θ + r sin 2 θ = r . Φ \Phi Φ は開集合 U = { ( r , θ ) ∣ r > 0 , 0 < θ < 2 π } U=\{(r,\theta)\mid r>0,\ 0<\theta<2\pi\} U = {( r , θ ) ∣ r > 0 , 0 < θ < 2 π } の上で単射で det J Φ = r > 0 \det J_\Phi=r>0 det J Φ = r > 0 、像は平面から半直線 { ( x , 0 ) ∣ x ≥ 0 } \{(x,0)\mid x\ge0\} {( x , 0 ) ∣ x ≥ 0 } を除いた開集合です。半径 a a a の円板 D a = { ( x , y ) ∣ x 2 + y 2 ≤ a 2 } D_a=\{(x,y)\mid x^2+y^2\le a^2\} D a = {( x , y ) ∣ x 2 + y 2 ≤ a 2 } に対しては E = [ 0 , a ] × [ 0 , 2 π ] E=[0,a]\times[0,2\pi] E = [ 0 , a ] × [ 0 , 2 π ] と取れば Φ ( E ) = D a \Phi(E)=D_a Φ ( E ) = D a で、単射性が崩れるのは r = 0 r=0 r = 0 の辺と θ = 0 , 2 π \theta=0,2\pi θ = 0 , 2 π の 2 辺だけ、いずれも面積 0 0 0 です。正確には、W = R 2 W=\mathbb{R}^2 W = R 2 、Z = ∂ E Z=\partial E Z = ∂ E (長方形の周、面積 0 0 0 )とすると i n t E ∖ Z = ( 0 , a ) × ( 0 , 2 π ) ⊂ U \mathrm{int}\,E\setminus Z=(0,a)\times(0,2\pi)\subset U int E ∖ Z = ( 0 , a ) × ( 0 , 2 π ) ⊂ U の上で Φ \Phi Φ は単射で det J Φ = r ≠ 0 \det J_\Phi=r\ne0 det J Φ = r = 0 なので、系 6.9 (定理 6.7 の系)が適用でき、連続関数 f f f について
∬ D a f ( x , y ) d x d y = ∫ 0 2 π ∫ 0 a f ( r cos θ , r sin θ ) r d r d θ \iint_{D_a} f(x,y)\,dx\,dy = \int_0^{2\pi}\!\!\int_0^{a} f(r\cos\theta,\ r\sin\theta)\;r\,dr\,d\theta ∬ D a f ( x , y ) d x d y = ∫ 0 2 π ∫ 0 a f ( r cos θ , r sin θ ) r d r d θ が使えます(右辺は長方形上の累次積分なので 系 4.3 で順序も自由です)。
確認 1(円の面積)。 f ≡ 1 f\equiv1 f ≡ 1 とすると
∫ 0 2 π ∫ 0 a r d r d θ = ∫ 0 2 π a 2 2 d θ = π a 2 . \int_0^{2\pi}\!\!\int_0^a r\,dr\,d\theta = \int_0^{2\pi}\frac{a^2}{2}\,d\theta = \pi a^2 . ∫ 0 2 π ∫ 0 a r d r d θ = ∫ 0 2 π 2 a 2 d θ = π a 2 . 確認 2(半球の体積)。 f ( x , y ) = a 2 − x 2 − y 2 f(x,y)=\sqrt{a^2-x^2-y^2} f ( x , y ) = a 2 − x 2 − y 2 とすると f ( Φ ( r , θ ) ) = a 2 − r 2 f(\Phi(r,\theta))=\sqrt{a^2-r^2} f ( Φ ( r , θ )) = a 2 − r 2 で
∫ 0 2 π ∫ 0 a a 2 − r 2 r d r d θ = 2 π [ − ( a 2 − r 2 ) 3 / 2 3 ] 0 a = 2 π ⋅ a 3 3 = 2 3 π a 3 . \int_0^{2\pi}\!\!\int_0^{a}\sqrt{a^2-r^2}\;r\,dr\,d\theta
= 2\pi\left[-\frac{(a^2-r^2)^{3/2}}{3}\right]_0^{a}
= 2\pi\cdot\frac{a^3}{3} = \frac{2}{3}\pi a^3 . ∫ 0 2 π ∫ 0 a a 2 − r 2 r d r d θ = 2 π [ − 3 ( a 2 − r 2 ) 3/2 ] 0 a = 2 π ⋅ 3 a 3 = 3 2 π a 3 . これは半径 a a a の球の体積 4 3 π a 3 \frac43\pi a^3 3 4 π a 3 のちょうど半分で、正しい値です。
因子 r r r を忘れるとどうなるかも見ておきます。∫ 0 2 π ∫ 0 a a 2 − r 2 d r d θ = 2 π ⋅ π a 2 4 \int_0^{2\pi}\int_0^a\sqrt{a^2-r^2}\,dr\,d\theta = 2\pi\cdot\frac{\pi a^2}{4} ∫ 0 2 π ∫ 0 a a 2 − r 2 d r d θ = 2 π ⋅ 4 π a 2 となり、半球の体積とは似ても似つきません。r r r は「外側の θ \theta θ 方向の弧が半径に比例して長くなる」ことを表す、なくてはならない因子です。
例 6.12 (ガウス積分 )
I = ∫ − ∞ ∞ e − x 2 d x I = \int_{-\infty}^{\infty} e^{-x^2}\,dx I = ∫ − ∞ ∞ e − x 2 d x を求めます。まず収束を確認します。∣ x ∣ ≥ 1 \lvert x\rvert\ge1 ∣ x ∣ ≥ 1 では x 2 ≥ ∣ x ∣ x^2\ge \lvert x\rvert x 2 ≥ ∣ x ∣ なので e − x 2 ≤ e − ∣ x ∣ e^{-x^2}\le e^{-\lvert x\rvert} e − x 2 ≤ e − ∣ x ∣ であり、∫ 1 ∞ e − x d x = e − 1 < ∞ \int_1^{\infty}e^{-x}dx = e^{-1}<\infty ∫ 1 ∞ e − x d x = e − 1 < ∞ ですから広義積分は収束します。
I R = ∫ − R R e − x 2 d x I_R=\int_{-R}^{R}e^{-x^2}dx I R = ∫ − R R e − x 2 d x とおきます。e − ( x 2 + y 2 ) = e − x 2 e − y 2 e^{-(x^2+y^2)}=e^{-x^2}e^{-y^2} e − ( x 2 + y 2 ) = e − x 2 e − y 2 は連続なので 系 4.3 が使え、正方形 S R = [ − R , R ] 2 S_R=[-R,R]^2 S R = [ − R , R ] 2 上で
∬ S R e − ( x 2 + y 2 ) d x d y = ∫ − R R ( ∫ − R R e − x 2 e − y 2 d x ) d y = ∫ − R R e − y 2 I R d y = I R 2 . \iint_{S_R} e^{-(x^2+y^2)}dx\,dy
= \int_{-R}^{R}\!\left(\int_{-R}^{R} e^{-x^2}e^{-y^2}\,dx\right)dy
= \int_{-R}^{R} e^{-y^2}I_R\,dy = I_R^{\,2}. ∬ S R e − ( x 2 + y 2 ) d x d y = ∫ − R R ( ∫ − R R e − x 2 e − y 2 d x ) d y = ∫ − R R e − y 2 I R d y = I R 2 . 次に円板で挟みます。D R ⊂ S R ⊂ D 2 R D_R\subset S_R\subset D_{\sqrt2 R} D R ⊂ S R ⊂ D 2 R (正方形の外接円の半径は 2 R \sqrt2R 2 R )であり、被積分関数は正なので、注意 3.7 の加法性と単調性から
∬ D R e − ( x 2 + y 2 ) d x d y ≤ I R 2 ≤ ∬ D 2 R e − ( x 2 + y 2 ) d x d y . \iint_{D_R} e^{-(x^2+y^2)}\,dx\,dy \ \le\ I_R^{\,2}\ \le\ \iint_{D_{\sqrt2 R}} e^{-(x^2+y^2)}\,dx\,dy . ∬ D R e − ( x 2 + y 2 ) d x d y ≤ I R 2 ≤ ∬ D 2 R e − ( x 2 + y 2 ) d x d y . 円板上の積分は 例 6.11 の極座標公式で計算できます。t = r 2 t=r^2 t = r 2 と置換して(d t = 2 r d r dt=2r\,dr d t = 2 r d r )
∬ D a e − ( x 2 + y 2 ) d x d y = ∫ 0 2 π ∫ 0 a e − r 2 r d r d θ = 2 π [ − e − r 2 2 ] 0 a = π ( 1 − e − a 2 ) . \iint_{D_a} e^{-(x^2+y^2)}dx\,dy = \int_0^{2\pi}\!\!\int_0^{a} e^{-r^2}\,r\,dr\,d\theta
= 2\pi\left[-\frac{e^{-r^2}}{2}\right]_0^{a} = \pi\bigl(1-e^{-a^2}\bigr). ∬ D a e − ( x 2 + y 2 ) d x d y = ∫ 0 2 π ∫ 0 a e − r 2 r d r d θ = 2 π [ − 2 e − r 2 ] 0 a = π ( 1 − e − a 2 ) . したがって
π ( 1 − e − R 2 ) ≤ I R 2 ≤ π ( 1 − e − 2 R 2 ) . \pi\bigl(1-e^{-R^2}\bigr)\ \le\ I_R^{\,2}\ \le\ \pi\bigl(1-e^{-2R^2}\bigr). π ( 1 − e − R 2 ) ≤ I R 2 ≤ π ( 1 − e − 2 R 2 ) . R → ∞ R\to\infty R → ∞ とすると両端はともに π \pi π に収束し、I R → I I_R\to I I R → I なので、はさみうちの原理から I 2 = π I^2=\pi I 2 = π 。e − x 2 > 0 e^{-x^2}>0 e − x 2 > 0 より I > 0 I>0 I > 0 ですから
∫ − ∞ ∞ e − x 2 d x = π . \int_{-\infty}^{\infty}e^{-x^2}\,dx = \sqrt{\pi}. ∫ − ∞ ∞ e − x 2 d x = π . 1 変数のままでは e − x 2 e^{-x^2} e − x 2 の原始関数が初等関数で書けないため手が出ませんが、2 乗して 2 次元にし、極座標に移ると r d r r\,dr r d r の r r r が e − r 2 e^{-r^2} e − r 2 の微分の因子をちょうど供給してくれます。ヤコビアンが計算を可能にした、いちばん有名な例です。
演習 7.1 易
∬ [ 0 , 1 ] × [ 0 , 2 ] ( x + y 2 ) d x d y \displaystyle\iint_{[0,1]\times[0,2]}(x+y^2)\,dx\,dy ∬ [ 0 , 1 ] × [ 0 , 2 ] ( x + y 2 ) d x d y を、2 通りの累次積分で計算し、値が一致することを確かめてください。
解答 被積分関数 x + y 2 x+y^2 x + y 2 は長方形 R = [ 0 , 1 ] × [ 0 , 2 ] R=[0,1]\times[0,2] R = [ 0 , 1 ] × [ 0 , 2 ] 上で連続なので、系 4.3 よりどちらの順序でもよく、値は重積分に一致します。
先に y y y で積分すると
∫ 0 1 ( ∫ 0 2 ( x + y 2 ) d y ) d x = ∫ 0 1 [ x y + y 3 3 ] y = 0 y = 2 d x = ∫ 0 1 ( 2 x + 8 3 ) d x = [ x 2 ] 0 1 + 8 3 = 1 + 8 3 = 11 3 . \int_0^1\!\left(\int_0^2 (x+y^2)\,dy\right)dx
= \int_0^1\left[xy+\frac{y^3}{3}\right]_{y=0}^{y=2}dx
= \int_0^1\left(2x+\frac83\right)dx
= \bigl[x^2\bigr]_0^1+\frac83 = 1+\frac83 = \frac{11}{3}. ∫ 0 1 ( ∫ 0 2 ( x + y 2 ) d y ) d x = ∫ 0 1 [ x y + 3 y 3 ] y = 0 y = 2 d x = ∫ 0 1 ( 2 x + 3 8 ) d x = [ x 2 ] 0 1 + 3 8 = 1 + 3 8 = 3 11 . 先に x x x で積分すると
∫ 0 2 ( ∫ 0 1 ( x + y 2 ) d x ) d y = ∫ 0 2 [ x 2 2 + x y 2 ] x = 0 x = 1 d y = ∫ 0 2 ( 1 2 + y 2 ) d y = 1 + 8 3 = 11 3 . \int_0^2\!\left(\int_0^1 (x+y^2)\,dx\right)dy
= \int_0^2\left[\frac{x^2}{2}+xy^2\right]_{x=0}^{x=1}dy
= \int_0^2\left(\frac12+y^2\right)dy
= 1+\frac83 = \frac{11}{3}. ∫ 0 2 ( ∫ 0 1 ( x + y 2 ) d x ) d y = ∫ 0 2 [ 2 x 2 + x y 2 ] x = 0 x = 1 d y = ∫ 0 2 ( 2 1 + y 2 ) d y = 1 + 3 8 = 3 11 . どちらも 11 3 \dfrac{11}{3} 3 11 です。
演習 7.2 標準
∫ 0 1 ( ∫ y 1 sin x x d x ) d y \int_0^1\!\left(\int_y^1 \frac{\sin x}{x}\,dx\right)dy ∫ 0 1 ( ∫ y 1 x sin x d x ) d y の値を求めてください。ただし x = 0 x=0 x = 0 では被積分関数の値を 1 1 1 と定めます(lim x → 0 sin x / x = 1 \lim_{x\to0}\sin x/x=1 lim x → 0 sin x / x = 1 なので、こうすると [ 0 , 1 ] [0,1] [ 0 , 1 ] 上連続になります)。
解答 g ( x ) = sin x x g(x)=\dfrac{\sin x}{x} g ( x ) = x sin x (x ≠ 0 x\ne0 x = 0 )、g ( 0 ) = 1 g(0)=1 g ( 0 ) = 1 とおくと g g g は [ 0 , 1 ] [0,1] [ 0 , 1 ] 上連続です。∫ y 1 g ( x ) d x \int_y^1 g(x)\,dx ∫ y 1 g ( x ) d x の原始関数は初等関数で書けないので、この順序では計算できません。積分順序を交換します。
積分領域は、外側が 0 ≤ y ≤ 1 0\le y\le 1 0 ≤ y ≤ 1 、内側が y ≤ x ≤ 1 y\le x\le 1 y ≤ x ≤ 1 ですから
D = { ( x , y ) ∣ 0 ≤ y ≤ 1 , y ≤ x ≤ 1 } = { ( x , y ) ∣ 0 ≤ y ≤ x ≤ 1 } D = \{(x,y)\mid 0\le y\le 1,\ y\le x\le 1\} = \{(x,y)\mid 0\le y\le x\le 1\} D = {( x , y ) ∣ 0 ≤ y ≤ 1 , y ≤ x ≤ 1 } = {( x , y ) ∣ 0 ≤ y ≤ x ≤ 1 } で、( 0 , 0 ) , ( 1 , 0 ) , ( 1 , 1 ) (0,0),(1,0),(1,1) ( 0 , 0 ) , ( 1 , 0 ) , ( 1 , 1 ) を頂点とする三角形です。これを縦線集合として書き直すと、各 x ∈ [ 0 , 1 ] x\in[0,1] x ∈ [ 0 , 1 ] に対し y y y は 0 0 0 から x x x まで動きます。g ( x ) g(x) g ( x ) は D D D 上連続なので 定理 5.6 を両方の表示に適用でき、
∫ 0 1 ( ∫ y 1 g ( x ) d x ) d y = ∬ D g ( x ) d x d y = ∫ 0 1 ( ∫ 0 x g ( x ) d y ) d x = ∫ 0 1 x g ( x ) d x . \int_0^1\!\left(\int_y^1 g(x)\,dx\right)dy
= \iint_D g(x)\,dx\,dy
= \int_0^1\!\left(\int_0^x g(x)\,dy\right)dx
= \int_0^1 x\,g(x)\,dx . ∫ 0 1 ( ∫ y 1 g ( x ) d x ) d y = ∬ D g ( x ) d x d y = ∫ 0 1 ( ∫ 0 x g ( x ) d y ) d x = ∫ 0 1 x g ( x ) d x . x g ( x ) = sin x x\,g(x)=\sin x x g ( x ) = sin x (x = 0 x=0 x = 0 でも両辺 0 0 0 )なので
∫ 0 1 sin x d x = [ − cos x ] 0 1 = 1 − cos 1. \int_0^1 \sin x\,dx = \bigl[-\cos x\bigr]_0^1 = 1-\cos 1 . ∫ 0 1 sin x d x = [ − cos x ] 0 1 = 1 − cos 1. 数値では 1 − cos 1 = 0.4597 … 1-\cos 1 = 0.4597\ldots 1 − cos 1 = 0.4597 … です。
演習 7.3 標準
a , b , c > 0 a,b,c>0 a , b , c > 0 とする。楕円体
{ ( x , y , z ) | x 2 a 2 + y 2 b 2 + z 2 c 2 ≤ 1 } \left\{(x,y,z)\ \middle|\ \frac{x^2}{a^2}+\frac{y^2}{b^2}+\frac{z^2}{c^2}\le 1\right\} { ( x , y , z ) a 2 x 2 + b 2 y 2 + c 2 z 2 ≤ 1 } の体積が 4 3 π a b c \dfrac43\pi abc 3 4 π ab c であることを、変数変換 Φ ( r , θ ) = ( a r cos θ , b r sin θ ) \Phi(r,\theta)=(ar\cos\theta,\ br\sin\theta) Φ ( r , θ ) = ( a r cos θ , b r sin θ ) を使って示してください。
解答 楕円体は、楕円板 D = { ( x , y ) | x 2 a 2 + y 2 b 2 ≤ 1 } D=\left\{(x,y)\ \middle|\ \frac{x^2}{a^2}+\frac{y^2}{b^2}\le1\right\} D = { ( x , y ) a 2 x 2 + b 2 y 2 ≤ 1 } の上で、グラフ z = h ( x , y ) z=h(x,y) z = h ( x , y ) と z = − h ( x , y ) z=-h(x,y) z = − h ( x , y ) に挟まれた立体です。ここで
h ( x , y ) = c 1 − x 2 a 2 − y 2 b 2 . h(x,y)=c\sqrt{1-\frac{x^2}{a^2}-\frac{y^2}{b^2}} . h ( x , y ) = c 1 − a 2 x 2 − b 2 y 2 . したがって体積は V = ∬ D 2 h ( x , y ) d x d y V=\displaystyle\iint_D 2h(x,y)\,dx\,dy V = ∬ D 2 h ( x , y ) d x d y です。
Φ ( r , θ ) = ( a r cos θ , b r sin θ ) \Phi(r,\theta)=(ar\cos\theta,\ br\sin\theta) Φ ( r , θ ) = ( a r cos θ , b r sin θ ) のヤコビ行列と行列式は
J Φ = ( a cos θ − a r sin θ b sin θ b r cos θ ) , det J Φ = a b r cos 2 θ + a b r sin 2 θ = a b r . J_\Phi=\begin{pmatrix} a\cos\theta & -ar\sin\theta\\ b\sin\theta & br\cos\theta\end{pmatrix},
\qquad
\det J_\Phi = abr\cos^2\theta+abr\sin^2\theta = abr . J Φ = ( a cos θ b sin θ − a r sin θ b r cos θ ) , det J Φ = ab r cos 2 θ + ab r sin 2 θ = ab r . E = [ 0 , 1 ] × [ 0 , 2 π ] E=[0,1]\times[0,2\pi] E = [ 0 , 1 ] × [ 0 , 2 π ] とすると Φ ( E ) = D \Phi(E)=D Φ ( E ) = D で、単射性が崩れるのは r = 0 r=0 r = 0 と θ = 0 , 2 π \theta=0,2\pi θ = 0 , 2 π の辺(面積 0 0 0 )だけです。Z = ∂ E Z=\partial E Z = ∂ E とすると i n t E ∖ Z = ( 0 , 1 ) × ( 0 , 2 π ) \mathrm{int}\,E\setminus Z=(0,1)\times(0,2\pi) int E ∖ Z = ( 0 , 1 ) × ( 0 , 2 π ) 上で Φ \Phi Φ は単射で det J Φ = a b r ≠ 0 \det J_\Phi=abr\ne0 det J Φ = ab r = 0 なので、系 6.9 により変数変換が使えます。また
h ( Φ ( r , θ ) ) = c 1 − r 2 cos 2 θ − r 2 sin 2 θ = c 1 − r 2 h(\Phi(r,\theta)) = c\sqrt{1-r^2\cos^2\theta-r^2\sin^2\theta}=c\sqrt{1-r^2} h ( Φ ( r , θ )) = c 1 − r 2 cos 2 θ − r 2 sin 2 θ = c 1 − r 2 なので、系 6.9 と 系 4.3 より
V = ∫ 0 2 π ∫ 0 1 2 c 1 − r 2 ⋅ a b r d r d θ = 2 a b c ⋅ 2 π ∫ 0 1 r 1 − r 2 d r . V = \int_0^{2\pi}\!\!\int_0^1 2c\sqrt{1-r^2}\cdot abr\,dr\,d\theta
= 2abc\cdot 2\pi\int_0^1 r\sqrt{1-r^2}\,dr . V = ∫ 0 2 π ∫ 0 1 2 c 1 − r 2 ⋅ ab r d r d θ = 2 ab c ⋅ 2 π ∫ 0 1 r 1 − r 2 d r . 最後の 1 変数積分は t = 1 − r 2 t=1-r^2 t = 1 − r 2 (d t = − 2 r d r dt=-2r\,dr d t = − 2 r d r )と置換して
∫ 0 1 r 1 − r 2 d r = [ − ( 1 − r 2 ) 3 / 2 3 ] 0 1 = 1 3 . \int_0^1 r\sqrt{1-r^2}\,dr = \left[-\frac{(1-r^2)^{3/2}}{3}\right]_0^1 = \frac13 . ∫ 0 1 r 1 − r 2 d r = [ − 3 ( 1 − r 2 ) 3/2 ] 0 1 = 3 1 . よって V = 4 π a b c ⋅ 1 3 = 4 3 π a b c V = 4\pi abc\cdot\dfrac13 = \dfrac43\pi abc V = 4 π ab c ⋅ 3 1 = 3 4 π ab c です。a = b = c a=b=c a = b = c とすれば球の体積 4 3 π a 3 \frac43\pi a^3 3 4 π a 3 に一致します。
演習 7.4 難
D D D を ( 0 , 0 ) , ( 1 , 0 ) , ( 0 , 1 ) (0,0),(1,0),(0,1) ( 0 , 0 ) , ( 1 , 0 ) , ( 0 , 1 ) を頂点とする三角形、すなわち D = { ( x , y ) ∣ x ≥ 0 , y ≥ 0 , x + y ≤ 1 } D=\{(x,y)\mid x\ge0,\ y\ge0,\ x+y\le1\} D = {( x , y ) ∣ x ≥ 0 , y ≥ 0 , x + y ≤ 1 } とする。
∬ D exp ( y − x y + x ) d x d y \iint_D \exp\left(\frac{y-x}{y+x}\right)dx\,dy ∬ D exp ( y + x y − x ) d x d y を求めてください(原点では被積分関数が定義されませんが、値を任意に 1 1 1 と定めます)。
解答 下準備。 ( x , y ) ∈ D (x,y)\in D ( x , y ) ∈ D かつ ( x , y ) ≠ ( 0 , 0 ) (x,y)\ne(0,0) ( x , y ) = ( 0 , 0 ) のとき x + y > 0 x+y>0 x + y > 0 で、∣ y − x ∣ ≤ x + y \lvert y-x\rvert\le x+y ∣ y − x ∣ ≤ x + y より指数の値は [ − 1 , 1 ] [-1,1] [ − 1 , 1 ] に入ります。よって被積分関数は [ e − 1 , e ] [e^{-1},e] [ e − 1 , e ] の範囲に収まり有界で、不連続点は原点と ∂ D \partial D ∂ D だけ、すなわち面積 0 0 0 の集合に含まれます。命題 5.4 より D D D 上可積分です。
変数変換。 u = y − x u=y-x u = y − x 、v = y + x v=y+x v = y + x と置きたいので、その逆写像
Φ ( u , v ) = ( v − u 2 , v + u 2 ) \Phi(u,v)=\left(\frac{v-u}{2},\ \frac{v+u}{2}\right) Φ ( u , v ) = ( 2 v − u , 2 v + u ) を考えます。Φ \Phi Φ は R 2 \mathbb{R}^2 R 2 全体の線形同型で
J Φ = ( − 1 2 1 2 1 2 1 2 ) , det J Φ = − 1 4 − 1 4 = − 1 2 , ∣ det J Φ ∣ = 1 2 . J_\Phi = \begin{pmatrix} -\tfrac12 & \tfrac12\\[2pt] \tfrac12 & \tfrac12\end{pmatrix},
\qquad
\det J_\Phi = -\frac14-\frac14 = -\frac12,
\qquad
\lvert\det J_\Phi\rvert = \frac12 . J Φ = ( − 2 1 2 1 2 1 2 1 ) , det J Φ = − 4 1 − 4 1 = − 2 1 , ∣ det J Φ ∣ = 2 1 . Φ ( u , v ) = ( x , y ) \Phi(u,v)=(x,y) Φ ( u , v ) = ( x , y ) とすると、x ≥ 0 ⟺ u ≤ v x\ge0\iff u\le v x ≥ 0 ⟺ u ≤ v 、y ≥ 0 ⟺ u ≥ − v y\ge0\iff u\ge -v y ≥ 0 ⟺ u ≥ − v 、x + y = v ≤ 1 x+y=v\le1 x + y = v ≤ 1 なので
E = { ( u , v ) ∣ 0 ≤ v ≤ 1 , − v ≤ u ≤ v } E=\{(u,v)\mid 0\le v\le 1,\ -v\le u\le v\} E = {( u , v ) ∣ 0 ≤ v ≤ 1 , − v ≤ u ≤ v } が Φ \Phi Φ によって D D D に全単射に写ります(E E E も ( 0 , 0 ) , ( 1 , 1 ) , ( − 1 , 1 ) (0,0),(1,1),(-1,1) ( 0 , 0 ) , ( 1 , 1 ) , ( − 1 , 1 ) を頂点とする三角形です)。
適用。 被積分関数 f f f は原点で連続でないので、η ∈ ( 0 , 1 ) \eta\in(0,1) η ∈ ( 0 , 1 ) を取り E η = { ( u , v ) ∈ E ∣ η ≤ v ≤ 1 } E_\eta=\{(u,v)\in E\mid \eta\le v\le1\} E η = {( u , v ) ∈ E ∣ η ≤ v ≤ 1 } に 定理 6.7 を適用します。E η E_\eta E η 上で f ( Φ ( u , v ) ) = e u / v f(\Phi(u,v))=e^{u/v} f ( Φ ( u , v )) = e u / v は連続です。E η E_\eta E η は縦線集合(v v v を「横軸」と見て、各 v v v に対し u u u が − v -v − v から v v v まで)なので 定理 5.6 が使え、
∬ Φ ( E η ) f d x d y = 1 2 ∫ η 1 ( ∫ − v v e u / v d u ) d v = 1 2 ∫ η 1 v [ e u / v ] u = − v u = v d v = e − e − 1 2 ∫ η 1 v d v . \iint_{\Phi(E_\eta)} f\,dx\,dy
= \frac12\int_\eta^1\!\left(\int_{-v}^{v} e^{u/v}\,du\right)dv
= \frac12\int_\eta^1 v\Bigl[e^{u/v}\Bigr]_{u=-v}^{u=v}dv
= \frac{e-e^{-1}}{2}\int_\eta^1 v\,dv . ∬ Φ ( E η ) f d x d y = 2 1 ∫ η 1 ( ∫ − v v e u / v d u ) d v = 2 1 ∫ η 1 v [ e u / v ] u = − v u = v d v = 2 e − e − 1 ∫ η 1 v d v . (内側で ∫ − v v e u / v d u = v ( e 1 − e − 1 ) \int_{-v}^{v}e^{u/v}du = v\,(e^{1}-e^{-1}) ∫ − v v e u / v d u = v ( e 1 − e − 1 ) を使いました。v > 0 v>0 v > 0 なので u ↦ e u / v u\mapsto e^{u/v} u ↦ e u / v の原始関数は v e u / v v\,e^{u/v} v e u / v です。)
一方 D ∖ Φ ( E η ) D\setminus\Phi(E_\eta) D ∖ Φ ( E η ) は Φ ( E ∩ { v ≤ η } ) \Phi(E\cap\{v\le\eta\}) Φ ( E ∩ { v ≤ η }) で、その面積は ∣ det J Φ ∣ ⋅ η 2 = η 2 / 2 \lvert\det J_\Phi\rvert\cdot\eta^2=\eta^2/2 ∣ det J Φ ∣ ⋅ η 2 = η 2 /2 (命題 6.5 、E ∩ { v ≤ η } E\cap\{v\le\eta\} E ∩ { v ≤ η } は面積 η 2 \eta^2 η 2 の三角形)。∣ f ∣ ≤ e \lvert f\rvert\le e ∣ f ∣ ≤ e なので、この部分の積分の寄与は e η 2 / 2 e\eta^2/2 e η 2 /2 以下で、η → 0 \eta\to0 η → 0 とすれば 0 0 0 に収束します。したがって
∬ D exp ( y − x y + x ) d x d y = e − e − 1 2 ∫ 0 1 v d v = e − e − 1 2 ⋅ 1 2 = e − e − 1 4 . \iint_D \exp\left(\frac{y-x}{y+x}\right)dx\,dy
= \frac{e-e^{-1}}{2}\int_0^1 v\,dv
= \frac{e-e^{-1}}{2}\cdot\frac12
= \frac{e-e^{-1}}{4}. ∬ D exp ( y + x y − x ) d x d y = 2 e − e − 1 ∫ 0 1 v d v = 2 e − e − 1 ⋅ 2 1 = 4 e − e − 1 . 数値では 0.5876 … 0.5876\ldots 0.5876 … です。y − x y-x y − x と y + x y+x y + x という「斜めの座標」に乗り換えることで、絡み合っていた変数が分離した点がこの問題の要点です。
杉浦光夫『解析入門 II』東京大学出版会、1985 — 重積分、フビニの定理、変数変換公式を扱う章。この記事で証明を省いた事実(リーマンの可積分条件、変数変換公式)の完全な証明があります。
高木貞治『解析概論』改訂第三版、岩波書店、1961 — 重積分と面積・体積を扱う章。古典的な記述で、累次積分の扱いが丁寧です。
W. Rudin, Principles of Mathematical Analysis , 3rd ed., McGraw-Hill, 1976 — Chapter 10。primitive mapping への分解による変数変換定理の証明。
M. Spivak, Calculus on Manifolds , W. A. Benjamin, 1965 — Chapter 3 (Integration)。フビニの定理と変数変換定理を、ルベーグ測度 0 0 0 の言葉で簡潔に扱っています。
G. Fubini, “Sugli integrali multipli”, Rendiconti della Reale Accademia dei Lincei 16 (1907) — 累次積分と重積分の関係を一般的に論じた原論文。
T. Tao, An Introduction to Measure Theory , American Mathematical Society, 2011 — ルベーグ積分の枠組みでのフビニ・トネリの定理。非有界関数や無限区間を扱いたいときはこちらへ。
ここまでの議論は次元を上げてもそのまま通用します。R 3 \mathbb{R}^3 R 3 の直方体を小直方体に分割してリーマン和を作れば 3 重積分が定義でき、定理 4.2 と同じ論法で「3 重積分 = 1 変数の積分 3 回」が示せます。変数変換公式も、命題 6.5 の ∣ det A ∣ \lvert\det A\rvert ∣ det A ∣ が「平行六面体の体積」を表すことに置き換えれば、同じ形
∭ Φ ( E ) f d x d y d z = ∭ E f ( Φ ( u , v , w ) ) ∣ det J Φ ∣ d u d v d w \iiint_{\Phi(E)} f\,dx\,dy\,dz = \iiint_{E} f(\Phi(u,v,w))\,\lvert\det J_\Phi\rvert\,du\,dv\,dw ∭ Φ ( E ) f d x d y d z = ∭ E f ( Φ ( u , v , w )) ∣ det J Φ ∣ d u d v d w
で成り立ちます。Appendix の 定理 6.7 と 系 6.9 の証明は、正方形を立方体に、L 2 L^2 L 2 を L 3 L^3 L 3 に、( 1 + η ) 2 (1+\eta)^2 ( 1 + η ) 2 を ( 1 + η ) 3 (1+\eta)^3 ( 1 + η ) 3 に替えればそのまま 3 次元で通ります。2 次元の証明から持ち込めない材料は 命題 6.5 の 3 次元版(一次変換 A A A は体積を ∣ det A ∣ \lvert\det A\rvert ∣ det A ∣ 倍する)だけで、この記事ではそれを証明していません。ルベーグ積分の枠での一般次元の証明は、ℝ^d の変数変換公式と極座標 の 定理 3.4[ℝ^d の変数変換公式と極座標] (線形写像)と 定理 4.5[ℝ^d の変数変換公式と極座標] (C 1 C^1 C 1 微分同相)にあり、例外集合を許す形は 系 5.2[ℝ^d の変数変換公式と極座標] 、球座標とその面積要素は 命題 6.10[ℝ^d の変数変換公式と極座標] です。以下の球座標の計算は、この 3 次元版の系を、E = [ 0 , a ] × [ 0 , π ] × [ 0 , 2 π ] E=[0,a]\times[0,\pi]\times[0,2\pi] E = [ 0 , a ] × [ 0 , π ] × [ 0 , 2 π ] 、例外集合 Z = ∂ E Z=\partial E Z = ∂ E に適用したものです(i n t E \mathrm{int}\,E int E の上では Φ \Phi Φ は単射で、ヤコビアンは下で計算するとおり r 2 sin θ > 0 r^2\sin\theta\gt0 r 2 sin θ > 0 )。
もっともよく使うのが球座標
Φ ( r , θ , φ ) = ( r sin θ cos φ , r sin θ sin φ , r cos θ ) ( r > 0 , 0 < θ < π , 0 < φ < 2 π ) \Phi(r,\theta,\varphi) = (r\sin\theta\cos\varphi,\ r\sin\theta\sin\varphi,\ r\cos\theta)
\qquad (r>0,\ 0<\theta<\pi,\ 0<\varphi<2\pi) Φ ( r , θ , φ ) = ( r sin θ cos φ , r sin θ sin φ , r cos θ ) ( r > 0 , 0 < θ < π , 0 < φ < 2 π )
です(θ \theta θ が天頂角、φ \varphi φ が方位角)。ヤコビ行列は
J Φ = ( sin θ cos φ r cos θ cos φ − r sin θ sin φ sin θ sin φ r cos θ sin φ r sin θ cos φ cos θ − r sin θ 0 ) J_\Phi = \begin{pmatrix}
\sin\theta\cos\varphi & r\cos\theta\cos\varphi & -r\sin\theta\sin\varphi\\
\sin\theta\sin\varphi & r\cos\theta\sin\varphi & r\sin\theta\cos\varphi\\
\cos\theta & -r\sin\theta & 0
\end{pmatrix} J Φ = sin θ cos φ sin θ sin φ cos θ r cos θ cos φ r cos θ sin φ − r sin θ − r sin θ sin φ r sin θ cos φ 0
で、第 3 行に沿って余因子展開します。第 ( 3 , 1 ) (3,1) ( 3 , 1 ) 余因子は
+ det ( r cos θ cos φ − r sin θ sin φ r cos θ sin φ r sin θ cos φ ) = r 2 sin θ cos θ ( cos 2 φ + sin 2 φ ) = r 2 sin θ cos θ , +\det\begin{pmatrix} r\cos\theta\cos\varphi & -r\sin\theta\sin\varphi\\ r\cos\theta\sin\varphi & r\sin\theta\cos\varphi\end{pmatrix}
= r^2\sin\theta\cos\theta(\cos^2\varphi+\sin^2\varphi) = r^2\sin\theta\cos\theta, + det ( r cos θ cos φ r cos θ sin φ − r sin θ sin φ r sin θ cos φ ) = r 2 sin θ cos θ ( cos 2 φ + sin 2 φ ) = r 2 sin θ cos θ ,
第 ( 3 , 2 ) (3,2) ( 3 , 2 ) 余因子は
− det ( sin θ cos φ − r sin θ sin φ sin θ sin φ r sin θ cos φ ) = − r sin 2 θ ( cos 2 φ + sin 2 φ ) = − r sin 2 θ -\det\begin{pmatrix}\sin\theta\cos\varphi & -r\sin\theta\sin\varphi\\ \sin\theta\sin\varphi & r\sin\theta\cos\varphi\end{pmatrix}
= -r\sin^2\theta(\cos^2\varphi+\sin^2\varphi) = -r\sin^2\theta − det ( sin θ cos φ sin θ sin φ − r sin θ sin φ r sin θ cos φ ) = − r sin 2 θ ( cos 2 φ + sin 2 φ ) = − r sin 2 θ
なので
det J Φ = cos θ ⋅ r 2 sin θ cos θ + ( − r sin θ ) ⋅ ( − r sin 2 θ ) = r 2 sin θ ( cos 2 θ + sin 2 θ ) = r 2 sin θ . \det J_\Phi = \cos\theta\cdot r^2\sin\theta\cos\theta + (-r\sin\theta)\cdot(-r\sin^2\theta)
= r^2\sin\theta\,(\cos^2\theta+\sin^2\theta) = r^2\sin\theta . det J Φ = cos θ ⋅ r 2 sin θ cos θ + ( − r sin θ ) ⋅ ( − r sin 2 θ ) = r 2 sin θ ( cos 2 θ + sin 2 θ ) = r 2 sin θ .
0 < θ < π 0<\theta<\pi 0 < θ < π では sin θ > 0 \sin\theta>0 sin θ > 0 なので絶対値も r 2 sin θ r^2\sin\theta r 2 sin θ です。これを使えば半径 a a a の球の体積は
∫ 0 2 π ∫ 0 π ∫ 0 a r 2 sin θ d r d θ d φ = 2 π ⋅ [ − cos θ ] 0 π ⋅ a 3 3 = 2 π ⋅ 2 ⋅ a 3 3 = 4 3 π a 3 \int_0^{2\pi}\!\!\int_0^{\pi}\!\!\int_0^{a} r^2\sin\theta\,dr\,d\theta\,d\varphi
= 2\pi\cdot\bigl[-\cos\theta\bigr]_0^{\pi}\cdot\frac{a^3}{3}
= 2\pi\cdot2\cdot\frac{a^3}{3} = \frac43\pi a^3 ∫ 0 2 π ∫ 0 π ∫ 0 a r 2 sin θ d r d θ d φ = 2 π ⋅ [ − cos θ ] 0 π ⋅ 3 a 3 = 2 π ⋅ 2 ⋅ 3 a 3 = 3 4 π a 3
と一息で求まります。例 6.11 の r r r が「弧長が半径に比例する」ことの反映だったのと同じく、r 2 sin θ r^2\sin\theta r 2 sin θ は「球面の面積要素が r 2 sin θ d θ d φ r^2\sin\theta\,d\theta\,d\varphi r 2 sin θ d θ d φ である」ことの反映です。
証明(定理 6.7) 記号を決めます。x = ( x 1 , x 2 ) \boldsymbol x=(x_1,x_2) x = ( x 1 , x 2 ) に対し ∣ x ∣ ∞ = max ( ∣ x 1 ∣ , ∣ x 2 ∣ ) |\boldsymbol x|_\infty=\max(|x_1|,|x_2|) ∣ x ∣ ∞ = max ( ∣ x 1 ∣ , ∣ x 2 ∣ ) 、2 次正方行列 B = ( b i j ) B=(b_{ij}) B = ( b ij ) に対し ∥ B ∥ = max i ( ∣ b i 1 ∣ + ∣ b i 2 ∣ ) \|B\|=\max_i(|b_{i1}|+|b_{i2}|) ∥ B ∥ = max i ( ∣ b i 1 ∣ + ∣ b i 2 ∣ ) とおくと ∣ B x ∣ ∞ ≤ ∥ B ∥ ∣ x ∣ ∞ |B\boldsymbol x|_\infty\le\|B\|\,|\boldsymbol x|_\infty ∣ B x ∣ ∞ ≤ ∥ B ∥ ∣ x ∣ ∞ です。Q ( a , r ) = { x : ∣ x − a ∣ ∞ ≤ r } Q(\boldsymbol a,r)=\{\boldsymbol x:|\boldsymbol x-\boldsymbol a|_\infty\le r\} Q ( a , r ) = { x : ∣ x − a ∣ ∞ ≤ r } は中心 a \boldsymbol a a 、一辺 2 r 2r 2 r の閉正方形です。積分の線形性・単調性(g ≤ h g\le h g ≤ h なら ∬ g ≤ ∬ h \iint g\le\iint h ∬ g ≤ ∬ h 。命題 3.6 (1)(2))と、平行移動で面積が変わらず h h h 倍の拡大で面積が h 2 h^2 h 2 倍になること(補題 6.1 )を使います。とくに Q ( a , r ) = a + r Q ( 0 , 1 ) Q(\boldsymbol a,r)=\boldsymbol a+r\,Q(\boldsymbol 0,1) Q ( a , r ) = a + r Q ( 0 , 1 ) なので ∣ Q ( a , r ) ∣ = r 2 ∣ Q ( 0 , 1 ) ∣ = 4 r 2 \lvert Q(\boldsymbol a,r)\rvert=r^2\lvert Q(\boldsymbol 0,1)\rvert=4r^2 ∣ Q ( a , r )∣ = r 2 ∣ Q ( 0 , 1 )∣ = 4 r 2 です。
第 0 段:面積 0 0 0 の集合。 有界集合 Z Z Z が面積 0 0 0 (定義 5.1 の意味で面積確定かつ ∣ Z ∣ = 0 |Z|=0 ∣ Z ∣ = 0 )であることは、任意の ε > 0 \varepsilon\gt0 ε > 0 に対し総面積 ε \varepsilon ε 未満の有限個の閉長方形で Z Z Z を覆えることと同値です。実際、前者なら 命題 3.3 により過剰和 S ( 1 Z , Δ ) < ε S(\mathbf 1_Z,\Delta)\lt\varepsilon S ( 1 Z , Δ ) < ε となる分割があり、sup 1 Z = 1 \sup\mathbf 1_Z=1 sup 1 Z = 1 となる小長方形(Z Z Z と交わるもの)が Z Z Z を覆います。後者なら各長方形を面積が 2 倍未満になる範囲で少し広げて内部で Z Z Z を覆い、それらの辺を分点に含む分割をとると、Z Z Z と交わる小長方形はどれかの広げた長方形に含まれるので S ( 1 Z , Δ ) < 2 ε S(\mathbf 1_Z,\Delta)\lt2\varepsilon S ( 1 Z , Δ ) < 2 ε 、s ( 1 Z , Δ ) ≥ 0 s(\mathbf 1_Z,\Delta)\ge0 s ( 1 Z , Δ ) ≥ 0 で、命題 3.3 から ∣ Z ∣ = 0 |Z|=0 ∣ Z ∣ = 0 です。
第 1 段:面積確定と境界。 有界集合 E E E が面積確定 ⟺ \iff ⟺ 境界 ∂ E \partial E ∂ E が面積 0 0 0 です。(⇐ \Leftarrow ⇐ ) 1 E \mathbf 1_E 1 E は ∂ E \partial E ∂ E の外で連続で、∂ E \partial E ∂ E は面積 0 0 0 なので 補題 5.3 によりその条件を満たし、命題 5.4 により可積分です。(⇒ \Rightarrow ⇒ ) E ⊂ R E\subset R E ⊂ R とし、S ( 1 E , Δ ) − s ( 1 E , Δ ) < ε S(\mathbf 1_E,\Delta)-s(\mathbf 1_E,\Delta)\lt\varepsilon S ( 1 E , Δ ) − s ( 1 E , Δ ) < ε となる分割 Δ \Delta Δ をとります。∂ E \partial E ∂ E の点が小長方形 R i j R_{ij} R ij の内部にあれば R i j R_{ij} R ij は E E E の点と E E E の外の点をともに含み、M i j − m i j = 1 M_{ij}-m_{ij}=1 M ij − m ij = 1 です。そうでない点は分割の線分(有限個、面積 0 0 0 )の上にあります。よって ∂ E \partial E ∂ E は総面積 ε \varepsilon ε 未満の長方形と面積 0 0 0 の集合で覆え、第 0 段より面積 0 0 0 です。
第 2 段:Φ \Phi Φ と Φ − 1 \Phi^{-1} Φ − 1 。 y 0 = Φ ( u 0 ) \boldsymbol y_0=\Phi(\boldsymbol u_0) y 0 = Φ ( u 0 ) とすると、逆関数定理(定理 4.1[逆関数定理・陰関数定理と条件付き極値] )により u 0 \boldsymbol u_0 u 0 の開近傍 W W W で Φ ∣ W \Phi|_W Φ ∣ W は開集合 Φ ( W ) \Phi(W) Φ ( W ) への全単射で逆写像は C 1 C^1 C 1 級です。Φ \Phi Φ は単射なので Φ − 1 \Phi^{-1} Φ − 1 は Φ ( W ) \Phi(W) Φ ( W ) 上でこの逆写像に一致します。よって Φ \Phi Φ は開集合を開集合に写し、Φ − 1 : V → U \Phi^{-1}\colon V\to U Φ − 1 : V → U は C 1 C^1 C 1 級で J Φ − 1 ( y ) = J Φ ( Φ − 1 ( y ) ) − 1 J_{\Phi^{-1}}(\boldsymbol y)=J_\Phi(\Phi^{-1}(\boldsymbol y))^{-1} J Φ − 1 ( y ) = J Φ ( Φ − 1 ( y ) ) − 1 、とくに det J Φ − 1 ≠ 0 \det J_{\Phi^{-1}}\ne0 det J Φ − 1 = 0 です。
第 3 段:コンパクト集合の近くでの一様な評価。 K ⊂ U K\subset U K ⊂ U を有界閉集合とします。U c U^c U c は閉で K K K と交わらないので、x ↦ inf z ∉ U ∣ x − z ∣ ∞ \boldsymbol x\mapsto\inf_{\boldsymbol z\notin U}|\boldsymbol x-\boldsymbol z|_\infty x ↦ inf z ∈ / U ∣ x − z ∣ ∞ の K K K 上の最小値は正で、K ′ = { x : inf k ∈ K ∣ x − k ∣ ∞ ≤ 2 r 0 } ⊂ U K'=\{\boldsymbol x:\inf_{\boldsymbol k\in K}|\boldsymbol x-\boldsymbol k|_\infty\le2r_0\}\subset U K ′ = { x : inf k ∈ K ∣ x − k ∣ ∞ ≤ 2 r 0 } ⊂ U となる r 0 > 0 r_0\gt0 r 0 > 0 があります。K ′ K' K ′ は有界閉集合なので、その上で J Φ J_\Phi J Φ の成分と det J Φ \det J_\Phi det J Φ は一様連続で、C = max K ′ ∥ J Φ − 1 ∥ C=\max_{K'}\|J_\Phi^{-1}\| C = max K ′ ∥ J Φ − 1 ∥ と D = max K ′ ∣ det J Φ ∣ D=\max_{K'}|\det J_\Phi| D = max K ′ ∣ det J Φ ∣ は有限です(J Φ − 1 J_\Phi^{-1} J Φ − 1 の成分は連続)。K K K と交わる一辺 2 r ≤ 2 r 0 2r\le2r_0 2 r ≤ 2 r 0 の正方形は K ′ K' K ′ に含まれます。
(i) 面積 0 0 0 の像。 L = max K ′ ∥ J Φ ∥ L=\max_{K'}\|J_\Phi\| L = max K ′ ∥ J Φ ∥ とおきます。Z ⊂ K Z\subset K Z ⊂ K が面積 0 0 0 なら、総面積 ε \varepsilon ε 未満の有限個の長方形 P l P_l P l (辺の長さ h l , k l h_l,k_l h l , k l )で覆います。一辺 2 r 2r 2 r の正方形の格子のうち P l P_l P l と交わるものは辺が h l + 4 r h_l+4r h l + 4 r 、k l + 4 r k_l+4r k l + 4 r の長方形に含まれるので、r ≤ r 0 r\le r_0 r ≤ r 0 を十分小さくとれば、⋃ l P l \bigcup_lP_l ⋃ l P l と交わる格子の正方形の総面積は ∑ l ( h l + 4 r ) ( k l + 4 r ) < 2 ε \sum_l(h_l+4r)(k_l+4r)\lt2\varepsilon ∑ l ( h l + 4 r ) ( k l + 4 r ) < 2 ε です。そのうち Z Z Z と交わるものだけを残すと Z Z Z を覆い、どれも K K K と交わるので K ′ K' K ′ に含まれます。そのような正方形 Q ( a , r ) ⊂ K ′ Q(\boldsymbol a,r)\subset K' Q ( a , r ) ⊂ K ′ の点 x \boldsymbol x x について、線分 [ a , x ] ⊂ Q [\boldsymbol a,\boldsymbol x]\subset Q [ a , x ] ⊂ Q 上で Φ ( x ) − Φ ( a ) = ∫ 0 1 J Φ ( a + t ( x − a ) ) ( x − a ) d t \Phi(\boldsymbol x)-\Phi(\boldsymbol a)=\int_0^1J_\Phi(\boldsymbol a+t(\boldsymbol x-\boldsymbol a))(\boldsymbol x-\boldsymbol a)\,dt Φ ( x ) − Φ ( a ) = ∫ 0 1 J Φ ( a + t ( x − a )) ( x − a ) d t なので ∣ Φ ( x ) − Φ ( a ) ∣ ∞ ≤ L r |\Phi(\boldsymbol x)-\Phi(\boldsymbol a)|_\infty\le Lr ∣Φ ( x ) − Φ ( a ) ∣ ∞ ≤ L r 、すなわち Φ ( Q ) ⊂ Q ( Φ ( a ) , L r ) \Phi(Q)\subset Q(\Phi(\boldsymbol a),Lr) Φ ( Q ) ⊂ Q ( Φ ( a ) , L r ) で、この正方形の面積は 4 L 2 r 2 = L 2 ∣ Q ∣ 4L^2r^2=L^2\lvert Q\rvert 4 L 2 r 2 = L 2 ∣ Q ∣ です(L = 0 L=0 L = 0 なら Φ ( Q ) \Phi(Q) Φ ( Q ) は 1 点で面積 0 0 0 )。よって Φ ( Z ) \Phi(Z) Φ ( Z ) は総面積 2 L 2 ε 2L^2\varepsilon 2 L 2 ε 未満の正方形で覆え、第 0 段より面積 0 0 0 です。
(ii) 小さな正方形の像。 η > 0 \eta\gt0 η > 0 に対し、∣ x − y ∣ ∞ ≤ 2 r 1 |\boldsymbol x-\boldsymbol y|_\infty\le2r_1 ∣ x − y ∣ ∞ ≤ 2 r 1 (x , y ∈ K ′ \boldsymbol x,\boldsymbol y\in K' x , y ∈ K ′ )なら ∥ J Φ ( x ) − J Φ ( y ) ∥ ≤ η / C \|J_\Phi(\boldsymbol x)-J_\Phi(\boldsymbol y)\|\le\eta/C ∥ J Φ ( x ) − J Φ ( y ) ∥ ≤ η / C となる r 1 ≤ r 0 r_1\le r_0 r 1 ≤ r 0 をとります。K K K と交わる Q = Q ( a , r ) Q=Q(\boldsymbol a,r) Q = Q ( a , r ) (r ≤ r 1 r\le r_1 r ≤ r 1 、a ∈ K ′ \boldsymbol a\in K' a ∈ K ′ )と A = J Φ ( a ) A=J_\Phi(\boldsymbol a) A = J Φ ( a ) について、x ∈ Q \boldsymbol x\in Q x ∈ Q なら
A − 1 ( Φ ( x ) − Φ ( a ) ) = ( x − a ) + ∫ 0 1 A − 1 ( J Φ ( a + t ( x − a ) ) − A ) ( x − a ) d t A^{-1}\bigl(\Phi(\boldsymbol x)-\Phi(\boldsymbol a)\bigr)=(\boldsymbol x-\boldsymbol a)+\int_0^1A^{-1}\bigl(J_\Phi(\boldsymbol a+t(\boldsymbol x-\boldsymbol a))-A\bigr)(\boldsymbol x-\boldsymbol a)\,dt A − 1 ( Φ ( x ) − Φ ( a ) ) = ( x − a ) + ∫ 0 1 A − 1 ( J Φ ( a + t ( x − a )) − A ) ( x − a ) d t で、右辺第 2 項の ∣ ⋅ ∣ ∞ |\cdot|_\infty ∣ ⋅ ∣ ∞ は C ⋅ ( η / C ) ⋅ r = η r C\cdot(\eta/C)\cdot r=\eta r C ⋅ ( η / C ) ⋅ r = η r 以下です。よって Φ ( Q ) ⊂ Φ ( a ) + A Q ( 0 , ( 1 + η ) r ) \Phi(Q)\subset\Phi(\boldsymbol a)+A\,Q(\boldsymbol 0,(1+\eta)r) Φ ( Q ) ⊂ Φ ( a ) + A Q ( 0 , ( 1 + η ) r ) です。Φ ( Q ) \Phi(Q) Φ ( Q ) は有界閉集合で、境界は第 2 段により Φ ( ∂ Q ) \Phi(\partial Q) Φ ( ∂ Q ) に含まれ(Φ \Phi Φ は Q Q Q の内部を Φ ( Q ) \Phi(Q) Φ ( Q ) の内部に写す)、(i)(K K K を有界閉集合 Q ⊂ U Q\subset U Q ⊂ U に取り替えて適用。∂ Q \partial Q ∂ Q は 4 本の線分なので面積 0 0 0 )により面積 0 0 0 なので、第 1 段より面積確定です。一方 Q ( 0 , ( 1 + η ) r ) Q(\boldsymbol 0,(1+\eta)r) Q ( 0 , ( 1 + η ) r ) は面積 ( 1 + η ) 2 ∣ Q ∣ (1+\eta)^2\lvert Q\rvert ( 1 + η ) 2 ∣ Q ∣ の正方形(補題 6.1 )で、命題 6.5 の 2 によりその A A A による像の面積は ∣ det A ∣ ( 1 + η ) 2 ∣ Q ∣ \lvert\det A\rvert(1+\eta)^2\lvert Q\rvert ∣ det A ∣ ( 1 + η ) 2 ∣ Q ∣ 、補題 6.1 (平行移動)により Φ ( a ) + A Q ( 0 , ( 1 + η ) r ) \Phi(\boldsymbol a)+A\,Q(\boldsymbol 0,(1+\eta)r) Φ ( a ) + A Q ( 0 , ( 1 + η ) r ) の面積も同じです。面積の単調性(命題 3.6 (2) を指示関数に使う)から
∣ Φ ( Q ) ∣ ≤ ∣ det J Φ ( a ) ∣ ( 1 + η ) 2 ∣ Q ∣ . |\Phi(Q)|\le|\det J_\Phi(\boldsymbol a)|\,(1+\eta)^2|Q| . ∣Φ ( Q ) ∣ ≤ ∣ det J Φ ( a ) ∣ ( 1 + η ) 2 ∣ Q ∣. 第 4 段:Φ ( E ) \Phi(E) Φ ( E ) は面積確定。 Φ ( E ) \Phi(E) Φ ( E ) は有界閉集合(連続写像によるコンパクト集合の像)です。Φ ( i n t E ) \Phi(\mathrm{int}\,E) Φ ( int E ) は開集合で Φ ( E ) \Phi(E) Φ ( E ) に含まれるので Φ ( E ) \Phi(E) Φ ( E ) の内部に含まれ、∂ Φ ( E ) = Φ ( E ) ∖ i n t Φ ( E ) ⊂ Φ ( E ∖ i n t E ) = Φ ( ∂ E ) \partial\Phi(E)=\Phi(E)\setminus\mathrm{int}\,\Phi(E)\subset\Phi(E\setminus\mathrm{int}\,E)=\Phi(\partial E) ∂ Φ ( E ) = Φ ( E ) ∖ int Φ ( E ) ⊂ Φ ( E ∖ int E ) = Φ ( ∂ E ) です。第 1 段より ∂ E \partial E ∂ E は面積 0 0 0 、第 3 段 (i)(K = E K=E K = E )より Φ ( ∂ E ) \Phi(\partial E) Φ ( ∂ E ) も面積 0 0 0 なので、第 1 段より Φ ( E ) \Phi(E) Φ ( E ) は面積確定です。同じ理由で、Φ ( E ) \Phi(E) Φ ( E ) 上の連続関数を 0 0 0 で延長した関数は ∂ Φ ( E ) \partial\Phi(E) ∂ Φ ( E ) の外で連続で、∂ Φ ( E ) \partial\Phi(E) ∂ Φ ( E ) は面積 0 0 0 (上で示しました)なので、両辺の積分は 補題 5.3 と 命題 5.4 により存在します。
第 5 段:f ≥ 0 f\ge0 f ≥ 0 なら左辺 ≤ \le ≤ 右辺。 g = ( f ∘ Φ ) ∣ det J Φ ∣ g=(f\circ\Phi)|\det J_\Phi| g = ( f ∘ Φ ) ∣ det J Φ ∣ を E E E 上の連続関数とし、g ~ \tilde g g ~ を E E E の外で 0 0 0 とした延長とします。K = E K=E K = E として第 3 段を使い、η > 0 \eta\gt0 η > 0 に対する r 1 r_1 r 1 をとります。∣ det J Φ ∣ |\det J_\Phi| ∣ det J Φ ∣ の一様連続性から、r 1 r_1 r 1 を小さくとり直して ∣ x − y ∣ ∞ ≤ 2 r 1 |\boldsymbol x-\boldsymbol y|_\infty\le2r_1 ∣ x − y ∣ ∞ ≤ 2 r 1 なら ∣ ∣ det J Φ ( x ) ∣ − ∣ det J Φ ( y ) ∣ ∣ ≤ η \bigl||\det J_\Phi(\boldsymbol x)|-|\det J_\Phi(\boldsymbol y)|\bigr|\le\eta ∣ det J Φ ( x ) ∣ − ∣ det J Φ ( y ) ∣ ≤ η としておきます。K ′ K' K ′ を含む正方形 R R R を、一辺 2 r 2r 2 r (r ≤ r 1 r\le r_1 r ≤ r 1 )の正方形の格子で分割し、その分割を Δ r \Delta_r Δ r 、E E E と交わる正方形を Q 1 , … , Q N Q_1,\dots,Q_N Q 1 , … , Q N (中心 a j \boldsymbol a_j a j )とします。
Φ ( E ) = ⋃ j Φ ( E ∩ Q j ) \Phi(E)=\bigcup_j\Phi(E\cap Q_j) Φ ( E ) = ⋃ j Φ ( E ∩ Q j ) で、各 Φ ( E ∩ Q j ) \Phi(E\cap Q_j) Φ ( E ∩ Q j ) は第 4 段と同じ理由で面積確定(∂ ( E ∩ Q j ) ⊂ ∂ E ∪ ∂ Q j \partial(E\cap Q_j)\subset\partial E\cup\partial Q_j ∂ ( E ∩ Q j ) ⊂ ∂ E ∪ ∂ Q j )です。f ≥ 0 f\ge0 f ≥ 0 なので 1 Φ ( E ) f ≤ ∑ j 1 Φ ( E ∩ Q j ) f \mathbf 1_{\Phi(E)}f\le\sum_j\mathbf 1_{\Phi(E\cap Q_j)}f 1 Φ ( E ) f ≤ ∑ j 1 Φ ( E ∩ Q j ) f で、積分の単調性と第 3 段 (ii) から
∬ Φ ( E ) f ≤ ∑ j ( max E ∩ Q j f ∘ Φ ) ∣ Φ ( Q j ) ∣ ≤ ( 1 + η ) 2 ∑ j ( max E ∩ Q j f ∘ Φ ) ∣ det J Φ ( a j ) ∣ ∣ Q j ∣ . \iint_{\Phi(E)}f\le\sum_j\Bigl(\max_{E\cap Q_j}f\circ\Phi\Bigr)|\Phi(Q_j)|\le(1+\eta)^2\sum_j\Bigl(\max_{E\cap Q_j}f\circ\Phi\Bigr)|\det J_\Phi(\boldsymbol a_j)|\,|Q_j| . ∬ Φ ( E ) f ≤ j ∑ ( E ∩ Q j max f ∘ Φ ) ∣Φ ( Q j ) ∣ ≤ ( 1 + η ) 2 j ∑ ( E ∩ Q j max f ∘ Φ ) ∣ det J Φ ( a j ) ∣ ∣ Q j ∣. E ∩ Q j E\cap Q_j E ∩ Q j の点 x \boldsymbol x x では ∣ det J Φ ( a j ) ∣ ≤ ∣ det J Φ ( x ) ∣ + η |\det J_\Phi(\boldsymbol a_j)|\le|\det J_\Phi(\boldsymbol x)|+\eta ∣ det J Φ ( a j ) ∣ ≤ ∣ det J Φ ( x ) ∣ + η なので、M f = max Φ ( E ) f M_f=\max_{\Phi(E)}f M f = max Φ ( E ) f とおくと ( max E ∩ Q j f ∘ Φ ) ∣ det J Φ ( a j ) ∣ ≤ sup Q j g ~ + M f η \bigl(\max_{E\cap Q_j}f\circ\Phi\bigr)|\det J_\Phi(\boldsymbol a_j)|\le\sup_{Q_j}\tilde g+M_f\eta ( max E ∩ Q j f ∘ Φ ) ∣ det J Φ ( a j ) ∣ ≤ sup Q j g ~ + M f η です(g ≥ 0 g\ge0 g ≥ 0 なので sup Q j g ~ = max E ∩ Q j g \sup_{Q_j}\tilde g=\max_{E\cap Q_j}g sup Q j g ~ = max E ∩ Q j g )。E E E と交わらない格子の正方形では g ~ = 0 \tilde g=0 g ~ = 0 なので、∑ j ( sup Q j g ~ ) ∣ Q j ∣ \sum_j(\sup_{Q_j}\tilde g)|Q_j| ∑ j ( sup Q j g ~ ) ∣ Q j ∣ は過剰和 S ( g ~ , Δ r ) S(\tilde g,\Delta_r) S ( g ~ , Δ r ) に等しく、
∬ Φ ( E ) f ≤ ( 1 + η ) 2 ( S ( g ~ , Δ r ) + M f η ∣ R ∣ ) . \iint_{\Phi(E)}f\le(1+\eta)^2\bigl(S(\tilde g,\Delta_r)+M_f\eta|R|\bigr). ∬ Φ ( E ) f ≤ ( 1 + η ) 2 ( S ( g ~ , Δ r ) + M f η ∣ R ∣ ) . g ~ \tilde g g ~ は可積分(第 4 段)なので、定義 3.2 と、過剰和が代表点を動かしたリーマン和の上限であること(命題 3.3 の証明の準備)から、r → 0 r\to0 r → 0 で S ( g ~ , Δ r ) → ∬ E g S(\tilde g,\Delta_r)\to\iint_E g S ( g ~ , Δ r ) → ∬ E g です。r → 0 r\to0 r → 0 ののち η → 0 \eta\to0 η → 0 とすれば ∬ Φ ( E ) f ≤ ∬ E g \iint_{\Phi(E)}f\le\iint_Eg ∬ Φ ( E ) f ≤ ∬ E g を得ます。
第 6 段:等号。 第 2 段より Φ − 1 : V → U \Phi^{-1}\colon V\to U Φ − 1 : V → U は定理の仮定を満たし、Φ ( E ) ⊂ V \Phi(E)\subset V Φ ( E ) ⊂ V は有界閉かつ面積確定(第 4 段)です。第 5 段を Φ − 1 \Phi^{-1} Φ − 1 、集合 Φ ( E ) \Phi(E) Φ ( E ) 、関数 g ≥ 0 g\ge0 g ≥ 0 に適用すると、g ( Φ − 1 ( y ) ) ∣ det J Φ − 1 ( y ) ∣ = f ( y ) g(\Phi^{-1}(\boldsymbol y))|\det J_{\Phi^{-1}}(\boldsymbol y)|=f(\boldsymbol y) g ( Φ − 1 ( y )) ∣ det J Φ − 1 ( y ) ∣ = f ( y ) なので
∬ E g ≤ ∬ Φ ( E ) g ∘ Φ − 1 ∣ det J Φ − 1 ∣ = ∬ Φ ( E ) f . \iint_Eg\le\iint_{\Phi(E)}g\circ\Phi^{-1}\,|\det J_{\Phi^{-1}}|=\iint_{\Phi(E)}f . ∬ E g ≤ ∬ Φ ( E ) g ∘ Φ − 1 ∣ det J Φ − 1 ∣ = ∬ Φ ( E ) f . 第 5 段と合わせて、f ≥ 0 f\ge0 f ≥ 0 のとき等号が成り立ちます。一般の連続関数 f f f は f = f + − f − f=f_+-f_- f = f + − f − (f ± = max ( ± f , 0 ) f_\pm=\max(\pm f,0) f ± = max ( ± f , 0 ) は連続で非負)と分けて、積分の線形性から等号を得ます。
∎
証明(系 6.9) 上の証明の記号と段を使います。まず、第 3 段 (i) の議論は Φ \Phi Φ の単射性も det J Φ ≠ 0 \det J_\Phi\ne0 det J Φ = 0 も使わず、有界閉集合 K K K を含む開集合の上で Φ \Phi Φ が C 1 C^1 C 1 級であることだけを使っています。そこで (i) を開集合 W W W と K = E K=E K = E に対して使い、その r 0 r_0 r 0 、K ′ ⊂ W K'\subset W K ′ ⊂ W 、L = max K ′ ∥ J Φ ∥ L=\max_{K'}\|J_\Phi\| L = max K ′ ∥ J Φ ∥ をとります。
定理を当てはめる開集合。 U 0 = i n t E ∖ Z U_0=\mathrm{int}\,E\setminus Z U 0 = int E ∖ Z は開集合で、仮定から Φ ∣ U 0 \Phi|_{U_0} Φ ∣ U 0 は単射、det J Φ ≠ 0 \det J_\Phi\ne0 det J Φ = 0 です。第 2 段の議論(逆関数定理)から V 0 = Φ ( U 0 ) V_0=\Phi(U_0) V 0 = Φ ( U 0 ) は開集合で、Φ ∣ U 0 : U 0 → V 0 \Phi|_{U_0}\colon U_0\to V_0 Φ ∣ U 0 : U 0 → V 0 は 定理 6.7 の仮定を満たします。Z ′ = Z ∪ ∂ E Z'=Z\cup\partial E Z ′ = Z ∪ ∂ E とおくと、E E E は閉なので ∂ E ⊂ E \partial E\subset E ∂ E ⊂ E で、E ∖ Z ′ = U 0 E\setminus Z'=U_0 E ∖ Z ′ = U 0 です。第 1 段より ∂ E \partial E ∂ E は面積 0 0 0 なので、Z ′ Z' Z ′ は面積 0 0 0 の有界閉集合です。
Φ ( E ) \Phi(E) Φ ( E ) は面積確定。 Φ ( E ) \Phi(E) Φ ( E ) は有界閉集合です。開集合 Φ ( U 0 ) \Phi(U_0) Φ ( U 0 ) は Φ ( E ) \Phi(E) Φ ( E ) に含まれるので Φ ( E ) \Phi(E) Φ ( E ) の内部に含まれ、∂ Φ ( E ) ⊂ Φ ( E ∖ U 0 ) = Φ ( Z ′ ) \partial\Phi(E)\subset\Phi(E\setminus U_0)=\Phi(Z') ∂ Φ ( E ) ⊂ Φ ( E ∖ U 0 ) = Φ ( Z ′ ) です。(i) より Φ ( Z ′ ) \Phi(Z') Φ ( Z ′ ) は面積 0 0 0 なので、第 1 段より Φ ( E ) \Phi(E) Φ ( E ) は面積確定で、第 4 段と同じ理由で両辺の積分が存在します。
例外集合を小さな正方形で取り除く。 ε > 0 \varepsilon\gt0 ε > 0 とします。(i) の証明と同じく、Z ′ Z' Z ′ と交わる一辺 2 r 2r 2 r の格子の正方形で総面積 2 ε 2\varepsilon 2 ε 未満のものをとり(r ≤ r 0 / 2 r\le r_0/2 r ≤ r 0 /2 )、各正方形を中心を保って一辺 2 r ′ 2r' 2 r ′ (r < r ′ ≤ r 0 r\lt r'\le r_0 r < r ′ ≤ r 0 、総面積が 4 ε 4\varepsilon 4 ε 未満に収まる範囲)に広げたものを P 1 , … , P N P_1,\dots,P_N P 1 , … , P N とします。各 P l P_l P l は Z ′ ⊂ E Z'\subset E Z ′ ⊂ E と交わり一辺が 2 r 0 2r_0 2 r 0 以下なので K ′ K' K ′ に含まれ、Z ′ Z' Z ′ は内部の和 O = ⋃ l i n t P l O=\bigcup_l\mathrm{int}\,P_l O = ⋃ l int P l に含まれます。E ε = E ∖ O E_\varepsilon=E\setminus O E ε = E ∖ O は有界閉集合で E ε ⊂ E ∖ Z ′ = U 0 E_\varepsilon\subset E\setminus Z'=U_0 E ε ⊂ E ∖ Z ′ = U 0 、境界は ∂ E ∪ ⋃ l ∂ P l \partial E\cup\bigcup_l\partial P_l ∂ E ∪ ⋃ l ∂ P l に含まれて面積 0 0 0 なので、第 1 段より面積確定です。定理 6.7 を Φ ∣ U 0 \Phi|_{U_0} Φ ∣ U 0 と E ε E_\varepsilon E ε に適用すると
∬ Φ ( E ε ) f = ∬ E ε g , g = ( f ∘ Φ ) ∣ det J Φ ∣ . \iint_{\Phi(E_\varepsilon)}f=\iint_{E_\varepsilon}g,\qquad g=(f\circ\Phi)\,|\det J_\Phi| . ∬ Φ ( E ε ) f = ∬ E ε g , g = ( f ∘ Φ ) ∣ det J Φ ∣. 差の評価。 M f = max Φ ( E ) ∣ f ∣ M_f=\max_{\Phi(E)}|f| M f = max Φ ( E ) ∣ f ∣ 、M g = max E ∣ g ∣ M_g=\max_E|g| M g = max E ∣ g ∣ とおきます。E ∖ E ε ⊂ ⋃ l P l E\setminus E_\varepsilon\subset\bigcup_lP_l E ∖ E ε ⊂ ⋃ l P l なので、積分の線形性と単調性から
∣ ∬ E g − ∬ E ε g ∣ = ∣ ∬ 1 E ∖ E ε g ∣ ≤ M g ∑ l ∣ P l ∣ < 4 M g ε . \Bigl|\iint_Eg-\iint_{E_\varepsilon}g\Bigr|=\Bigl|\iint\mathbf 1_{E\setminus E_\varepsilon}\,g\Bigr|\le M_g\sum_l|P_l|\lt4M_g\varepsilon . ∬ E g − ∬ E ε g = ∬ 1 E ∖ E ε g ≤ M g l ∑ ∣ P l ∣ < 4 M g ε . 一方 Φ ( E ) ∖ Φ ( E ε ) ⊂ Φ ( E ∖ E ε ) ⊂ ⋃ l Φ ( P l ) \Phi(E)\setminus\Phi(E_\varepsilon)\subset\Phi(E\setminus E_\varepsilon)\subset\bigcup_l\Phi(P_l) Φ ( E ) ∖ Φ ( E ε ) ⊂ Φ ( E ∖ E ε ) ⊂ ⋃ l Φ ( P l ) で、(i) の証明のとおり Φ ( P l ) \Phi(P_l) Φ ( P l ) は面積 L 2 ∣ P l ∣ L^2|P_l| L 2 ∣ P l ∣ の正方形に含まれるので、同様に
∣ ∬ Φ ( E ) f − ∬ Φ ( E ε ) f ∣ ≤ M f L 2 ∑ l ∣ P l ∣ < 4 M f L 2 ε . \Bigl|\iint_{\Phi(E)}f-\iint_{\Phi(E_\varepsilon)}f\Bigr|\le M_fL^2\sum_l|P_l|\lt4M_fL^2\varepsilon . ∬ Φ ( E ) f − ∬ Φ ( E ε ) f ≤ M f L 2 l ∑ ∣ P l ∣ < 4 M f L 2 ε . 3 つの式を合わせると ∣ ∬ Φ ( E ) f − ∬ E g ∣ < 4 ( M f L 2 + M g ) ε \bigl|\iint_{\Phi(E)}f-\iint_Eg\bigr|\lt4(M_fL^2+M_g)\varepsilon ∬ Φ ( E ) f − ∬ E g < 4 ( M f L 2 + M g ) ε で、ε > 0 \varepsilon\gt0 ε > 0 は任意なので等号が成り立ちます。
∎ 執筆: Claude Code(Opus 5.5)、Claude Code(Sonnet 5.5)・監査 8 回(Claude Code)・最終 2026-10-12