円筒座標（空間極座標）を活用した3重積分の計算

円筒座標を直交座標へ変換する微分同相写像

空間\(\mathbb{R} ^{3}\)上に存在する点\(P\)の直交座標が\(\left( x,y,z\right) \in \mathbb{R} ^{3}\)であり、円筒座標（空間極座標）が\(\left( r,\theta,z\right) \in \mathbb{R} ^{3}\)であるものとします。このとき、以下の関係\begin{equation*}\left(
\begin{array}{c}
x \\
y \\
z\end{array}\right) =\left(
\begin{array}{c}
r\cos \left( \theta \right) \\
r\sin \left( \theta \right) \\
z\end{array}\right)
\end{equation*}が成り立ちます。

直交座標と円座標の間には1対1の関係が成立しません。そこで、両者の間に1対1の関係を成立させるために、円筒座標がとり得る値の範囲を、\begin{equation*}
\left( r,\theta \right) \in \left( 0,+\infty \right) \times \left( 0,2\pi
\right) \times \left( -\infty ,+\infty \right)
\end{equation*}に制限します。その上で、円筒座標を直交座標へ変換するベクトル値関数\begin{equation*}
\boldsymbol{T}:\mathbb{R} ^{3}\subset \left( 0,+\infty \right) \times \left( 0,2\pi \right) \times
\left( -\infty ,+\infty \right) \rightarrow \mathbb{R} ^{3}
\end{equation*}を定義します。つまり、この関数\(\boldsymbol{T}\)はそれぞれの円筒座標\(\left( r,\theta ,z\right) \in \left( 0,+\infty \right) \times \left( 0,2\pi\right) \times \left( -\infty ,+\infty \right) \)に対して、それに対応する直交座標\begin{equation*}\boldsymbol{T}\left( r,\theta ,z\right) =\left(
\begin{array}{c}
T_{1}\left( r,\theta ,z\right) \\
T_{2}\left( r,\theta ,z\right) \\
T_{3}\left( r,\theta ,z\right)
\end{array}\right) =\left(
\begin{array}{c}
r\cos \left( \theta \right) \\
r\sin \left( \theta \right) \\
z\end{array}\right)
\end{equation*}を値として定めます。円筒座標がとり得る値の範囲を制限しているため、直交座標において点が\(xz\)平面上の\(x\geq 0\)となる領域に位置する可能性はなく、直交座標がとり得る値の範囲、すなわち関数\(\boldsymbol{T}\)の値域は、\begin{equation*}R\left( \boldsymbol{T}\right) =\mathbb{R} ^{3}\backslash \left\{ \left( x,0,z\right) \in \mathbb{R} ^{3}\ |\ x\geq 0\right\}
\end{equation*}となります。そこで、関数\(\boldsymbol{T}\)の終集合を値域に制限して、\begin{equation*}\boldsymbol{T}:\mathbb{R} ^{3}\supset \left( 0,+\infty \right) \times \left( 0,2\pi \right) \times
\left( -\infty ,+\infty \right) \rightarrow R\left( \boldsymbol{T}\right)
\end{equation*}とすることにより\(\boldsymbol{T}\)は全単射になるため、逆関数\begin{equation*}\boldsymbol{T}^{-1}:\mathbb{R} ^{3}\supset R\left( \boldsymbol{T}\right) \rightarrow \left( 0,+\infty
\right) \times \left( 0,2\pi \right) \times \left( -\infty ,+\infty \right)
\end{equation*}の存在を保証できます。この逆関数\(\boldsymbol{T}^{-1}\)はそれぞれの直交座標\(\left( x,y,z\right) \in R\left( \boldsymbol{T}\right) \)に対して、それに対応する円筒座標\begin{equation*}\boldsymbol{T}^{-1}\left( x,y,z\right) =\left(
\begin{array}{c}
T_{1}^{-1}\left( x,y,z\right) \\
T_{2}^{-1}\left( x,y,z\right) \\
T_{3}^{-1}\left( x,y,z\right)
\end{array}\right)
\end{equation*}を値として定めます。

ベクトル値関数\(\boldsymbol{T}:\mathbb{R} ^{2}\supset \left( 0,+\infty \right) \times \left( 0,2\pi \right) \times\left( -\infty ,+\infty \right) \rightarrow R\left( \boldsymbol{T}\right) \)は偏微分可能であり、ヤコビ行列関数\(J_{\boldsymbol{T}}:\mathbb{R} ^{3}\supset \left( 0,+\infty \right) \times \left( 0,2\pi \right) \times\left( -\infty ,+\infty \right) \rightarrow R\left( \boldsymbol{T}\right) \)はそれぞれの\(\left( r,\theta ,z\right)\in \left( 0,+\infty \right) \times \left( 0,2\pi \right) \times \left( -\infty ,+\infty \right) \)に対して、\begin{eqnarray*}J_{\boldsymbol{T}}\left( r,\theta ,z\right) &=&\begin{pmatrix}
\frac{\partial T_{1}\left( r,\theta ,z\right) }{\partial r} & \frac{\partial
T_{1}\left( r,\theta ,z\right) }{\partial \theta } & \frac{\partial
T_{1}\left( r,\theta ,z\right) }{\partial z} \\
\frac{\partial T_{2}\left( r,\theta ,z\right) }{\partial r} & \frac{\partial
T_{2}\left( r,\theta ,z\right) }{\partial \theta } & \frac{\partial
T_{2}\left( r,\theta ,z\right) }{\partial z} \\
\frac{\partial T_{3}\left( r,\theta ,z\right) }{\partial r} & \frac{\partial
T_{3}\left( r,\theta ,z\right) }{\partial \theta } & \frac{\partial
T_{3}\left( r,\theta ,z\right) }{\partial z}\end{pmatrix}
\\
&=&\begin{pmatrix}
\frac{\partial r\cos \left( \theta \right) }{\partial r} & \frac{\partial
r\cos \left( \theta \right) }{\partial \theta } & \frac{\partial r\cos
\left( \theta \right) }{\partial z} \\
\frac{\partial r\sin \left( \theta \right) }{\partial r} & \frac{\partial
r\sin \left( \theta \right) }{\partial \theta } & \frac{\partial r\sin
\left( \theta \right) }{\partial z} \\
\frac{\partial z}{\partial r} & \frac{\partial z}{\partial \theta } & \frac{\partial z}{\partial z}\end{pmatrix}
\\
&=&\begin{pmatrix}
\cos \left( \theta \right) & -r\sin \left( \theta \right) & 0 \\
\sin \left( \theta \right) & r\cos \left( \theta \right) & 0 \\
0 & 0 & 1\end{pmatrix}\end{eqnarray*}を定めますが、これは連続関数です。以上より、\(\boldsymbol{T}\)は\(C^{1}\)級であることが明らかになりました。さらに、ヤコビアンは、\begin{eqnarray*}\det \left( J_{\boldsymbol{T}}\left( r,\theta ,z\right) \right)
&=&\left\vert
\begin{pmatrix}
\cos \left( \theta \right) & -r\sin \left( \theta \right) & 0 \\
\sin \left( \theta \right) & r\cos \left( \theta \right) & 0 \\
0 & 0 & 1\end{pmatrix}\right\vert \\
&=&r\cos ^{2}\theta +r\sin ^{2}\theta \\
&=&r \\
&>&0\quad \because r\in \left( 0,+\infty \right)
\end{eqnarray*}を満たします。したがって、逆関数定理より、逆関数\(\boldsymbol{T}^{-1}\)もまた\(C^{1}\)級になることが保証されます。

これまでの議論の結論をまとめます。

円筒座標（空間極座標）のもとでの3重積分

円筒座標を直交座標へ変換するベクトル値関数\begin{equation*}
\boldsymbol{T}:\mathbb{R} ^{3}\supset \left( 0,+\infty \right) \times \left( 0,2\pi \right) \times
\left( -\infty ,+\infty \right) \rightarrow R\left( \boldsymbol{T}\right)
\end{equation*}は\(C^{1}\)級の微分同相写像であることが明らかになりました。\(\boldsymbol{T}\)の定義域の部分集合\begin{equation*}D\subset \left( 0,+\infty \right) \times \left( 0,2\pi \right) \times \left(
-\infty ,+\infty \right)
\end{equation*}を選んだ上で\(\boldsymbol{T}\)の定義域を\(D\)に制限すると、それにあわせて\(\boldsymbol{T}\)の値域は、\begin{eqnarray*}\boldsymbol{T}\left( D\right) &=&\left\{ \boldsymbol{T}\left( r,\theta
,z\right) \in \mathbb{R} ^{3}\ |\ \left( r,\theta ,z\right) \in D\right\} \\
&=&\left\{ \left(
\begin{array}{c}
r\cos \left( \theta \right) \\
r\sin \left( \theta \right) \\
z\end{array}\right) \in \mathbb{R} ^{3}\ |\ \left( r,\theta ,z\right) \in D\right\}
\end{eqnarray*}になります。関数\(\boldsymbol{T}\)の定義域と終集合をともに制限して、\begin{equation*}\boldsymbol{T}:\mathbb{R} ^{3}\supset D\rightarrow \boldsymbol{T}\left( D\right)
\end{equation*}とします。その上で、その値域\(\boldsymbol{T}\left( D\right) \)上に定義された3変数関数の実数値関数\begin{equation*}f:\mathbb{R} ^{3}\supset \boldsymbol{T}\left( D\right) \rightarrow \mathbb{R} \end{equation*}が連続である場合には、多重積分の変数置換公式より、\begin{eqnarray*}
\int \int \int_{\boldsymbol{T}\left( D\right) }f\left( x,y,z\right) dxdydz
&=&\int \int \int_{D}f\left( \boldsymbol{T}\left( r,\theta ,z\right) \right)
\left\vert \det \left( J_{\boldsymbol{T}}\left( r,\theta ,z\right) \right)
\right\vert drd\theta dz \\
&=&\int \int \int_{D}f\left( r\cos \left( \theta \right) ,r\sin \left(
\theta \right) ,z\right) rdrd\theta dz
\end{eqnarray*}を得ます。

これまでは円筒座標の動径\(r\)と偏角\(\theta \)が非ゼロである状況を想定しましたが、\(r\)と\(\theta \)の少なくとも一方が\(0\)をとり得る状況を想定しても、すなわち積分範囲\(D\)に\(xy\)平面上の\(x\geq 0\)となる領域が加わる状況を想定しても3重積分の結果は変わらないため、これまでの議論より以下を得ます。

この命題はどのような意味において有用なのでしょうか。空間\(\mathbb{R} ^{3}\)上の領域\(\boldsymbol{T}\left( D\right) \)に定義された3変数の連続関数\(f\left( x,y,z\right) \)を直交座標のもとで3重積分しようとしている状況、すなわち、\begin{equation}\int \int \int_{\boldsymbol{T}\left( D\right) }f\left( x,y,z\right) dxdydz
\quad \cdots (1)
\end{equation}を特定しようとしている状況を想定します。これを計算するのが困難である場合、まずは、積分範囲と被積分関数を直交座標から円筒座標へそれぞれ変換します。つまり、\begin{eqnarray*}
\boldsymbol{T}\left( D\right) &\rightarrow &D \\
f\left( x,y,z\right) &\rightarrow &f\left( r\cos \left( \theta \right)
,r\sin \left( \theta \right) ,z\right) r
\end{eqnarray*}とするということです。変換後の関数\(f\left(r\cos \left( \theta \right) ,r\sin \left( \theta \right) ,z\right) r\)を変換後の領域\(D\)上において3重積分すれば、\begin{equation}\int \int \int_{D}f\left( r\cos \left( \theta \right) ,r\sin \left( \theta
\right) ,z\right) rdrd\theta dz \quad \cdots (2)
\end{equation}が得られますが、先の命題より、これは\(\left( 1\right) \)と一致することが保証されます。\(\left(1\right) \)を計算するよりも\(\left( 2\right) \)の計算のほうが容易である場合、この手法は有用です。

演習問題