多重積分

概要

多重積分とは、関数を多変数で積分することである。

多重積分の種類

多重積分の種類は、以下の4種類に大別される。

• 積分領域が定数のみで決まり、被積分関数が変数分離できる場合

  積分領域が定数のみで決まるということは、$\int _{b}^{a}\int _{d}^{c}f(x,y)dxdy$ のように個々の変数の積分範囲が定数で表されることを指す。

• 積分領域が定数のみで決まり、被積分関数が変数分離できない場合

  被積分関数が変数分離できるということは、$f(x,y)=a(x)b(y)$ のように被積分関数が一変数関数の積で表すことができることを意味する。

• 積分領域が変数に依存し、変数変換する必要がない場合

  積分領域が変数に依存するということは、$\int _{b}^{a}\int _{h(x)}^{g(x)}f(x,y)dxdy$ のように、ある変数の積分範囲が別の変数の関数で表されることを意味する。

• 積分領域が変数に依存し、変数変換する必要がある場合

  変数変換する必要があるということは、与えられた座標系では積分領域が複雑で、変数変換しないと解析的に積分できないことを指す。

積分領域が定数のみで決まり、被積分関数が変数分離できる場合

積分領域が定数のみで決まり、被積分関数が変数分離できる場合、すなわち、
${\begin{aligned}\int _{a_{1}}^{b_{1}}dx_{1}\int _{a_{2}}^{b_{2}}dx_{2}\cdots \int _{a_{n}}^{b_{n}}dx_{n}f(x_{1},x_{2},\cdots ,x_{n})=\int _{a_{1}}^{b_{1}}dx_{1}\int _{a_{2}}^{b_{2}}dx_{2}\cdots \int _{a_{n}}^{b_{n}}dx_{n}\prod _{k=1}^{n}g_{k}(x_{k})\end{aligned}}$ となる場合を考える。

この場合、各1変数関数g_k(x_k)は他の変数の積分には寄与しないため、最終的に次式となり、1変数での積分の積の形に帰着する。
${\begin{aligned}&\int _{a_{1}}^{b_{1}}dx_{1}\int _{a_{2}}^{b_{2}}dx_{2}\cdots \int _{a_{n}}^{b_{n}}dx_{n}\prod _{k=1}^{n}g_{k}(x_{k})\\=&\left(\int _{a_{1}}^{b_{1}}g_{1}(x_{1})dx_{1}\right)\left(\int _{a_{2}}^{b_{2}}g_{2}(x_{2})dx_{2}\right)\cdots \left(\int _{a_{n}}^{b_{n}}g_{n}(x_{n})dx_{n}\right)\\=&\prod _{k=1}^{n}\int _{a_{k}}^{b_{k}}g_{k}(x_{k})dx_{k}\end{aligned}}$

例題. 次の定積分を求めよ。
(1) 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \int_{0}^{1}dx \int_{0}^{\pi/2}dy \, e^{2x} \cos y}

(2) 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \int_{0}^{1}dx \int_{0}^{1}dy \, xye^{-x^{2}-y^{2}}}

例題(1)、(2)ともに積分領域が定数のみで決まり、被積分関数が変数分離できるため、それぞれ独立した1変数での積分の積で表される。

(1)
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{align} \int_{0}^{1} \int_{0}^{\pi / 2} e^{2x} \cos y \, dxdy &= \left ( \int_{0}^{1}e^{2x} dx \right ) \left ( \int_{0}^{\pi/2} \cos y dy \right ) \\ &= \left [ \frac{e^{2x}}{2} \right ]_{0}^{1} \biggl [ \sin y \biggl ]_{0}^{\pi / 2} \\ &= \frac{e^2 - 1}{2} \times 1 \\ &= {\frac{e^{2} - 1}{2}} \end{align} }

(2)
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{align} \int_{0}^{1} \int_{0}^{1} xye^{-x^{2}-y^{2}} \, dxdy &= \left ( \int_{0}^{1} xe^{-x^{2}} dx \right ) \left ( \int_{0}^{1} ye^{-y^{2}} dy \right ) \\ &= \left [ - \frac{e^{-x^{2}}}{2} \right ]^{1}_{0} \left [ - \frac{e^{-y^{2}}}{2} \right ]^{1}_{0} \\ &= \frac{1}{4} \left ( \frac{1}{e} - 1 \right ) \left ( \frac{1}{e} - 1 \right ) \\ &= \frac{1}{4} \left ( \frac{1}{e} - 1 \right )^{2} \end{align} }

積分領域が定数のみで決まり、被積分関数が変数分離できない場合

ある変数で積分する際は他の変数を固定(定数とみなして)して積分して、次の変数で積分する際も同様に他の変数を固定して積分・・・を繰り返す。
積分変数を逐次変えて積分を進めていくことから、このような積分を逐次積分と呼ぶ。

例えば、被積分関数が2変数関数f(x, y)の場合、先にyを固定してxで積分して、次にyで積分するという手順となる。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \int_{b_{y}}^{a_{y}} \left \{ \int_{b_{x}}^{a_{x}} f(x, y) \, dx \right \} dy}

積分変数を選ぶ順番、例えば、上記の場合では、先にxで積分するかyで積分するかは物理数学の範囲では特に気にしなくてよい。
もし、何らかの指定がある場合は、その順序で計算する。

例題. 次の重積分を求めよ。

(1) 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \int_{0}^{\pi / 2} dx \int_{0}^{\pi} \sin (x+y) \, dy}

(2) 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \int_{\pi}^{2 \pi} dx \int_{0}^{1} xy \sin(xy) \, dy}

(1)
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{align} \int_{0}^{\pi / 2} dx \int_{0}^{\pi} \sin(x+y) \, dy &= \int_{0}^{\pi / 2} dx \biggl[- \cos(x + y) \biggr ]_{0}^{\pi} \\ &= \int_{0}^{\pi / 2} \left \{ - \cos(x + \pi) + \cos x \right \} \, dx \\ &= \int_{0}^{\pi / 2} 2 \cos x \, dx \\ &= 2 \biggl [ \sin x \biggr ]_{0}^{\pi / 2} \\ &= 2 \end{align} }

(2)
$f(y)=xy,\quad g(y)=-\cos(xy)$として部分積分を用いる。
${\begin{aligned}\int _{\pi }^{2\pi }dx\int _{0}^{1}xy\sin(xy)\,dy&=\int _{\pi }^{2\pi }dx\left\{{\biggl [}-xy\cos(xy){\biggr ]}_{0}^{1}-\int _{0}^{1}\left\{-x\cos(xy)\right\}dy\right\}\\&=\int _{\pi }^{2\pi }(\sin(x)-x\cos(x))\,dx\\&=\int _{\pi }^{2\pi }\sin(x)\,dx-\int _{\pi }^{2\pi }x\cos(x)\,dx\end{aligned}}$
次に、$f(x)=x,\quad g(x)=\sin(x)$として部分積分を用いる。</math>
${\begin{aligned}&\int _{\pi }^{2\pi }\sin(x)\,dx-\int _{\pi }^{2\pi }x\cos(x)\,dx\\&=\int _{\pi }^{2\pi }\sin(x)dx\,-\left\{{\biggl [}x\sin(x){\biggr ]}_{\pi }^{2\pi }-\int _{\pi }^{2\pi }\sin(x)\,dx\right\}\\&=\int _{\pi }^{2\pi }\sin(x)\,dx-2\\&={\biggl [}-\cos(x){\biggr ]}_{\pi }^{2\pi }-2\\&=-2-2\\&=-4\end{aligned}}$

積分領域が変数に依存し、変数変換する必要がない場合

与えられた座標系でそのまま積分できる(変数変換する必要がない)場合から扱う。

このパターンの最大の特徴は、積分する順番が決まっていることであり、最初に積分領域が変数に依存する積分から実行する。

例えば、2変数関数f(x, y)を$\int _{b}^{a}dx\int _{h(x)}^{g(x)}f(x,y)\,dy$ のように積分する場合、yの積分範囲がxに依存するため、先にyで積分した後にxで積分するという順番をとる。

この順番になる理由は、yでの積分結果がxの関数になるためであり、これを先に計算しないとxでの積分の被積分関数が決定しないからである。
そのため、被積分関数が変数分離できる場合でも、同時並行で各変数での積分を実行することはできない。

積分領域が変数に依存する積分から実行して、他の変数の積分の被積分関数を決定させる必要がある。

例題. 次の重積分を求めよ。

(1) $\int _{0}^{1}\int _{0}^{-x^{2}+1}xe^{-y}\,dxdy$
(2) $\int _{0}^{1}int_{x^{2}}^{x}(x+y)^{2}\,dxdy$

(1)
積分領域は右図のようになる。 被積分関数が変数分離できる形であるが、yの積分がxに依存しているため、まず、yで積分する。
${\begin{aligned}&\int _{0}^{1}\int _{0}^{-x^{2}+1}xe^{-y}\,dxdy\\=&\int _{0}^{1}x\,dx{\biggl [}-e^{-y}{\biggr ]}_{0}^{-x^{2}+1}\\=&\int _{0}^{1}x(-e^{x^{2}-1}+1)\,dx\end{aligned}}$

次に、xで積分する。
${\begin{aligned}\int _{0}^{1}x(-e^{x^{2}-1}+1)\,dx&=\int _{0}^{1}x\,dx-\int _{0}^{1}xe^{x^{2}-1}\,dx\\&=\left[{\frac {1}{2}}x^{2}\right]_{0}^{1}-\left[{\frac {1}{2}}e^{x^{2}-1}\right]_{0}^{1}\\&=\left[{\frac {1}{2}}x^{2}\right]_{0}^{1}-\int _{0}^{1}xe^{x^{2}-1}\,dx\\&={\frac {1}{2}}-\int _{0}^{1}xe^{x^{2}-1}\,dx\end{aligned}}$

ここで、第2項において、$u=x^{2}-1$とおいて変数変換を行う。
$u=x^{2}-1\quad {\mbox{よ り }}\quad {\frac {du}{dx}}=2x\quad {\frac {du}{2}}=xdx$ となる。
また、$x=0\to u=-1,\quad x=1\to u=0$ となる。

上記を代入すると、以下のようになる。
${\begin{aligned}{\frac {1}{2}}-\int _{0}^{1}xe^{x^{2}-1}\,dx&={\frac {1}{2}}-\int _{-1}^{0}{\frac {e^{u}}{2}}\,du\\&={\frac {1}{2}}-{\biggl [}{\frac {e^{u}}{2}}{\biggl ]}_{-1}^{0}\\&={\frac {1}{2}}-\left({\frac {1}{2}}-{\frac {1}{2e}}\right)\\&={\frac {1}{2e}}\end{aligned}}$

(2)
まず、yでの積分がxに依存しているため、yで積分する。
${\begin{aligned}\int _{0}^{1}int_{x^{2}}^{x}(x+y)^{2}\,dxdy&=\int _{0}^{1}dx{\biggl [}{\frac {(x+y)^{3}}{3}}{\biggr ]}_{x^{2}}^{x}\\&={\frac {1}{3}}\int _{0}^{1}\left\{8x^{3}-(x+x^{2})^{3}\right\}\,dx\\&={\frac {1}{3}}\int _{0}^{1}(-x^{6}-3x^{5}-3x^{4}+7x^{3})\,dx\end{aligned}}$

次に、xで積分する。
${\begin{aligned}{\frac {1}{3}}\int _{0}^{1}(-x^{6}-3x^{5}-3x^{4}+7x^{3})\,dx&={\frac {1}{3}}\left[-{\frac {1}{7}}x^{7}-{\frac {1}{2}}x^{6}-{\frac {3}{5}}x^{5}+{\frac {7}{4}}x^{4}\right]_{0}^{1}\\&={\frac {1}{3}}\left(-{\frac {1}{7}}-{\frac {1}{2}}-{\frac {3}{5}}+{\frac {7}{4}}\right)\\&={\frac {71}{420}}\end{aligned}}$

積分領域が変数に依存し、変数変換する必要がある場合

上記までのパターンでは、積分領域が比較的単純だった。
しかし、積分領域が複雑になる場合、計算が煩雑になったり、そもそも与えられた座標系では解析的に積分できない場合がある。

その場合には、座標変換によって積分変数の変換を行い、積分領域を単純化した上で積分計算を進めていく。

積分変数の変換は、次のように行う。
ある多変数関数f(x1, ..., x_n)について、領域Dにわたって積分することを考える。
$\int _{D}f(x_{1},\cdots ,x_{n})\,dx_{1}\cdots dx_{n}$
ここで、$(x_{1},\cdots ,x_{n})\to (\alpha _{1},\cdots ,\alpha _{n})$の変数変換を行う場合、変数変換後の領域をD'とする時、
上記の積分は、次式のようになる。
${\begin{aligned}&\int _{D}f(x_{1},\cdots ,x_{n})\,dx_{1}\cdots dx_{n}\\=&\int _{D'}\left|{\frac {\partial (x_{1},\cdots ,x_{n})}{\partial (\alpha _{1},\cdots ,\alpha _{n})}}\right|d\alpha _{1}\cdots d\alpha _{n}\quad f(x_{1}(\alpha _{1},\cdots ,\alpha _{n}),\cdots ,x_{n}(\alpha _{1},\cdots ,\alpha _{n}))\end{aligned}}$

ここで、$\left|{\frac {\partial (x_{1},\cdots ,x_{n})}{\partial (\alpha _{1},\cdots ,\alpha _{n})}}\right|$ は、ヤコビアンの絶対値である。
すなわち、積分変数を変換したのち、被積分関数にヤコビアンの絶対値を乗算すればよい。

例題. 次の重積分を求めよ。

(1) $f(x,y)={\frac {a}{(x^{2}+y^{2}+a^{2})^{\frac {3}{2}}}},\qquad D:R_{1}^{2}\leq x^{2}+y^{2}\leq R_{2}^{2}$
(2) $f(x,y,z)=2z^{3}(2x^{2}-y^{2}),\qquad D:x^{2}+y^{2}+z^{2}\leq a^{2},\,z\geq 0$

(1)
与えられた座標系では、積分領域Dは下図左のようにドーナツ型になる。

まず、極座標系(r, θ)に座標変換する。
${\begin{aligned}{\begin{cases}x=r\cos \theta \\y=r\sin \theta \end{cases}}\qquad (R_{1}\leq r\leq R_{2},\,0\leq \theta \leq 2\pi )\end{aligned}}$

この時、新たな積分領域D'は下図右のようになる。


また、ヤコビアンは、${\frac {\partial (x,y)}{\partial (r,\theta )}}=r$ となるため、求める多重積分は次式となる。
${\begin{aligned}&\iint _{D}{\frac {a}{(x^{2}+y^{2}+a^{2})^{\frac {3}{2}}}}\,dxdy\\&=\iint _{D'}{\frac {a}{(r^{2}\cos ^{2}\theta +r^{2}\sin ^{2}\theta +a^{2})^{\frac {3}{2}}}}\left|{\frac {\partial (x,y)}{\partial (r,\theta )}}\right|\,drd\theta \\&=\int _{0}^{2\pi }\,d\theta \int _{R_{1}}^{R_{2}}{\frac {ar}{(r^{2}+a^{2})^{3/2}}}\,dr\\&=2\pi a\left[-{\frac {1}{(r^{2}+a^{2})^{\frac {1}{2}}}}\right]_{R_{1}}^{R_{2}}\\&=2\pi a\left({\frac {1}{\sqrt {R_{1}^{2}-a^{2}}}}-{\frac {1}{\sqrt {R_{2}^{2}-a^{2}}}}\right)\end{aligned}}$

(2)
与えられた座標系では、積分領域Dは下図のように半球型になる。
そのため、極座標系$(r,\theta ,\phi )$ に座標変換する。

${\begin{aligned}{\begin{cases}x=r\sin \theta \cos \phi \\y=r\sin \theta \sin \phi \\z=r\cos \theta \end{cases}}\qquad \left(0\leq r\leq a,\,0\leq \theta \leq {\frac {\pi }{2}},\,0\leq \phi \leq 2\pi \right)\end{aligned}}$

また、ヤコビアンは、${\frac {\partial (x,y,z)}{\partial (r,\theta ,\phi )}}=r^{2}\sin \theta$ となる。

したがって、求める多重積分は次式となる。
${\begin{aligned}&\iint _{D}2z^{3}(2x^{2}-y^{2})\,dxdydz\\&=\iint _{D'}2r^{3}\cos ^{3}\theta (2r^{2}\sin ^{2}\theta \cos ^{2}\phi -r^{2}\sin ^{2}\theta \sin ^{2}\phi )\,\left|{\frac {\partial (x,y,z)}{\partial (r,\theta ,\phi )}}\right|drd\theta d\phi \\&=2\int _{0}^{a}dr\int _{0}^{\pi /2}d\theta \int _{0}^{2\pi }r^{7}\sin ^{3}\theta \cos ^{3}\theta (2\cos ^{2}\phi -\sin ^{2}\phi )\,d\phi \\&=2\left(\int _{0}^{a}r^{7}dr\,\right)\left(\int _{0}^{\pi /2}\sin ^{3}\theta \cos ^{3}\theta \,d\theta \right)\left\{\int _{0}^{2\pi }(2\cos ^{2}\phi -\sin ^{2}\phi )\,d\phi \right\}\\&=\left(\int _{0}^{a}r^{7}\,dr\right)\left\{\int _{0}^{\pi /2}\sin ^{3}\theta (1-\sin ^{2}\theta )\cos \theta \,d\theta \right\}\left\{\int _{0}^{2\pi }(3\cos ^{2}\phi -1)\,d\phi \right\}\\&=\left(\int _{0}^{a}r^{7}\,dr\right)\left\{\int _{0}^{\pi /2}\cos \theta (\sin ^{3}\theta -\sin ^{5}\theta )\,d\theta \right\}\left\{\int _{0}^{2\pi }\left({\frac {3}{2}}\cos 2\phi +{\frac {1}{2}}\right)\,d\phi \right\}\\&=2\left[{\frac {r^{8}}{8}}\right]_{0}^{a}\left[{\frac {\sin ^{4}\theta }{4}}-{\frac {\sin ^{6}\theta }{6}}\right]_{0}^{\pi /2}{\frac {1}{2}}\left[{\frac {3\sin 2\phi }{2}}+\phi \right]_{0}^{2\pi }\\&={\frac {a^{8}}{8}}\left({\frac {1}{4}}-{\frac {1}{6}}\right)2\pi \\&={\frac {\pi a^{8}}{48}}\end{aligned}}$