解析学
重積分
2次元以上の領域にわたる積分で、体積や高次元の総量を計算する。
直観面積の短冊から体積の角柱へ
一変数の積分 ∫abf(x)dx は区間を幅 Δx の細い区間に分割し、長方形の面積を足し合わせる。二重積分 ∬Df(x,y)dA はまったく同じ発想を二次元領域 D に拡張する。D を面積 ΔA=ΔxΔy の微小な長方形で敷き詰め、各タイルの上に高さ f(xi,yj) の細い角柱を立てて格子を細かくしていく。角柱の体積を足し合わせると、曲面 z=f(x,y) の下の符号付き体積、あるいは面密度 f(x,y) をもつ薄板の全質量が得られる。
平面領域 D 上の曲面 z=f(x,y): 二重積分 ∬Df(x,y)dA は曲面下の鉛直な角柱の体積を累積し、フビニの定理は断面ごとに積分することでそれを計算する。大学厳密な定義:二重積分・三重積分、フビニの累次積分、ヤコビアン
定義: 二重および三重リーマン積分
有界な可測領域 D 上の有界関数 f:D⊆R2→R に対し、二重積分 ∬Df(x,y)dA は分割の細かさを 0 に近づけたときのリーマン和 ∑i,jf(xi∗,yj∗)ΔxiΔyj の極限である。同様に、立体領域 E⊆R3 上の三重積分 ∭Ef(x,y,z)dV は体積 ΔV=ΔxΔyΔz の微小直方体についての和の極限として定義される。
∬Df(x,y)dA=∫ab(∫g1(x)g2(x)f(x,y)dy)dx=∫cd(∫h1(y)h2(y)f(x,y)dx)dy 積分領域が円形、円柱形、あるいは球形であるとき、滑らかな座標変換 Φ(u,v)=(x(u,v),y(u,v)) によって境界の積分範囲が簡単になる。写像 Φ により、微小な長方形 du×dv は引き伸ばされて平行四辺形に移り、その面積はヤコビ行列式の絶対値 ∣detDΦ∣=∂(u,v)∂(x,y) 倍される。極座標 (x,y)=(rcosθ,rsinθ) ではこの因子は r であり、球面座標 (x,y,z)=(ρsinφcosθ,ρsinφsinθ,ρcosφ) では ρ2sinφ となる。
∬Φ(U)f(x,y)dxdy=∬Uf(Φ(u,v))∣detDΦ(u,v)∣dudv,∫−∞∞e−x2dx=π 2次元・3次元の座標系とヤコビアン面積・体積要素| 座標系 | 変換式 | 面積・体積要素 |
|---|
| 2次元直交座標 | (x,y) | dA=dxdy |
| 2次元極座標 | x=rcosθ,y=rsinθ | dA=rdrdθ |
| 3次元円柱座標 | x=rcosθ,y=rsinθ,z=z | dV=rdrdθdz |
| 3次元球面座標 | x=ρsinφcosθ,y=ρsinφsinθ,z=ρcosφ | dV=ρ2sinφdρdφdθ |
大学基本定理:フビニの定理とヤコビアン変数変換公式
f(x,y) が長方形 R=[a,b]×[c,d] 上で連続ならば、二重積分は二つの一変数累次積分のいずれにも等しい: ∬Rf(x,y)dA=∫ab(∫cdf(x,y)dy)dx=∫cd(∫abf(x,y)dx)dy。より一般に、縦線単純領域 D={(x,y):a≤x≤b,g1(x)≤y≤g2(x)} 上では ∬Df(x,y)dA=∫ab∫g1(x)g2(x)f(x,y)dydx が成り立つ。
なぜ正しいのか?
格子上の微小な角柱を足し合わせて体積を求めるとき、先に y 方向の列を足して断面積 A(x) を求めてから A(x) を x で積分しても、先に y 軸に垂直な断面で切ってから積分しても結果は同じになる。
証明
[a,b] を幅 Δx の m 個の小区間 [xi−1,xi] に、[c,d] を幅 Δy の n 個の小区間 [yj−1,yj] に分割する。各 x を固定するごとに断面積分 A(x)=∫cdf(x,y)dy を定める。積分の平均値の定理により、各区間 [yj−1,yj] においてある点 yij∗∈[yj−1,yj] が存在し ∫yj−1yjf(xi,y)dy=f(xi,yij∗)Δy を満たす。
j=1,…,n にわたって足し合わせると A(xi)=∑j=1nf(xi,yij∗)Δy となる。これに Δx を掛けて i=1,…,m について和をとると ∑i=1mA(xi)Δx=∑i=1m∑j=1nf(xi,yij∗)ΔxΔy を得る。f はコンパクトな長方形 R 上で一様連続であるから、Δx,Δy→0 とすると左辺は ∫abA(x)dx=∫ab(∫cdf(x,y)dy)dx に収束し、右辺は ∬Rf(x,y)dA に収束する。x と y を入れ替えて同様の議論を行えば逆順の累次積分との等号も示される。
Φ:U→Φ(U)⊆R2 をヤコビ行列 DΦ(u,v)=(xuyuxvyv) をもつ C1 級微分同相写像とする。Φ(U) 上の任意の可積分関数 f に対して ∬Φ(U)f(x,y)dxdy=∬Uf(x(u,v),y(u,v))∣detDΦ(u,v)∣dudv が成り立つ。ただし detDΦ=xuyv−xvyu である。
なぜ正しいのか?
一変数の置換積分で dx=g′(u)du が長さを伸縮させるのと同様に、∣detDΦ(u,v)∣ は微小な (u,v) 長方形が (x,y) 平面の平行四辺形へ写されるときの局所的な面積拡大率を表す。
証明
U 内の微小な長方形 Rij=[ui,ui+Δu]×[vj,vj+Δv] を考える。(ui,vj) の周りでの一階テイラー展開により、辺 (Δu,0) と (0,Δv) は近似的に接ベクトル a=Φu(ui,vj)Δu=(xu,yu)Δu と b=Φv(ui,vj)Δv=(xv,yv)Δv に写される。R2 において a と b が張る平行四辺形の面積は、高次の誤差 o(ΔuΔv) を除いて ∣xuyv−xvyu∣ΔuΔv=∣detDΦ(ui,vj)∣ΔuΔv に等しい。
ΔAij≈∣detDΦ(ui,vj)∣ΔuΔv をリーマン和 ∑i,jf(Φ(ui,vj))ΔAij に代入すると ∑i,jf(Φ(ui,vj))∣detDΦ(ui,vj)∣ΔuΔv が得られる。特に極座標 x=rcosθ、y=rsinθ では detDΦ=(cosθ)(rcosθ)−(−rsinθ)(sinθ)=r(cos2θ+sin2θ)=r となり、dxdy=rdrdθ が導かれる。
大学実世界での応用と具体例
重積分は物理学、確率論、工学のあらゆる場面に現れる。力学では、三重積分によって全質量 M=∭Eρ(x,y,z)dV、重心、および回転機械の慣性モーメント Iz=∭E(x2+y2)ρdV を計算する。確率・統計学では、正規分布の正規化定数がガウス積分 ∫−∞∞e−x2dx=π に基づいており、これは二乗して極座標の二重積分に直すことで証明される。また天体物理学や電磁気学では、球面座標の体積要素 ρ2sinφdρdφdθ を用いて恒星や惑星の重力場や電荷分布を積分する。
例: 内側の積分に初等的な原始関数がないときの積分順序の交換
積分領域を図示し、フビニの定理によって積分の順序を交換することで、累次積分 I=∫01∫y1ex2dxdy の値を求めよ。
解答
ex2 は初等的な原始関数をもたないため、内側の積分 ∫y1ex2dx をそのまま閉じた形で計算することはできない。しかし、不等式 0≤y≤1 かつ y≤x≤1 は三角形領域 D={(x,y):0≤x≤1,0≤y≤x} を表している。
フビニの定理により、先に縦方向に切ることで I=∫01(∫0xex2dy)dx=∫01xex2dx と書き換えられる。ここで u=x2、du=2xdx と置換すれば I=[21ex2]01=2e−1 が得られる。
例: 極座標を用いたガウス積分の計算
I2 を R2 上の二重積分として表し極座標に変換することで、ガウス積分の公式 I=∫−∞∞e−x2dx=π を証明せよ。
解答
ダミー変数 x と y を用いて I の二つの独立な積を書くと、I2=(∫−∞∞e−x2dx)(∫−∞∞e−y2dy)=∬R2e−(x2+y2)dxdy となる。
極座標 x=rcosθ、y=rsinθ (0≤r<∞、0≤θ≤2π) に変換すると、ヤコビアン因子より dxdy=rdrdθ である。したがって I2=∫02π∫0∞e−r2rdrdθ=2π[−21e−r2]0∞=2π⋅21=π となる。e−x2>0 より正の平方根をとれば I=π が得られる。
長方形領域 R=[0,2]×[0,1] 上の二重積分 ∬R(2x+4y)dA の値を求めよ。
累次積分 ∫02∫x24f(x,y)dydx の積分順序を交換すると、どの累次積分になるか。
3次元球面座標 (x,y,z)=(ρsinφcosθ,ρsinφsinθ,ρcosφ) における体積要素 dV はどれか。
ガウス積分 ∫−∞∞e−x2dx=π を用いると、2次元の確率積分 ∬R2e−(x2+y2)/2dxdy の値はいくらか。