============================================================ rotの座標変換 ============================================================ ベクトル解析で出てくるrot（ローテーション、回転）は、 他のデカルト座標系では果たして本当にベクトルとして振舞うのかを 調べてみました。 前の記事は、 divの座標変換不変性_ です。 次の記事は、 gradの座標変換_ です。 rotについて基本的なこと ============================= ベクトル $\bm{A}$ に対して、 $\textrm{rot}$ は以下のように 表わされます。デカルト座標系 $S$ の基底ベクトルを $\bm{e}_1,\bm{e}_2,\bm{e}_3$ とすると [*]_ 、 .. [*] : この記事では添え字について、 $1,2,3$ は、それぞれ $x,y,z$ を表すものとします。 \textrm{rot}\bm{A} &= \bigl( \frac{\partial A_3}{\partial 2} - \frac{\partial A_2}{\partial 3} \bigr) \bm{e}_1 \\ &+\bigl( \frac{\partial A_1}{\partial 3} - \frac{\partial A_3}{\partial 1} \bigr) \bm{e}_2 \\ &+\bigl( \frac{\partial A_2}{\partial 1} - \frac{\partial A_1}{\partial 2} \bigr) \bm{e}_3 \tag{##} この式は、エディントンのイプシロンと呼ばれる次の記号を使うと簡単に書けます。 \varepsilon_{ijk} = \begin{cases} 1 & (\textrm{even}) \\ -1 & (\textrm{odd}) \\ 0 & (\textrm{the \ others}) \end{cases} \tag{##} ただし、添え字が偶置換の時を $(\textrm{even})$ 、奇置換の時を $(\textrm{odd})$ それ以外の３つの添え字の中に同じ数が入る時を $(\textrm{the \ others})$ で表わしました。 大抵、この記号はアインシュタインの縮約という「同じ添え字について、和を取ってΣ記号を省略する」ルール [*]_ と一緒に用いられます。 .. [*] : たとえば、 $\bm{A}\cdot \bm{B}= \sum_i A_i B_i = A_i B_i$ のように三次元ベクトルの内積を表します。 更に、例えば $x$ に関する偏微分を、 $\partial_1$ などと書くことにします。 すると、式 $(1)$ は、 \textrm{rot}\bm{A} = \varepsilon_{ijk} \partial_j A_k \bm{e}_i \tag{##} となります。簡潔ですね。 これから、座標の変換を議論します。そこで、 divの座標変換不変性_ で用いたのと同じ変換行列 $U$ を用います。 つまり、座標系 $S$ のベクトル $A_i$ と座標系 $S^\prime$ のベクトル $A_i^\prime$ について、 実直交行列 $U$ が存在して、次のように表わされます。 \begin{pmatrix} A_1 \\ A_2 \\ A_3 \end{pmatrix} &= U \begin{pmatrix} A_1^\prime \\ A_2^\prime \\ A_3^\prime \end{pmatrix} \\ &= \begin{pmatrix} u_{11} & u_{12} & u_{13} \\ u_{21} & u_{22} & u_{23} \\ u_{31} & u_{32} & u_{33} \end{pmatrix} \begin{pmatrix} A_1^\prime \\ A_2^\prime \\ A_3^\prime \end{pmatrix} \tag{##} \begin{pmatrix} \partial_1 \\ \partial_1 \\ \partial_1 \end{pmatrix} &= \begin{pmatrix} u_{11} & u_{12} & u_{13} \\ u_{21} & u_{22} & u_{23} \\ u_{31} & u_{32} & u_{33} \end{pmatrix} \begin{pmatrix} \partial_1^\prime \\ \partial_2^\prime \\ \partial_3^\prime \end{pmatrix} \tag{##} ここで、例えば式 $(4)$ を書き直すと、 \bm{A}= A_i \bm{e}_i = u_{ih} A_h^\prime \bm{e}_i^\prime \tag{##} 成分を抜き出すと、 A_i= u_{ij} A_j^\prime となります。 これが、ベクトルの（偏微分も一緒）変換則です。 すると、今示したい変換は、 \textrm{rot}\bm{A} &= (\varepsilon_{ijk} \partial_j A_k) \bm{e}_i \\ &= u_{ih} (\varepsilon_{hlm} \partial_l^\prime A_m^\prime) \bm{e}_i^\prime \tag{##} となります。これは、位置ベクトルの「成分」の変換と同様に振舞います [*]_ 。 .. [*] 「成分」とカギ括弧をつけたのは、基底ベクトルとは異なる変換をされるからです。 詳しくは、 ベクトルの基底の変換_ をご覧ください。 本題 ================ それでは準備が整いましたので、式 $(7)$ の左辺を計算してみましょう。 \textrm{rot}\bm{A} &= (\varepsilon_{ijk} \partial_j A_k) \bm{e}_i \\ &= (\varepsilon_{ijk}u_{jl}\partial_l^\prime u_{km} A_{m}^\prime)\bm{e}_i^\prime \tag{##} これから先は、成分を抜き出して考えていきます。 例えば、第一成分（ $x$ 方向の成分）を書き出すと、 (\textrm{rot} \bm{A})_1 &= \varepsilon_{123} \partial_2 A_3 +\varepsilon_{132} \partial_3 A_2 \\ &= \varepsilon_{123} (u_{21}\partial_1^\prime +u_{22}\partial_2^\prime +u_{23}\partial_3^\prime) (u_{31}A_1^\prime+u_{32}A_2^\prime+u_{33}A_3^\prime) \\ &+\varepsilon_{132} (u_{31}\partial_1^\prime +u_{32}\partial_2^\prime +u_{33}\partial_3^\prime) (u_{21}A_1^\prime+u_{22}A_2^\prime+u_{23}A_3^\prime) \\ &= (u_{21}u_{32}-u_{31}u_{22})\partial_1^\prime A_2^\prime + (u_{21}u_{33}-u_{31}u_{23})\partial_1^\prime A_3^\prime \\ &+(u_{22}u_{31}-u_{32}u_{21})\partial_2^\prime A_1^\prime + (u_{22}u_{33}-u_{32}u_{23})\partial_2^\prime A_3^\prime \\ &+(u_{23}u_{31}-u_{33}u_{21})\partial_3^\prime A_1^\prime + (u_{23}u_{32}-u_{33}u_{22})\partial_3^\prime A_2^\prime \\ &= (u_{22}u_{33}-u_{32}u_{23})(\partial_2^\prime A_3^\prime-\partial_3^\prime A_2^\prime) \\ &+(u_{23}u_{31}-u_{33}u_{21})(\partial_3^\prime A_1^\prime-\partial_1^\prime A_3^\prime) \\ &+(u_{21}u_{32}-u_{31}u_{22})(\partial_1^\prime A_2^\prime-\partial_2^\prime A_1^\prime) \tag{##} 同様に第二、第三成分も書くと、 (\textrm{rot} \bm{A})_2 &= (u_{32}u_{13}-u_{12}u_{33})(\partial_2^\prime A_3^\prime-\partial_3^\prime A_2^\prime) \\ &+(u_{33}u_{11}-u_{13}u_{31})(\partial_3^\prime A_1^\prime-\partial_1^\prime A_3^\prime) \\ &+(u_{31}u_{12}-u_{11}u_{32})(\partial_1^\prime A_2^\prime-\partial_2^\prime A_1^\prime) \tag{##} (\textrm{rot} \bm{A})_3 &= (u_{12}u_{23}-u_{22}u_{13})(\partial_2^\prime A_3^\prime-\partial_3^\prime A_2^\prime) \\ &+(u_{13}u_{21}-u_{23}u_{11})(\partial_3^\prime A_1^\prime-\partial_1^\prime A_3^\prime) \\ &+(u_{11}u_{22}-u_{21}u_{12})(\partial_1^\prime A_2^\prime-\partial_2^\prime A_1^\prime) \tag{##} となります。 ここで、行列 $U$ の余因子行列 $V=\Tilde{U}$ を考えます。これは、逆行列 $U^{-1}$ の $\det U$ 倍であり、 以下のような形をしています。（ $\det U=1$ です。） V &= \begin{pmatrix} v_{11} & v_{12} & v_{13} \\ v_{21} & v_{22} & v_{23} \\ v_{31} & v_{32} & v_{33} \end{pmatrix} \\ &= \Tilde{U} \\ &= \begin{pmatrix} \begin{vmatrix} u_{22} & u_{23} \\ u_{32} & u_{33} \end{vmatrix} & -\begin{vmatrix} u_{12} & u_{13} \\ u_{32} & u_{33} \end{vmatrix} & \begin{vmatrix} u_{12} & u_{13} \\ u_{22} & u_{23} \end{vmatrix} \\ -\begin{vmatrix} u_{21} & u_{23} \\ u_{31} & u_{33} \end{vmatrix} & \begin{vmatrix} u_{11} & u_{13} \\ u_{31} & u_{33} \end{vmatrix} & -\begin{vmatrix} u_{11} & u_{13} \\ u_{21} & u_{23} \end{vmatrix} \\ \begin{vmatrix} u_{21} & u_{22} \\ u_{31} & u_{32} \end{vmatrix} & -\begin{vmatrix} u_{11} & u_{12} \\ u_{31} & u_{32} \end{vmatrix} & \begin{vmatrix} u_{11} & u_{12} \\ u_{21} & u_{22} \end{vmatrix} \end{pmatrix} \\ &= (\det U)U^{-1} \\ &= U^{-1} \tag{##} .. image :: chromel-rotAnotherCoordinates-01-t.png 今、行列 $U$ は実直交行列なので、 $U^{-1}=U^T$ （ $^T$ は転置を表す。）となっています。 よって、 V &= \begin{pmatrix} \begin{vmatrix} u_{22} & u_{23} \\ u_{32} & u_{33} \end{vmatrix} & -\begin{vmatrix} u_{12} & u_{13} \\ u_{32} & u_{33} \end{vmatrix} & \begin{vmatrix} u_{12} & u_{13} \\ u_{22} & u_{23} \end{vmatrix} \\ -\begin{vmatrix} u_{21} & u_{23} \\ u_{31} & u_{33} \end{vmatrix} & \begin{vmatrix} u_{11} & u_{13} \\ u_{31} & u_{33} \end{vmatrix} & -\begin{vmatrix} u_{11} & u_{13} \\ u_{21} & u_{23} \end{vmatrix} \\ \begin{vmatrix} u_{21} & u_{22} \\ u_{31} & u_{32} \end{vmatrix} & -\begin{vmatrix} u_{11} & u_{12} \\ u_{31} & u_{32} \end{vmatrix} & \begin{vmatrix} u_{11} & u_{12} \\ u_{21} & u_{22} \end{vmatrix} \end{pmatrix} \\ &= U^{-1} \\ &= U^T \\ &= \begin{pmatrix} u_{11} & u_{21} & u_{31} \\ u_{12} & u_{22} & u_{32} \\ u_{13} & u_{23} & u_{33} \end{pmatrix} \tag{##} となります。 式 $(9)$ に式 $(13)$ を代入すると、見事に、 (\textrm{rot} \bm{A})_1 &= (u_{22}u_{33}-u_{32}u_{23})(\partial_2^\prime A_3^\prime-\partial_3^\prime A_2^\prime) \\ &+(u_{23}u_{31}-u_{33}u_{21})(\partial_3^\prime A_1^\prime-\partial_1^\prime A_3^\prime) \\ &+(u_{21}u_{32}-u_{31}u_{22})(\partial_1^\prime A_2^\prime-\partial_2^\prime A_1^\prime) \\ &= u_{11}(\partial_2^\prime A_3^\prime-\partial_3^\prime A_2^\prime) \\ &+ u_{12}(\partial_3^\prime A_1^\prime-\partial_1^\prime A_3^\prime) \\ &+ u_{13}(\partial_1^\prime A_2^\prime-\partial_2^\prime A_1^\prime) \tag{##} (\textrm{rot} \bm{A})_2 &= (u_{32}u_{13}-u_{12}u_{33})(\partial_2^\prime A_3^\prime-\partial_3^\prime A_2^\prime) \\ &+(u_{33}u_{11}-u_{13}u_{31})(\partial_3^\prime A_1^\prime-\partial_1^\prime A_3^\prime) \\ &+(u_{31}u_{12}-u_{11}u_{32})(\partial_1^\prime A_2^\prime-\partial_2^\prime A_1^\prime) \\ &= u_{21}(\partial_2^\prime A_3^\prime-\partial_3^\prime A_2^\prime) \\ &+ u_{22}(\partial_3^\prime A_1^\prime-\partial_1^\prime A_3^\prime) \\ &+ u_{23}(\partial_1^\prime A_2^\prime-\partial_2^\prime A_1^\prime) \tag{##} (\textrm{rot} \bm{A})_3 &= (u_{12}u_{23}-u_{22}u_{13})(\partial_2^\prime A_3^\prime-\partial_3^\prime A_2^\prime) \\ &+(u_{13}u_{21}-u_{23}u_{11})(\partial_3^\prime A_1^\prime-\partial_1^\prime A_3^\prime) \\ &+(u_{11}u_{22}-u_{21}u_{12})(\partial_1^\prime A_2^\prime-\partial_2^\prime A_1^\prime) \\ &= u_{31}(\partial_2^\prime A_3^\prime-\partial_3^\prime A_2^\prime) \\ &+ u_{32}(\partial_3^\prime A_1^\prime-\partial_1^\prime A_3^\prime) \\ &+ u_{33}(\partial_1^\prime A_2^\prime-\partial_2^\prime A_1^\prime) \tag{##} が成り立ちます。 よって無事、式 $(7)$ 、 \textrm{rot}\bm{A} &= (\varepsilon_{ijk} \partial_j A_k) \bm{e}_i \\ &= u_{ih} (\varepsilon_{hlm} \partial_l^\prime A_m^\prime) \bm{e}_i^\prime \tag{7} が示せました。 それでは、今日はこの辺で。 続きは こちら_ 。 .. _ベクトルの基底の変換: http://hooktail.sub.jp/vectoranalysis/vectorTransform/ .. _divの座標変換不変性: http://hooktail.sub.jp/vectoranalysis/divAnotherCoordinates/ .. _gradの座標変換: http://hooktail.sub.jp/vectoranalysis/gradAnotherCoordinates/ .. _こちら: http://hooktail.sub.jp/vectoranalysis/gradAnotherCoordinates/ @@author:クロメル@@ @@accept:2010-01-10@@ @@category:ベクトル解析@@ @@id:rotAnotherCoordinates@@