学びのヒント by SLA

東北大生のための「学びのヒント」をSLAがお届けします。
理解が深まったり、学びがもっと面白くなる、そんな情報を発信していきます。

2019/09/02 学習ポイント

ヤコビアンを用いた座標変換

前置き

ヤコビアン J とは、N 個の変数 p_{1},p_{2}, \dots ,p_{N}N 個の変数 q_{1},q_{2}, \dots ,q_{N} に依存している時、つまり

(1)   \begin{equation*} p_{i} = p_{i}(q_{1}, \dots , q_{N}) \end{equation*}

である時、

(2)   \begin{equation*} J \equiv \frac{\partial(p_{1}, p_{2} \dots , p_{N})}{\partial(q_{1}, q_{2} \dots , q_{N})} \equiv \mathrm{det} \begin{bmatrix} \frac{\partial p_{1}}{\partial q_{1}} & \cdots & \frac{\partial p_{N}}{\partial q_{1}} \\ \hdotsfor{3} \\ \frac{\partial p_{1}}{\partial q_{N}} & \dots & \frac{\partial p_{N}}{\partial q_{N}} \end{bmatrix} \end{equation*}

で定義される行列式です。(2)式の真ん中の辺はただの記法で、どの変数たちとどの変数たちが対応しているかを分かりやすくするためのものです。
このヤコビアンが、実際に物理数学でどのように利用されているか見てみましょう。

 

ヤコビアンの利用方法

皆さんは、重積分を計算する際の変数変換(例えば x,y,z の直交座標から r,\theta,\varphi の極座標への変換)が面倒だと思ったことはありませんか。特に習いたての頃は、いちいち図を描いて微小面積要素

(3)   \begin{equation*} dxdy = rdrd\theta \end{equation*}

及び微小体積要素

(4)   \begin{equation*} dxdydz = r^{2} \sin \theta dr d\theta d\varphi \end{equation*}

を求めていた(いる)のではないでしょうか。

この求め方は座標変換の基本で大事な考え方ですが、複雑な座標変換や変数の数が4つ以上になると対応できなくなってしまいます。しかし、ヤコビアンを用いることで、面倒な座標変換を図に頼らず、機械的にできるようになります(極座標で図を描いて微小体積要素を求められない人は、そちらに慣れてからヤコビアンを使うことをお勧めします。)

さて、ヤコビアンの利用方法ですが、ある座標系 p_{1}, \dots ,p_{N} とある座標系 q_{1}, \dots ,q_{N} が(1)式の関係にあるとき、p_{1}, \dots ,p_{N} での重積分

    \[ \int \cdots \int f(p_{1}, \dots p_{N}) dp_{1} \cdots dp_{N} \]

は、(2)で定義したヤコビアン J を用いて、

(5)   \begin{equation*} \int \cdots \int f(p_{1}, \dots p_{N}) dp_{1} \cdots dp_{N} = \int \cdots \int f(p_{1} \{q_{i}\}, \dots p_{N} \{q_{i}\}) J dq_{1} \cdots dq_{N} \end{equation*}

と表せます。ただし略記法として p_{j} \{q_{i}\}=p_{j} (q_1,q_2,⋯,q_N ) を用いた。(なぜこのように表せるのかはここでは議論しません。興味のある人は、多変数関数の微分積分についての数学の専門書を読んでみてください。)

 さて、これを実際の重積分での座標変換、直交座標系から極座標系への変換に適用してみましょう。まず、2次元では(3)より

(6)   \begin{equation*} \int \int f(x, y) dxdy = \int \int f(r, \theta) rdrd\theta \end{equation*}

と変換できることが分かっています。では、ヤコビアンを利用してこれを求めてみましょう。この場合のヤコビアンは(2)から、

    \[ J = \frac{\partial (x, y)}{\partial (t, \theta)} = \mathrm{det} \begin{bmatrix} \frac{\partial x}{\partial r} & \frac{\partial y}{\partial r} \\ \frac{\partial x}{\partial \theta} & \frac{\partial y}{\partial \theta} \end{bmatrix}\]

となって、

    \[ x = r \cos \theta, y = r \sin \theta \]

から最右辺を計算すると、

    \[ J = r \]

となります。これを(5)に適用すると、

    \[ \int \int f(x, y) dxdy = \int \int f(r, \theta) Jdrd\theta = \int \int f(r, \theta) rdrd\theta \]

となって、(6)と一致していることが分かります。

 次に、3次元において直交座標系から極座標系への変換を行ってみましょう。やることは2次元と一緒です。まず、ヤコビアンを求めます。(4)より、

(7)   \begin{equation*} \int \int \int f(x, y, z) dxdydz = \int \int \int f(r, \theta, \varphi) r^{2} \sin \theta dr d\theta d\varphi \end{equation*}

であることが分かっています。では、ヤコビアンを用いてこの表式を求めましょう。(2)より、この場合のヤコビアンは、

    \[ J = \frac{\partial (x, y, z)}{\partial (r, \theta, \varphi)} = \mathrm{det} \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 \varphi} & \frac{\partial y}{\partial \varphi} & \frac{\partial z}{\partial \varphi} \\ \end{bmatrix} \]

となります。変数の数が多くなると、必然的に行列式の計算が大変になるので、線形代数学の勉強をしっかりしておきましょう。

    \[ x = r \sin \theta \cos \varphi, y = r \sin \theta \sin \varphi, z = r \cos \theta \]

から最右辺を計算すると、

    \[ J = r^{2} \sin \theta \]

となります。これを(5)に適用すると、

    \begin{align*} \int \int \int f(x, y, z)dxdydz &= \int \int \int f(r, \theta, \varphi) Jdr d\theta d\varphi \\ &= \int \int \int f(r, \theta, \varphi) r^{2}\sin \theta dr d\theta d\varphi \end{align*}

となって、(7)と一致していることが分かります。

 このように、ヤコビアンを用いることで、重積分の変数変換を機械的に行うことができるようになり、覚えておくと便利です。ぜひ、いろいろな座標系に適用して、使いこなせるようになってください。

作成者:SLA