【力学の道具箱】ベクトルの成分表示と座標変換

力学

力や速度など、大きさだけでなく向きも問題となる量はベクトルと呼ばれます。ベクトルは図としては矢印で表示できますが、計算などでベクトルを数値的に扱いたい場合には、ベクトルの成分表示を利用します。

ベクトルの成分表示

下図のように、ベクトル \(\vec{a}\) を成分表示する場合には、\(\vec{a}\) を平行移動して基準となる座標系 \(O-xyz\) の原点に始点が来るように置きます。このときの終点の座標 \((a_x, a_y, a_z)\) がベクトル \(\vec{a}\) の成分と呼ばれるものです。[1]

座標系の \(x\)軸、\(y\)軸、\(z\) 軸方向の単位ベクトル(大きさ1のベクトル)をそれぞれ \(\vec{e}_x\)、\(\vec{e}_y\)、\(\vec{e}_z\) とすれば、

$$\vec{a}=a_x\vec{e}_x+a_y\vec{e}_y+a_z\vec{e}_z$$

が成り立ちます。上式の右辺の各項は \(\vec{a}\) の成分ベクトルと呼ばれます。[1]

ベクトル \(\vec{a}\) の成分表示は、次のように成分を列行列にして表します。

$$\mathbf{a}=\begin{bmatrix} a_x \\ a_y \\ a_z \end{bmatrix}$$

\(\vec{a}\) を幾何ベクトル、\(\mathbf{a}\) を代数ベクトルと呼ぶこともあります。[2]

ベクトル成分表示の演算

任意のベクトル \(\vec{a}\)、\(\vec{b}\) の成分表示をそれぞれ \(\mathbf{a}\)、\(\mathbf{b}\) とすると、幾何ベクトルの演算は、以下のように代数ベクトルの演算に置き換えられます。[2]

  • 和: \(\vec{a}+\vec{b} \quad \longrightarrow \quad \mathbf{a}+\mathbf{b}\)
  • 内積: \(\vec{a}\cdot\vec{b} \quad \longrightarrow \quad \mathbf{a}^T\mathbf{b}\)
  • 絶対値: \(|\vec{a}| \quad \longrightarrow \quad \sqrt{\mathbf{a}^T\mathbf{a}}\)
  • 外積: \(\vec{a}\times\vec{b} \quad \longrightarrow \quad \tilde{\mathbf{a}}\mathbf{b}\)

ここで、上付きの \(T\) は行列の転置を表し、\(\tilde{\mathbf{a}}\) は \(\mathbf{a}\) のチルダ行列と呼ばれる歪対称行列で次のように定義されます。

$$\mathbf{a}=\begin{bmatrix} a_x \\ a_y \\ a_z \end{bmatrix} \quad \Rightarrow \quad \tilde{\mathbf{a}}=\begin{bmatrix} 0 & -a_z & a_y \\ a_z & 0 & -a_x \\ -a_y & a_x & 0 \end{bmatrix}$$

ベクトル成分表示の座標変換

ベクトルの成分表示は座標系に依存します。下図のように、座標系 \(O-xyz\) と平行でない座標系 \(O^{\prime}-x^{\prime}y^{\prime}z^{\prime}\) を考えます。(一般性を失うことなく、原点 \(O\) と \(O^\prime\) が一致したケースで考えることができます。)

ベクトル \(\vec{a}\) の座標系 \(O-xyz\) と座標系 \(O^{\prime}-x^{\prime}y^{\prime}z^{\prime}\) による成分表示をそれぞれ

$$\mathbf{a}=\begin{bmatrix} a_x \\ a_y \\ a_z \end{bmatrix}, \quad \mathbf{a}^\prime=\begin{bmatrix} a_x^\prime \\ a_y^\prime \\ a_z^\prime \end{bmatrix}$$

とすると、一般に \(\mathbf{a}\ne\mathbf{a}^\prime\) です。

では、\(\mathbf{a}\) と \(\mathbf{a}^\prime\) の関係はどのように記述できるでしょうか。

上図に示されているように、\(x^\prime\)軸、\(y^\prime\)軸、\(z^\prime\)軸方向の単位ベクトルの座標系 \(O-xyz\) による成分表示をそれぞれ \(\mathbf{n}\)、\(\mathbf{t}\)、\(\mathbf{b}\) とすると、次の関係式が成り立ちます。

$$\begin{align} \mathbf{a} &= a_x^\prime\mathbf{n}+a_y^\prime\mathbf{t}+a_z^\prime\mathbf{b} \\ &= \begin{bmatrix} \mathbf{n} & \mathbf{t} & \mathbf{b} \end{bmatrix}\begin{bmatrix} a_x^\prime \\ a_x^\prime \\ a_x^\prime \end{bmatrix} \end{align}$$

ここで、

$$\mathbf{R}=\begin{bmatrix} \mathbf{n} & \mathbf{t} & \mathbf{b} \end{bmatrix}$$

と置くと

$$\mathbf{a}=\mathbf{R}\mathbf{a}^\prime \tag{1}$$

を得ます。行列 \(\mathbf{R}\) は座標系 \(O-xyz\) に対する座標系 \(O^{\prime}-x^{\prime}y^{\prime}z^{\prime}\) の回転行列と呼ばれるものです。すなわち、回転行列を用いれば、ベクトルのある座標系による成分表示を別の座標系による成分表示に変換できるというわけです。[3]

回転行列は正規直交性という性質を持っており

$$\mathbf{R}^T=\mathbf{R}^{-1} \tag{2}$$

を満たすため、式(1)と(2)より、次の逆変換の関係も得られます。

$$\mathbf{a}^\prime=\mathbf{R}^T\mathbf{a}$$

行列 \(\mathbf{R}\) と \(\mathbf{R}^T\) は座標変換行列とも呼ばれます。[2][3]

回転行列については、以下の記事で詳しく解説していますのでよろしければご覧ください。

回転行列(方向余弦行列)とは?定義・性質・3つの物理的な意味・公式のまとめ

参考文献

[1] Spiegel, M. R; Liu, J. Mathematical Handbook of Formulas and Tables, 2nd edition, Schaum’s Outline Series, McGraw-Hill, 1999.

[2] 日本機械学会編. マルチボディダイナミクス(1)基礎理論, コンピュータダイナミクスシリーズ3, コロナ社, 2006.

[3] Asada, H.; Slotine, J.-J. E.  Robot Analysis and Control, John Wiley & Sons, 1986.

コメント