Loading [MathJax]/jax/output/HTML-CSS/jax.js
Skip to main content
Library homepage
 

Text Color

Text Size

 

Margin Size

 

Font Type

Enable Dyslexic Font
Physics LibreTexts

3.7: Rotation of Axes, Three Dimensions. Eulerian Angles

( \newcommand{\kernel}{\mathrm{null}\,}\)

We now consider two sets of orthogonal axes Ox, Oy, Oz and Ox, Oy, Oz in three-dimensional space and inclined to each other. A point in space can be described by its coordinates (x,y,z) with respect to one basis set or (x,y,z) with respect to the other. What is the relation between the coordinates (x,y,z) and the coordinates (x,y,z)?

We first need to describe exactly how the primed axes are inclined with respect to the unprimed axes. In the figure below are shown the axes Ox, Oy and Oz. Also shown are the axes Ox and Oz; the axis Oy is directed behind the plane of the paper and is not drawn. The orientation of the primed axes with respect to the unprimed axes is described by three angles θ, ϕ and ψ, known as the Eulerian angles, and they are shown in figure III.19.

The precise definitions of the three angles can be understood by three consecutive rotations, illustrated in figures III.20,21,22.

Figure 3.19.png
FIGURE III.19

First, a rotation through ϕ counterclockwise around the Oz axis to form a set of intermediate axes Ox1, Oy1, Oz1, as shown in figure III.20. The Oz and Oz1 axes are identical. Part (b) shows the rotation as seen when looking directly down the Oz (or Oz1 ) axis.

Figure 3.20.png
FIGURE III.20

The relation between the (x,y,z) and (x1, y1, z1) coordinates is

(x1y1)=(cosϕsinϕsinϕcosϕ)(xy).

Next, a rotation through θ counterclockwise around the Ox1 axis to form a set of axes Ox2, Oy2, Oz2. The Ox1 and Ox2 axes are identical (Figure III.21). Part (b) of the figure shows the rotation as seen when looking directly towards the origin along the Ox1 (or Ox2 ) axis.

Figure 3.21.png
FIGURE III.21

The relation between the (x1,y1,z1) and (x2,y2,z2) coordinates is

(y2z2)=(cosθsinθsinθcosθ)(y1z1).

Lastly, a rotation through ψ counterclockwise around the Oz2 axis to form the set of axes Ox,Oy,Oz (figure III.22). The Oz2 and Oz axes are identical. Part (b) of the figure shows the rotation as seen when looking directly down the Oz2 (or Oz) axis.

Figure 3.22.png
FIGURE III.22

The relation between the (x2,y2,z2) and (x,y,z) coordinates is

(xyz)=(cosψsinψ1sinψcosψ0000)(x2y2z2).

Thus we have for the relations between (x,y,z) and (x,y,z)

(xyz)=(cosψsinψ0sinψcosψ0001)(1000cosθsinθ0sinθcosθ)(cosϕsinϕ0sinϕcosϕ0001)(xyz).

On multiplication of these matrices, we obtain

(xyz)=(cosψcosϕcosθsinϕsinψcosψsinϕ+cosθcosϕsinψsinψsinθsinψcosϕcosθsinϕcosψsinψsinϕ+cosθcosϕcosψcosψsinθsinθsinϕsinθcosϕcosθ)(xyz).

The inverse of this may be found, as in the two-dimensional case, either by solving these three Equations for x, y and z (which would be rather tedious); or by interchanging the primed and unprimed quantities and reversing the order and signs of all operations (replace ψ with ϕ, θ with θ, and ϕ with ψ) which is less tedious; or by recognizing that the determinant of the matrix is unity and therefore its reciprocal is its transpose, which is hardly tedious at all. The reader should verify that the determinant of the matrix is unity by multiplying it out and making use of trigonometric identities. The reason that the determinant must be unity, however, and that the rotation matrix must be orthogonal, is that rotation of axes cannot change the magnitude of a vector.

Each element of the matrix is the cosine of the angle between an axis in one basis set and an axis in the other basis set. For example, the second element in the first row is the cosine of the angles between Ox and Oy. The first element of the third row is the cosine of the angles between Oz and Ox. The matrix can be referred to as the matrix of direction cosines between the axes of one basis set and the axes of the other basis set, and the relations between the coordinates can be written

(xyz)=(c11c12c13c21c22c23c31c32c33)(xyz)

R=CR.

You will note the similarity of the forms of the direction cosines to the cosine formula for the solution of a spherical triangle, and indeed the direction cosines can all be derived by drawing and solving the relevant spherical triangles. You might (or might not!)enjoy trying to do this.

The matrix C of the direction cosines is orthogonal, and the properties of an orthogonal matrix are as follows. The reader should verify this using the formulas for the direction cosines in terms of the Eulerian angles. The properties also apply, of course, although more trivially, to the rotation matrix in two dimensions.

  1. det C=±1 and (det C=1 implies that the two basis sets are of opposite chirality or "handedness"; that is, if one basis set is right-handed, the other is left-handed.)
  2. The sum of the squares of the elements in any row or any column is unity. This merely means that the magnitudes of unit orthogonal vectors are indeed unity.
  3. The sum of the products of corresponding elements in any two rows or any two columns is zero. This is merely a reflection of the fact that the scalar o dot product of any two unit orthogonal vectors is zero.
  4. Every element is equal to its own cofactor. This a reflection of the fact that the vector or cross product of any two unit orthogonal vectors in cyclic order is equal to the third.
  5. C1=˜C, or the reciprocal of an orthogonal matrix is equal to its transpose.

The first four properties above can be (and should be) used in a numerical case to verify that the matrix is indeed orthogonal, and they can be used for detecting and for correcting mistakes.

For example, the following matrix is supposed to be orthogonal, but there are, in fact, two mistakes in it. Using properties (b) and (c) above, locate and correct the mistakes. (It will become clear when you do this why verification of property (b) alone is not sufficient.) When you have corrected the matrix, see if you can find the Eulerian angles θ, ϕ and ψ without ambiguity of quadrant. As a hint, start at the bottom right hand side of the matrix and note, from the way in which the Eulerian angles are set up, that θ must be between 0 and 180, so that there is no ambiguity of quadrant. The other two angles, however, can lie between 0 and 360 and must be determined by examining the signs of their sines and cosines. When you have calculated the Eulerian angles, a further useful exercise would be to prepare a drawing showing the orientation of the primed axes with respect to the unprimed axes.

(+0.075 284 882 70.518 674 468 2+0.851 650 739 60.553 110 473 20.732 363 000 8+0.397 131 261 90.829 699 337 5+0.442 158 963 20.342 020 143 3)

Note, as a matter of good computational practice, that the numbers are written in groups of three separated by half spaces after the decimal point, all numbers, positive and negative, are signed, and leading zeroes are not omitted.


This page titled 3.7: Rotation of Axes, Three Dimensions. Eulerian Angles is shared under a CC BY-NC 4.0 license and was authored, remixed, and/or curated by Jeremy Tatum via source content that was edited to the style and standards of the LibreTexts platform.

Support Center

How can we help?