Dual quaternion

From David's Wiki
\( \newcommand{\P}[]{\unicode{xB6}} \newcommand{\AA}[]{\unicode{x212B}} \newcommand{\empty}[]{\emptyset} \newcommand{\O}[]{\emptyset} \newcommand{\Alpha}[]{Α} \newcommand{\Beta}[]{Β} \newcommand{\Epsilon}[]{Ε} \newcommand{\Iota}[]{Ι} \newcommand{\Kappa}[]{Κ} \newcommand{\Rho}[]{Ρ} \newcommand{\Tau}[]{Τ} \newcommand{\Zeta}[]{Ζ} \newcommand{\Mu}[]{\unicode{x039C}} \newcommand{\Chi}[]{Χ} \newcommand{\Eta}[]{\unicode{x0397}} \newcommand{\Nu}[]{\unicode{x039D}} \newcommand{\Omicron}[]{\unicode{x039F}} \DeclareMathOperator{\sgn}{sgn} \def\oiint{\mathop{\vcenter{\mathchoice{\huge\unicode{x222F}\,}{\unicode{x222F}}{\unicode{x222F}}{\unicode{x222F}}}\,}\nolimits} \def\oiiint{\mathop{\vcenter{\mathchoice{\huge\unicode{x2230}\,}{\unicode{x2230}}{\unicode{x2230}}{\unicode{x2230}}}\,}\nolimits} \)

Dual quaternions are an 8-dimensional number system (i.e. isomorphic to \(\displaystyle \mathbb{R}^8\)) which can be used to jointly represent rotations and translations in 3D space. They can be used in place of the standard \(\displaystyle 4 \times 4\) homogeneous transformation matrices.

Background

A dual quaternion can be written as \(\displaystyle \mathbf{q} = \mathbf{q}_r + \mathbf{q}_d \varepsilon\).
Here, \(\displaystyle \varepsilon^2=0\).

Scalar Multiplication

\(\displaystyle s\mathbf{q} = s\mathbf{q}_r + s \mathbf{q}_d \varepsilon\)

Addition

\(\displaystyle \mathbf{q}_1 + \mathbf{q}_2 = \mathbf{q}_{r1} +\mathbf{q}_{r2} + (\mathbf{q}_{d1} + \mathbf{q}_{d2}) \varepsilon\)

Multiplication

\(\displaystyle \mathbf{q}_1 \mathbf{q}_2 = \mathbf{q}_{r1} \mathbf{q}_{r2} + (\mathbf{q}_{r1}\mathbf{q}_{d2} + \mathbf{q}_{d1} \mathbf{q}_{r2})\varepsilon\).

Conjugate

\(\displaystyle \mathbf{q}^* = \mathbf{q}_{r}^* + \mathbf{q}_{d}^*\varepsilon\)

Magnitude

\(\displaystyle \Vert \mathbf{q} \Vert = \mathbf{q}\mathbf{q}^*\)

Rotations and Translations

A translation is represented as:
\(\displaystyle \mathbf{q}_t = [1,0,0,0][0, \frac{t_x}{2}, \frac{t_y}{2}, \frac{t_z}{2}] = 1 + \frac{\varepsilon}{2}\mathbf{t}\)

A rotation is represented as:
\(\displaystyle \mathbf{q}_r = [\cos(\frac{\theta}{2}), \sin(\frac{\theta}{2})n_x, \sin(\frac{\theta}{2})n_y, \sin(\frac{\theta}{2})n_z][0,0,0,0] = \cos(\frac{\theta}{2}) + \sin(\frac{\theta}{2}) \mathbf{n}\)

These can be combined as \(\displaystyle \mathbf{q} = \mathbf{q}_t * \mathbf{q}_r = \mathbf{q}_r + \frac{\varepsilon}{2}\mathbf{t}\mathbf{q}_r\).
Applying the transformation to a point \(\displaystyle \mathbf{v} \in \mathbb{R}^3\) is:

\(\displaystyle \mathbf{p}' = \mathbf{q}(1 + \varepsilon\mathbf{v})\mathbf{q}^*\)

Resources