Quaternion
|
- For other uses, see Quaternion (disambiguation).
In mathematics, the quaternions are a non-commutative extension of the complex numbers. They were first described by Sir William Rowan Hamilton of Ireland in 1843. At first, the quaternions were regarded as pathological, because they disobeyed the commutative law ab = ba. However, today they find many uses in both theoretical and applied mathematics.
In modern language, the quaternions form a 4-dimensional normed division algebra over the real numbers. The algebra of quaternions is often denoted by H (for Hamilton), or in blackboard bold by <math>\mathbb H<math>.
Contents |
Definition
While the complex numbers are obtained by adding the element i to the real numbers which satisfies i2 = −1, the quaternions are obtained by adding the elements i, j and k to the real numbers which satisfy the following relations.
- <math>i^2 = j^2 = k^2 = ijk = -1<math>
- <math>\begin{matrix}
ij & = & k, & & & & ji & = & -k, \\ jk & = & i, & & & & kj & = & -i, \\ ki & = & j, & & & & ik & = & -j. \end{matrix}<math>
Every quaternion is a real linear combination of the basis quaternions 1, i, j, and k, i.e. every quaternion is uniquely expressible in the form a + bi + cj + dk where a, b, c, and d are real numbers. In other words, as a vector space over the real numbers, the set H of all quaternions has dimension 4, whereas the complex number plane has dimension 2. Addition of quaternions is accomplished by adding corresponding coefficients, as with the complex numbers. By linearity, multiplication of quaternions is completely determined by the multiplication table above for the basis quaternions. Under this multiplication, the basis quaternions, with their negatives, form the quaternion group of order 8, Q8.The scalar part of the quaternion is a while the remainder is the vector part. Thus a vector in the context of quaternions has zero for scalar part.
Example
Let
- <math>\begin{matrix}
x & = & 3 + i \\ y & = & 5i + j - 2k \end{matrix}<math> Then
- <math>\begin{matrix}
x + y & = & 3 + 6i + j - 2k \\ \\ xy & = & (3 + i)(5i + j - 2k) \\ & = & 15i + 3j - 6k + 5i^2 + ij - 2ik \\ & = & 15i + 3j - 6k - 5 + k + 2j \\ & = & -5 + 15i + 5j - 5k \\ \\ yx & = & (5i + j - 2k)(3 + i) \\ & = & 15i + 5i^2 + 3j + ji - 6k - 2ki \\ & = & 15i - 5 + 3j - k - 6k - 2j \\ & = & -5 + 15i + j - 7k \end{matrix}<math>
Arithmetic
Unlike real or complex numbers, multiplication of quaternions is not commutative: e.g.
- <math>\begin{matrix}
ij & = & k \\ ji & = & -k \\ jk & = & i \\ kj & = & -i \\ ki & = & j \\ ik & = & -j \\ \end{matrix}<math>
The quaternions are an example of a division ring, an algebraic structure similar to a field except for commutativity of multiplication. In particular, multiplication is still associative and every non-zero element has a unique inverse.
Quaternions form a 4-dimensional associative algebra over the reals (in fact a division algebra) and contain the complex numbers, but they do not form an associative algebra over the complex numbers. The quaternions, along with the complex and real numbers, are the only finite-dimensional associative division algebras over the field of real numbers. The non-commutativity of multiplication has some unexpected consequences, among them that polynomial equations over the quaternions can have more distinct solutions than the degree of the polynomial.
The equation <math>z^2 + 1 = 0<math>, for instance, has the infinitely-many quaternion solutions <math>z = bi + cj + dk<math> with <math>b^2 + c^2 + d^2 = 1<math>.
The conjugate <math>z^*<math> of the quaternion <math>z = a + bi + cj + dk<math> is defined as
- <math>z^* = a - bi - cj - dk \,<math>
and the absolute value of z is the non-negative real number defined by
- <math>|z| = \sqrt{zz^*} = \sqrt{a^2+b^2+c^2+d^2}. \,<math>
Note that <math>(wz)^* = z^* w^*<math>, which is not in general equal to <math>w^* z^*<math>. The multiplicative inverse of the non-zero quaternion z can be conveniently computed as z−1 = z* / |z|2.
By using the distance function d(z, w) = |z − w|, the quaternions form a metric space (isometric to the usual Euclidean metric on R4) and the arithmetic operations are continuous. We also have |zw| = |z| |w| for all quaternions z and w. Using the absolute value as norm, the quaternions form a real Banach algebra.
Fundamental formula
The set of equations
- <math> i^2 = j^2 = k^2 = i j k = -1 <math>
is the fundamental formula for quaternion multiplication. The multiplication table of basis quaternions is easily derivable from it. Since
- <math> i j k = -1 <math>
then, right-multiplying both sides by k,
- <math> i j k k = -k, <math>
- <math> i j (-1) = -k, <math>
- <math> i j = k. <math>
Alternatively, left-multiplying both sides by i,
- <math> i i j k = -i, <math>
- <math> (-1) j k = -i, <math>
- <math> j k = i. <math>
This last equation can consequently be left-multiplied on both sides by j,
- <math> j j k = j i, <math>
- <math> -k = j i, <math>
and continuing in this fashion the rest of the multiplication table is forthwith derived.
Profile
The set of quaternions that square to -1 is the set of vectors of absolute value 1, that is
- <math>\left\{ q : q ^2 = -1 \right\} = \left\{ q : q^* = -q\ \mbox{and}\ q q^* = 1 \right\} = S^2 \, <math>
From this set equality one can view H as the union of complex planes sharing the same real line and taking an imaginary unit from the set. Furthermore, the unit sphere in H, the 3-sphere, is formed by the collection of unit circles in these complex planes. As the general point on a circle is
- <math>e^{ai}= \cos {a} + i \sin {a}\,<math>
(known as Euler's formula), the general point on the 3-sphere is ear where r is a unit vector of H, <math>r \in S^2<math>.
Group rotation
As is explained in more detail in quaternions and spatial rotation, the multiplicative group of non-zero quaternions acts by conjugation on the copy of R3 consisting of quaternions with real part equal to zero. The conjugation by a unit quaternion (a quaternion of absolute value 1) with real part cos(t) is a rotation by an angle 2t, the axis of the rotation being the direction of the imaginary part. The advantages of Quaternions are:
- Non singular representation (compared with Euler angles for example)
- More compact (and faster) than matrices
- Pairs of unit quaternions can represent a rotation in 4d space.
The set of all unit quaternions forms a 3-dimensional sphere S3 and a group (a Lie group) under multiplication. S3 is the double cover of the group SO(3,R) of real orthogonal 3×3 matrices of determinant 1 since two unit quaternions correspond to every rotation under the above correspondence. The group S3 is isomorphic to SU(2), the group of complex unitary 2×2 matrices of determinant 1. Let A be the set of quaternions of the form a + bi + cj + dk where a, b, c and d are either all integers or all rational numbers with odd numerator and denominator 2. The set A is a ring and a lattice. There are 24 unit quaternions in this ring, and they are the vertices of a 24-cell regular polytope with Schläfli symbol {3,4,3}.
Representing quaternions by matrices
There are at least two ways of representing quaternions as matrices, in such a way that quaternion addition and multiplication correspond to matrix addition and matrix multiplication (i.e., quaternion-matrix homomorphisms). One is to use 2×2 complex matrices, and the other is to use 4×4 real matrices.
In the first way, the quaternion a + bi + cj + dk is represented as
- <math>\begin{pmatrix} a-di & -b+ci \\ b+ci & \;\; a+di \end{pmatrix}<math>
This representation has several nice properties.
- All complex numbers (c = d = 0) correspond to matrices with only real entries.
- The square of the absolute value of a quaternion is the same as the determinant of the corresponding matrix.
- The conjugate of a quaternion corresponds to the conjugate transpose of the matrix.
- Restricted to unit quaternions, this representation provides the isomorphism between S3 and SU(2). The latter group is important in quantum mechanics when dealing with spin; see also Pauli matrices.
In the second way, the quaternion a + bi + cj + dk is represented as
- <math>\begin{pmatrix}
\;\; a & -b & \;\; d & -c \\ \;\; b & \;\; a & -c & -d \\ -d & \;\; c & \;\; a & -b \\ \;\; c & \;\; d & \;\; b & \;\; a
\end{pmatrix}<math>
In this representation, the conjugate of a quaternion corresponds to the transpose of the matrix.
Quaternion operations
Quaternion operations have extended applications in electrodynamics and general relativity. The use of quaternions can replace tensors in representation. It is sometimes easier to use quaternions with complex elements, leading to a form that is not a division algebra. However, the same operations can be performed using a combination of conjugate operations. Only quaternions with real elements will be discussed here. The discussion will involve describing quaternions in two forms. One as a combination of a vector and a scalar, and the other as a combination of the two constructors and the bivector (i, j, and k).
Define two quaternions:
- <math>q = a + \vec{u} = a + bi + cj + dk<math>
- <math>p = t + \vec{v} = t + xi + yj + zk<math>
where <math>\vec{u}<math> represents the vector (b, c, d), and <math>\vec{v}<math> represents the vector (x, y, z).
Addition, products, and general functions
- Quaternion addition
- p + q :
Like complex numbers, vectors, and matrices, the addition of two quaternions is equivalent to summing the elements together:
- <math>p + q = a + t + \vec{u} + \vec{v} = (a + t) + (b + x)i + (c + y)j + (d + z)k<math>
Addition follows all of the commutativity and associativity rules of real and complex number.
- Quaternion multiplication
- pq :
The usual non-commutative multiplication between two quaternions is termed the Grassmann product. This product has been described briefly above. The complete form is described below:
- <math>pq = at - \vec{u}\cdot\vec{v} + a\vec{v} + t\vec{u} + \vec{v}\times\vec{u}<math>
- <math>pq = (at - bx - cy - dz) + (bt + ax + cz - dy)i + (ct + ay + dx - bz)j + (dt + az + by - cx)k<math>
Due to the non-commutative nature of the quaternion multiplication, pq is not equivalent to qp. The Grassmann product is useful to describe many other algebraic functions. The vector portion of the multiplication of qp follows:
- <math>qp = at - \vec{u}\cdot\vec{v} + a\vec{u} + t\vec{v} - \vec{v}\times\vec{u}<math>
- Quaternion dot-product
- p · q :
The dot-product is also referred to as the Euclidean inner product, and is equivalent to a 4-vector dot product. The dot product is the sum of the quantity of each element of p multiplied by each element of q. It is a commutative product between quaternions, and returns a scalar quantity.
- <math>p \cdot q = at + \vec{u}\cdot\vec{v} = at + bx + cy + dz<math>
The dot-product can be rewritten using the Grassmann product:
- <math>p \cdot q = \frac{p^*q + q^*p}{2}<math>
This product is useful to isolate an element from a quaternion. For instance, the i term can be pulled out from p:
- <math>p \cdot i = b<math>
- Quaternion outer-product
- Outer(p,q) :
The Euclidean outer-product is not used often; however, it is mentioned as a pair with the inner-product because of the similarity in the Grassmann product form:
- <math>\operatorname{Outer}(p,q) = \frac{p^*q - q^*p}{2}<math>
- <math>\operatorname{Outer}(p,q) = a\vec{u} - t\vec{v} - \vec{v}\times\vec{u}<math>
- <math>\operatorname{Inner}(p,q) = (ax - bt - cz + dy)i + (ay - ct - dx + bz)j + (az - dt - by + cx)k<math>
- Quaternion even-product
- Even(p,q) :
The even-product of quaternions is also not widely used, but it mentioned due to the similarity between it and the odd product. It is the purely symmetric product; therefore, it is completely commutative.
- <math>\operatorname{Even}(p,q) = \frac{pq + qp}{2}<math>
- <math>\operatorname{Even}(p,q) = at - \vec{u}\cdot\vec{v} + a\vec{v} + t\vec{u}<math>
- <math>\operatorname{Outer}(p,q) = (at - bx - cy - dz) + (ax + bt)i + (ay + ct)j + (az + dt)k<math>
- Quaternion cross-product
- p × q
The cross-product of quaternions is also known as the odd-product. It is equivalent to the vector cross-product, and returns a vector quantity only:
- <math>p \times q = \frac{pq - qp}{2}<math>
- <math>p \times q = \vec{u}\times\vec{v}<math>
- <math>p \times q = (cz - dy)i + (dx - bz)j + (by - xc)k<math>
- Quaternion reciprocal
- p−1
The inverse of a quaternion is defined in a way that p−1p = 1. It is defined above in the definition section, under properties (note the difference in variable notation). It is formed the same way that the complex inverse is found:
- <math>p^{-1} = \frac{p^*}{p\cdot p}<math>
The dot product of a quaternion is a scalar. The division of a quaternion by a scalar is equivalent to multiplication by the scalar inverse, such that each element of the quaternion is divided by the divisor.
- Quaternion division
- p−1q
The non-commutativity of quaternions allows for two divisions of numbers p−1q and qp−1. This means that the notation of q/p cannot be used unless p is a scalar only.
- Quaternion scalar
- Scalar(p) :
The scalar of a quaternion can be isolated in the same way that was described earlier with the dot-product:
<math>1\cdot p = \frac{p + p^*}{2} = a<math>
- Quaternion vector
- Vector(p) :
The vector of a quaternion can be isolated using the outer-product in the same way the inner product is used to isolate the scalar:
<math>\operatorname{Outer}(1, p) = \frac{p - p^*}{2} = \vec{u} = bi + cj + dk<math>
- Quaternion modulus
- <math>|p|<math> :
The absolute value of a quaternion is the scalar quantity that determines the length of the quaternion from the origin.
- <math>|p| = \sqrt{p \cdot p} = \sqrt{p^*p} = \sqrt{a^2 + b^2 + c^2 + d^2}<math>
- Quaternion sign
- sgn(p) :
The sign of a complex number finds the complex number of the same direction found on the unit circle. The quaternion sign also produces the unit quaternion:
- <math>\sgn(p) = \frac{p}{|p|}<math>
- Quaternion argument
- arg(p) :
The argument finds the angle of the 4-vector quaternion from the unit scalar (i.e. 1). This returns a scalar angle.
- <math>\arg(p) = \arccos\left(\frac{\operatorname{Scalar}(p)}{|p|}\right)<math>
Exponentials and logarithms
Exponential and logarithmic functions can be defined, because quaternions have a division algebra.
- Natural exponential: <math>\exp(p) = \exp(a)(\cos(|\vec{u}|) + \sgn(\vec{u})\sin(|\vec{u}|))<math>
- Natural logarithm: <math>\ln(p) = \ln(|p|) + \sgn(\vec{u})\arg(p)<math>
- Power: <math>p^q = \exp(\ln(p)q)<math>
Trigonometry
- Sine: <math>\sin(p) = \sin(a)\cosh(|\vec{u}|) + \cos(a)\sgn(\vec{u})\sinh(|\vec{u}|)<math>
- Cosine: <math>\cos(p) = \cos(a)\cosh(|\vec{u}|) - \sin(a)\sgn(\vec{u})\sinh(|\vec{u}|)<math>
- Tangent: <math>\tan(p) = \frac{\sin(p)}{\cos(p)}<math>
Hyperbolic
- Hyperbolic sine: <math>\sinh(p) = \sinh(a)\cos(|\vec{u}|) + \cosh(a)\sgn(|\vec{u}|)\sin(|\vec{u}|)<math>
- Hyperbolic cosine: <math>\cosh(p) = \cosh(a)\cos(|\vec{u}|) + \sinh(a)\sgn(|\vec{u}|)\sin(|\vec{u}|)<math>
- Hyperbolic tangent: <math>\tanh(p) = \frac{\sinh(p)}{\cosh(p)}<math>
Inverse hyperbolic functions
- Inverse hyperbolic sine: <math>\operatorname{arcsinh}(p) = \ln(p + \sqrt{p^2 + 1})<math>
- Inverse hyperbolic cosine: <math>\operatorname{arccosh}(p) = \ln(p + \sqrt{p^2 - 1})<math>
- Inverse hyperbolic tangent: <math>\operatorname{arctanh}(p) = \frac{\ln(1+p)-\ln(1-p)}{2}<math>
Inverse trigonometric functions
This was listed last, because the inverse hyperbolic functions needed to be defined over quaternions first.
- Inverse sine: <math>\arcsin(p) = -\sgn(\vec{u})\operatorname{arcsinh}(p \sgn(\vec{u}))<math>
- Inverse cosine: <math>\arccos(p) = -\sgn(\vec{u})\operatorname{arccosh}(p)<math>
- Inverse tangent: <math>\arctan(p) = -\sgn(\vec{u})\operatorname{arctanh}(p \sgn(\vec{u}))<math>
Construction of quaternions from complex numbers
According to the Cayley-Dickson construction, a quaternion is an ordered pair of complex numbers. Letting j be a new root of −1, different from both i and −i, and given u and v are a pair of complex numbers, then
- <math> q = u + j v <math>
is a quaternion.
If <math> u = a + i b <math> and <math> v = c + i d <math> then
- <math> q = a + i b + j c + j i d <math>.
Moreover, let
- <math> j i = - i j <math>,
so that
- <math> q = a + i b + j c + i j (-d) <math>,
and also let the product of quaternions be associative.
With these rules, we can now derive the multiplication table for i, j and i j, the imaginary components of a quaternion:
- <math> i i = -1, <math>
- <math> i j = (i j), <math>
- <math> i (i j) = (i i) j = -j, <math>
- <math> j i = - (i j), <math>
- <math> j j = -1, <math>
- <math> j (i j) = - j (j i) = - (j j) i = i, <math>
- <math> (i j) i = - (j i) i = -j (i i) = j, <math>
- <math> (i j) j = i (j j) = -i, <math>
- <math> (i j) (i j) = -(i j) (j i) = -i (j j) i = i i = -1. <math>
Notice how the dyad i j behaves just like the k in the definition.
For any complex number v = c + i d, its product with j has the following property:
- <math> j v = v^* j <math>
since
- <math> j v = j c + j i d = j c - (i j) d = (c - i d) j = v^* j <math>.
Let p be the quaternion with complex components w and z:
- <math> p = w + j z <math>.
Then the product q p is
- <math> q p = (u + j v) (w + j z) = u w + u j z + j v w + j v j z <math>
- <math> = u w + j u^* z + j v w + j j v^* z <math>
- <math> = (u w - v^* z) + j (u^* z + v w). <math>
Since the product of complex numbers is commutative, we have
- <math> (u + j v) (w + j z) = (u w - z v^*) + j (u^* z + w v) <math>
which is precisely how quaternion multiplication is defined by the Cayley-Dickson construction.
Note that if u = a + i b, v = c + i d, and p = a + i b + j c + k d then p′s construction from u and v is rather
- <math> p = u + v j = u + j v^* <math>.
Generalizations
If F is any field with characteristic different from 2, and a and b are elements of F, one may define a four-dimensional unitary associative algebra over F by using two generators i and j and the relations i2 = a, j2 = b and ij = −ji. These algebras are either isomorphic to the algebra of 2×2 matrices over F, or they are division algebras over F. They are called quaternion algebras.
History
Quaternions were introduced by Sir William Rowan Hamilton of Ireland in 1843. Hamilton was looking for ways of extending complex numbers (which can be viewed as points on a plane) to higher spatial dimensions. He could not do so for 3 dimensions, but 4 dimensions produce quaternions. According to the story Hamilton told, on October 16, he was out walking along the Royal Canal in Dublin with his wife when the solution in the form of the equation
- <math>i^2 = j^2 = k^2 = ijk = -1\,<math>
suddenly occurred to him; Hamilton then promptly carved this equation into the side of the nearby Brougham Bridge (now called Broom Bridge). This involved abandoning the commutative law, a radical step for the time. Vector algebra and matrices were still in the future.
Not only this, but Hamilton had in a sense invented the cross and dot products of vector algebra. Hamilton also described a quaternion as an ordered quadruple (4-tuple) of real numbers, and described the first coordinate as the 'scalar' part, and the remaining three as the 'vector' part. If two quaternions with zero scalar parts are multiplied, the scalar part of the product is the negative of the dot product of the vector parts, while the vector part of the product is the cross product. But the significance of these was still to be discovered. Hamilton proceeded to popularize quaternions with several books, the last of which, Elements of Quaternions, had 800 pages and was published shortly after his death.
Use controversy
Even by this time there was controversy about the use of quaternions. Some of Hamilton's supporters vociferously opposed the growing fields of vector algebra and vector calculus (developed by Oliver Heaviside and Willard Gibbs among others), maintaining that quaternions provided a superior notation. While this is debatable in three dimensions, quaternions cannot be directly applied in higher dimensions (though extensions like octonions and Clifford algebras may be more applicable). Vector notation had nearly universally replaced quaternions in science and engineering by the mid-20th century.
Some early formulations of Maxwell's equations used a quaternion-based notation (although Maxwell's original formulation simply used 20 equations in 20 variables), but it proved unpopular compared to the vector-based notation of Heaviside. (All of these formulations were mathematically equivalent.)
Recent years
Quaternions are often used in computer graphics (and associated geometric analysis) to represent rotations (see quaternions and spatial rotation) and orientations of objects in 3d space. They are smaller than other representations such as matrices, and operations on them such as composition can be computed more efficiently. Quaternions also see use in control theory, signal processing, attitude control, physics, and orbital mechanics, mainly for representing rotations/orientations in three dimensions. For example, it is common for spacecraft attitude-control systems to be commanded in terms of quaternions, which are also used to telemeter their current attitude. The rationale is that combining many quaternion transformations is more numerically stable than combining many matrix transformations, avoiding such phenomena as gimbal lock.
Since 1989, the National University of Ireland, Maynooth has organized a pilgrimage, where mathematicians (including Murray Gell-Mann in 2002 and Andrew Wiles in 2003) take a walk from Dunsink observatory to the Royal Canal bridge where, unfortunately, no trace of Hamilton's carving remains.
See also
- coquaternion
- associative algebra
- complex number
- division algebra
- hypercomplex number
- octonion
- quaternions and spatial rotation
- biquaternion
- hyperbolic quaternion
- tesseract
- Hurwitz quaternion
External links and resources
- Doing Physics with Quaternions (http://world.std.com/~sweetser/quaternions/qindex/qindex.html)
- Quaternions for Computer Graphics and Mechanics (Gernot Hoffman) (http://www.fho-emden.de/~hoffmann/quater12012002.pdf)
- Quaternion Calculator (http://theworld.com/~sweetser/java/qcalc/qcalc.html) [Java]
- The Physical Heritage of Sir W. R. Hamilton (http://arxiv.org/pdf/math-ph/0201058) (PDF)
- Crowe, Michael J. (1967). A History of Vector Analysis: The Evolution of the Idea of a Vectorial System University of Notre Dame Press. Surveys the major and minor vector systems of the 19th century (Hamilton, Möbius, Bellavitis, Grassmann, Tait, Peirce, Maxwell, Clifford, Gibbs, Heaviside). The competition between quaternions and other systems is a major theme.
- Kuipers, Jack (2002). Quaternions and Rotation Sequences: A Primer With Applications to Orbits, Aerospace, and Virtual Reality (Reprint edition). Princeton University Press. ISBN 0691102988
- "Quaternion (http://31.1911encyclopedia.org/Q/QU/QUATERNIONS.htm)". 1911 encyclopedia.
- Tait, Peter Guthrie, "Quaternion (http://www.ugcs.caltech.edu/~presto/papers/Quaternions-Britannica.ps)". M.A. Sec. R.S.E. Encyclopaedia Britannica, Ninth Edition, 1886, Vol. XX, pp. 160-164. (PostScript file)
- Hamilton’s Research on Quaternions (http://www.maths.tcd.ie/pub/HistMath/People/Hamilton/Quaternions.html)
- The Matrix and Quaternions FAQ (http://www.j3d.org/matrix_faq/matrfaq_latest.html) (Deprecated)
Topics in mathematics related to quantity | |
Numbers | Natural numbers | Integers | Rational numbers | Constructible numbers | Algebraic numbers | Computable numbers | Real numbers | Complex numbers | Split-complex numbers | Bicomplex numbers | Hypercomplex numbers | Quaternions | Octonions | Sedenions | Superreal numbers | Hyperreal numbers | Surreal numbers | Nominal numbers | Ordinal numbers | Cardinal numbers | p-adic numbers | Integer sequences | Mathematical constants | Large numbers | Infinity |
de:Quaternion es:Cuaterniones fr:Quaternion ko:사원수 is:Fertala it:Quaternione lt:Kvarternionas nl:Quaternion ja:四元数 pl:Kwaterniony ru:Кватернион sl:Kvaternion sv:Kvaternion zh:四元數