#rst2hooktail_source
============================================================
四次元空間中のベクトルの回転
============================================================
この記事では、Joh氏の ベクトルの回転_ を拡張して四次元空間中での回転を考えます。
特殊な場合
=======================
四次元空間中の回転では、回転軸が二つに
なり、 $\bm{\omega}_1,\bm{\omega}_2$ となります。
簡単のため、それらは正規直交化されているものとします。
一番明らかな回転は、回転軸が
<tex>
\bm{\omega}_1=\begin{pmatrix}
\sin \theta \cos \phi \\
\sin \theta \sin \phi \\
\cos \theta \\
0
\end{pmatrix}, \ \ \ \ \
\bm{\omega}_2=\begin{pmatrix}
0 \\
0 \\
0 \\
1
\end{pmatrix} \tag{##}
</tex>
の時です。これは、四番目の次元の方向 $\bm{\omega}_2$ にまず直交し、
二次元球面の極座標表示の示す方向に対しても直交する回転です。
つまり、先に言いましたが、四番目の次元方向 $\bm{\omega}_2$ に直交していても、それに直交する単位ベクトル $\bm{\omega}_1$ をとれることが分かるでしょう。
一般化
=======================
Joh氏の ベクトルの回転_ を一般化します。
拡張すべき式は、
<tex>
\bm{r}^\prime = (\bm{n} \cdot \bm{r})\bm{n} + [\bm{r}-(\bm{n} \cdot \bm{r})\bm{n}]\cos \eta - (\bm{r} \times \bm{n})\sin \eta \tag{##}
</tex>
です。四次元の任意のベクトル $\bm{r}$ を
回転軸
<tex>
\bm{\omega}_1 =
\begin{pmatrix}
\omega_{11} \\
\omega_{12} \\
\omega_{13} \\
\omega_{14}
\end{pmatrix}
, \ \ \
\bm{\omega}_2 =
\begin{pmatrix}
\omega_{21} \\
\omega_{22} \\
\omega_{23} \\
\omega_{24}
\end{pmatrix} \tag{##}
</tex>
の周りに角度 $ \eta $ だけ回転した
ベクトルを $ \bm{r}^\prime $ とするのです。
注意深く考えると、回転軸の方向にはいくら回転しても、値が変化することはありません。
よって、その方向の自由度を奪って回転面内のベクトルにしてしまいましょう。それには、
<tex>
\bm{\alpha}_1 = \bm{r} - (\bm{\omega}_1 \cdot \bm{r}) \bm{\omega}_1 - (\bm{\omega}_2 \cdot \bm{r}) \bm{\omega}_2
\tag{##}
</tex>
を考えれば、回転面内のベクトルで回転の始点となるベクトル
<tex>
\bm{\alpha}_1=
\begin{pmatrix}
\alpha_{11}\\
\alpha_{12}\\
\alpha_{13}\\
\alpha_{14}
\end{pmatrix} \tag{##}
</tex>
ができます。後は、回転面内で $\bm{\alpha}_1$ に直交するベクトルを $\bm{\alpha}_2$ と置いて、 $ \bm{alpha}_1 \cos \eta \pm \bm{alpha}_2 \sin \eta $ としてやれば良さそうですね。 [*]_ この $\bm{\alpha}_2$ は、回転面内と言うところがポイントです。つまり、 $\bm{\alpha}_2$ は、 $\bm{\omega}_1$ , $\bm{\omega}_2$ , $\bm{\alpha}_1$ の全てに直交しているのです。
ができます。後は、回転面内で $\bm{\alpha}_1$ に直交するベクトルを $\bm{\alpha}_2$ と置いて、 $ \bm{\alpha}_1 \cos \eta \pm \bm{\alpha}_2 \sin \eta $ としてやれば良さそうですね。 [*]_ この $\bm{\alpha}_2$ は、回転面内と言うところがポイントです。つまり、 $\bm{\alpha}_2$ は、 $\bm{\omega}_1$ , $\bm{\omega}_2$ , $\bm{\alpha}_1$ の全てに直交しているのです。
これは次の様に求まります。 $x,y,z,w$ 軸方向の単位ベクトルを $\bm{i},\bm{j},\bm{k},\bm{l}$ とすると、行列式の記法で、
.. [*] 三次元では三つのベクトルに対して右手系と左手系が定義されていますが、四次元では標準の様なものが無いので、 $\pm$ を決定できません。
<tex>
\bm{\alpha}_2 = \begin{vmatrix}
\bm{i} & \bm{j} & \bm{k} & \bm{l} \\
\omega_{11} & \omega_{12} & \omega_{13} & \omega_{14} \\
\omega_{21} & \omega_{22} & \omega_{23} & \omega_{24} \\
\alpha_{11} & \alpha_{12} & \alpha_{13} & \alpha_{14}
\end{vmatrix}
\tag{##}
</tex>
となります。そんなに怖いものではありません。例えば、 $\bm{\alpha}_2$ と、 $\bm{\omega}_1$ の内積をとると、式 $(5)$ の第一行目が $\bm{\omega}_1$ に置き換わります。つまり、同じ行が第二行目にあるので式の値はゼロとなり、直交していることが分かります。これで、無事 $\bm{\alpha}_1$ と $\bm{\alpha}_2$ が求まったので、回転が記述できます。
それは、
<tex>
\bm{r}^\prime &= (\bm{\omega}_1 \cdot \bm{r}) \bm{\omega}_1 + (\bm{\omega}_2 \cdot \bm{r}) \bm{\omega}_2 +
[ \bm{r}- (\bm{\omega}_1 \cdot \bm{r}) \bm{\omega}_1 - (\bm{\omega}_2 \cdot \bm{r}) \bm{\omega}_2 ] \cos \eta \pm
\begin{vmatrix}
\bm{i} & \bm{j} & \bm{k} & \bm{l} \\
\omega_{11} & \omega_{12} & \omega_{13} & \omega_{14} \\
\omega_{21} & \omega_{22} & \omega_{23} & \omega_{24} \\
\alpha_{11} & \alpha_{12} & \alpha_{13} & \alpha_{14}
\end{vmatrix} \sin \eta \\
&= (\bm{\omega}_1 \cdot \bm{r}) \bm{\omega}_1 + (\bm{\omega}_2 \cdot \bm{r}) \bm{\omega}_2 + \bm{\alpha}_1 \cos \eta
\pm \bm{\alpha}_2 \sin \eta
\tag{##}
</tex>
となります。
四次元以上への拡張は容易ですね。五次元なら、回転の軸を正規直交した $\bm{\omega}_1$ , $\bm{\omega}_2$ , $\bm{\omega}_3$ にすればいいのです。行列式は、これプラス $\bm{\alpha}_1$ を用いる５×５の行列式になります。
今日はここまで、お疲れさまでした。
.. _ベクトルの回転: http://hooktail.sub.jp/vectoranalysis/vectorRot/
@@author:クロメル@@
@@accept:2013-04-03@@
@@category:ベクトル解析@@
@@id:vectorRotDim4@@
記事ソース/四次元空間中のベクトルの回転 の変更点
