Dual quaternion

Last updated
Plaque on Broom bridge (Dublin) commemorating Hamilton's invention of quaternions Quaternion Plaque on Broom Bridge.jpg
Plaque on Broom bridge (Dublin) commemorating Hamilton's invention of quaternions

In mathematics, the dual quaternions are an 8-dimensional real algebra isomorphic to the tensor product of the quaternions and the dual numbers. Thus, they may be constructed in the same way as the quaternions, except using dual numbers instead of real numbers as coefficients. A dual quaternion can be represented in the form A + εB, where A and B are ordinary quaternions and ε is the dual unit, which satisfies ε2 = 0 and commutes with every element of the algebra. Unlike quaternions, the dual quaternions do not form a division algebra.

Contents

In mechanics, the dual quaternions are applied as a number system to represent rigid transformations in three dimensions. [1] Since the space of dual quaternions is 8-dimensional and a rigid transformation has six real degrees of freedom, three for translations and three for rotations, dual quaternions obeying two algebraic constraints are used in this application. Since unit quaternions are subject to two algebraic constraints, unit quaternions are standard to represent rigid transformations. [2]

Similar to the way that rotations in 3D space can be represented by quaternions of unit length, rigid motions in 3D space can be represented by dual quaternions of unit length. This fact is used in theoretical kinematics (see McCarthy [3] ), and in applications to 3D computer graphics, [4] robotics [5] [6] and computer vision. [7] Polynomials with coefficients given by (non-zero real norm) dual quaternions have also been used in the context of mechanical linkages design. [8] [9]

History

W. R. Hamilton introduced quaternions [10] [11] in 1843, and by 1873 W. K. Clifford obtained a broad generalization of these numbers that he called biquaternions, [12] [13] which is an example of what is now called a Clifford algebra. [3]

In 1898 Alexander McAulay used Ω with Ω2 = 0 to generate the dual quaternion algebra. [14] However, his terminology of "octonions" did not stick as today's octonions are another algebra.

In Russia, Aleksandr Kotelnikov [15] developed dual vectors and dual quaternions for use in the study of mechanics.

In 1891 Eduard Study realized that this associative algebra was ideal for describing the group of motions of three-dimensional space. He further developed the idea in Geometrie der Dynamen in 1901. [16] B. L. van der Waerden called the structure "Study biquaternions", one of three eight-dimensional algebras referred to as biquaternions.

Formulas

In order to describe operations with dual quaternions, it is helpful to first consider quaternions. [17]

A quaternion is a linear combination of the basis elements 1, i, j, and k. Hamilton's product rule for i, j, and k is often written as

Compute i ( i j k ) = −j k = −i, to obtain j k = i, and ( i j k ) k = −i j = −k or i j = k. Now because j ( j k ) = j i = −k, we see that this product yields i j = −j i, which links quaternions to the properties of determinants.

A convenient way to work with the quaternion product is to write a quaternion as the sum of a scalar and a vector (strictly speaking a bivector), that is A = a0 + A, where a0 is a real number and A = A1i + A2j + A3k is a three dimensional vector. The vector dot and cross operations can now be used to define the quaternion product of A = a0 + A and C = c0 + C as

A dual quaternion is usually described as a quaternion with dual numbers as coefficients. A dual number is an ordered pair â = ( a, b ). Two dual numbers add componentwise and multiply by the rule â ĉ = ( a, b ) ( c, d ) = (a c, a d + b c). Dual numbers are often written in the form â = a + εb, where ε is the dual unit that commutes with i, j, k and has the property ε2 = 0.

The result is that a dual quaternion can be written as an ordered pair of quaternions ( A, B ). Two dual quaternions add componentwise and multiply by the rule,

It is convenient to write a dual quaternion as the sum of a dual scalar and a dual vector, Â = â0 + A, where â0 = ( a, b ) and A = ( A, B ) is the dual vector that defines a screw. This notation allows us to write the product of two dual quaternions as

Addition

The addition of dual quaternions is defined componentwise so that given,

and

then

Multiplication

Multiplication of two dual quaternion follows from the multiplication rules for the quaternion units i, j, k and commutative multiplication by the dual unit ε. In particular, given

and

then

Notice that there is no BD term, because the definition of dual numbers requires that ε2 = 0.

This gives us the multiplication table (note the multiplication order is row times column):

Multiplication table for dual quaternion units
(Row x Column)1ijkεεiεjεk
11ijkεεiεjεk
ii−1kjεi−εεk−εj
jjk−1iεj−εk−εεi
kkji−1εkεj−εi−ε
εεεiεjεk0000
εiεi−εεk−εj0000
εjεj−εk−εεi0000
εkεkεj−εi−ε0000

Conjugate

The conjugate of a dual quaternion is the extension of the conjugate of a quaternion, that is

As with quaternions, the conjugate of the product of dual quaternions, Ĝ = ÂĈ, is the product of their conjugates in reverse order,

It is useful to introduce the functions Sc(∗) and Vec(∗) that select the scalar and vector parts of a quaternion, or the dual scalar and dual vector parts of a dual quaternion. In particular, if  = â0 + A, then

This allows the definition of the conjugate of  as

or,

The product of a dual quaternion with its conjugate yields

This is a dual scalar which is the magnitude squared of the dual quaternion.

Dual number conjugate

A second type of conjugate of a dual quaternion is given by taking the dual number conjugate, given by

The quaternion and dual number conjugates can be combined into a third form of conjugate given by

In the context of dual quaternions, the term "conjugate" can be used to mean the quaternion conjugate, dual number conjugate, or both.

Norm

The norm of a dual quaternion |Â| is computed using the conjugate to compute |Â| = Â Â*. This is a dual number called the magnitude of the dual quaternion. Dual quaternions with |Â| = 1 are unit dual quaternions.

Dual quaternions of magnitude 1 are used to represent spatial Euclidean displacements. Notice that the requirement that  Â* = 1, introduces two algebraic constraints on the components of Â, that is

Inverse

If p + ε q is a dual quaternion, and p is not zero, then the inverse dual quaternion is given by

p−1 (1 − ε qp−1).

Thus the elements of the subspace { ε q : q ∈ H } do not have inverses. This subspace is called an ideal in ring theory. It happens to be the unique maximal ideal of the ring of dual numbers.

The group of units of the dual number ring then consists of numbers not in the ideal. The dual numbers form a local ring since there is a unique maximal ideal. The group of units is a Lie group and can be studied using the exponential mapping. Dual quaternions have been used to exhibit transformations in the Euclidean group. A typical element can be written as a screw transformation.

Dual quaternions and spatial displacements

A benefit of the dual quaternion formulation of the composition of two spatial displacements DB = ([RB], b) and DA = ([RA],a) is that the resulting dual quaternion yields directly the screw axis and dual angle of the composite displacement DC = DBDA.

In general, the dual quaternion associated with a spatial displacement D = ([A], d) is constructed from its screw axis S = (S, V) and the dual angle (φ, d) where φ is the rotation about and d the slide along this axis, which defines the displacement D. The associated dual quaternion is given by,

Let the composition of the displacement DB with DA be the displacement DC = DBDA. The screw axis and dual angle of DC is obtained from the product of the dual quaternions of DA and DB, given by

That is, the composite displacement DC=DBDA has the associated dual quaternion given by

Expand this product in order to obtain

Divide both sides of this equation by the identity

to obtain

This is Rodrigues' formula for the screw axis of a composite displacement defined in terms of the screw axes of the two displacements. He derived this formula in 1840. [18]

The three screw axes A, B, and C form a spatial triangle and the dual angles at these vertices between the common normals that form the sides of this triangle are directly related to the dual angles of the three spatial displacements.

Matrix form of dual quaternion multiplication

The matrix representation of the quaternion product is convenient for programming quaternion computations using matrix algebra, which is true for dual quaternion operations as well.

The quaternion product AC is a linear transformation by the operator A of the components of the quaternion C, therefore there is a matrix representation of A operating on the vector formed from the components of C.

Assemble the components of the quaternion C = c0 + C into the array C = (C1, C2, C3, c0). Notice that the components of the vector part of the quaternion are listed first and the scalar is listed last. This is an arbitrary choice, but once this convention is selected we must abide by it.

The quaternion product AC can now be represented as the matrix product

The product AC can also be viewed as an operation by C on the components of A, in which case we have

The dual quaternion product ÂĈ = (A, B)(C, D) = (AC, AD+BC) can be formulated as a matrix operation as follows. Assemble the components of Ĉ into the eight dimensional array Ĉ = (C1, C2, C3, c0, D1, D2, D3, d0), then ÂĈ is given by the 8x8 matrix product

As we saw for quaternions, the product ÂĈ can be viewed as the operation of Ĉ on the coordinate vector Â, which means ÂĈ can also be formulated as,

More on spatial displacements

The dual quaternion of a displacement D=([A], d) can be constructed from the quaternion S=cos(φ/2) + sin(φ/2)S that defines the rotation [A] and the vector quaternion constructed from the translation vector d, given by D = d1i + d2j + d3k. Using this notation, the dual quaternion for the displacement D=([A], d) is given by

Let the Plücker coordinates of a line in the direction x through a point p in a moving body and its coordinates in the fixed frame which is in the direction X through the point P be given by,

Then the dual quaternion of the displacement of this body transforms Plücker coordinates in the moving frame to Plücker coordinates in the fixed frame by the formula

Using the matrix form of the dual quaternion product this becomes,

This calculation is easily managed using matrix operations.

Dual quaternions and 4×4 homogeneous transforms

It might be helpful, especially in rigid body motion, to represent unit dual quaternions as homogeneous matrices. As given above a dual quaternion can be written as: where r and d are both quaternions. The r quaternion is known as the real or rotational part and the quaternion is known as the dual or displacement part.

The rotation part can be given by

where is the angle of rotation about the direction given by unit vector . The displacement part can be written as

.

The dual-quaternion equivalent of a 3D-vector is

and its transformation by is given by [19]

.

These dual quaternions (or actually their transformations on 3D-vectors) can be represented by the homogeneous transformation matrix

where the 3×3 orthogonal matrix is given by

For the 3D-vector

the transformation by T is given by

Connection to Clifford algebras

Besides being the tensor product of two Clifford algebras, the quaternions and the dual numbers, the dual quaternions have two other formulations in terms of Clifford algebras.

First, dual quaternions are isomorphic to the Clifford algebra generated by 3 anticommuting elements , , with and . If we define and , then the relations defining the dual quaternions are implied by these and vice versa. Second, the dual quaternions are isomorphic to the even part of the Clifford algebra generated by 4 anticommuting elements with

For details, see Clifford algebras: dual quaternions.

Eponyms

Since both Eduard Study and William Kingdon Clifford used and wrote about dual quaternions, at times authors refer to dual quaternions as "Study biquaternions" or "Clifford biquaternions". The latter eponym has also been used to refer to split-biquaternions. Read the article by Joe Rooney linked below for view of a supporter of W.K. Clifford's claim. Since the claims of Clifford and Study are in contention, it is convenient to use the current designation dual quaternion to avoid conflict.

See also

Related Research Articles

<span class="mw-page-title-main">Pauli matrices</span> Matrices important in quantum mechanics and the study of spin

In mathematical physics and mathematics, the Pauli matrices are a set of three 2 × 2 complex matrices that are Hermitian, involutory and unitary. Usually indicated by the Greek letter sigma, they are occasionally denoted by tau when used in connection with isospin symmetries.

In mathematics, de Moivre's formula states that for any real number x and integer n it holds that

<span class="mw-page-title-main">Cross product</span> Mathematical operation on vectors in 3D space

In mathematics, the cross product or vector product is a binary operation on two vectors in a three-dimensional oriented Euclidean vector space, and is denoted by the symbol . Given two linearly independent vectors a and b, the cross product, a × b, is a vector that is perpendicular to both a and b, and thus normal to the plane containing them. It has many applications in mathematics, physics, engineering, and computer programming. It should not be confused with the dot product.

<span class="mw-page-title-main">Unit vector</span> Vector of length one

In mathematics, a unit vector in a normed vector space is a vector of length 1. A unit vector is often denoted by a lowercase letter with a circumflex, or "hat", as in .

In mechanics and geometry, the 3D rotation group, often denoted SO(3), is the group of all rotations about the origin of three-dimensional Euclidean space under the operation of composition.

In continuum mechanics, the infinitesimal strain theory is a mathematical approach to the description of the deformation of a solid body in which the displacements of the material particles are assumed to be much smaller than any relevant dimension of the body; so that its geometry and the constitutive properties of the material at each point of space can be assumed to be unchanged by the deformation.

Unit quaternions, known as versors, provide a convenient mathematical notation for representing spatial orientations and rotations of elements in three dimensional space. Specifically, they encode information about an axis-angle rotation about an arbitrary axis. Rotation and orientation quaternions have applications in computer graphics, computer vision, robotics, navigation, molecular dynamics, flight dynamics, orbital mechanics of satellites, and crystallographic texture analysis.

In the mathematical field of differential geometry, a metric tensor is an additional structure on a manifold M that allows defining distances and angles, just as the inner product on a Euclidean space allows defining distances and angles there. More precisely, a metric tensor at a point p of M is a bilinear form defined on the tangent space at p, and a metric field on M consists of a metric tensor at each point p of M that varies smoothly with p.

<span class="mw-page-title-main">Hooke's law</span> Physical law: force needed to deform a spring scales linearly with distance

In physics, Hooke's law is an empirical law which states that the force needed to extend or compress a spring by some distance scales linearly with respect to that distance—that is, Fs = kx, where k is a constant factor characteristic of the spring, and x is small compared to the total possible deformation of the spring. The law is named after 17th-century British physicist Robert Hooke. He first stated the law in 1676 as a Latin anagram. He published the solution of his anagram in 1678 as: ut tensio, sic vis. Hooke states in the 1678 work that he was aware of the law since 1660.

Linear elasticity is a mathematical model of how solid objects deform and become internally stressed due to prescribed loading conditions. It is a simplification of the more general nonlinear theory of elasticity and a branch of continuum mechanics.

In linear algebra, a rotation matrix is a transformation matrix that is used to perform a rotation in Euclidean space. For example, using the convention below, the matrix

Screw theory is the algebraic calculation of pairs of vectors, such as angular and linear velocity, or forces and moments, that arise in the kinematics and dynamics of rigid bodies.

In statistics, ordinary least squares (OLS) is a type of linear least squares method for choosing the unknown parameters in a linear regression model by the principle of least squares: minimizing the sum of the squares of the differences between the observed dependent variable in the input dataset and the output of the (linear) function of the independent variable.

<span class="mw-page-title-main">Transverse isotropy</span>

A transversely isotropic material is one with physical properties that are symmetric about an axis that is normal to a plane of isotropy. This transverse plane has infinite planes of symmetry and thus, within this plane, the material properties are the same in all directions. Hence, such materials are also known as "polar anisotropic" materials. In geophysics, vertically transverse isotropy (VTI) is also known as radial anisotropy.

<span class="mw-page-title-main">Cartesian tensor</span>

In geometry and linear algebra, a Cartesian tensor uses an orthonormal basis to represent a tensor in a Euclidean space in the form of components. Converting a tensor's components from one such basis to another is done through an orthogonal transformation.

In the study of differential equations, the Ritz method is a direct method to find an approximate solution for boundary value problems. The method is named after Walther Ritz. Some alternative formulations include the Rayleigh–Ritz method and the Ritz-Galerkin method.

In geometry, various formalisms exist to express a rotation in three dimensions as a mathematical transformation. In physics, this concept is applied to classical mechanics where rotational kinematics is the science of quantitative description of a purely rotational motion. The orientation of an object at a given instant is described with the same tools, as it is defined as an imaginary rotation from a reference placement in space, rather than an actually observed rotation from a previous placement in space.

<span class="mw-page-title-main">Composite laminate</span> Assembly of layers of fibrous composite materials

In materials science, a composite laminate is an assembly of layers of fibrous composite materials which can be joined to provide required engineering properties, including in-plane stiffness, bending stiffness, strength, and coefficient of thermal expansion.

The following are important identities in vector algebra. Identities that involve the magnitude of a vector , or the dot product of two vectors A·B, apply to vectors in any dimension. Identities that use the cross product A×B are defined only in three dimensions. Most of these relations can be dated to Josiah Willard Gibbs, founder of vector calculus, if not earlier.

In this article, we discuss certain applications of the dual quaternion algebra to 2D geometry. At this present time, the article is focused on a 4-dimensional subalgebra of the dual quaternions which we will call the planar quaternions.

References

Notes

  1. A.T. Yang, Application of Quaternion Algebra and Dual Numbers to the Analysis of Spatial Mechanisms, Ph.D thesis, Columbia University, 1963.
  2. Valverde, Alfredo; Tsiotras, Panagiotis (2018). "Dual Quaternion Framework for Modeling of Spacecraft-Mounted Multibody Robotic Systems". Frontiers in Robotics and AI. 5. doi: 10.3389/frobt.2018.00128 . ISSN   2296-9144. PMC   7805728 .
  3. 1 2 McCarthy, J. M. (1990). An Introduction to Theoretical Kinematics. MIT Press. pp. 62–5. ISBN   9780262132527.
  4. Kenwright, Ben. "Dual-Quaternions: From Classical Mechanics to Computer Graphics and Beyond" (PDF). Retrieved December 24, 2022.
  5. Figueredo, L.F.C.; Adorno, B.V.; Ishihara, J.Y.; Borges, G.A. (2013). "Robust kinematic control of manipulator robots using dual quaternion representation". 2013 IEEE International Conference on Robotics and Automation. pp. 1949–1955. doi:10.1109/ICRA.2013.6630836. ISBN   978-1-4673-5643-5. S2CID   531000.
  6. Vilhena Adorno, Bruno (2017). Robot Kinematic Modeling and Control Based on Dual Quaternion Algebra — Part I: Fundamentals.
  7. A. Torsello, E. Rodolà and A. Albarelli, Multiview Registration via Graph Diffusion of Dual Quaternions, Proc. of the XXIV IEEE Conference on Computer Vision and Pattern Recognition, pp. 2441-2448, June 2011.
  8. Li, Zijia; Schröcker, Hans-Peter; Scharler, Daniel F. (2022-09-07). "A Complete Characterization of Bounded Motion Polynomials Admitting a Factorization with Linear Factors". arXiv: 2209.02306 [math.RA].
  9. Huczala, D.; Siegele, J.; Thimm, D.; Pfurner, M.; Schröcker, H.-P. (2024). Rational Linkages: From Poses to 3D-printed Prototypes. Advances in Robot Kinematics 2024. arXiv: 2403.00558 .
  10. W. R. Hamilton, "On quaternions, or on a new system of imaginaries in algebra," Phil. Mag. 18, installments July 1844 – April 1850, ed. by D. E. Wilkins (2000)
  11. W. R. Hamilton, Elements of Quaternions, Longmans, Green & Co., London, 1866
  12. W. K. Clifford, "Preliminary sketch of bi-quaternions, Proc. London Math. Soc. Vol. 4 (1873) pp. 381–395
  13. W. K. Clifford, Mathematical Papers, (ed. R. Tucker), London: Macmillan, 1882.
  14. Alexander McAulay (1898) Octonions: a development of Clifford's Biquaternions, link from Internet Archive
  15. A. P. Kotelnikov (1895) Screw calculus and some applications to geometry and mechanics, Annal. Imp. Univ. Kazan
  16. Eduard Study (1901) Geometrie der Dynamen, Teubner, Leipzig
  17. O. Bottema and B. Roth, Theoretical Kinematics, North Holland Publ. Co., 1979
  18. Rodrigues, O. (1840), Des lois géométriques qui régissent les déplacements d’un système solide dans l’espace, et la variation des coordonnées provenant de ses déplacements considérés indépendamment des causes qui peuvent les produire, Journal de Mathématiques Pures et Appliquées de Liouville 5, 380–440.
  19. Dual Quaternions for Rigid Transformation Blending, p. 4.

Sources

Further reading