主题
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 参数, 欧拉参数, 欧拉旋转定理, 无穷小旋转, 四元数, 旋转, 旋转公式, 旋转矩阵

使用 Wolfram|Alpha 探索

参考文献

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.

在 Wolfram|Alpha 上被引用

欧拉角

引用为

Weisstein, Eric W. “欧拉角。” 来自 MathWorld——Wolfram Web 资源。 https://mathworld.net.cn/EulerAngles.html

主题分类