quaternions and Euler angles

Quaternions and Euler Angles

Euler Angles and quaternions are two common descriptions of three dimensional orientation.

Euler angles are intuitive because they build an orientation from a sequence of elementary rotations. In aerospace applications, yaw, pitch, and roll are especially familiar. Quaternions are less intuitive geometrically, but they are compact, compose efficiently, propagate naturally from angular velocity, and avoid the coordinate singularity of an Euler angle chart.

The main difficulty in converting between them is not the algebra. It is convention management.

This article uses the PhysicsLibrary passive frame convention and develops the intrinsic moving axis 3-2-1 yaw, pitch, roll sequence.

Let

The intrinsic sequence is:

1. rotate the frame by yaw ψ about the original 3 axis;
2. rotate by pitch 𝜃 about the new moving 2 axis;
3. rotate by roll ϕ about the final moving 1 axis.

Under the PhysicsLibrary passive convention, the resulting frame transformation from A coordinates into B coordinates is

(1)

The matching unit quaternion is

(2)

The quaternion and DCM therefore use the same written order.

1 Convention declaration

PhysicsLibrary uses Hamilton multiplication,

with reversed products changing sign.

Quaternion components are displayed scalar first:

(3)

The canonical passive attitude quaternion maps coordinate components according to

(4)

For a positive frame rotation through angle α about a positive unit axis e_i,

(5)

The corresponding elementary passive DCM is denoted C_i(α).

The reverse frame quaternion and DCM satisfy

(6)

(7)

2 Euler angle families

Three angle orientation sequences are commonly divided into two families.

1. Proper Euler sequences repeat one axis. Examples include 3-1-3 and 1-2-1.
2. Tait Bryan sequences use three different axes. Examples include 3-2-1 and 1-2-3.

The aerospace yaw, pitch, roll sequence treated here is a Tait Bryan sequence.

There are six Tait Bryan sequences and six proper Euler sequences. Conversion formulas depend on which one is selected. A formula for 3-2-1 must not be applied to 1-2-3 data merely because the same three angles appear.

3 Intrinsic and extrinsic language

An intrinsic sequence uses moving axes. After the first elementary rotation, the next rotation axis moves with the frame.

An extrinsic sequence uses axes fixed in the reference frame.

The intrinsic sequence

with angles

is geometrically equivalent to an extrinsic sequence about fixed axes

with angles

This equivalence does not mean finite rotations commute. It means the same final orientation can be described from two different axis viewpoints.

PhysicsLibrary uses intrinsic moving axis language for the aerospace sequence unless an article explicitly states otherwise.

4 Elementary passive axis quaternions

For roll,

(8)

For pitch,

(9)

For yaw,

(10)

The intrinsic 3-2-1 frame sequence is therefore represented by

(11)

Although yaw occurs first in the moving axis geometric construction, its elementary quaternion appears at the right of the final passive product. This is the same order as the associated DCM product.

5 Step by step geometry

The following illustration shows the intrinsic moving axis 3-2-1 sequence.

Figure. The intrinsic 3-2-1 yaw, pitch, roll frame sequence. Step 1 applies a positive frame yaw ψ about the original 3 axis. Step 2 applies a positive frame pitch 𝜃 about the intermediate moving 2 axis. Step 3 applies a positive frame roll ϕ about the final moving 1 axis. Under the PhysicsLibrary passive convention,

and

6 Elementary passive rotation matrices

The passive elementary 1 axis matrix is

(12)

The passive elementary 2 axis matrix is

(13)

The passive elementary 3 axis matrix is

(14)

Each is the transpose of the familiar positive active elementary rotation matrix for the same axis and angle.

7 The passive -- DCM

Introduce

and

Multiplying

gives

(15)

This matrix maps coordinates from the initial frame A into the final frame B.

8 Quaternion formula for the passive -- sequence

Introduce the half angle abbreviations

(16)

(17)

(18)

For compactness in the multiplication below, write

and similarly for 𝜃 and ψ.

Then

and

Multiplying

gives

(19)

with

(20)

(21)

(22)

(23)

The scalar component is the same as in the corresponding active quaternion. The three vector components have the opposite sign because the passive frame quaternion is the conjugate of the active rotor representing the same final orientation.

9 Relation to the active -- rotor

Let

be the Hamilton quaternion that actively rotates a physical vector through the same geometric orientation while the coordinate frame remains fixed.

Then

(24)

Thus if

then

(25)

Similarly,

(26)

This is why many published Euler to quaternion formulas have the same scalar term but the opposite three vector signs from the PhysicsLibrary formulas.

10 Cross check against the quaternion generated DCM

The preceding PhysicsLibrary article defines

Substituting equations (21) through (24) into the quaternion generated DCM must reproduce equation (16).

This provides an important implementation check because the Euler product, quaternion product, and matrix product are three independently useful ways of describing the same passive frame transformation:

(27)

11 Recovering roll, pitch, and yaw from the passive DCM

For the principal nonsingular branch, choose

(28)

From equation (16),

Therefore

(29)

When cos 𝜃≠0,

and

Thus

(30)

Likewise,

and

so

(31)

These formulas are often the clearest route from a quaternion to Euler angles: first generate the passive DCM using Q09, then extract the angles from the declared DCM sequence.

12 Recovering roll, pitch, and yaw directly from a quaternion

Using the passive quaternion generated DCM,

and

Therefore the principal 3-2-1 angles can be recovered from the PhysicsLibrary passive quaternion by

(32)

(33)

(34)

These are the passive PhysicsLibrary versions of the familiar active yaw, pitch, roll extraction formulas.

13 Check: pure positive yaw

Let

The passive quaternion formula reduces to

and

Therefore

(35)

The matching DCM is

(36)

This agrees with the positive 90^∘ passive z frame rotation used as a convention diagnostic in Q07 through Q09.

14 Check: pure positive pitch

Let

Then

(37)

The DCM becomes

(38)

For 𝜃 = 90^∘,

15 Check: pure positive roll

Let

Then

(39)

The DCM becomes

(40)

For ϕ = 90^∘,

16 A nontrivial numerical example

Take

The passive quaternion is obtained from equations (21) through (24). Numerically,

Substitution gives approximately

(41)

Converting this quaternion to a DCM and then applying equations (30) through (32) recovers the original roll, pitch, and yaw values to numerical roundoff.

This round trip is an important software verification case because all three angles are nonzero.

17 Gimbal lock

Euler angles are local coordinates on the orientation manifold rather than a global one to one orientation representation.

For the intrinsic 3-2-1 sequence, the singularity occurs when

(42)

At these pitch angles,

The formulas

and

lose independent information because both numerator and denominator pairs collapse.

For

the matrix depends on the combination

For

it depends on the combination

Thus yaw and roll cannot be determined independently at the singular configuration.

The physical orientation remains perfectly well defined. Only the chosen Euler coordinate chart becomes singular.

Unit quaternions do not have this gimbal lock singularity.

18 Near the singularity

Near

small attitude perturbations can produce large numerical changes in the individual yaw and roll angles.

For estimation, simulation, and control, it is therefore usually preferable to propagate and update attitude with quaternions or DCMs and convert to Euler angles primarily for display, operator interfaces, or applications in which the sequence is physically meaningful.

A robust software conversion should also clamp the argument of arcsin into the interval [−1, 1] when small floating point errors place it just outside that range.

19 Quaternion sign and Euler extraction

Because

and

represent the same orientation, they must produce the same Euler angles on the same extraction branch.

Every expression in equations (33) through (35) is quadratic in quaternion components. Replacing all four components by their negatives leaves those expressions unchanged.

Therefore

(43)

up to the ordinary nonuniqueness and branch choices of Euler Angle coordinates.

20 Euler angle nonuniqueness

Even away from gimbal lock, Euler angle triples are not globally unique.

Angles are periodic. Adding 2π to an elementary rotation angle does not change the physical orientation.

There are also alternative triples associated with different branches of the inverse trigonometric functions.

For the principal 3-2-1 representation, PhysicsLibrary normally chooses

(44)

(45)

(46)

with an explicit singularity policy at 𝜃 = ±π∕2.

21 Implementation guidance

When converting between quaternions and Euler angles in software:

1. Declare the Euler sequence explicitly. Do not write only “yaw pitch roll” without defining the product.
2. Declare whether the Euler rotations are intrinsic moving axis rotations or extrinsic fixed axis rotations.
3. Declare the quaternion frame direction. PhysicsLibrary uses ^Bq _A for A → B coordinates.
4. Declare Hamilton multiplication separately from scalar first or scalar last storage.
5. Verify the three pure axis cases.
6. Verify at least one case with all three Euler angles nonzero.
7. Verify quaternion to DCM to Euler and Euler to quaternion to DCM round trips.
8. Treat gimbal lock explicitly rather than allowing division or atan2(0, 0) behavior to define the result accidentally.
9. Compare physical orientation rather than raw Euler triples when testing for equality.

22 Common pitfalls

1. Using the old active PhysicsLibrary formulas.

   The current house convention uses passive frame quaternions. A positive elementary frame rotation has a negative quaternion vector part.

2. Using the active product order for the passive intrinsic sequence.

   PhysicsLibrary uses

   
3. Using the wrong DCM order.

   The matching matrix is

   
4. Mixing intrinsic and extrinsic descriptions.

   An intrinsic 3-2-1 sequence is equivalent to an extrinsic 1-2-3 description with reversed angle order, not to the same verbal sequence about fixed axes.

5. Copying quaternion extraction formulas from an active source.

   The PhysicsLibrary passive formulas have different signs in the terms that are linear in q_w.

6. Confusing component storage with quaternion semantics.

   Scalar first versus scalar last storage does not determine Euler sequence, frame direction, or active versus passive interpretation.

7. Ignoring gimbal lock.

   At

   

   yaw and roll are not independently recoverable.

8. Assuming Euler angles are unique.

   Different triples can represent the same physical orientation.

9. Comparing and as different attitudes.

   They produce the same DCM and the same Euler orientation.

23 Relationship to adjacent PhysicsLibrary entries

The preceding article, quaternions and direction cosine matrices, establishes

and gives the passive quaternion generated matrix.

The present article combines that DCM relation with an intrinsic moving axis 3-2-1 sequence to derive yaw, pitch, roll conversion in both directions.

The separate examples entry, Quaternions and Euler Angles: Examples, Exercises, and Solutions, should use the same passive formulas developed here.

The next main article, quaternion kinematics and angular velocity, derives the differential equations that propagate the passive frame quaternion from body or reference resolved angular velocity.

24 Sources and convention notes

Euler angle conversion formulas are particularly sensitive to frame direction, sequence order, intrinsic versus extrinsic language, and active versus passive interpretation. A formula should therefore be accompanied by its defining matrix or quaternion product.

Henderson provides an important aerospace reference for Euler angle, quaternion, and transformation matrix relationships. Moore provides an openly licensed modern treatment of reference frame orientation. Sommer and coauthors provide a useful convention analysis for Hamilton versus flipped quaternion multiplication and passive frame transformations.

References

License

Unless otherwise noted, this PhysicsLibrary entry is intended for release under the Creative Commons Attribution ShareAlike 4.0 International license.