■3次元の回転(その13)

複素数球面に作用する1次分数変換の性質を用いて、ケーリーの公式を導き、以下の式を与える

複素数平面CとそのコピーC'から{0}除いた集合上の点zとz'を

z'=-1/z

という関係でもって同一視する。こうして作られた集合を複素数球面(リーマン球面)とよぶ=CU{∞}

x'=-x/(x^2+y^2)

y'=y/(x^2+y^2)

x=-x'/(x'^2+y'^2)

y=y'/(x'^2+y'^2)

一次分数変換

w=φ(z)=(αz+β)/(γz+δ)

は複素数球面(リーマン球面)上の円を円に写す等角写像である。(これはz平面上の円か直線を意味するものとする)

極射影ではP(v1,v2,v3),z=x+yiとおくと

x=v1/(1+v3), y=v2/(1+v3)

z=v1/(1+v3)+iv2/(1+v3)

x^2+y^2=(v1^2+v2^2)/(1+v3)^2=(1-v3^2)/(1+v3)^2=(1-v3)/(1+v3)

v1=2x/(1+x^2+y^2)

v2=2y/(1+x^2+y^2)

v3=(1-x^2-y^2)/(1+x^2+y^2)

===================================

【1】ケーリーの公式

k=(l,m,n),l^2+m^2+n^2=1を軸とする角θの回転は1次分数変換

φ(z)={(cosθ/2+insinθ/2)z+(msinθ/2-ilsinθ/2)}/{(-msinθ/2-ilsinθ/2)z+(cosθ/2-insinθ/2)}

と書かれる

===================================

a=cosθ/2,b=nsinθ/2

c=-msinθ/2,d=-lsinθ/2とおくと、a^2+b^2+c^2+d^2=1となり、

ケーリーの公式は

φ(z)={(a+bi)z+(-c+di)}/{(c+di)z+(a-bi)}

と書かれる。

U=[(a+bi),(-c+di)]

[(-c+di),(a-bi)]

det(U)=a^2+b^2+c^2+d^2=1

は2次特殊ユニタリー行列と呼ばれる。

U=[0,1]

[-1,0]

も2次特殊ユニタリー行列である

===================================

a,b,c,dを用いて、3次直交行列を書くと

A(1,1)=a^2-b^2-c^2+d^2

A(1,2)=2(cd-ab)

A(1,3)=-2(ac+bd)

A(2,1)=2(ab+cd)

A(2,2)=a^2-b^2+c^2-d^2

A(2,3)=2(ad-bc)

A(3,1)=2(ac-bd)

A(3,2)=-2(ad+bc)

A(3,3)=a^2+b^2-c^2-d^2

l,m,n,θを用いて、3次直交行列を書くと最終的に次が得られる。

A(1,1)=l^2+(1-l^2)cosθ

A(1,2)=lm-lmcosθ-nsinθ

A(1,3)=ln-lncosθ+msinθ

A(2,1)=lm-lmcosθ+nsinθ

A(2,2)=m^2+(1-m^2)cosθ

A(2,3)=mn-mncosθ+msinθ

A(3,1)=ln-lncosθ

A(3,2)=mn-mncosθ-lsinθ

A(3,3)=n^2+(1-n^2)cosθ

===================================

【3】単位四元数と3次元の回転

 x,y,z軸の周りの回転では使いにくい.そこで,単位ベクトル

  n=(α,β,γ)

を回転軸とし,その周りに正の回転方向にθだけ回転する回転行列はα,β,γは方向余弦で,α^2+β^2+γ^2=1を満たすものとして

  R(1,1)=α^2(1-cosθ)+cosθ

  R(2,2)=β^2(1-cosθ)+cosθ

  R(3,3)=γ^2(1-cosθ)+cosθ

  R(1,2)=αβ(1-cosθ)+γsinθ

  R(2,1)=αβ(1-cosθ)-γsinθ

  R(1,3)=αγ(1-cosθ)-βsinθ

  R(3,1)=αγ(1-cosθ)+βsinθ

  R(2,3)=βγ(1-cosθ)+αsinθ

  R(3,2)=βγ(1-cosθ)-αsinθ

で表される.

===================================

これが先にあって、後知恵で直交行列Kを見つけると

k(1,1)=α^2/(1+γ)-1

k(1,2)=αβ/(1+γ)-1

k(1,3)=α

k(2,1)=αβ/(1+γ)-1

k(2,2)=β^2/(1+γ)-1

k(2,3)=β

k(3,1)=α

k(3,2)=β

k(3,3)=γ

===================================

     [ cosθ,sinθ,0]

  R3 =[−sinθ,cosθ,0]

     [   0,   0, 1]

として,R=KR3K^-1のようなものを考えている.

===================================