Conversion between quaternions and Euler angles: Difference between revisions

Revision as of 02:11, 20 September 2023

Latest revision as of 23:19, 13 August 2025

(13 intermediate revisions by 9 users not shown)

Line 1:

{{Short description|Mathematical strategy}}

{{cleanup rewrite|date=January 2024}}

[[Rotation formalisms in three dimensions|Spatial rotations in three dimensions]] can be [[Coordinate system|parametrized]] using both [[Euler angles]] and [[Quaternions and spatial rotation|unit quaternions]]. This article explains how to convert between the two representations. Actually this simple use of "quaternions" was first presented by [[Leonhard Euler|Euler]] some seventy years earlier than [[William Rowan Hamilton|Hamilton]] to solve the problem of [[magic square]]s. For this reason the dynamics community commonly refers to quaternions in this application as "Euler parameters".

==Definition==

There are [[Quaternions_and_spatial_rotation#Alternative_conventions|two representations]] of quaternions. This article uses the more popular Hamilton.

There are two representations of quaternions. Hamilton (where w is the first component) and [[JPL]] (where w is the last component).<ref>W. G. Breckenridge, "Quaternions proposed standard conventions," NASA Jet Propulsion Laboratory, Technical Report, Oct. 1979.</ref> This article uses Hamilton for some formulas. A unit [[quaternion]] can be described as:

:<math>\mathbf{q} = \begin{bmatrix} q_w & q_x & q_y & q_z \end{bmatrix}^T</math>

A quaternion has 4 real values: {{mvar|q_w}} (the real part or the scalar part) and {{mvar|q_x q_y q_z}} (the imaginary part).

:<math>|\mathbf{q}|^2 = q_w^2 + q_x^2 + q_y^2 + q_z^2 = 1</math>

Defining the [[Quaternion#Conjugation,_the_norm,_and_reciprocal|norm of the quaternion]] as follows:

:<math display=block>|\mathbf{lVert q}|^2 \rVert = \sqrt{\,q_w^2 + q_x^2 + q_y^2 + q_z^2 = 1~}</math>

A ''unit quaternion'' satisfies:

<math display=block>\lVert q \rVert = 1</math>

We can associate a [[quaternion]] with a rotation around an axis by the following expression

Line 13 ⟶ 20:

:<math>\mathbf{q}_y = \sin(\alpha/2)\cos(\beta_y)</math>

:<math>\mathbf{q}_z = \sin(\alpha/2)\cos(\beta_z)</math>

where α is a simple rotation angle (the value in radians of the [[angle of rotation]]) and cos(β_''x''), cos(β_''y'') and cos(β_''z'') are the "[[direction cosine]]s" of the angles between the three coordinate axes and the axis of rotation. ([[Euler's rotation theorem|Euler's Rotation Theorem]]).

==Intuition==

Line 63 ⟶ 70:

== Rotation matrices ==

The [[orthogonal matrix]] (post-multiplying a [[Row and column vectors|column vector]]) corresponding to a clockwise/[[Right-hand rule|left-handed]] (looking along positive axis to origin) rotation by the unit [[quaternion]] <math>q=q_w+iq_x+jq_y+kq_z</math> is given by the [[Quaternions and spatial rotation#Quaternion-derived rotation matrix|inhomogeneous expression]]:

:<math>R = \begin{bmatrix}

Line 166 ⟶ 173:

};

// This is not in game format, it is in mathematical format.
Quaternion ToQuaternion(double roll, double pitch, double yaw) // roll (x), pitch (y), yaw (z), angles are in radians

{

// Abbreviations for the various angular functions

Line 179 ⟶ 187:

Quaternion q;

q.w = cr * cp * cy + sr * sp * sy;

q.x = sr * cp * cy +- cr * sp * sy;

q.y = cr * sp * cy -+ sr * cp * sy;

q.z = cr * cp * sy - sr * sp * cy;

Line 259 ⟶ 267:

:<math>\vec{t} = 2\vec{q} \times \vec{v}</math>

:<math>\vec{v}^{\,\prime} = \vec{v} + q_w \vec{t} + \vec{q} \times \vec{t}</math>

where <math>\times</math> indicates a three-dimensional vector [[cross product]]. This involves fewer multiplications and is therefore computationally faster. Numerical tests indicate this latter approach may be up to 30% <ref>{{cite journal |pmc=4435132|year=2015|last1=Janota|first1=A|title=Improving the Precision and Speed of Euler Angles Computation from Low-Cost Rotation Sensor Data|journal=Sensors|volume=15|issue=3|pages=7016–7039|last2=Šimák|first2=V|last3=Nemec|first3=D|last4=Hrbček|first4=J|doi=10.3390/s150307016|pmid=25806874|bibcode=2015Senso..15.7016J |doi-access=free}}</ref> faster than the original for vector rotation.

=== Proof ===