Appendix/Mathematics/EulerAngles

From jetwiki
Jump to navigation Jump to search


Euler Angles

Here we will follow quite closely the online class notes prepared by Professor Mona Berciu, Department of Physics & Astronomy, University of British Columbia.

Basic Relations

In terms of the unit vectors of the (X,Y,Z) Cartesian coordinate system shown in the left panel of Figure 1, we can uniquely specify the (red) vector, A, by the expression,

A=eX(AX)+eY(AY)+eZ(AZ)

where the coefficient triplet, (AX,AY,AZ), give the length (and ± direction) of each vector component. Alternatively, in terms of the unit vectors of the (x1,x2,x3) Cartesian coordinate system shown in the right panel of Figure 1, the same (red) vector is specified by the expression,

A=e1(A1)+e2(A2)+e3(A3).

Figure 1

Berciu Figure 1a

Berciu Figure 1b

Now, it is clear from the nature of unit vectors and vector dot-products that the value of the coefficient, A1 — which explicitly appears in the second of these two expressions — may be obtained from the dot product, e1A. The same must be true if we insert, for A, the first of the two expressions; that is to say,

A1=e1A

=

e1eX(AX)+e1eY(AY)+e1eZ(AZ).

Analogously, we can write,

A2=e2A

=

e2eX(AX)+e2eY(AY)+e2eZ(AZ),

A3=e3A

=

e3eX(AX)+e3eY(AY)+e3eZ(AZ).

This set of three relations therefore provides a generic way to express one set of projections in terms of the other. When written in matrix form, the set becomes,

[A1A2A3]

=

[e1eXe1eYe1eZe2eXe2eYe2eZe3eXe3eYe3eZ][AXAYAZ]=R^[AXAYAZ].

Alternatively, we could "dot" the unit-vector triplet (eX,eY,eZ) into the vector, A, in which case it would be easy to demonstrate that mapping the other direction is accomplished via the relation,

[AXAYAZ]

=

R^1[A1A2A3],

where,

R^1

[eXe1eXe2eXe3eYe1eYe2eYe3eZe1eZe2eZe3].

A Sequence of Rotations

It is quite generally true that we can transition/map/migrate from one set of orthogonal unit vectors — such as the (eX,eY,eZ) inertial/laboratory frame illustrated by the black arrows in the left panel of Figure 1 — to any other set of orthogonal unit vectors — such as the (e1,e2,e3) body frame illustrated by the black arrows in the right panel of Figure 1 — by carrying out three rotations. A fairly standard sequence of rotations is depicted in the three panels of Figure 2:

  1. Rotate the triplet of unit vectors about the Z (i.e., eZ) axis by an angle, ϕ. The result is the green coordinate system labeled, (1,2,3). Note: (a) The 3' axis remains aligned with the inertial-frame Z-axis; and (b) we will refer to the 1' axis as the line of nodes.
  2. Rotate the triplet of unit vectors about the line of nodes by an angle, θ. The result is the light-blue coordinate system labeled, (1", 2", 3"). Note that the 1" axis is aligned with the 1' axis (line of nodes).
Figure 2

Berciu Figure 2a

Berciu Figure 2b

Berciu Figure 2c

Rotation #1

Rotation #2

Rotation #3

Rotation Matrix

From the last row of the column labeled "Proper Euler angles" in Wikipedia's discussion of the rotation matrix, we find,

Z1X2Z3

=

[(c1c3c2s1s3)(c1s3c2c3s1)(s1s2)(c3s1+c1c2s3)(c1c2c3s1s3)(c1s2)(s2s3)(c3s2)(c2)]

The equivalent expression can be found in Professor Berciu's online class notes; it reads,

R^(ϕ,θ,ψ)

=

[(cosϕcosψcosθsinϕsinψ)(cosϕsinψcosθcosψsinϕ)(sinϕsinθ)(cosψsinϕ+cosϕcosθsinψ)(cosϕcosθcosψsinϕsinψ)(cosϕsinθ)(sinθsinψ)(cosψsinθ)(cosθ)]

See Also


 

Tiled Menu

Appendices: | VisTrailsEquations | VisTrailsVariables | References | Ramblings | VisTrailsImages | myphys.lsu | ADS |