Quaternion

The quaternions , noted \ mathbb H, are a type of numbers hypercomplexes, constituting an extension of the complex numbers, extension similar to that which had led real numbers \ mathbb R to the complex numbers \ mathbb C.

Origins and principles

History

The quaternions “were discovered” by William Rowan Hamilton in 1843 starting from work of Carl Friedrich Gauss and with the previous century Leonhard Euler. He then studied the geometrical interpretation of arithmetic of complex numbers in the plan and sought to obtain results similar in space to three dimensions.

After years of research on the construction of a algebra with “triplets” of three real numbers, it butted against the multiplication, and in particular the conservation of the standard S (Georg Ferdinand Frobenius showed in 1877 that such a multiplication of triplets was impossible to define).

It had then the idea to use “quadruplets” by employing a additional Dimension. According to its dires, it went, on October 16th, 1843, along the royal channel, with his wife when suddenly to mind the solution in the form came to him from the relations: i^2 = j^2 = k^2 = ijk = -1 \, . It then engraved promptly these relations with a knife in a stone of the bridge of Brougham (now called Broom Bridge) with Dublin. This inscription, unfortunately erased by time, was replaced by a Plaque.jpg with the memory of Sir William Rowan Hamilton.

The theory was generalized, other whole like the Octonion S discovered thereafter. An element of a whole of this nature was described as Nombre hypercomplexe until the First World War. These units are now regarded as examples of semi-simple algebras. The Théorème of Artin-Wedderburn provides a method of construction generic, it is based on the theory of the Représentations of a group finished. The construction of the quaternions is given in the article Représentations of the group of the quaternions. It corresponds to the single faithful simple algebra of the representation of the Groupe of quaternions on the body of the real numbers.

Principle

Hamilton described a quaternion like quadruplet real numbers, the first element being a “scalar”, and the three elements remaining forming a “vector”, or “imaginary pure”.

It could thus define a multiplication with the good properties. This one can be summarized with this multiplication table:

Any quaternion H can be regarded as a linear combination of the four quaternions " unités" 1 , I , J , and K :

H = has \ cdot 1 + B \ cdot I + C \ cdot J + D \ cdot K \,
(where has , B , C , D is real numbers).

H can be also written: H = Z + z'·J (with Z and z' of the complex numbers of the form has + B·I )

The real numbers has , B , C and D is characteristic of H : there exists only one way of writing H in this form, and any quaternion comprising these same 4 characteristics is necessarily equal to H (the reciprocal one is true).

has is called the component real or scalar of H , while B , C and D is the complex components of H . It is also said that has is the scalar of H and that the triplet { B , C , D } or I + C \ cdot J + D \ cdot K \, is the vector of H (or its vectorial part).

This discovery involved the abandonment of the exclusive use of the commutative laws, a radical projection for the time. The Vector S and the matrix S still formed part of the future, but Hamilton to some extent had just introduced the vector Product and the scalar Produit vectors.

Mathematical properties

Relative to the other classes of numbers

The algebra of the quaternions is not any more commutative, but partially anticommutative: 1 · I = I · 1 = I but I · J = K and J · I = - K .
  • One carries out a rotation around the axis Y followed by a rotation around axis X:

  • the two cubes underwent same rotations, but in a different order. The end result is different, which expresses in a graphic way not-commutation of rotations.
|}

Summon

The sum of 2 quaternions Q_1 \, and Q_2 \, , noted Q_1 + Q_2 \, is defined as follows:

if Q_1 = has + B I + C J + D K \, and

if Q_2 = a'+ b' I + it J + of K \, , then:

Q_1 + Q_2 = (a+a') + (b+b') I + (c+c') J + (d+d') K \,

The sum is commutative and associative.

Opposed

The quaternions Q_1 \, and Q_2 \, are known as opposite if their amount is null:

Q_1 + Q_2 = 0 + 0 I + 0 J + 0 K \,

In this case, one writes:

Q_2 = - Q_1 \,

Product

The product of 2 quaternions Q_1 and Q_2, noted Q_1 \ cdot Q_2 is defined as follows:

if Q_1 = has + B I + C J + D K \, and

if Q_2 = a'+ b' I + it J + of K \, , then:

Q_1 \ cdot Q_2 = aa'-bb'-cc'-dd' + (ab'+ba'+cd'-dc') I + (ac'+ca'+db'-bd') J + (ad'+da'+bc'-cb') K \,

which can be still written:
Q_1 \ cdot Q_2 = aa'- (bb'+cc'+dd') + (ab'+ba') I + (ac'+ca') J + (ad'+da') K + (bi+cj+dk) \ wedge (b' i+c' j+d' K) \, , Q_1 \ cdot Q_2 = aa'- (\ vec V_1 \ bullet \ vec V_2) + (ab'+ba') I + (ac'+ca') J + (ad'+da') K + (\ vec V_1 \ wedge \ vec V_2) \, ,

the last formula using at the same time the produces scalar (symbol \ bullet \, ) and the vector product (symbol \ wedge \, ) of the vectorial components \ vec V_1 = (bi+cj+dk) \, and \ vec V_2 = (b' i+c' j+d' K) \, of the two quaternions.

The product is associative but, as one mentioned above, it is not (except exceptions) commutative!

The product is distributive compared to the addition:

\ begin {matrix} Q_1 \ (Q_2 + Q_3) &=& Q_1 Q_2 + Q_1 Q_3 \ \ \ \ mbox {and} \ \ (Q_1 + Q_2) \ Q_3 &=& Q_1 Q_3 + Q_2 Q_3 \ \ \ \ \ \ \ \ end {matrix} \,

Scalar

A scalar \ lambda \, can be regarded as a quaternion whose 3 complex components are null (just as a real number can be regarded as a complex number of which the imaginary part is null). One can thus define the sum and the product of a scalar and a quaternion. In this particular case, the product is commutative:

\ lambda \ Q = Q \ \ lambda \,

Conjugation

One defines the Conjugué (noted * ) of the quaternion: Q = has + B I + C J + D K \, of components has , B , C , D by:

Q^* = has - B I - C J - D K \,

The product of Q \, by its combined Q^* \, gives:

Q \ cdot Q^* = (+ B has I + C J + D K). (- B I - C J has - D K) = a^2 + b^2 + c^2 + d^2 \,
who is the square of the standard \|Q \|\, of Q \, .

Remarks

  • Attention this conjugation is not a Morphisme, it is a antimorphism, considering that \ forall (has, b) \ in \ mathbb H^2, (a*b)^*= (b^*) * (a^*) .
  • the invariants numbers by this conjugation, such as a*=a, are the realities number.

Normalizes and unit quaternions

When the standard (a^2 + b^2 + c^2 + d^2) ^ {\ frac {1} {2}} \, of a quaternion Q = has + B I + C J + D K \, is worth 1, it is said that the quaternion is normalized or that it is about a unit quaternion .

We will see below that one can establish a kind of correspondence between a unit quaternion and a vectorial rotation in the Euclidean space of dimension 3, and that this characteristic allows a very simple representation of the product of two vectorial rotations.

Opposite

If the quaternion Q \, is not null, it has single opposite,
Q^ {- 1} = \ frac {1} {a^2+b^2+c^2+d^2} \ Q^* \,
.

Division

The product being noncommutative, one can define two ways of dividing the quaternion P \, by the quaternion Q \, (not no one):

première way: P \ cdot Q^ {- 1} \,
seconde way: Q^ {- 1} \ cdot P \,

Combined of a reverse, combined sum and product of two quaternions

The equalities easily are shown:

\ begin {matrix} (Q^*) ^* &=& Q \ \ (Q^ {- 1}) ^* &=& \ frac {Q} {\|Q \|^2} \ \ (Q^*) ^ {- 1} &=& \ frac {Q} {\|Q \|^2} \ \ (Q^ {- 1}) ^ {- 1} &=& Q \ \ (Q_1 + Q_2) ^* &=& Q^*_1 + Q^*_2 \ \ (Q_1 \ cdot Q_2) ^* &=& Q^*_2 \ cdot Q^*_1 \ \ (Q_1 \ cdot Q_2) ^ {- 1} &=& Q^ {- 1} _2 \ cdot Q^ {- 1} _1 \ end {matrix} \,

The notation (has, V)

The quaternion Q = has \ cdot 1 + B \ cdot I + C \ cdot J + D \ cdot K \, peut to be broken up (and in a single way) into a formed couple of reality a \, and vector \ vec V of \ mathbb R^3 whose coordinates are ( B , C , D ).

One writes: Q = (has \, \ \ vec V) \, .

This notation makes it possible to define the sum and the product in the following way:

It also makes it possible to redefine or define the 3 following concepts:

  • combined the Q^* = (has \, \ - \ vec V) \, of Q \, ,
  • the scalar product of two quaternions: Q_1 \ bullet Q_2 = (a_1 \, \ \ vec V_1) \ bullet (a_2 \, \ \ vec V_2) = a_1 \ cdot a_2 + \ vec V_1 \ bullet \ vec V_2
from where one deduces:
  • the standard of a quaternion: \|Q \| = \ sqrt {Q \ bullet Q} = \ sqrt {Q.Q^*} = (Q.Q^*) ^ \ frac {1} {2} = (a^2 + \ vec V \ bullet \ vec V) ^ \ frac {1} {2} = (a^2+ \|\ vec V \|^2) ^ \ frac {1} {2} \,

foot-note : le produces scalar definite above is commutative and it is thus of course différent product of quaternions defined higher.

That is to say now a quaternion Q = (has \, \ \ vec V) \, unspecified; let us note q = \|Q \|\, and v = \|\ vec V \|\, . If reality v \, positive is not null, reality q \, is not it either and one can thus always write:

Q = Q \ cdot \ left (\ frac {has} {Q} \, \ \ frac {1} {Q} \ cdot \ vec V \ right) = Q \ cdot \ left (\ frac {has} {Q} \, \ \ frac {v} {Q} \ frac {1} {v} \ cdot \ vec V \ right) \,

However \ frac {1} {v} \ cdot \ vec V is a normalized vector and one can write: q^2 = a^2 + v^2 \, , or: \ left (\ frac {has} {Q} \ right) ^2 + \ left (\ frac {v} {Q} \ right) ^2= 1 .

It results from it that there exists:

  • an angle \ varphi \, (whose cosine and sine are worth \ frac {respectively has} {Q} \, and \ frac {v} {Q} \, ) and
  • a vector normalized \ vec U = \ frac {1} {v} \ cdot \ vec V

who are such as one can write the quaternion Q \, (of vector \ vec V not no one) in the form:

This way of writing a quaternion is important: the terms of the couple, q \ cos \ varphi \, and q \ sin \ varphi \ cdot \ vec U \, , are indeed respectively the scalar product and the vector product of two vectors \ vec orthogonal V_1 and \ vec V_2 with \ vec V, these 2 vectors forming between them an angle equal to \ varphi \, . And this writing makes it possible to build the multiplication of the quaternions thanks to the composition of the similarities of \ mathbb R ³ as one can see it while clicking here

Similarities of space and quaternions

To demystify the quaternions, we will make a small instructive turning by the elementary geometry and in particular by the similarities in space. A Similitude in {\ mathbb R} ^3 is entirely defined by triple given:

  • of an axis of rotation well directed (an unit vector U),
  • of a Angle 2φ defined in 2kπ near and
  • of a report/ratio of Homothety K, a reality strictly positif.
    the effect of a similarity on all the vectors can be regarded coarsely as a screwing with expansion .

Voyage and ways

More precisely, the image of transformed of a vector V (of which the origin is supposed to be located on the axis U) is obtained initially by a multiplication (homothety) of this vector per K, followed by a rotation of angle 2φ around the axis of rotation (one could also start with rotation and make it follow homothety, but it would be necessary to modify a little the explanations which will follow…). This rotation makes turn of an angle 2φ the end of vector Kv on a circle (C) centered on the axis and located in a plan perpendicular to U. However on this circle, there are two ways of carrying out the way: either by using an arc, or by using its complementary, these arcs not being able unfortunately to be distinguished by only 2φ measurement + 2kπ.
It is precisely this difficulty which the concept of quaternion makes it possible to solve. Schematically, one can say that a quaternion, it is as a similarity which could distinguish the 2 ways that can borrow associated rotation.

In the everyday life, so for a voyage between two localities L and L, you have a priori two possible ways, the distinction between these ways can be made by indicating two intermediate site-stages S and S. And while speaking about the way S and the way S, you will imply arrival and the starting localities L and L.

By preserving this analogy, it is necessary for us thus to define two intermediate points out of the two arcs of the way.

Halfway

The points located halfway are perfect for this mission. Indeed, if I divide the angle of vectors 2φ + 2kπ by 2, I obtain two distinct angles φ + 2kπ and φ-π + 2kπ. However, if I use the rotation of axis U and angle φ + 2kπ, I definite a site-stage different from that which I obtain with rotation φ-π + 2kπ. Thus to the similarity sim (U, 2φ, K), it corresponds two distinct ways which are represented by the two distinct quaternions quat (U, φ, K) and quat (U, −π+φ, K).

the formalism

The triplet (U, φ, K) can be written in an equivalent way in the shape of the couple (kcos (φ), ksin (φ) ∙U) of the notation (has, V). And by using vectors has and B orthogonal with U suitable, it is easy to show that this couple takes the form (a.b, a^b). Thus, our site-stages enable us to return to very simple operations on vectors. And as these operations are rich remarkable properties, one can define (as one saw above) a multiplication and an addition of the quaternions. You can “see” these two operations on the quaternions here: http://www.alcys.com

To announce that a promising way of research can be consulted on the site. One precisely defines in it the similarity in space in three dimensions by a Bivecteur which is with the couple of vectors what the vector is with the couple of points. The law of composition introduced into these bivecteurs is indeed noncommutative, and the restriction of this unit on the plan is the whole of the complexes.

Double product of quaternions

Just as one can calculate a double vector product, it is possible to calculate a double product of quaternions.

Correspondence between unit quaternion and vectorial rotation

One can show that transformed the \ vec V' = \ mathbf R_ {\ left \ vec NR \ right} (\ vec V) \, of any vector \ vec V \, unspecified (of the Euclidean space of dimension 3) in the rotation \ mathbf R \ left NR \ right of angle 2 \, \ varphi \, and of axis \ vec NR \, ( \ vec NR \, being a normalized vector) can be calculated thanks to the product of following quaternions:

where (\ cos \ varphi, \ \ sin \ varphi \ \ vec NR) and (\ cos \ varphi, \ - \ sin \ varphi \ \ vec NR) are two combined unit quaternions and where (0, \ \ vec V) and (0, \ \ vec V') are quaternions whose scalar component is null.

One can also write this transformed with the notation Q = has \ cdot 1 + B \ cdot I + C \ cdot J + D \ cdot K \, . If rotation is around an axis directed according to the vector \ vec V \, of coordinates (X, there, Z) (the vector being normalized) and of angle \ varphi, the associated quaternion is worth:

Composition of vectorial rotations and produced quaternions

The preceding property justifies the fact that one has habit to say, but in a not very rigorous way, that the quaternion (\ cos \ varphi, \ \ sin \ varphi \ \ vec NR) represents rotation \ mathbf R \ left NR \ right.

By using the same approximate language, one can say that the composition of two successive rotations \ mathbf R_1 puis \ mathbf R_2 is a rotation \ mathbf R which are represented by the quaternion Q = Q_2 \ cdot Q_1 \, , the quaternions Q_1 \, and Q_2 \, being the respective representatives rotations \ mathbf R_1 and \ mathbf R_2.

Let us show it!

While posing: \ vec V' = \, \ mathbf R_1 \, (\ vec V) \, , then \ vec V = \, \ mathbf R_2 \, (\ vec V') \, , the formula framed above gives us, written in a condensed way, the 2 equalities:
(0, \ \ vec V') = Q_1 \ cdot (0, \ \ vec V) \ cdot Q^*_1 et
(0, \ \ vec V
) = Q_2 \ cdot (0, \ \ vec V') \ cdot Q^*_2, which can thus be still written:

(0, \ \ vec V ) = Q_2 \ cdot \ left (0, \ \ vec V) \ cdot Q^*_1 \ right \ cdot Q^*_2 or, if one takes account of the associativeness of the product of quaternions:
(0, \ \ vec V
) = (Q_2 \ cdot Q_1) \ cdot (0, \ \ vec V) \ cdot (Q^*_1 \ cdot Q^*_2) , or:
(0, \ \ vec V ) = (Q_2 \ cdot Q_1) \ cdot (0, \ \ vec V) \ cdot (Q_2 \ cdot Q_1) ^*, by taking account of the value of combined of two quaternions.

What establishes the property announced for the composition of two rotations and that we will write:

\ Bigg (0, \ \ mathbf R_ {\ left \ vec N_2 \ right} \ left (\ mathrm R_ {\ left \ vec N_1 \ right} (\ vec V) \ right) \ Bigg) = (\ cos \ varphi_2, \ \ sin \ varphi_2 \ \ vec N_2) \ cdot (\ cos \ varphi_1, \ \ sin \ varphi_1 \ \ vec N_1) \ cdot (0, \ \ vec V) \ cdot (\ cos \ varphi_1, \ - \ sin \ varphi_1 \ \ vec N_1) \ cdot (\ cos \ varphi_2, \ - \ sin \ varphi_2 \ \ vec N_2)

Matric notations

Just as it is possible to put in correspondence the complex number z = has + I B \, with the matrix: \ begin {bmatrix} has & - B \ \ B & has \ end {bmatrix} \, , it is possible to make correspond the quaternion Q = has + B I + C J + D K \, with the following complex matrix:

or with the following real matrix:

There exist several matric representations of a quaternion. The preceding matrix is one. That which follows is more often used. Thus, the real matrix created starting from a quaternion is written this way (by keeping q=a+ib+jc+kd):

If the unit quaternion represents a rotation since the origin, one can represent it using a matrix 3x3

With these equivalences, the sum and the product of two quaternions correspond respectively to the sum and the product of the matrices which correspond to them.

Note::

The matrix complexes \ begin {bmatrix} help & - b+ic \ \ b+ic & a+id \ end {bmatrix} \, can be still written in the form:

where 4 matrices: E = \ begin {bmatrix} 1&0 \ \ 0&1 \ end {bmatrix} , I = \ begin {bmatrix} 0&-1 \ \ 1&0 \ end {bmatrix} , J = \ begin {bmatrix} 0&i \ \ i&0 \ end {bmatrix} and K = \ begin {bmatrix} - i&0 \ \ 0&i \ end {bmatrix} is the complex matrices which correspond to four quaternion-units 1, I , J and K evoked in the first definition of the quaternions.

Applications

Whereas that is debatable in dimensions three, the quaternions cannot be employed in other dimensions (although extensions like those of the Biquaternion S and the algebras of Clifford are usable). In any event, the concept of vector had almost universally replaced that of the quaternions in science and technology in the middle of the 20th century.

Today, the quaternions find their place in Infographie, Control theory, in the Treatment of the signal, in the ordering of movement and the orbital Mécanique, mainly to represent the Rotation S and the orientations in dimension three. For example, it is frequent that the control devices of displacement of a spaceship are governed in terms of quaternions. The reason is that to carry out much operations on the quaternions is numerically more stable than to carry out many operations on the matrices.

Interpolation of rotations

If one takes 2 rotations of space r_a and r_b, the linear interpolation of these rotations is in general not a rotation. To be able to interpolate, it is necessary is

In the last case, 2 rotations are represented by 2 quaternions q_a and q_b on the sphere unit S_3, and the interpolation corresponds to the Géodésique between these 2 points

See too

Related articles

External bonds

  • Definition

  • On the quaternions

Sources

Random links:Lev Pram | Lucio Fulci | Anaximene of Lampsaque | Lake Junín | Mathieu Career | L'Israël_Putnam