多重積分

概要

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

多重積分の種類

多重積分の種類は、以下の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) $\int _{0}^{1}dx\int _{0}^{\pi /2}dy\,e^{2x}\cos y$
(2) $\int _{0}^{1}dx\int _{0}^{1}dy\,xye^{-x^{2}-y^{2}}$

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

(1)
${\begin{aligned}\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 ydy\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{aligned}}$

(2)
${\begin{aligned}\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]_{0}^{1}\left[-{\frac {e^{-y^{2}}}{2}}\right]_{0}^{1}\\&={\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{aligned}}$

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

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

例えば、被積分関数が2変数関数f(x, y)の場合、先にyを固定してxで積分して、次にyで積分するという手順となる。
$\int _{b_{y}}^{a_{y}}\left\{\int _{b_{x}}^{a_{x}}f(x,y)\,dx\right\}dy$

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

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

(1) $\int _{0}^{\pi /2}dx\int _{0}^{\pi }\sin(x+y)\,dy$
(2) $\int _{\pi }^{2\pi }dx\int _{0}^{1}xy\sin(xy)\,dy$

(1)
${\begin{aligned}\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{aligned}}$

(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で積分する。
構文解析に失敗 (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\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}-\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$とおいて変数変換を行う。
構文解析に失敗 (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle 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'は下図右のようになる。


また、ヤコビアンは、構文解析に失敗 (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\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 )$ に座標変換する。

構文解析に失敗 (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\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}}$

多重積分の極座標変換

軸対称や球対称の関数を積分する際に用いられる極座標による積分である。
$\int _{-\infty }^{+\infty }\int _{-\infty }^{+\infty }{e^{-(x^{2}+y^{2})}\,dxdy}\qquad D=\{(x,y)\,|\,0\leq x^{2}+y^{2}\leq +\infty ,\,-\infty <x<+\infty ,\,-\infty <y<+\infty \}$

積分範囲の中に円 $0\leq x^{2}+y^{2}\leq +\infty$ が含まれている。
このように、円が含まれている時は極座標変換 $x=r\cos(\theta ),\,y=r\sin(\theta ),\quad (0\leq r,\,\,0\leq \theta \leq 2\pi )$ とおく。
（ $r$ は半径のため必ず $0$ 以上、 $\theta$ は最大でも1周分なので $0\leq \theta \leq 2\pi$ または 構文解析に失敗 (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle -\pi \leq \theta \leq \pi } の制約がかかる点に注意する）

ここで変数変換をおこなったのでヤコビアンを計算して $dxdy$ と $drd\theta$ の関係式を求める必要がある。
${\begin{aligned}J=&\left|{\begin{array}{ccc}{\frac {\partial x}{\partial r}}&{\frac {\partial x}{\partial \theta }}\\{\frac {\partial y}{\partial r}}&{\frac {\partial y}{\partial \theta }}\end{array}}\right|\\=&\left|{\begin{array}{ccc}\cos \theta &-r\sin \theta \\\sin \theta &r\cos \theta \end{array}}\right|\\=&r\left(\cos ^{2}\theta +\sin ^{2}\theta \right)\\=&r\end{aligned}}$

したがって、 構文解析に失敗 (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle dxdy=r\ drd\theta } となる。

積分範囲は $0\leq r^{2}<\infty ,\ \ \ -1\leq r\cos(\theta )\leq 1,\ \ \ -1\leq r\sin(\theta )\leq 1$ となり、
$r\geq 0$ のため、 構文解析に失敗 (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle 0\leqq r\leqq \infty ,\ \ \ -1\leq \cos(\theta )\leq 1,\ \ \ -1\leq \sin(\theta )\leq 1} となる。

$\cos(\theta )\geq 0$ かつ $\sin(\theta )\geq 0$ を満たすような $\theta$ は $0\leq \theta \leq 2\pi$ なので、
変換後の積分範囲 構文解析に失敗 (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle D'} は $D'=\{(r,\theta )\,|\,0\leq r\leq \infty ,\,\,0\leq \theta \leq 2\pi \}$ の形に変換でき、2重積分を計算することができる。

構文解析に失敗 (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\begin{aligned}\int _{-\infty }^{+\infty }\int _{-\infty }^{+\infty }{e^{-(x^{2}+y^{2})}\,dxdy}&=\int _{0}^{2\pi }\int _{0}^{\infty }{e^{-(r^{2}\cos ^{2}(\theta )+r^{2}\cos ^{2}(\theta ))}\,rdrd\theta }\\&=\int _{0}^{2\pi }\int _{0}^{\infty }{e^{-r^{2}}\,rdrd\theta }\end{aligned}}}

ここで、更に変数変換を行う。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle t = r^{2}} とおくと、 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle dt = 2r dr} となり、変数 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle t} の積分範囲は、 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle 0 \le t < \infty} となる。

構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{align} \int_{0}^{2 \pi} \int_{0}^{\infty} {e^{- r^{2}} \, r dr d \theta} &= \frac{1}{2} \int_{0}^{2 \pi} \int_{0}^{\infty} {e^{-t} \, dt d \theta} \\ &= \pi \int_{0}^{\infty} {e^{-t} \, dt} \\ &= \pi \Big[ - e^{-t} \Big]_{0}^{\infty} \\ &= \pi \, (0 + 1) \\ &= \pi \end{align} }