\( \require{cancel} \)
earth
湘南理工学舎
[戻る]
2024/08/24   
2020/02/24

 楽しく学ぶ…物理数学

 座標の回転と座標変換

(rotation of coodinate and coordinate transformation)

 --目 次--
回転移動と回転行列
回転座標の座標変換(1)
回転座標の座標変換(2)

1.回転移動と回転行列
 2次元空間内(※)におけるx y 直交座標(デカルト座標)の点の回転移動と回転行列(rotation matrix)について説明します。
注:(※)ユークリッド空間ともいいます、ユークリッド空間は簡単にいうと直交座標における幾何であり、相対論での曲がった空間と区別するときなどに使います。
 【fig1】
rotation matrix
上図のx y 直交座標において原点O を中心として点 P(x,y) を角度θ だけ回転させたときの点 P'(x',y') とすると下式が成り立つ。
\((x,y)=( r\ cosα,r\ sinα )\)

\((x',y')= (r\ cos(α+θ),r\ sin(α+θ) )\)


以下は加法定理を使い式展開する。
\(x'=r\ cos(α+θ)=r\ (cosα\ cosθ -sin α\ sinθ)\)
\(\quad =r\ cosα\ cosθ -r\ sinα\ sinθ \)

\(\underline{x'= x\ cosθ-y\ sinθ}\)


\(y'=r\ sin(α+θ)=r\ (sinα\ cosθ +cosα\ sinθ)\)
\(\quad =r\ sinα\ cosθ + r\ cosα\ sinθ \) \(= y\ cosθ+x\ sinθ \)

\( \underline{y'= x\ sinθ+y\ cosθ} \)

上式を以下の行列の積の形に変形します。
\(\begin{bmatrix}x'= ax + cy \\ y'=bx + dy \end{bmatrix}\) \(=\begin{bmatrix} a & c \\ b & d \end{bmatrix}\) \(\begin{bmatrix} x \\ y \end{bmatrix}\)
行列の積の【参照先】

すなわち:

\( \begin{bmatrix} x’ \\ y’ \end{bmatrix}\) \(=\begin{bmatrix} cosθ & -sinθ \\ sinθ& cosθ\end{bmatrix}\) \( \begin{bmatrix} x \\ y \end{bmatrix}\)


この式の中の
\(R(θ)=\begin{bmatrix} cosθ & -sinθ \\ sinθ& cosθ\end{bmatrix}\)\(:❶\)

回転行列といいます。 (反時計回り(CCW)が正)

時計回転り(cw)は負、下式です。
\(R(-θ)=\begin{bmatrix} cos(-θ) & -sin(-θ) \\ sin(-θ)& cos(-θ)\end{bmatrix}\)

\( \quad =\begin{bmatrix} cosθ & sin(θ) \\ -sin(θ)& cos(θ)\end{bmatrix}\)\(:❷\)


2.回転座標の座標変換(1)
  ここでの目的は xy座標での点P(x,y) を新しい座標(x',y') で表す…すなわち座標変換を行うことです。
座標変換とはある座標上の点の位置を別の座標で見るとどういった位置になるのか。
具体的に新座標での位置(x',y') が 「元の座標(x,y) と変換係数(ここでは座標変換行列)を介したもの」で 表すことです。
ここでは点の位置ですが, ベクトルでも同じです。
座標変換では上記で求めた回転行列の \(R(θ)\) と\(R(-θ)\) を座標変換行列として使います。
原点が同じで、元の xy 座標より角度 θ を CCW方向に回転した新座標です。
新座標での点 P' は元の座標の点 P と合じ点にいます。それが fig2 です。
座標表示は 新座標= P'(x', y')、 元の座標= P ( x, y) です。
この状態(fig2)から(fig3)を経由して, 新座標を角度θ、逆回転(CW方向)させ元の座標と同じすれば(fig4)、
P と P’ の角度は θ となります、注意すべきは、逆転させたので、回転行列は\(R(-θ)\)となります。

回転座標の場合は上記の回転行列\(R(-θ)\) を座標変換行列として使います。
すなわち目的の式は:

\( \begin{bmatrix} x’ \\ y’ \end{bmatrix}\) \( =\begin{bmatrix} cosθ & sinθ \\ -sinθ& cosθ\end{bmatrix}\) \( \begin{bmatrix} x \\ y \end{bmatrix}\)

\(=R(-θ)\) \( \begin{bmatrix} x \\ y \end{bmatrix}\)\(:❸\)

これを展開すると:
\(x'=x\ cosθ +y\ sinθ\)

\(y'=-x\ sin θ +y\ cosθ\)

となります。
目的であった、新座標での成分が元の座標成分で表示できました。

rotation matrix
  【fig2】
rotation matrix
  【fig3】

rotation matrix
【fig4】
補足すると、上記のことは座標変換しても、変換係数を介しての変換なので、2次元空間でのベクトル、点の位置などの情報量は保存されます。

3.回転座標の座標変換(2)
 上記と異なる方法により式❸を導出してみます。
rotation matrix
【fig5】
rotation matrix
【fig6】

 上図(fig5,6)を参照し, 下式を見て下さい。
\(x=rcosα\)\(\ ,\ y=rsinα\) \(:ⓐ\)

\(x'=rcos(α-β)\) \(\ ,\ y'=rsin(α-β)\)

xy座標とx'y'座標での同一点に対しての座標変換です。
同一点での(x, y) と(x' y') の関係を求めるため,fig5 のx'y'座標(青色線)をfig6 のように水平の位置におきました。
このことは x'y'座標(青色線)を時計方向にβ分回転したことです。
これが加法定理の角が\((α-β)\)である由来で, 結果として以下のように, 直接, 逆行列\(R(-θ)\) を求まるイメージです。

\(x'=rcos(α-β)\) \(=r(cosα cosβ+sinα sinβ)\)
\(= \ul{rcosα}\ cosβ+ \ul{rsinα}\ sinβ \)
上式に\(ⓐ\) を代入, これにより\(x'\) と\(x, y\) の関係式になる
\(=\ul{x}\ cosβ+ \ul{y}\ sinβ\)\(:ⓑ\)

同様にして
\(y'=rsin(α-β)\) \(=r(sinα\ cosβ-cosα\ sinβ)\)
\(=-x\ sinβ + y\ cosβ\) \(:©\)

これを行列式を使い,式ⓑ,©をまとめると:
\( \begin{bmatrix} x’ \\ y’ \end{bmatrix}\) \( =\begin{bmatrix} cosθ & sinθ \\ -sinθ& cosθ\end{bmatrix}\) \( \begin{bmatrix} x \\ y \end{bmatrix}\)

\(=R(-θ)\) \( \begin{bmatrix} x \\ y \end{bmatrix}\)

上式の\(❸\) と同じになりました。

…以上

coffe

[コーヒーブレイク/閑話]…お疲れさまでした