主題
Search

尤拉角


EulerAngles

根據 尤拉旋轉定理,任何 旋轉 都可以用三個 來描述。如果 旋轉旋轉矩陣 DCB 表示,那麼一般 旋轉 A 可以寫成

 A=BCD.
(1)

給出三個旋轉矩陣的三個角稱為尤拉角。尤拉角有幾種約定,取決於旋轉所繞的軸。將矩陣 A 寫成

 A=[a_(11) a_(12) a_(13); a_(21) a_(22) a_(23); a_(31) a_(32) a_(33)].
(2)

上面所示的所謂“x-約定”是最常見的定義。在這個約定中,由尤拉角 (phi,theta,psi) 給出的旋轉,其中

1. 第一次旋轉是繞 z旋轉角度 phi,使用 D

2. 第二次旋轉是繞先前的 x(現在是 x^')旋轉角度 theta in [0,pi]theta in [0,pi] 在 [0,pi] 範圍內,使用 C,以及

3. 第三次旋轉是繞先前的 z(現在是 z^')旋轉角度 psi,使用 B

然而,請注意,角的幾種符號約定是通用的。Goldstein (1980, pp. 145-148) 和 Landau and Lifschitz (1976) 使用 (phi,theta,psi),Tuma (1974) 說 (psi,theta,phi) 在航空工程中用於分析航天器(但聲稱 (phi,theta,psi) 用於分析陀螺運動),而 Bate et al. (1971) 使用 (Omega,i,omega)。Goldstein 評論說,歐洲大陸的作者通常使用 (psi,theta,phi),並警告說偶爾也會使用左手座標系 (Osgood 1937, Margenau and Murphy 1956-64)。Varshalovich (1988, pp. 21-23) 使用符號 (alpha,beta,gamma)(alpha^',beta^',gamma^') 來表示尤拉角,並給出了三種不同的角度約定,沒有一種對應於 x-約定。

在這裡,使用符號 (phi,theta,psi),這是一種可以在 6 之前的 Wolfram 語言 版本中使用的約定,如RotationMatrix3D[phi, theta, psi](可以在載入後執行Geometry`Rotations`)以及RotateShape[g, phi, theta, psi](可以在載入後執行Geometry`Shapes`)。在 x-約定中,分量旋轉由下式給出

D=[cosphi sinphi 0; -sinphi cosphi 0; 0 0 1]
(3)
C=[1 0 0; 0 costheta sintheta; 0 -sintheta costheta]
(4)
B=[cospsi sinpsi 0; -sinpsi cospsi 0; 0 0 1],
(5)

所以

a_(11)=cospsicosphi-costhetasinphisinpsi
(6)
a_(12)=cospsisinphi+costhetacosphisinpsi
(7)
a_(13)=sinpsisintheta
(8)
a_(21)=-sinpsicosphi-costhetasinphicospsi
(9)
a_(22)=-sinpsisinphi+costhetacosphicospsi
(10)
a_(23)=cospsisintheta
(11)
a_(31)=sinthetasinphi
(12)
a_(32)=-sinthetacosphi
(13)
a_(33)=costheta
(14)

為了獲得 角速度 omega 在本體座標系中的分量,請注意對於 矩陣

 A=[A_1 A_2 A_3],
(15)

這是成立的

[a_(11) a_(12) a_(13); a_(21) a_(22) a_(23); a_(31) a_(32) a_(33)][omega_x; omega_y; omega_z]=[a_(11)omega_x+a_(12)omega_y+a_(13)omega_z; a_(21)omega_x+a_(22)omega_y+a_(23)omega_z; a_(31)omega_x+a_(32)omega_y+a_(33)omega_z]
(16)
=A_1omega_x+A_2omega_y+A_3omega_z.
(17)

現在,omega_z 對應於繞 phi 軸的旋轉,所以檢視 Aomegaomega_z 分量,

omega_phi=A_3omega_z
(18)
=[sinpsisintheta; cospsisintheta; costheta]phi^..
(19)

交點線對應於繞 xi 軸旋轉 theta,所以檢視 Bomegaomega_xi 分量,

omega_theta=B_1omega_xi
(20)
=B_1theta^.
(21)
=[cospsi; -sinpsi; 0]theta^..
(22)

類似地,要找到繞剩餘軸旋轉 psi,請檢視 Bomegaomega_psi 分量,

omega_psi=B_3omega_psi
(23)
=B_3psi^.
(24)
=[0; 0; 1]psi^..
(25)

結合這些部分得到

 omega=[sinpsisinthetaphi^.+cospsitheta^.; cospsisinthetaphi^.-sinpsitheta^.; costhetaphi^.+psi^..]
(26)

更多細節,請參見 Goldstein (1980, p. 176) 和 Landau and Lifschitz (1976, p. 111)。

x-約定尤拉角由 Cayley-Klein 引數 表示為

phi=-2iln[+/-(alpha^(1/2)gamma^(1/4))/(beta^(1/4)(1+betagamma)^(1/4))],-2iln[+/-(ialpha^(1/2)gamma^(1/4))/(beta^(1/4)(1+betagamma)^(1/4))]
(27)
psi=-2iln[+/-(alpha^(1/2)beta^(1/4))/(gamma^(1/4)(1+betagamma)^(1/4))],-2iln[+/-(ialpha^(1/2)beta^(1/4))/(gamma^(1/4)(1+betagamma)^(1/4))]
(28)
theta=+/-2cos^(-1)(+/-sqrt(1+betagamma)).
(29)

在“y-約定”中,

phi_x=phi_y+1/2pi
(30)
psi_x=psi_y-1/2pi.
(31)

因此,

sinphi_x=cosphi_y
(32)
cosphi_x=-sinphi_y
(33)
sinpsi_x=-cospsi_y
(34)
cospsi_x=sinpsi_y,
(35)

給出旋轉矩陣

D=[-sinphi cosphi 0; -cosphi -sinphi 0; 0 0 1]
(36)
C=[1 0 0; 0 costheta sintheta; 0 -sintheta costheta]
(37)
B=[sinpsi -cospsi 0; cospsi sinpsi 0; 0 0 1]
(38)

並且 A 由下式給出

a_(11)=-sinpsisinphi+costhetacosphicospsi
(39)
a_(12)=sinpsicosphi+costhetasinphicospsi
(40)
a_(13)=-cospsisintheta
(41)
a_(21)=-cospsisinphi-costhetacosphisinpsi
(42)
a_(22)=cospsicosphi-costhetasinphisinpsi
(43)
a_(23)=sinpsisintheta
(44)
a_(31)=sinthetacosphi
(45)
a_(32)=sinthetasinphi
(46)
a_(33)=costheta.
(47)

在 “xyz (俯仰-滾轉-偏航) 約定”中,theta 是俯仰角,psi 是滾轉角,phi 是偏航角。

D=[cosphi sinphi 0; -sinphi cosphi 0; 0 0 1]
(48)
C=[costheta 0 -sintheta; 0 1 0; sintheta 0 costheta]
(49)
B=[1 0 0; 0 cospsi sinpsi; 0 -sinpsi cospsi]
(50)

並且 A 由下式給出

a_(11)=costhetacosphi
(51)
a_(12)=costhetasinphi
(52)
a_(13)=-sintheta
(53)
a_(21)=sinpsisinthetacosphi-cospsisinphi
(54)
a_(22)=sinpsisinthetasinphi+cospsicosphi
(55)
a_(23)=costhetasinpsi
(56)
a_(31)=cospsisinthetacosphi+sinpsisinphi
(57)
a_(32)=cospsisinthetasinphi-sinpsicosphi
(58)
a_(33)=costhetacospsi.
(59)

有時使用一組引數代替角,即 尤拉引數 e_0e_1e_2e_3,定義為

e_0=cos(phi/2)
(60)
e=[e_1; e_2; e_3]=n^^sin(phi/2).
(61)

使用 尤拉引數(它們是四元數),可以由下式描述任意旋轉矩陣

a_(11)=e_0^2+e_1^2-e_2^2-e_3^2
(62)
a_(12)=2(e_1e_2+e_0e_3)
(63)
a_(13)=2(e_1e_3-e_0e_2)
(64)
a_(21)=2(e_1e_2-e_0e_3)
(65)
a_(22)=e_0^2-e_1^2+e_2^2-e_3^2
(66)
a_(23)=2(e_2e_3+e_0e_1)
(67)
a_(31)=2(e_1e_3+e_0e_2)
(68)
a_(32)=2(e_2e_3-e_0e_1)
(69)
a_(33)=e_0^2-e_1^2-e_2^2+e_3^2
(70)

(Goldstein 1980, p. 153)。

如果已知兩組 n 個點 x_ix_i^' 的座標,一組相對於另一組旋轉,那麼可以使用 最小二乘擬合 以直接的方式獲得尤拉旋轉矩陣。將點寫成向量陣列,因此

 [x_1^' ... x_n^']=A[x_1 ... x_n].
(71)

將向量陣列寫成矩陣得到

 X^'=AX
(72)
 X^'X^(T)=AXX^(T),
(73)

並且求解 A 得到

 A=X^'X^(T)(XX^(T))^(-1).
(74)

然而,我們想要的是角 thetaphipsi,而不是它們在矩陣 A 中包含的組合。因此,將 3×3 矩陣

 A=[f_1(theta,phi,psi) f_2(theta,phi,psi) f_3(theta,phi,psi); f_4(theta,phi,psi) f_5(theta,phi,psi) f_6(theta,phi,psi); f_7(theta,phi,psi) f_7(theta,phi,psi) f_9(theta,phi,psi)]
(75)

作為 1×9 向量

 f=[f_1(theta,phi,psi); |; f_9(theta,phi,psi)].
(76)

現在設定矩陣

 [(partialf_1)/(partialtheta)|_(theta_i,phi_i,psi_i) (partialf_1)/(partialphi)|_(theta_i,phi_i,psi_i) (partialf_1)/(partialpsi)|_(theta_i,phi_i,psi_i); | | |; (partialf_9)/(partialtheta)|_(theta_i,phi_i,psi_i) (partialf_9)/(partialphi)|_(theta_i,phi_i,psi_i) (partialf_9)/(partialpsi)|_(theta_i,phi_i,psi_i)][dtheta; dphi; dpsi]=df.
(77)

然後使用 非線性最小二乘擬合 給出收斂到 (theta,phi,psi) 的解。


另請參閱

Cayley-Klein 引數, 尤拉引數, 尤拉旋轉定理, 無窮小旋轉, 四元數, 旋轉, 旋轉公式, 旋轉矩陣

使用 探索

參考文獻

Arfken, G. 物理學家數學方法,第 3 版。 Orlando, FL: Academic Press, pp. 198-200, 1985.Bate, R. R.; Mueller, D. D.; and White, J. E. 天體動力學基礎。 New York: Dover, 1971.Goldstein, H. "The Euler Angles" and "Euler Angles in Alternate Conventions." §4-4 and Appendix B in 經典力學,第 2 版。 Reading, MA: Addison-Wesley, pp. 143-148 and 606-610, 1980.Kraus, M. "LiveGraphics3D Example: Euler Angles." http://wwwvis.informatik.uni-stuttgart.de/~kraus/LiveGraphics3D/examples/Euler.html.Landau, L. D. and Lifschitz, E. M. 力學,第 3 版。 Oxford, England: Pergamon Press, 1976.Margenau, H. and Murphy, G. M. 物理和化學數學,2 卷。 Princeton, NJ: Van Nostrand, 1956-64.Osgood, W. F. 力學。 New York: Macmillan, 1937.Tuma, J. J. 動力學。 New York: Quantum Publishers, 1974.Varshalovich, D. A.; Moskalev, A. N.; and Khersonskii, V. K. "Description of Rotation in Terms of the Euler Angles." §1.4.1 in 角動量量子理論。 Singapore: World Scientific, pp. 21-23, 1988.

在 上被引用

尤拉角

引用為

Weisstein, Eric W. “尤拉角。” 來自 —— 資源。 https://mathworld.tw/EulerAngles.html

主題分類