倭算数理研究所

科学・数学・学習関連の記事を、「倭マン日記」とは別に書いていくのだ!

極座標のヤコビ行列とヤコビアン : n次元

いろいろな次元で極座標のヤコビ行列とヤコビアンを求めるシリーズ(目次)。 今回は { n } 次元。 証明はともかく、結果(体積要素)は後で使うので把握しておいてネ。

{ n } 次元極座標

{ n } 次元の極座標は次のように定義されます:

  { \displaystyle \begin{align*}
  &\begin{cases}
      x_1 &= r \cos\theta_1 \\
      x_2 &= r \sin\theta_1 \cos\theta_2 \\
      x_3  &= r\sin\theta_1 \sin\theta_2 \cos\theta_3 \\
              &\vdots \\
      x_{n-1} &= r \sin\theta_1 \cdots \sin\theta_{n-2} \cos \theta_{n-1} \\
      x_n &= r \sin\theta_1 \cdots \sin\theta_{n-2} \sin \theta_{n-1}
  \end{cases} &
  \begin{pmatrix}
      0 \le r < \infty \\
      0 \le \theta_i \le \pi \\
      0 \le \theta_{n-1} \le 2\pi
  \end{pmatrix}
\end{align*}}

ただし { i = 1,\,2,\,\cdots ,\,n-2 }

{ n } 次元極座標の体積要素

{ n }次元極座標ヤコビアンと体積要素は以下のようになります:

  { \displaystyle \begin{align*}
  &\left|\frac{\partial (x_1,\,\cdots,\,x_n)}{\partial (r,\,\theta_1,\,\cdots,\,\theta_{n-1})}\right|
      = r^{n-1} \prod_{i=1}^{n-1}\left(\sin\theta_i\right)^{n-i-1} \\[4mm]
  &dx_1 \cdots dx_n
      = r^{n-1} \left\{\prod_{i=1}^{n-1}\left(\sin\theta_i\right)^{n-i-1} \right\} dr d\theta_1 \cdots d\theta_{n-1}
\end{align*}}

{ n } 次元極座標のヤコビ行列は、書き下すのがあまりにも面倒なのでスルーします。

証明

体積要素の表式はヤコビアンのものを書き換えただけなので、ヤコビアンの表式の証明のみをします。 証明は数学的帰納法でやりましょう。

{ n=2 } の場合を証明する
{ n=2 } のとき、ヤコビアンの表式の右辺は

  { \displaystyle \begin{align*}
  r^{2-1} \prod_{i=1}^{2-1}\left(\sin\theta_i\right)^{2-i-1} = r
\end{align*}}

となり、「極座標のヤコビ行列とヤコビアン : 2次元」で示した2次元極座標の表式と一致します。

{ n-1 } 次元の場合成り立つと仮定して { n } 次元の場合を証明する
以下、{ n \ge 3 } とします。 { n } 次元の直交座標から極座標への変数変換を考えるのですが、この変換に関してまず直交座標から { x_1 } 座標を変換しない円柱座標に変換した上で、そこから極座標に変換するとしましょう。

直交座標 → 円柱座標
まずは直交座標から円柱座標への変換。 { x_1 } 座標を変換しない円柱座標では { x_2,\,x_3,\,\cdots,\,x_n }{ n-1 } 次元極座標へ変換されます。 { n-1 } 次元極座標の動径座標を { \rho } とすると

  { \displaystyle \begin{align*}
  (x_1,\,\cdots,\,x_n) \mapsto (x_1,\,\rho,\,\theta_2,\,\cdots,\,\theta_{n-1})
\end{align*}}

ここで

  { \displaystyle \begin{align*}
  &\begin{cases}
      x_1 &= x_1 \\
      x_2 &= \rho \cos\theta_2 \\
      x_3 &= \rho \sin\theta_2 \cos\theta_3 \\
             &\vdots \\
      x_{n-1} &= \rho \sin\theta_2 \cdots \sin\theta_{n-2} \cos \theta_{n-1} \\
      x_n &= \rho  \sin\theta_2 \cdots \sin\theta_{n-2} \sin \theta_{n-1}
  \end{cases} &
  \begin{pmatrix}
      0 \le \rho < \infty \\
      0 \le \theta_i \le \pi \\
      0 \le \theta_{n-1} \le 2\pi
  \end{pmatrix}
\end{align*}}

となります。 この変換のヤコビアン{ n-1 } 次元極座標ヤコビアンに等しく、それは仮定より

  { \displaystyle \begin{align*}
  \left|\frac{\partial (x_1,\,\cdots,\,x_n)}{\partial (x_1,\,\rho,\,\theta_2,\,\cdots,\,\theta_{n-1})}\right|
      &= \left|\frac{\partial (x_2,\,\cdots,\,x_n)}{\partial (\rho,\,\theta_2,\,\cdots,\,\theta_{n-1})}\right| \\
      &= \rho^{n-2} \prod_{i=2}^{n-1}\left(\sin\theta_i\right)^{n-i-1}
\end{align*}}

となります。

円柱座標 → 極座標
次は円柱座標から極座標への変換。 これは、角度座標 { \theta_2, \cdots, \theta_{n-1} } に関しては変換せず、{ r,\,\rho } に関しては2次元の極座標へ変換します:

  { \displaystyle \begin{align*}
  (x_1,\,\rho,\,\theta_2,\,\cdots,\,\theta_{n-1}) \mapsto (r,\,\theta_1,\,\theta_2\cdots,\,\theta_{n-1})
\end{align*}}

ここで

  { \displaystyle \begin{align*}
  &\begin{cases}
      x_1 &= r \cos\theta_1 \\
      \rho &= r \sin\theta_1
  \end{cases} &
  \begin{pmatrix}
      0 \le r < \infty \\
      0 \le \theta_1 \le \pi
  \end{pmatrix}
\end{align*}}

です。 角度座標 { \theta_2, \cdots, \theta_{n-1} } は変換に関与しないので省略。 この変換のヤコビアンは、2次元極座標の場合と同じく

  { \displaystyle \begin{align*}
  \left|\frac{\partial (x_1,\,\rho,\,\theta_2,\,\cdots,\,\theta_{n-1})}{\partial (r,\,\theta_1,\,\cdots,\,\theta_{n-1})}\right|
      &= \left|\frac{\partial (x_1,\,\rho)}{\partial (r,\,\theta_1)}\right| \\
      &= r
\end{align*}}

となります。

直交座標 → 極座標
以上の結果を踏まえて直交座標から極座標への変換。 合成変換のヤコビ行列はそれぞれの変換のヤコビ行列の積に等しく、行列の積の行列式は各行列の行列式の積に等しいので*1、合成変換のヤコビアンは各変換のヤコビアンの積に等しくなります。 したがって

  { \displaystyle \begin{align*}
  \left|\frac{\partial (x_1,\,\cdots,\,x_n)}{\partial (r,\,\theta_1,\,\cdots,\,\theta_{n-1})}\right|
      &= \left|\frac{\partial (x_1,\,\cdots,\,x_n)}{\partial (x_1,\,\rho,\,\theta_2,\,\cdots,\,\theta_{n-1})}\right|
          \left|\frac{\partial (x_1,\,\rho,\,\theta_2,\,\cdots,\,\theta_{n-1})}
          {\partial (r,\,\theta_1,\,\cdots,\,\theta_{n-1})}\right| \\
      &= \left\{\rho^{n-2} \prod_{i=2}^{n-1}\left(\sin\theta_i\right)^{n-i-1}\right\} r \\
      &= r \left(r\sin\theta_1\right)^{n-2} \prod_{i=2}^{n-1}\left(\sin\theta_i\right)^{n-i-1} \qquad
        (\because \rho = r\sin\theta_1) \\
      &= r^{n-1} \prod_{i=1}^{n-1}\left(\sin\theta_i\right)^{n-i-1} 
\end{align*}}

よって { n } 次元の場合も成立する。

*1:要は、行列 { A,\,B } に対して、{ |AB| = |A||B| } が成り立つってこと。