行列 4). 2次元と3次元の座標変換と回転行列 オイラー角の入り口まで

はじめに

回転行列の見かたについては、前回の記事にて

この記事ではベクトルの内積（投影長）を切り口に行列の積を組み合わせた軸変換とその見かた、そこから回転行列への展開までを書いた。

今回は回転行列の中身を再度。方向余弦を使って導く。

まずは、図的に理解しやすい平面の回転行列（2×2 行列)からはじめて、3次元の回転行列（ 3×3 行列 ) へと話をすすめる。

ま、とりあえず3次元の回転行列までおさえておけば、大概なんとかなる
（オイラー角とかとか）。。

二次元(平面)での回転行列

まずは方向余弦を使った座標変換から
（２ｘ２行列）：二次元(平面)

使用する座標軸

共通の原点Oをもつ、二つの座標系 O-X₀Y₀ と O-X₁Y₁ を考える。
座標軸は右手系、角度はX軸をY軸に向かう方向に回転させる向きを正）とする。ここで

• X₀軸、 Y₀軸方向の単位ベクトルをそれぞれ \(\overrightarrow{i_0} \)、\(\overrightarrow{j_0} \)
• X₁軸、 Y₁軸方向の単位ベクトルをそれぞれ \(\overrightarrow{i_1} \)、\(\overrightarrow{j_1} \)

とする。

各軸の方向余弦

O-X₀Y₀から見たX₁Y₁ 各軸の方向余弦

二つの座標系 O-X₀Y₀ と O-X₁Y₁ の各軸についてお互いのなす角を

• X₁軸に対して、X₀軸 、Y₀軸がなす角を \( \alpha_0 \)、\( \beta_0 \)
• Y₁軸に対してX₀軸 、Y₀軸がなす角を \( \alpha_1 \)、\( \beta_1 \)

（ここで、X₀軸 ,Y₀軸は直行していることから \( \alpha_n \)+\( \beta_n \) =\(\displaystyle\frac{\pi}{2} ~~~(n=0,1)\) ）

とする。

また、このO-X₀Y₀ の座標軸の方向ベクトル \(\overrightarrow{i_0}=\left ( \begin{array}{c} 1\\ 0 \end{array} \right ), \overrightarrow{j_0}=\left ( \begin{array}{c} 0\\ 1 \end{array} \right ) \)を使えば、
\(\overrightarrow{i_1} \)、\(\overrightarrow{j_1} \) はそれぞれ、

\( \begin{align} \overrightarrow{i_1}&= \cos\alpha_0\cdot \overrightarrow{i_0} + \cos\beta_0\cdot \overrightarrow{j_0} =\cos\alpha_0\left ( \begin{array}{c} 1\\ 0 \end{array} \right ) +\cos\beta_0\left ( \begin{array}{c} 0\\ 1 \end{array} \right ) \\ &=\left ( \begin{array}{c} \cos\alpha_0\\ \cos\beta_0 \end{array} \right ) ・・・① \end{align}\) \(\\[4pt]\) \( \begin{align} \overrightarrow{j_1}&= \cos\alpha_1\cdot \overrightarrow{i_0} + \cos\beta_1\cdot \overrightarrow{j_0} =\cos\alpha_1\left ( \begin{array}{c} 1\\ 0 \end{array} \right ) +\cos\beta_1\left ( \begin{array}{c} 0\\ 1 \end{array} \right ) \\& =\left ( \begin{array}{c} \cos\alpha_1\\\cos\beta_1 \end{array} \right ) ・・・② \end{align}\)

この余弦(cos)を使った方向ベクトル表示は、方向余弦とよばれる

• X₁軸の方向余弦: \( \overrightarrow{i_1}=\left ( \begin{array}{c} \cos\alpha_0\\\cos\beta_0 \end{array} \right ) \) \(\\[4pt]\)
• Y₁軸の方向余弦: \(\overrightarrow{j_1}=\left ( \begin{array}{c} \cos\alpha_1\\\cos\beta_1 \end{array} \right ) \)

\(\overrightarrow{i_1}\)、\(\overrightarrow{j_1}\) は、定義から原点Oを始点にもち常に大きさ1のベクトル。角度は任意にとる事ができるのでベクトルの終点の軌跡は円を描く。

（円の公式に当てはめれば）以下も成立

cos²α₀ + cos²β₀ = 1、かつ、cos²α₁ + cos²β₁ = 1 ・・・③

O-X₁Y₁から見たX₀Y₀ 各軸の方向余弦

また逆に X₁軸＆Y₁軸 側から、X₀軸＆Y₀軸を表すと、

• X₀軸に対してX₁軸 、Y₁軸がなす角 ：α₀、α₁
• Y₀軸に対してX₁軸 、Y₁軸がなす角 ：β₀、β1

となるため、①②式同様に

\(\overrightarrow{i_0}= \cos\alpha_0\cdot \overrightarrow{i_1} + \cos\alpha_1\cdot \overrightarrow{j_1} =\left ( \begin{array}{c} \cos\alpha_0 \\ \cos\alpha_1 \end{array} \right ) ・・・④ \) \(\\[2pt]\) \(\overrightarrow{j_0}= \cos\beta_0 \cdot \overrightarrow{i_1} + \cos\beta_1\cdot \overrightarrow{j_1} =\left ( \begin{array}{c} \cos\beta_0 \\\cos\beta_1 \end{array} \right ) ・・・⑤ \)

が成立。つまり

• \(\left ( \begin{array}{c} \cos\alpha_0 \\ \cos\alpha_1 \end{array} \right ) \) : \(\overrightarrow{i_0} \)の方向余弦（X₀軸）\(\\[2pt] \)
• \(\left ( \begin{array}{c} \cos\beta_0 \\\cos\beta_1 \end{array} \right ) \) : \(\overrightarrow{j_0} \)の方向余弦（Y₀軸）

また \(\overrightarrow{i_0} \)、\(\overrightarrow{j_0} \)の大きさは、定義から常に１,よってその方向余弦も③と同じ性質をもつ

cos²α₀ + cos²α₁ = 1、 cos²β₀ + cos²β₁ = 1 ・・・⑥

行列化: 2ｘ2行列の座標軸変換

さて、次に平面上の任意の点Aのベクトル OA→ を考える。

その成分は、それぞれの座標系にて

• O-X₀Y₀ 上のOA→成分: (a_x0 a_y0)
• O-X₁Y₁ 上のOA→成分: (a_x1 a_y1)

とする。

(同じ点Aを表す) OA→ を O-X₀Y₀ 、O-X₁Y₁上でそれぞれ表すと、

OA→ = a_x0・ i₀→ + a_y0・ j₀→ ・・・⑦
OA→ = a_x1・ i₁→ + a_y1・ j₁→ ・・・⑧ である。

さて、⑦式 に i₀→、j₀→ （④⑤式）を代入する

OA→ = a_x0・ i₀→ + a_y0・ j₀→
= a_x0 (cos α₀・ i₁→ + cos α₁・j₁→)
+ a_y0 (cos β₀・ i₁→ + cos β₁・j₁→)

並び替えて

= (a_x0・cos α₀ + a_y0・cos β₀) i₁→
+ (a_x0・cos α₁ + a_y0・cos β₁) j₁→ ・・・⑨

同じOA→を表している事から、⑧式＝⑨式。

つまり、

a_x1・ i₁→ + a_y1・ j₁→
= (a_x0・cos α₀ + a_y0・cos β₀) i₁→
+ (a_x0・cos α₁ + a_y0・cos β₁) j₁→

\( \overrightarrow{i_1}≠0 、\overrightarrow{j_1}≠0 \) より両辺の係数は等しい。よって

a_x1 = a_x0・cos α₀ + a_y0・cos β₀
a_y1 = a_x0・cos α₁ + a_y0・cos β₁

である事がわかる。これを行列で表わして

\(\left ( \begin{array}{c} a_{x1} \\ a_{y1} \end{array} \right ) =\left ( \begin{array}{cc} \cos\alpha_0&\cos\beta_0 \\ \cos\alpha_ 1 & \cos\beta_1 \end{array} \right ) \left ( \begin{array}{c} a_{x0} \\ a_{y0} \end{array} \right ) ・・・⑩ \)

⑩式の 行列部分

\(R_\alpha =\left ( \begin{array}{cc} \cos\alpha_0&\cos\beta_0 \\ \cos\alpha_ 1 & \cos\beta_1 \end{array} \right ) ・・・⑪ \)

これが、二次元の座標変換 ( <- 回転行列）。