Web Analytics Made Easy - StatCounter

工業大学生ももやまのうさぎ塾

うさぎでもわかるをモットーに大学レベルの数学・情報科目をわかりやすく解説! 数式が読み込まれない場合は1回再読み込みしてみてください。

うさぎでもわかる解析 Part24 変数変換を用いた2重積分の求め方・ヤコビアン

こんにちは、ももやまです。

今回は置換積分の2重積分バージョンに進化したものを紹介していきたいと思います。

 

 

前回の記事(Part23 2重積分の基礎・積分範囲の交換)はこちら!

www.momoyama-usagi.com

(2重積分がうーんな人は上の記事で復習してからこちらの記事を読むことをおすすめします。)

1.ヤコビアンとは

領域  D が\[
D = \{ (x,y) \mid a \leqq x \leqq b, \ c \leqq y \leqq d \ \}
\]の形になっていない場合、2重積分の値を求めることができません。

このような場合、\[
D' = \{ (p,q) \mid a' \leqq p \leqq b', \ c' \leqq q \leqq d' \ \}
\]の形に置き換えることで2重積分の値を求めていきます。

(置換積分の2重積分バージョンだと思ってくれればOK!)

 

具体的に1題練習してみましょう。

例題

つぎの重積分\[
\iint_D (x-2y) e^{x+y} \  dxdy  \\  D = \{ (x,y) \mid 1 \leqq x - 2y \leqq 3, \ 0 \leqq x + y \leqq 1 \ \}
\]を計算しなさい。

解説

範囲の形を変えるため、 x - 2y = p,  x + y = q とおきます。

ここで1変数関数の置換積分を思い出してください。1変数関数のときは、\[
dx = a \ dt
\]のような  dx,  dt の関係式を微分することにより求めていましたね。

2変数関数も同様に  dxdy dpdq の関係式を求めなければなりません。

2変数関数の積分範囲は領域(平面)で表されていましたね。この積分領域の面積の変化率 dxdy dpdq の関係式になります。

変換における面積変化の具合\[
dxdy = a \ dpdq
\]の関係式を求めるために  x,  y p,  q を用いてどう表されるかを計算します。\[
p + 2q = (x-2y) + 2(x+y) = 3x \\ -p + q = -(x-2y)+(x+y) = 3y
\]となるので、\[
\left\{ \begin{array}{l} x = \ \ \ \frac{1}{3} p + \frac{2}{3} q \\ y =- \frac{1}{3} p + \frac{1}{3} q \end{array}\right. 
\]と変換できます。

 

ここで線形代数を思い出してみましょう。

2次の行列式は面積のイメージでしたね。なのでそれぞれの変数を偏微分し、行列式で計算してあげれば  dxdy dpdq の変化率が求められそうですね。

これをヤコビアン、もしくはヤコビ行列式といい、下のような数式で表されます。

 

2重積分における変数変換・ヤコビアン

2重積分の領域を  x = ap + bq ,  y = cp + dq と置き換える場合のヤコビアン  J は\[
J =  \left| \begin{array}{ccc} \frac{\partial x}{\partial p} & \frac{\partial x}{\partial q} \\ \frac{\partial y}{\partial p} & \frac{\partial y}{\partial q} \end{array} \right| = \left| \begin{array}{ccc} a & b \\ c & d \end{array} \right| 
\]の絶対値  |J| となる。このときの変数変換における  dxdy dpdq の関係は\[
dxdy = |J| dpdq
\]となる。(面積の変化率なので絶対値がつくことに要注意!)

ヤコビアンは領域  D' D に戻した場合の面積の変化率を表しています*1。そのため、必ず正となることに注意してください。

 

行列式についてもし分からないことがあればこちらの記事をご覧ください!

www.momoyama-usagi.com

(2次のサラスの公式は覚えていてほしいです…)

 

例題の場合のヤコビアン  J は\[\begin{align*}
J = & \left| \begin{array}{ccc} \frac{\partial x}{\partial p} & \frac{\partial x}{\partial q} \\ \frac{\partial y}{\partial p} & \frac{\partial y}{\partial q} \end{array} \right|
\\ = & \left| \begin{array}{ccc} \frac{1}{3} & \frac{2}{3} \\ - \frac{1}{3} & \frac{1}{3} \end{array} \right| 
\\ = & \ \frac{1}{9} \left| \begin{array}{ccc} 1 & 2 \\ -1 & 1 \end{array} \right| 
\\ = & \ \frac{1}{3}
\end{align*}\]と求められるので、\[
dxdy = \frac{1}{3} \ dpdq, \\
D'= \{ (p,q) \mid 1 \leqq p \leqq 3, \ 0 \leqq q \leqq 1 \ \}
\]となる。

 

変換前の領域  D と変換後の領域  D' を図示すると下のようになります。

f:id:momoyama1192:20190923085139g:plain

※変換前  D と変換後  D' では領域の面積が変わってしまいますね(今回は3倍)、なので拡大率の逆数を積分結果に掛ける*2(今回は1/3倍)することで積分計算結果を変わらないようにしていると頭の中では思ってください。

 

あとは前回の2重積分と同じように計算するだけ。\[\begin{align*} &
\iint_D (x-2y) e^{x+y} \  dxdy 
\\ = & \iint_{D'} \frac{1}{3} u e^v \  dpdq
\\ = & \frac{1}{3} \cdot \int^3_1 p \ dp \cdot \int^1_0 e^q \ dq 
\end{align*} \]を計算すればよい。

 

それぞれの積分の値は\[\begin{align*}
\int^3_1 p \ dp = & \left[ \frac{1}{2} p^2 \right]^3_1
\\ = & \ \frac{1}{2} (3^2 - 1^2)
\\ = & \ 4
\end{align*} \]\[\begin{align*}
\int^1_0 e^q \ dq = & \left[ e^q \right]^1_0
\\ = & \ e - 1
\end{align*} \]となるので、\[\begin{align*} &
\iint_D (x-2y) e^{x+y} \  dxdy 
\\ = & \ \frac{1}{3} \cdot 4 \cdot (e-1)
\\ = & \ \frac{4}{3} (e-1)
\end{align*} \]となる。

 

2.練習問題

では1問だけですが練習してみましょう。

練習1

つぎの重積分\[
\iint_D x^2 \ dxdy  \\  D = \{ (x,y) \mid 0 \leqq x + y \leqq 1, \ 0 \leqq x - y \leqq 1 \ \}
\]を計算しなさい。

 

 

3.練習問題の答え

解答1

 x + y = p,  x - y = q とおく。すると、\[
2x = (x+y)+(x-y) = p + q \\
2y = (x+y)-(x-y) = p - q
\]となるので\[
\left\{ \begin{array}{l} x = \frac{1}{2} p + \frac{1}{2} q \\ y = \frac{1}{2} p - \frac{1}{2} q \end{array}\right. 
\]ヤコビアン  J は\[\begin{align*}
J & = \left| \begin{array}{ccc} \frac{1}{2} & \frac{1}{2} \\ \frac{1}{2} & - \frac{1}{2} \end{array} \right| 
\\ & = \ \frac{1}{4} \left| \begin{array}{ccc} 1 & 1 \\ 1 & -1 \end{array} \right|
\\ & = -\frac{1}{2}
\end{align*} \]の絶対値となる。よって変数変換は\[
dxdy = \frac{1}{2} \ dpdq, \\
D'= \{ (p,q) \mid 0 \leqq p \leqq 1, \ 0 \leqq q \leqq 1 \ \}
\]となる。(絶対値忘れに要注意!)

よって\[\begin{align*} &
\iint_D x^2 \  dxdy 
\\ = & \iint_{D'} \frac{1}{2} (\frac{1}{2} p+ \frac{1}{2} q)^2 \  dpdq
\\ = & \frac{1}{8} \int^1_0 \left( \int^1_0 (p+q)^2 \ dp \right) \ dq 
\end{align*} \]を計算すればよい。

 \[\begin{align*}
\int^1_0 (p+q)^2 \ dp & = \left[ \frac{1}{3} (p+q)^3 \right]^{1}_0
\\ & = \frac{1}{3} (q+1)^3 - \frac{1}{3} q^3
\\ & = q^2 + q + \frac{1}{3}
\end{align*} \]より、\[\begin{align*} &
\frac{1}{8} \int^1_0 \left( \int^1_0 (p+q)^2 \ dp \right) \ dq 
\\ = & \ \frac{1}{8} \int^1_0 q^2 + q + \frac{1}{3} \ dq
\\ = & \ \frac{1}{8} \left[ \frac{1}{3} q^3 + \frac{1}{2} q^2 + \frac{1}{3} q \right]^{1}_0
\\ = & \ \frac{1}{8} \left( \frac{1}{3} + \frac{1}{2} + \frac{1}{3} \right)
\\ = & \ \frac{1}{8} \cdot \frac{7}{6}
\\ = & \ \frac{7}{48}
\end{align*} \]となる。

 

4.さいごに

今回は2重積分における置換積分(変数変換を用いて2重積分を解く)方法についてまとめました。

ヤコビアンは置換後の領域  D' から置換前の領域  D における面積の変化率と頭の片隅にいれておきましょう。

また、この分野は数検1級でも頻出するので数検取りたい人もこの積分法はマスターしましょう!

 

次回は変数変換の中でも特別な極座標を用いた変数変換についてまとめていきたいと思います。

*1:例えば領域  D から領域  D' で面積が2倍になった場合、そのままだと積分結果が2倍になってしまうため、 D' から  D へ戻した場合の面積の変化率(ヤコビアン)を掛けることで積分結果が変わらないようにする。詳しくは図示しているのでそちらで確認してください。

*2:言い換えると変換後の領域  D' から変換前  D に戻す際の面積の変化率を掛ける