Euler angles: Difference between revisions
imported>Paul Wormer No edit summary |
mNo edit summary |
||
(19 intermediate revisions by 2 users not shown) | |||
Line 1: | Line 1: | ||
{{subpages}} | {{subpages}} | ||
{{Image|Euler angles.png|right|550px|Figure 1. Euler angles. From left to right: initial configuration, after rotation over angle α, after rotation over angle β, and after rotation over angle γ.}} | {{Image|Euler angles.png|right|550px|Figure 1. Euler angles. From left to right: initial configuration, after rotation over angle α, after rotation over angle β, and after rotation over angle γ.}} | ||
Different conventions are in use: a rotation can be active (the object is rotated, the system of axes is fixed in space), or passive (the object is fixed in space, the axes are rotated). | In [[physics]], [[mathematics]], and [[engineering]], '''[[Euler angles]]''' are three rotation angles, often denoted by 0 ≤ α ≤ 2π, 0 ≤ β ≤ π, and 0 ≤ γ ≤ 2π, although the notation φ, θ, ψ is also common. Any rotation of a 3-dimensional object can be performed by three consecutive rotations over the three Euler angles. | ||
Different conventions are in use: a rotation can be active (the object is rotated, the system of axes is fixed in space), or passive (the object is fixed in space, the axes are rotated). | |||
Also the choice of rotation axes may vary; an active convention common in quantum mechanical applications is the ''z-y′-z′'' convention. Attach a system of [[Cartesian coordinates|Cartesian coordinate]] axes to the body that is to be rotated (the coordinate frame is fixed to the body and is rotated simultaneously with it); in the figure the body-fixed frame is shown in red and labeled by lowercase letters. First rotate around ''z'', then around the new body-fixed ''y''-axis, ''y′'', and finally around ''z′''. Another convention often used is the ''z-x′-z′'' convention, where instead of over the new ''y''-axis the second rotation is over the new ''x''-axis. Also the ''z-y-x'' convention is used (and will be discussed [[Euler_angles#Theorem|below]]). | |||
The [[right-hand screw rule]] is practically always followed: the rotation axis is a directed line and a positive rotation is as a cork screw driven into the positive direction of the axis. In older literature ''left-handed'' Cartesian coordinate frames appear sometimes, but in modern literature right-handed frames are used exclusively. | |||
The | Euler angles are used in many different branches of physics and engineering. The present article is written from the point of view of molecular physics, where the objects to be rotated are [[molecules]] and applications are often of quantum mechanical nature. | ||
The [[ | The angles are named after the 18th century mathematician [[Leonhard Euler]] who introduced in 1765 two of the three for an axially symmetric body where the third angle, γ, does not play a role.<ref>[http://www.17centurymaths.com/contents/euler/mechvol3/tmvol1ch9tr.pdf Translation by Ian Bruce] of L. Euler, ''Theoria Motus Corporum Solidorum Seu Rigidorum'' (Theory of the motion of solid or rigid bodies), Rostock (1765), pdf page 11. Later Euler returned to the angles and gave an alternative derivation, see [http://www.17centurymaths.com/contents/euler/e478tr.pdf. Translation by Johan Sten] of ''Formulae generales pro translatione quacunque corporum rigidorum'' (General formulas for the translation of arbitrary rigid bodies), Novi Commentarii academiae scientiarum Petropolitanae, vol. '''20''', (1776), pp. 189-207</ref> | ||
==Geometric discussion== | ==Geometric discussion== | ||
In Figure 1 the space-fixed (laboratory) axes are labeled by capital X, Y, and Z and are shown in black. The body to be rotated is not shown, but a system of axes fixed to it is shown in red. One may use any convenient orthonormal frame as a body-fixed frame. Often the body-fixed axes are [[principal axes]], that means that they are [[eigenvectors]] of the [[inertia tensor]] of the body. Also symmetry axes, when present, may be used. When the body has symmetry axes, the principal axes often coincide with these. | In Figure 1 the space-fixed (laboratory) axes are labeled by capital X, Y, and Z and are shown in black. The body to be rotated is not shown, but a system of axes fixed to it is shown in red. One may use any convenient orthonormal frame as a body-fixed frame. Often the body-fixed axes are [[principal axes]], that means that they are [[eigenvectors]] of the [[inertia tensor]] of the body. Also symmetry axes, when present, may be used. When the body has symmetry axes, the principal axes often coincide with these. | ||
{{Image|Rot z Rot y.png|right|300px|Figure 2. Rotation of '''r''' to '''r'''′. On the left around ''z''-axis over α (φ increases), on the right around ''y''-axis over β (φ decreases). Both rotation axes point toward the reader.}} | {{Image|Rot z Rot y.png|right|300px|Figure 2. Rotation of '''r''' to '''r'''′. On the left around ''z''-axis over α (φ increases), on the right around ''y''-axis over β (φ decreases). Both rotation axes point toward the reader.}} | ||
The ''z-y′-z′'' convention will be followed. Initially, the two frames coincide, and the path to a final arbitrary orientation of the body—and its frame—is depicted on Figure 1. The first rotation is around the ''z''-axis, which coincides with the Z-axis. The ''x''- and ''y''-axis move in a plane perpendicular to the ''z''-axis over an angle α. The second rotation is in a plane through the origin perpendicular to the ''y′''-axis. The angle is β. The present convention has the | The ''z-y′-z′'' convention will be followed. Initially, the two frames coincide, and the path to a final arbitrary orientation of the body—and its frame—is depicted on Figure 1. The first rotation is around the ''z''-axis, which coincides with the Z-axis. The ''x''- and ''y''-axis move in a plane perpendicular to the ''z''-axis over an angle α. The second rotation is in a plane through the origin perpendicular to the ''y′''-axis. The angle is β. The present convention has the practical advantage that the ''z′''-axis has the usual [[spherical polar coordinates]] α ≡ φ (longitude angle) and β ≡ θ (colatitude angle) with respect to the space-fixed frame.<ref>In the ''z''-''x''′-''z''′ convention the first two Euler angles are not equal to spherical polar angles, in consequence the (''m'', ''m''′) [[Wigner D-matrix]]-element carries the complex phase exp[''i''π(''m''−''m''′)/2]. This phase is absent in the ''z''-''y''′-''z''′ convention </ref> The final rotation is in a plane perpendicular to the ''z′''-axis over an angle γ. From geometric considerations follows that any orientation of the body-fixed frame in space may be obtained. | ||
Write <math>\mathbf{R}(\varphi, \hat{n})</math> for the [[rotation matrix]] that describes a rotation around the unit vector <math>\hat{n}</math> over an angle <math>\varphi</math>. | Write <math>\mathbf{R}(\varphi, \hat{n})</math> for the [[rotation matrix]] that describes a rotation around the unit vector <math>\hat{n}</math> over an angle <math>\varphi</math>. | ||
Line 47: | Line 52: | ||
\end{align} | \end{align} | ||
</math> | </math> | ||
Note that the third column contains the Cartesian coordinates with respect to the space-fixed frame of <math>\hat{e}_{z''}</math> expressed in sines and cosines of spherical polar angles. The first and second column contain by definition expressions for the Cartesian coordinates of <math>\hat{e}_{x'''}</math> and <math>\hat{e}_{y''}</math>, respectively, but evidently these are ''not'' solely in terms of spherical polar angles, γ also enters. | |||
Before proving the first equality in the above equation (reversal of order), we derive the matrix for a rotation around the ''z''-axis, see the left drawing in Figure 2. The rotated vector has components | Before proving the first equality in the above equation (reversal of order), we derive the matrix for a rotation around the ''z''-axis, see the left drawing in Figure 2. The rotated vector has components | ||
Line 134: | Line 140: | ||
==Algebraic treatment== | ==Algebraic treatment== | ||
In the proof that any rotation can be written as three consecutive rotations, an appeal was made to the geometric insight of the reader. The same result can be proved more rigorously by algebraic means. To that end the notation is somewhat shortened: | In the proof that ''any'' rotation can be written as three consecutive rotations, an appeal was made to the geometric insight of the reader. The same result can be proved more rigorously by algebraic means. To that end the notation is somewhat shortened: | ||
:<math> | :<math> | ||
\mathbf{R}_z (\omega) \equiv \mathbf{R}(\omega, \hat{e}_{z}) \quad\hbox{and}\quad | \mathbf{R}_z (\omega) \equiv \mathbf{R}(\omega, \hat{e}_{z}) \quad\hbox{and}\quad | ||
Line 151: | Line 157: | ||
the ''Euler z-y-z parametrization''. | the ''Euler z-y-z parametrization''. | ||
====Proof==== | ====Proof==== | ||
First the Euler ''z-y-x''-parametrization will be proved by an algorithm for the factorization of a given matrix '''R''' ≡ ('''r'''<sub>1</sub>, '''r'''<sub>2</sub>, '''r'''<sub>3</sub>). | First the Euler ''z-y-x''-parametrization will be proved by an algorithm for the factorization of a given matrix '''R''' ≡ ('''r'''<sub>1</sub>, '''r'''<sub>2</sub>, '''r'''<sub>3</sub>). Second the ''z-y-z'' parametrization will be proved; this parametrization is—as shown above—equivalent to the ''z′-y′-z'' parametrization with angles in reverse order. | ||
'':A Fortran subroutine based on the algorithm is given on the [[Euler_angles/Code|code page]].'' | |||
To prove the ''z-y-x'' parametrization we consider the matrix product | |||
:<math> | :<math> | ||
\mathbf{R}_z (\omega_3 ) \, \mathbf{R}_y (\omega_2 ) = | \mathbf{R}_z (\omega_3 ) \, \mathbf{R}_y (\omega_2 ) = | ||
Line 230: | Line 239: | ||
\end{align} | \end{align} | ||
</math> | </math> | ||
This | This gives the required ''z-y-x'' factorization of the arbitrary proper orthogonal matrix '''R'''. | ||
The Euler ''z-y-z'' parametrization is obtained by a small modification of the previous proof. | The proof of the Euler ''z-y-z'' parametrization is obtained by a small modification of the previous proof. We start by retrieving the spherical polar coordinates <math>\omega_2\;</math> and | ||
<math>\omega_2\;</math> and | <math>\omega_3\;</math> of the unit vector <math>\mathbf{r}_3 = \mathbf{a}_3 </math>, the third column [the rightmost multiplication by '''R'''<sub>''z''</sub>(ω<sub>1</sub>) does not affect '''r'''<sub>3</sub>]. Then consider | ||
<math>\omega_3\;</math> | |||
:<math> | :<math> | ||
( \mathbf{r}_1, \; \mathbf{r}_2 ) = (\mathbf{a}_1, \; \mathbf{a}_2 ) | ( \mathbf{r}_1, \; \mathbf{r}_2 ) = (\mathbf{a}_1, \; \mathbf{a}_2 ) | ||
Line 255: | Line 262: | ||
\mathbf{R}_y (\omega_2 ) \, \mathbf{R}_z (\omega_1 ) \; , | \mathbf{R}_y (\omega_2 ) \, \mathbf{R}_z (\omega_1 ) \; , | ||
</math> | </math> | ||
which proves the Euler ''z-y-z'' parametrization. Clearly, this factorization is equal to the one given | which proves the Euler ''z-y-z'' parametrization. Clearly, this factorization is equal to the one given in the previous section, with | ||
:<math> | :<math> | ||
\omega_3 | \omega_3 \equiv \alpha,\quad \omega_2 \equiv \beta, \quad \omega_1 \equiv \gamma. | ||
</math> | </math> | ||
==Note== | |||
<references />[[Category:Suggestion Bot Tag]] |
Latest revision as of 06:00, 14 August 2024
In physics, mathematics, and engineering, Euler angles are three rotation angles, often denoted by 0 ≤ α ≤ 2π, 0 ≤ β ≤ π, and 0 ≤ γ ≤ 2π, although the notation φ, θ, ψ is also common. Any rotation of a 3-dimensional object can be performed by three consecutive rotations over the three Euler angles.
Different conventions are in use: a rotation can be active (the object is rotated, the system of axes is fixed in space), or passive (the object is fixed in space, the axes are rotated).
Also the choice of rotation axes may vary; an active convention common in quantum mechanical applications is the z-y′-z′ convention. Attach a system of Cartesian coordinate axes to the body that is to be rotated (the coordinate frame is fixed to the body and is rotated simultaneously with it); in the figure the body-fixed frame is shown in red and labeled by lowercase letters. First rotate around z, then around the new body-fixed y-axis, y′, and finally around z′. Another convention often used is the z-x′-z′ convention, where instead of over the new y-axis the second rotation is over the new x-axis. Also the z-y-x convention is used (and will be discussed below).
The right-hand screw rule is practically always followed: the rotation axis is a directed line and a positive rotation is as a cork screw driven into the positive direction of the axis. In older literature left-handed Cartesian coordinate frames appear sometimes, but in modern literature right-handed frames are used exclusively.
Euler angles are used in many different branches of physics and engineering. The present article is written from the point of view of molecular physics, where the objects to be rotated are molecules and applications are often of quantum mechanical nature.
The angles are named after the 18th century mathematician Leonhard Euler who introduced in 1765 two of the three for an axially symmetric body where the third angle, γ, does not play a role.[1]
Geometric discussion
In Figure 1 the space-fixed (laboratory) axes are labeled by capital X, Y, and Z and are shown in black. The body to be rotated is not shown, but a system of axes fixed to it is shown in red. One may use any convenient orthonormal frame as a body-fixed frame. Often the body-fixed axes are principal axes, that means that they are eigenvectors of the inertia tensor of the body. Also symmetry axes, when present, may be used. When the body has symmetry axes, the principal axes often coincide with these.
The z-y′-z′ convention will be followed. Initially, the two frames coincide, and the path to a final arbitrary orientation of the body—and its frame—is depicted on Figure 1. The first rotation is around the z-axis, which coincides with the Z-axis. The x- and y-axis move in a plane perpendicular to the z-axis over an angle α. The second rotation is in a plane through the origin perpendicular to the y′-axis. The angle is β. The present convention has the practical advantage that the z′-axis has the usual spherical polar coordinates α ≡ φ (longitude angle) and β ≡ θ (colatitude angle) with respect to the space-fixed frame.[2] The final rotation is in a plane perpendicular to the z′-axis over an angle γ. From geometric considerations follows that any orientation of the body-fixed frame in space may be obtained.
Write Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{R}(\varphi, \hat{n})} for the rotation matrix that describes a rotation around the unit vector Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \hat{n}} over an angle Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \varphi} . Clearly the three consecutive Euler rotations correspond to rotations around the unit vectors along the body-fixed axes z, y′, and z′ over angles α β, and γ, respectively. Because a matrix acts on a column vector to its right, the order in the matrix product is as in the leftmost term in the following equation. It will be shown that the corresponding matrix product can be written in reverse order (but around fixed, unprimed, axes z, y, z), that is,
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{align} \mathbf{R}(\gamma, \hat{e}_{z'}) \mathbf{R}(\beta, \hat{e}_{y'}) \mathbf{R}(\alpha, \hat{e}_{z})&= \mathbf{R}(\alpha, \hat{e}_{z }) \mathbf{R}(\beta, \hat{e}_{y }) \mathbf{R}(\gamma, \hat{e}_{z}) \\ &= \begin{pmatrix} \cos\alpha & -\sin \alpha & 0 \\ \sin\alpha & \cos \alpha & 0 \\ 0 & 0 & 1 \\ \end{pmatrix} \begin{pmatrix} \cos\beta & 0 & \sin\beta \\ 0 & 1 & 0 \\ -\sin\beta & 0 & \cos\beta \\ \end{pmatrix} \begin{pmatrix} \cos\gamma & -\sin\gamma & 0 \\ \sin\gamma & \cos\gamma & 0 \\ 0 & 0 & 1 \\ \end{pmatrix} \\ &= \begin{pmatrix} \cos\alpha\cos\beta\cos\gamma-\sin\alpha\cos\gamma \;&\; -\cos\alpha\cos\beta\sin\gamma-\sin\alpha\cos\gamma \;&\;\cos\alpha\sin\beta \\ \sin\alpha\cos\beta\cos\gamma+\cos\alpha\sin\gamma \;&\; -\sin\alpha\cos\beta\sin\gamma+\cos\alpha\cos\gamma \;&\; \sin\alpha\sin\beta \\ -\sin\beta\cos\gamma \;&\; \sin\beta\sin\gamma \;&\; \cos\beta \end{pmatrix} \end{align} }
Note that the third column contains the Cartesian coordinates with respect to the space-fixed frame of Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \hat{e}_{z''}} expressed in sines and cosines of spherical polar angles. The first and second column contain by definition expressions for the Cartesian coordinates of Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \hat{e}_{x'''}} and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \hat{e}_{y''}} , respectively, but evidently these are not solely in terms of spherical polar angles, γ also enters.
Before proving the first equality in the above equation (reversal of order), we derive the matrix for a rotation around the z-axis, see the left drawing in Figure 2. The rotated vector has components
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{pmatrix} \cos(\alpha + \phi) \\ \sin(\alpha + \phi) \\ 0 \end{pmatrix} = \begin{pmatrix} \cos\alpha\cos\phi -\sin\alpha\sin\phi \\ \sin\alpha\cos\phi + \cos\alpha\sin\phi \\ 0 \end{pmatrix} = \begin{pmatrix} \cos\alpha & -\sin\alpha & 0 \\ \sin\alpha & \cos\alpha & 0\\ 0 & 0 & 1 \\ \end{pmatrix} \begin{pmatrix} \cos\phi \\ \sin\phi \\ 0 \\ \end{pmatrix} \equiv \mathbf{R}(\alpha, \hat{e}_z) \begin{pmatrix} \cos\phi \\ \sin\phi \\ 0 \end{pmatrix} . }
We used here the relations well-known from trigonometry for the sine and cosine of a sum angle. The derivation of the matrix for a rotation around the y-axis proceeds along the same lines. Note, however, that the angle of a vector with the x-axis decreases by a rotation around the positive y-axis (see right-hand drawing in Figure 2).
To prove the first equality (reversal of the order in the angles), a property of rotation matrices is used. A rotation (orthogonal 3×3) matrix A, transforming a rotation axis, gives rise to the following similarity equation,
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{R}(\varphi, \mathbf{A}\hat{n}) = \mathbf{A}\mathbf{R}(\varphi, \hat{n})\mathbf{A}^{\mathrm{T}}, }
where the superscript T indicates the transpose of the matrix. For rotation matrices the transposed matrix is equal to the inverse of the matrix. From this similarity relation follows that
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{R}(\gamma, \hat{e}_{z'}) = \mathbf{R}(\beta, \hat{e}_{y'}) \mathbf{R}(\gamma, \hat{e}_{z}) \mathbf{R}(\beta, \hat{e}_{y'})^{\mathrm{T}}, }
so that
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{R}(\gamma, \hat{e}_{z'}) \mathbf{R}(\beta, \hat{e}_{y'}) \mathbf{R}(\alpha, \hat{e}_{z})= \mathbf{R}(\beta, \hat{e}_{y'}) \mathbf{R}(\gamma, \hat{e}_{z}) \mathbf{R}(\beta, \hat{e}_{y'})^{\mathrm{T}} \; \mathbf{R}(\beta, \hat{e}_{y'}) \mathbf{R}(\alpha, \hat{e}_{z}) = \mathbf{R}(\beta, \hat{e}_{y'}) \mathbf{R}(\gamma, \hat{e}_{z}) \mathbf{R}(\alpha, \hat{e}_{z}). }
Also
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{R}(\beta, \hat{e}_{y'}) = \mathbf{R}(\alpha, \hat{e}_{z})\mathbf{R}(\beta, \hat{e}_{y})\mathbf{R}(\alpha, \hat{e}_{z})^{\mathrm{T}}, }
so that
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{R}(\beta, \hat{e}_{y'}) \mathbf{R}(\gamma, \hat{e}_{z}) \mathbf{R}(\alpha, \hat{e}_{z}) = \mathbf{R}(\alpha, \hat{e}_{z})\mathbf{R}(\beta, \hat{e}_{y}) \mathbf{R}(\alpha, \hat{e}_{z})^{\mathrm{T}} \; \mathbf{R}(\gamma, \hat{e}_{z}) \mathbf{R}(\alpha, \hat{e}_{z}) = \mathbf{R}(\alpha, \hat{e}_{z})\mathbf{R}(\beta, \hat{e}_{y}) \mathbf{R}(\gamma, \hat{e}_{z}), }
where it is used that rotations around the same axis commute, that is,
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{R}(\alpha, \hat{e}_{z})^{\mathrm{T}} \mathbf{R}(\gamma, \hat{e}_{z}) \mathbf{R}(\alpha, \hat{e}_{z}) = \mathbf{R}(\gamma, \hat{e}_{z}) \mathbf{R}(\alpha, \hat{e}_{z})^{\mathrm{T}} \mathbf{R}(\alpha, \hat{e}_{z}) = \mathbf{R}(\gamma, \hat{e}_{z}) }
and the required result is proved.
Algebraic treatment
In the proof that any rotation can be written as three consecutive rotations, an appeal was made to the geometric insight of the reader. The same result can be proved more rigorously by algebraic means. To that end the notation is somewhat shortened:
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{R}_z (\omega) \equiv \mathbf{R}(\omega, \hat{e}_{z}) \quad\hbox{and}\quad \mathbf{R}_y (\omega )\equiv \mathbf{R}(\omega, \hat{e}_{y}). }
Theorem
A proper rotation matrix R can be factorized thus
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{R} = \mathbf{R}_z (\omega_3 ) \; \mathbf{R}_y (\omega_2 ) \; \mathbf{R}_x (\omega_1 ) }
which is referred to as the Euler z-y-x parametrization, or also as
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{R} = \mathbf{R}_z (\alpha) \; \mathbf{R}_y (\beta ) \; \mathbf{R}_z (\gamma ) \quad }
the Euler z-y-z parametrization.
Proof
First the Euler z-y-x-parametrization will be proved by an algorithm for the factorization of a given matrix R ≡ (r1, r2, r3). Second the z-y-z parametrization will be proved; this parametrization is—as shown above—equivalent to the z′-y′-z parametrization with angles in reverse order.
:A Fortran subroutine based on the algorithm is given on the code page.
To prove the z-y-x parametrization we consider the matrix product
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{R}_z (\omega_3 ) \, \mathbf{R}_y (\omega_2 ) = \begin{pmatrix} \cos \omega_3 \cos \omega_2 & -\sin \omega_3 & \cos \omega_3 \sin \omega_2 \\ \sin \omega_3 \cos \omega_2 & \cos \omega_3 & \sin \omega_3 \sin \omega_2 \\ -\sin \omega_2 & 0 & \cos \omega_2 \\ \end{pmatrix} \equiv (\mathbf{a}_1 , \mathbf{a}_2 , \mathbf{a}_3 ) . }
The columns of the matrix product are for ease of reference designated by a1, a2, and a3. Note that the multiplication by
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{R}_x(\omega_1) \equiv \begin{pmatrix} 1 & 0 & 0 \\ 0 & \cos\omega_1 & -\sin\omega_1 \\ 0 & \sin\omega_1 & \cos\omega_1 \end{pmatrix} }
on the right does not affect the first column, so that a1 = r1 (the first column of R). Solve Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega_2\;} and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega_3\;} from the first column of R (which is known),
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{a}_1 = \begin{pmatrix} \cos \omega_3 \; \cos \omega_2 \\ \sin \omega_3 \; \cos \omega_2 \\ -\sin \omega_2 \\ \end{pmatrix} = \begin{pmatrix} R_{11} \\ R_{21} \\ R_{31} \\ \end{pmatrix} \equiv \mathbf{r}_1 . }
This is possible. First solve Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega_2\;} for Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle -\pi/2 \leq \omega_2 \leq \pi/2} from
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sin \omega_2 = - R_{31}. \, }
Then solve Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega_3\;} for Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle 0 \leq \omega_3 \leq 2 \pi} from the two equations:
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \cos \omega_3 = {R_{11} \over \cos \omega_2},\qquad \sin \omega_3 = {R_{21} \over \cos \omega_2} . }
The angles Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega_2\;} and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega_3\;} determine fully the vectors a2 and a3.
Since a1, a2 and a3 are the columns of a proper rotation matrix they form an orthonormal right-handed system. The plane spanned by a2 and a3 is orthogonal to Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{a}_1 \equiv \mathbf{r}_1} and hence the plane contains Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{r}_2} and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{r}_3} . Thus the latter two vectors are a linear combination of the first two,
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle ( \mathbf{r}_2 , \mathbf{r}_3 ) = (\mathbf{a}_2 , \mathbf{a}_3 ) \begin{pmatrix} \cos \omega_1 & -\sin \omega_1 \\ \sin \omega_1 & \cos \omega_1 \\ \end{pmatrix} . }
Since Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{r}_2,\; \mathbf{a}_2,\; \mathbf{a}_3} are known orthonormal vectors, we can compute
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{align} \mathbf{a}_2 \cdot \mathbf{r}_2 =& \cos \omega_1 \\ \mathbf{a}_3 \cdot \mathbf{r}_2 =& \sin \omega_1. \end{align} }
These equations give Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega_1\;} with Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle 0 \leq \omega_1 \leq 2 \pi} .
The angle ω1 gives the matrix Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{R}_x(\omega_1)} with
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{align} \mathbf{R} \equiv ( \mathbf{r}_1 , \mathbf{r}_2 , \mathbf{r}_3 ) = ( \mathbf{r}_1 , \mathbf{a}_2 , \mathbf{a}_3 ) \mathbf{R}_x (\omega_1 ) = (\mathbf{a}_1, \mathbf{a}_2, \mathbf{a}_3)\mathbf{R}_x (\omega_1 ) = \mathbf{R}_z (\omega_3 ) \, \mathbf{R}_y (\omega_2 ) \, \mathbf{R}_x (\omega_1 ) . \end{align} }
This gives the required z-y-x factorization of the arbitrary proper orthogonal matrix R.
The proof of the Euler z-y-z parametrization is obtained by a small modification of the previous proof. We start by retrieving the spherical polar coordinates Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega_2\;} and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega_3\;} of the unit vector Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{r}_3 = \mathbf{a}_3 } , the third column [the rightmost multiplication by Rz(ω1) does not affect r3]. Then consider
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle ( \mathbf{r}_1, \; \mathbf{r}_2 ) = (\mathbf{a}_1, \; \mathbf{a}_2 ) \begin{pmatrix} \cos \omega_1 & -\sin \omega_1 \\ \sin \omega_1 & \cos \omega_1 \\ \end{pmatrix} }
or, Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{a}_1 \cdot \mathbf{r}_1 = \cos \omega_1 \; , \quad \mathbf{a}_2 \cdot \mathbf{r}_1 = \sin \omega_1 . } The equation for R can be written as
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle ( \mathbf{r}_1 , \mathbf{r}_2 , \mathbf{r}_3 ) = ( \mathbf{a}_1 , \mathbf{a}_2 , \mathbf{r}_3 ) \, \mathbf{R}_z (\omega_1 ) = \mathbf{R}_z (\omega_3 ) \, \mathbf{R}_y (\omega_2 ) \, \mathbf{R}_z (\omega_1 ) \; , }
which proves the Euler z-y-z parametrization. Clearly, this factorization is equal to the one given in the previous section, with
- Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega_3 \equiv \alpha,\quad \omega_2 \equiv \beta, \quad \omega_1 \equiv \gamma. }
Note
- ↑ Translation by Ian Bruce of L. Euler, Theoria Motus Corporum Solidorum Seu Rigidorum (Theory of the motion of solid or rigid bodies), Rostock (1765), pdf page 11. Later Euler returned to the angles and gave an alternative derivation, see Translation by Johan Sten of Formulae generales pro translatione quacunque corporum rigidorum (General formulas for the translation of arbitrary rigid bodies), Novi Commentarii academiae scientiarum Petropolitanae, vol. 20, (1776), pp. 189-207
- ↑ In the z-x′-z′ convention the first two Euler angles are not equal to spherical polar angles, in consequence the (m, m′) Wigner D-matrix-element carries the complex phase exp[iπ(m−m′)/2]. This phase is absent in the z-y′-z′ convention