うさぎでもわかる線形代数 第14羽 回転変換(回転行列)・対称変換

スポンサードリンク

こんにちは、ももやまです。

今回は回転変換(回転行列)・対称変換についてのまとめです。

変換の中でも、原点中心に回転させる(回転変換)、原点を通るある直線 l と対称移動させる変換(対称変換)の表現行列の作り方、実際に座標を回転変換、対称変換させる方法のまとめとなっております。

前回の記事はこちら!
線形写像における核空間・像空間と線形写像における全射・単射についてです。

www.momoyama-usagi.com

スポンサードリンク

1.回転変換(2次元・平面)

ある座標における点を原点を中心に反時計回りに θ 回転させた点を求めるような線形変換 f の表現行列 A を考えてみましょう。

表現行列を求める際には、2次元標準基底e1=(10),   e2=(01)の行き先を考えればOK。

単位円を考えると、x2+y2=cos2θ+sin2θ=1と変えられ、単位円上の点Aを (cosθ,sinθ) とすることができます。
θ は点A (1,0) を 0 とし、反時計回り方向を正とする)

e1θ=0 なので、θ 回転させた f(e1) は、f(e1)=(cosθsinθ)となります。

同様に e2θ=90 なので、θ 回転させた f(e2) は、f(e2)=(cos(θ+90)sin(θ+90))=(sinθcosθ)となります*1

f:id:momoyama1192:20190821120309g:plain

よって2次元の回転変換を表す f の表現行列 R(θ) は、R(θ)=(f(e1),f(e2))=(cosθsinθsinθcosθ)となります。

回転変換を表す行列(2次元回転行列)

2次元ベクトル空間において、ある点を原点中心の反時計回りに θ 移す線形変換の表現行列(回転行列)R(θ) は、R(θ)=(cosθsinθsinθcosθ)となる。

回転行列は直交行列となっています。

回転行列を用いることで、元の座標を (x,y) 、原点中心に θ 回転させた座標を (x,y) とすると、(xy)=(cosθsinθsinθcosθ)(xy)=(xcosθysinθxsinθ+ycosθ)と表すことができます。

スポンサードリンク

2.回転変換(3次元回転行列)

先ほど説明した2次元の回転変換を3次元の場合に拡張して考えてみましょう。
2次元の場合と同様に求めることができます。

しかし、3次元の場合は回転する軸(x 軸回り、y 軸回り、z 軸回りの3つ)によって合計3パターンの回転行列が存在します。

いずれの3パターンの場合でも表現行列を求める際には、3次元標準基底e1=(100),   e2=(010),   e3=(001)の行き先を考えればOKです。

(1) x軸回りの回転

まずは x 軸回りの回転変換 f の表現行列 Rx(θ) を求めてみましょう。

x 軸回りの回転なので x 軸方向の単位ベクトルを表す e1変換の前後で向きが変化しませんね。なので、f(e1)=(100)となります。

残りの e2, e3 はそのまま考えるのは少しむずかしいので、x の正の方向からみた yz 座標平面を書いて変換の前後でどのようにベクトルが変わるのかを考えることをおすすめします*2

f:id:momoyama1192:20190821171325g:plain

上のように図示すると、2次元の場合と同じような要領でf(e2)=(0cosθsinθ)f(e3)=(0sinθcosx)を求めることができます。

よって、 x 軸回りに θ 回転する変換 f の表現行列 Rx(θ) は、Rx(θ)=(f(e1),f(e2),f(e3)=(1000cosθsinθ0sinθcosθ)と求められます。

(2) y軸回りの回転

y 軸回りの回転変換 f の表現行列 Ry(θ) も(1)と同様に求められます。

y 軸回りの回転なので y 軸方向の単位ベクトルを表す e2 は変換の前後で向きが変化しませんね。なので、f(e2)=(010)となります。

f:id:momoyama1192:20190821171330g:plain

残り2つのベクトルは(1)と同じように y の正の方向から見た zx 平面を書くことで動きがわかりやすくなります。

z, x の順番に注意してください。

よって、残りの2つのベクトルもf(e1)=(cosθ0sinθ)   f(e3)=(sinθ0cosθ)と求められます。

よって、 y 軸回りに θ 回転する変換 f の表現行列 Ry(θ) は、Ry(θ)=(f(e1),f(e2),f(e3)=(cosθ0sinθ010sinθ0cosθ)と求められます。

(3) z軸回りの回転

z 軸回りの回転変換 f の表現行列 Rz(θ) も(1)と同様に求められます。

z 軸回りの回転なので z 軸方向の単位ベクトルを表す e2 は変化なし、よってf(e3)=(001)となります。

残り2つのベクトルも(1),(2)と同じように x の正の方向から見た xy 平面を書くとすぐわかります。

f:id:momoyama1192:20190821171337g:plain

残りの2つのベクトルはf(e1)=(cosθsinθ0)   f(e2)=(sinθcosθ0)と求められます。

よって、 z 軸回りに θ 回転する変換 f の表現行列 Rz(θ) は、Rz(θ)=(f(e1),f(e2),f(e3)=(cosθsinθ0sinθcosθ0001)と求められます。

3次元回転変換を表す行列(3次元回転行列)

3次元ベクトル空間において、ある点をそれぞれの軸まわりの反時計回りに θ 移す線形変換の表現行列(回転行列) Rx, Ry, Rz は、

x 軸回りの回転を表す回転行列 Rx Rx(θ)=(1000cosθsinθ0sinθcosθ)

y 軸回りの回転を表す回転行列 Rx Ry(θ)=(cosθ0sinθ010sinθ0cosθ)

z 軸回りの回転を表す回転行列 Rz Rz(θ)=(cosθsinθ0sinθcosθ0001)

3次元の場合でも回転行列は直交行列となります。

(4) ロドリゲスの回転公式(おまけ)

任意の軸回りの回転変換 f を表す表現行列を参考までに紹介しておきます。

任意の軸 n=(nx,ny,xz)回りの回転を表す回転行列は、

Rn(θ)=(nx2(1cosθ)+cosθnxny(1cosθ)nzsinθnxnz(1cosθ)+nysinθnynx(1cosθ)+nzsinθny2(1cosθ)+cosθnynz(1cosθ)nxsinθnznx(1cosθ)nysinθnzny(1cosθ)+nxsinθnz2(1cosθ)+cosθ)

と表せます(ロドリゲスの回転公式、ロドリゲスの定理と呼ばれます)。

スポンサードリンク

3.対称変換

原点を通るような直線 l : y=ax に関する対称移動を表す線形変換 f の表現行列を考えてみましょう。

直線 l の傾き a は、x 軸と l のなす角 θ を用いて、a=tanθと表すことができます。つまり、y=(tanθ)xに関する対称移動を求めればいいことになります。

今まで通り2次元標準基底e1=(10),   e2=(01)の行き先を考えることで表現行列を求めていきます。

e1e2 を対称移動させた様子を図示すると、下のようになります。

[修正] β=α2θ と書かれている部分は、正しくは β=αθ です。

f:id:momoyama1192:20190821171346g:plain

e1 から f(e1) への変換は比較的わかりやすいと思います。f(e1)=(cos2θsin2θ)となります。

e2 から e2 への変換が少しわかりにくいかもしれません。

図のように α, β を取ると、e2 から f(e2) に変換する際に 2α 時計回りに回転していることがわかると思います。

0 以降の回転している部分を β とすると、α=90θ なので、β=αθ=(90θ)θ=902θと計算できます。

時計回りに β 回転しているということは、反時計回りには β 回転しているので、e2 は、f(e2)=(cos(2θ90)sin(2θ90))=(sin2θcos2θ)と求められます*3

よって直線 l y=(tanθ)x に関する対称移動を表す線形変換 f の表現行列 S(θ) は、S(θ)=(f(e1),f(e2))=(cos2θsin2θsin2θcos2θ)となります。

ある直線に関する対称変換(対称変換行列)

2次元ベクトル空間において、ある点を直線 l y=(tanθ)x に関する対称移動する線形変換の表現行列(対称変換行列)S(θ) は、S(θ)=(cos2θsin2θsin2θcos2θ)となる。

対称変換行列も回転行列と同様に直交行列となります。

対称変換行列を用いることで、元の座標を (x,y) 、直線 l y=(tanθ)x に関する対称移動後の座標を (x,y) とすると、(xy)=(cos2θsin2θsin2θcos2θ)(xy)=(xcos2θ+ysin2θxsin2θycos2θ)と表すことができます。

4.練習問題

では、簡単にですが、練習をしてみましょう。

練習1

2次元ベクトル空間において、つぎの(1)~(5)が示す対称移動の線形変換行列を求め、①~⑥の中から選びなさい。

(1) x 軸に関する対称移動
(2) y 軸に関する対称移動
(3) 原点に関する対称移動
(4) y=x に関する対称移動
(5) y=x に関する対称移動

f:id:momoyama1192:20190827003415g:plain

練習2

2次元ベクトル空間上における点A (6,4) を(1),(2)の方法で移動させた点の座標を求めなさい。

(1) 原点中心に 30 回転させた点Bの座標
(2) 直線 y=2x について対称移動させた点Cの座標

5.練習問題の答え

解答1

(1) x 軸に関する対称移動

正解:①

θ=0 のときの直線 y=(tanθ)x に関する対称移動と同じ。

なので、S(θ)=(cos0sin0sin0cos0)=(1001)となる。

(2) y 軸に関する対称移動

正解:②

θ=90 のときの直線 y=(tanθ)x に関する対称移動と同じ。

なので、S(90)=(cos180sin180sin180cos180)=(1001)となる。

(3) 原点に関する対称移動

正解:③

やり方は2パターン

[パターン1:合成変換を利用]

原点に関する対称移動は、x 軸に関する対称移動 → y 軸に関する対称移動の2つの操作で実現できる。

よって、原点に関する対称変換行列は y 軸に関する対称移動と x 軸に関する対称移動の積で求められる*4

よって、(1001)(1001)=(1001)と求められる。

[パターン2:回転行列を利用]

原点に関する対称移動は、原点中心に 180 回転させた移動と同じ。

よって、R(180)=(cos180sin180sin180cos180)=(1001)と求められる。

(4) y=x に関する対称移動

正解:⑤

θ=45 のときの直線 y=(tanθ)x に関する対称移動と同じ。

なので、S(45)=(cos90sin90sin90cos90)=(0110)となる。

(5) y=x に関する対称移動

正解:⑥

θ=45 のときの直線 y=(tanθ)x に関する対称移動と同じ。

なので、S(45)=(cos90sin90sin90cos90)=(0110)となる。

解答2

(1)

原点中心に 30 回転させる表現行列 R(30) は、R(30)=(cos30sin30sin30cos30)=(32121232)となる。

よって座標 B は、(32121232)(64)=(33223+3)となる。

(2)

直線 l の傾きが2のときの x 軸と l がなす角度は、tanθ=2となるので、cos2θ=1tan2θ1+tan2θ=35sin2θ=2tanθ1+tan2θ=45となる*5

よって、(2)の変換を表す行列はS=(35454535)=15(3443)となる。

よって座標Cは、15(3443)(64)=15(236)となる。

6.さいごに

今回は行列の2次元・3次元の回転変換、対称変換についてまとめました。

練習問題では3次元の回転変換の問題を用意していませんが、もし必要な方は個人で3次元の回転変換の問題を解いていきましょう。

ちなみに、2次の直交行列は、回転を表す行列R(θ)=(cosθsinθsinθcosθ)もしくは、対称変換を表す行列 S(θ/2)=(cosθsinθsinθcosθ)2θθ に変えている)

の2パターンしかないことを頭に入れておくとどこかで役にたつかもしれません。

*1:わからなければ加法定理かなんかでバラしてゴリ押ししましょう。cos(θ+90)=cosθcos90sinθsin90=sinθsin(θ+90)=sinθcos90+sin90cosθ=cosθ 

*2:x 軸回りの回転では、 x 軸方向の成分が変化しないため、x 軸を考える必要はない。

*3:わからなければ加法定理かなんかでバラしてゴリ押ししましょう。cos(2θ90)=cos2θcos90+sin2θsin90=sin2θsin(2θ90)=sin2θcos90sin90cos2θ=cos2θ 

*4:計算順序に要注意(今回の場合はどちらの操作を先にやっても答えが変わらないため、計算順序によって答えがわからないが、基本的に操作順序を間違えると答えが変わる)。後に行う操作が左側、先に行う操作が右側になる。

*5:ここの計算は tanx2=t と置き換える積分を解くときに cosx, sinxt で置き換えるために必要な変形と一緒。詳しくはこちらから。

関連広告・スポンサードリンク

おすすめの記事