Dual quaternion

From David's Wiki
Revision as of 21:38, 15 October 2020 by David (talk | contribs) (→‎Background)
\( \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 \epsilon\).
Here, \(\displaystyle \epsilon^2=0\).

Multiplication is:
\(\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})\epsilon\).

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{\epsilon}{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 \times \mathbf{q}_r = \mathbf{q}_r + \frac{\epsilon}{2}\mathbf{t}\mathbf{q}_r\).
Applying the transformation to a point \(\displaystyle \mathbf{p}\) is:

\(\displaystyle \mathbf{p}' = \mathbf{q}\mathbf{p}\mathbf{q}^*\)

Resources