Help us understand the problem. What is going on with this article?

ベクトルの極座標変換

More than 3 years have passed since last update.

※ 本記事のきっかけは、「よく分かる量子力学(前野昌弘)」の第12章「12.1. 3次元極座標のシュレディンガー方程式」に出てきた極座標変換の式をもっと機械的に導出できないのかなと思ったことです。

通常のデカルト座標 $(x, y, z)$ と極座標 $(r, \theta, \phi)$ は次式の関係で結ばれている。

\begin{align}
x &= r\sin\theta\cos\phi \\
y &= r\sin\theta\sin\phi \\
z &= r\cos\theta \\
\end{align}

本記事の目的は、デカルト座標系のベクトル $\vec{A} = (A_x, A_y, A_z)$ を極座標系のベクトル $(A_r, A_\theta, A_\phi)$ に変換する式を導出することである。

多様体の理論で出てくる接ベクトルの表記法を用いると、ベクトル $A$ は次のように表記できる。

\begin{align}
A &= A_x\frac{\partial}{\partial x} + A_y\frac{\partial}{\partial y} + A_z\frac{\partial}{\partial z} \\
  &= A'_r\frac{\partial}{\partial r} + A'_\theta\frac{\partial}{\partial\theta} + A'_\phi\frac{\partial}{\partial\phi}
\end{align}

ここで極座標のベクトル成分にダッシュ($'$)を付けたのは、極座標ベクトルの基底

\frac{\partial}{\partial r}, \frac{\partial}{\partial\theta}, \frac{\partial}{\partial\phi}

が正規直交基底になっていないためである。上記の基底ベクトルは直交はしているものの、デカルト座標系におけるノルムは1ではない。そのため、後のほうでこれらの基底ベクトルを単位ベクトルに正規化する操作を行う。

(実は今回ハマったのはこの正規化であった。この程度の座標変換式なんて簡単に導出できるだろうと高をくくって計算を始めたのだが、教科書の式と答えが合わなくて随分焦った。単位ベクトルに正規化されていないことに気づいて漸く教科書の式と一致する式を導出することが出来た。)

ではこれらの基底ベクトルの変換式を計算しよう。

\begin{align}
\begin{pmatrix}
\frac{\partial}{\partial r} \\
\frac{\partial}{\partial\theta} \\
\frac{\partial}{\partial\phi}
\end{pmatrix}
&=
\begin{bmatrix}
\frac{\partial x}{\partial r} &\frac{\partial y}{\partial r} &\frac{\partial z}{\partial r} \\
\frac{\partial x}{\partial\theta} &\frac{\partial y}{\partial\theta} &\frac{\partial z}{\partial\theta} \\
\frac{\partial x}{\partial\phi} &\frac{\partial y}{\partial\phi} &\frac{\partial z}{\partial\phi} \\
\end{bmatrix}
\begin{pmatrix}
\frac{\partial}{\partial x} \\
\frac{\partial}{\partial y} \\
\frac{\partial}{\partial z}
\end{pmatrix} \\
&=
\begin{bmatrix}
\sin\theta\cos\phi & \sin\theta\sin\phi & \cos\theta \\
r\cos\theta\cos\phi & r\cos\theta\sin\phi & -r\sin\theta \\
-r\sin\theta\sin\phi & r\sin\theta\cos\phi & 0
\end{bmatrix}
\begin{pmatrix}
\frac{\partial}{\partial x} \\
\frac{\partial}{\partial y} \\
\frac{\partial}{\partial z}
\end{pmatrix}
\end{align}

次に極座標の基底ベクトルを正規化する。そのために基底ベクトルのデカルト座標系でのノルムを計算しよう。

\begin{align}
\|\frac{\partial}{\partial r}\|^2 &= (\sin\theta\cos\phi)^2 + (\sin\theta\sin\phi)^2 + \cos^2\theta &=& 1 \\
\|\frac{\partial}{\partial\theta}\|^2 &= (r\cos\theta\cos\phi)^2 + (r\cos\theta\sin\phi)^2 + (-r\sin\theta)^2 &=& r^2 \\
\|\frac{\partial}{\partial\phi}\|^2 &= (-r\sin\theta\sin\phi)^2 + (r\sin\theta\cos\phi)^2 &=& (r\sin\theta)^2 \\
\end{align}

以上の計算により、正規化された極座標の基底ベクトルは次式となることが分かる。

\frac{\partial}{\partial r}, \frac{1}{r}\frac{\partial}{\partial\theta}, \frac{1}{r\sin\theta}\frac{\partial}{\partial\phi}

従って正規化された基底ベクトルの座標変換式は次式となる。

\begin{pmatrix}
\frac{\partial}{\partial r} \\
\frac{1}{r}\frac{\partial}{\partial\theta} \\
\frac{1}{r\sin\theta}\frac{\partial}{\partial\phi}
\end{pmatrix}
=
\begin{bmatrix}
\sin\theta\cos\phi & \sin\theta\sin\phi & \cos\theta \\
\cos\theta\cos\phi & \cos\theta\sin\phi & -\sin\theta \\
-\sin\phi & \cos\phi & 0
\end{bmatrix}
\begin{pmatrix}
\frac{\partial}{\partial x} \\
\frac{\partial}{\partial y} \\
\frac{\partial}{\partial z}
\end{pmatrix}

上式の行列は直交行列であるから、逆行列は単に転置することで得られる。従って上式を逆に解くと次式となる。

\begin{pmatrix}
\frac{\partial}{\partial x} \\
\frac{\partial}{\partial y} \\
\frac{\partial}{\partial z}
\end{pmatrix}
=
\begin{bmatrix}
\sin\theta\cos\phi & \cos\theta\sin\phi & -\sin\theta \\
\sin\theta\sin\phi & \cos\theta\sin\phi & \cos\phi \\
\cos\theta & -\sin\theta & 0
\end{bmatrix}
\begin{pmatrix}
\frac{\partial}{\partial r} \\
\frac{1}{r}\frac{\partial}{\partial\theta} \\
\frac{1}{r\sin\theta}\frac{\partial}{\partial\phi}
\end{pmatrix}

この行列は教科書の (12.11) 式に一致している。

さて、ベクトル $A$ はそれぞれの基底を用いて次のように表記できる。

\begin{align}
A &= A_x\frac{\partial}{\partial x} + A_y\frac{\partial}{\partial y} + A_z\frac{\partial}{\partial z} \\
  &= A_r\frac{\partial}{\partial r} + A_\theta\frac{1}{r}\frac{\partial}{\partial\theta} + A_\phi\frac{1}{r\sin\theta}\frac{\partial}{\partial\phi}
\end{align}

これと上で求めた基底ベクトルの変換式を合わせると、次式の座標変換式が得られる。

\begin{pmatrix} A_r \\ A_\theta \\ A_\phi \end{pmatrix}
=
\begin{bmatrix}
\sin\theta\cos\phi & \sin\theta\sin\phi & \cos\theta \\
\cos\theta\cos\phi & \cos\theta\sin\phi & -\sin\theta \\
-\sin\phi & \cos\phi & 0
\end{bmatrix}
\begin{pmatrix} A_x \\ A_y \\ A_z \end{pmatrix}

以上で教科書の (12.15) 式に一致する式が導出できた。

Why not register and get more from Qiita?
  1. We will deliver articles that match you
    By following users and tags, you can catch up information on technical fields that you are interested in as a whole
  2. you can read useful information later efficiently
    By "stocking" the articles you like, you can search right away
Comments
No comments
Sign up for free and join this conversation.
If you already have a Qiita account
Why do not you register as a user and use Qiita more conveniently?
You need to log in to use this function. Qiita can be used more conveniently after logging in.
You seem to be reading articles frequently this month. Qiita can be used more conveniently after logging in.
  1. We will deliver articles that match you
    By following users and tags, you can catch up information on technical fields that you are interested in as a whole
  2. you can read useful information later efficiently
    By "stocking" the articles you like, you can search right away
ユーザーは見つかりませんでした