Quaternions et rotations dans l’espace

4 Pages·2013·0.239 MB·French

by  Alain Prouté

Quaternions et Rotations dans R3
Alain Prouté
Le «produit scalaire canonique» sur Rn est celui qui fait de la base canonique une base orthonormée. La norme associée
est appelée «norme euclidienne» et sera celle que nous utiliserons. L’ensemble des vecteurs de norme 1 dans Rn+1 est
noté Sn et appelé la «sphère de dimension n». C’est un espace topologique compact (car fermé borné dans Rn+1).
Rappelons que le corps C des nombres complexes peut être vu comme l’ensemble des matrices de similitude réelles,
c’est-à-dire des matrices de la forme
(cid:18) (cid:19)
a −b
b a
où a et b sont réels (une telle matrice s’identifie au nombre complexe noté a+ib). La conjugaison complexe x 7→ x est
donnée par la transposition des matrices.
Tout espace vectoriel complexe a un espace vectoriel réel sous-jacent. Si E est un espace vectoriel complexe de dimension
n, l’espace réel sous-jacent est de dimension 2n, plus précisément, si (e ,...,e ) est une base de E, (e ,ie ,...,e ,ie )
1 n 1 1 n n
est une base de l’espace réel sous-jacent. Si (x ,...,x ) sont les coordonnées (complexes) d’un vecteur x dans la base
1 n
canonique (e ,...,e ) de Cn, le carré de la norme de ce vecteur dans l’espace réel sous-jacent est donnée par
1 n
kxk2 =x x +···+x x
1 1 n n
puisqu’en posant x =a +ib , on a x x =a2 +b2.
k k k k k k k
D’une manière analogue, on peut considérer l’ensemble H des matrices de la forme
(cid:18) (cid:19)
a −b
q =
b a
où a et b sont maintenant des nombres complexes. De telles matrices sont appelées des «quaternions». Un produit de
deux quaternions est un quaternion :
(cid:18) (cid:19)(cid:18) (cid:19) (cid:18) (cid:19)
a −b c −d ac−bd −ad−bc
=
b a d c ad+bc ac−bd
De plus, la matrice identité est un quaternion. Le produit des quaternions étant le produit des matrices, les quaternions
forment une algèbre associative unitaire sur R.
(cid:18) (cid:19)
a −b
Le quaternion q = est dit «imaginaire pur» (ou plus simplement «pur») si le nombre complexe a est un
b a
imaginaire pur. H est bien sûr un espace vectoriel réel de dimension 4 avec pour «base canonique» les matrices(1)
(cid:18) (cid:19) (cid:18) (cid:19) (cid:18) (cid:19) (cid:18) (cid:19)
1 0 i 0 0 i 0 −1
1= i= j = k =
0 1 0 −i i 0 1 0
Lesquaternionsimaginairespurssontbiensûrlescombinaisonslinéaire(àcoefficientsréels)dei,j etk.Toutquaternion
q est d’une manière unique somme d’un «quaternion réel» (c’est-à-dire un multiple réel de 1 ∈ H), qui est sa «partie
réelle»,etd’unquaternionimaginairepur(quiestsa«partieimaginaire»).Le«conjugué»qdeqestobtenuenconjugant
lecomplexeaetenchangeantlesignedeb,cequirevientàconjuguerlesdeuxcomplexesaetbetàtransposerlamatrice.
Pour tout quaternion q, le produit qq est un quaternion réel. En effet, on a
(cid:18) (cid:19)(cid:18) (cid:19) (cid:18) (cid:19)
a −b a b aa+bb 0
=
b a −b a 0 aa+bb
et c’est le produit de 1 ∈ H par le carré de la norme de q dans l’espace vectoriel réel sous-jacent (dans lequel (1,i,j,k)
est la base canonique). Comme le conjugué de qq0 est clairement q0q, il est immédiat qu’on a kqkkq0k = kqq0k pour tous
quaternions q et q0.
1. Apparemment, il n’y a pas ce consensus sur la question de savoir qui est précisément i, j ou k. L’exercice 1 montre que cela n’a pas
d’importancepourvuquei2=j2=k2=ijk=−1.
Pour simplifier l’écriture, on pourra noter [[a,b]] le quaternion
(cid:18) (cid:19)
a −b
b a
On a donc [[a,b]][[c,d]]=[[ac−bd,bc+ad]]. Le conjugué de [[a,b]] est [[a,−b]].
(cid:43) Exercice 1. (a) Vérifier que H est un corps non commutatif.
(b)Montrerquel’ensembledesracinesquaternioniquesdel’équationq2+1=0estl’ensembledesquaternionsimaginaires
de norme 1, que i, j et k en font partie, et que cet ensemble est, comme sous-espace topologique de H, homéomorphe à
S2.
(c)Montrerqueleplan(vectorielréel)engendréparlevecteur1etl’unequelconquedesracinesdel’équationprécédente
est un sous-corps de H isomorphe à C.
(d) Soient i0, j0 et k0 trois quaternions tels que i02 =j02 =k02 =i0j0k0 =−1. Montrer que l’application linéaire réelle qui
envoie 1 sur 1, i sur i0, j sur j0 et k sur k0 est un automorphisme du corps H.
L’ensemble des matrices (carrées) réelles orthogonales(2) à n lignes (et donc n colonnes) est un groupe, noté O(n) et
appelé«groupeorthogonal»(pourladimensionn).SO(n)estlesous-groupedeO(n)desmatricededéterminantpositif
(et donc égal à +1). Il est appelé le «groupe spécial orthogonal» (pour la dimension n). Ces groupes sont des groupes
topologiques (les opérations du groupe sont continues) et ils sont compacts (fermés bornés dans l’espace des matrices
réelles n×n).
(cid:43) Exercice 2. (a) Montrer que l’ensemble des éléments de SO(n+1) qui laissent fixe un vecteur non nul donné de
Rn+1 est un sous-groupe isomorphe à SO(n).
(b) Montrer que deux tels sous-groupes sont conjugués (c’est-à-dire qu’il existe un automorphisme intérieur x7→axa−1
de SO(n+1) par lequel l’un des sous-groupes est l’image de l’autre).
Il est facile de vérifier que les quaternions i, j et k satisfont les équations :
i2 =j2 =k2 =ijk =−1
De ces règles on déduit que ij =−ji=k, jk =−kj =i et ki=−ik =j. On a ainsi une table de multiplication complète
sur les vecteurs de base, ce qui montre que les équations ci-dessus définissent le produit des quaternions sans ambiguïté.
Soit un quaternion q =a+bi+cj+dk. On a <(q)=a (partie réelle de q) et =(q)=bi+cj+dk (partie imaginaire de q).
L’ensembledesquaternionspurss’identifieàR3 (onidentifie(i,j,k)àlabasecanoniquedeR3).Soitq0 =a0+b0i+c0j+d0k
un second quaternion. Le produit qq0 se décompose en les termes suivants :
aa0 ab0i ac0j ad0k
ba0i −bb0 bc0k −bd0j
ca0j −cb0k −cc0 cd0i
da0k db0j −dc0i −dd0
En regroupant astucieusement ces termes, on obtient, en notant u et u0 les parties imaginaires de q et q0 :
qq0 =aa0−u.u0+au0+a0u+u∧u0
où le point désigne le produit scalaire et où ∧ désigne le produit vectoriel dans R3. Noter que la partie réelle de qq0 est
aa0−u.u0, c’est-à-dire le produit des parties réelles diminué du produit scalaire des parties imaginaires, et que la partie
imaginaire de qq0 est au0+a0u+u∧u0, c’est-à-dire la somme des produits «croisés» (partie réelle par partie imaginaire)
augmentée du produit vectoriel des parties imaginaire. La principale différence avec le produit des nombres conplexes est
donc la présence du produit vectoriel u∧u0.
En particulier, si u et u0 sont des quaternions purs, on a uu0 =u∧u0−u.u0. On fera attention au fait que les quaternions
sontsujetsàtroismultiplicationsdifférentes(uu0,u.u0 etu∧u0 ci-dessus).Commecelavanousservirplusloin,onvatout
2. C’est-à-direlesmatricesM tellesquetMM =I,oùI estlamatriceidentitéetoùM 7→tM estl’opérationdetranspositiondesmatrices.
de suite calculer le produit (a+u)x(a−u) où a est réel, et où u et x sont purs. On a d’abord (a+u)x=ax+u∧x−u.x,
puis
(a+u)x(a−u) = a2x+a(u∧x)−a(u.x)−axu−(u∧x)u+(u.x)u
= a2x+a(u∧x)−a(u.x)−a(x∧u)+a(x.u)−(u∧x)u+(u.x)u
= a2x+2a(u∧x)+2(u.x)u−((u∧x)u+(u.x)u)
Par ailleurs, (u∧x)u+(u.x)u=−(x∧u)u+(x.u)u=−xuu=xuu=kuk2x. On a donc finalement :
(a+u)x(a−u)=(a2−kuk2)x+2a(u∧x)+2(u.x)u
L’ensemble des quaternions de norme 1 n’est autre que la sphère S3 usuelle. Comme le produit de deux quaternions de
norme 1 est encore un quaternion de norme 1, comme 1 est de norme 1 et comme l’inverse d’un quaternion de norme 1
est un quaternion de norme 1, on voit que S3 est un groupe.(3)
Soit maintenant q un quaternion de norme 1, et soit x un quaternion quelconque. On a kqxk = kqkkxk = kxk. Ainsi,
la multiplication par un quaternion de norme 1 (aussi bien à gauche qu’à droite) préserve la norme, donc le produit
scalaire. C’est donc une isométrie de R4. Un produit d’isométries étant une isométrie, l’application x 7→ qxq est encore
une isométrie. De plus cette isométrie laisse 1 invariant, donc laisse globalement invariant le sous-espace orthogonal à 1,
c’est-à-dire le sous-espace R3 des quaternions purs. On voit qu’on vient de définir une application :
S3 π (cid:47)(cid:47)O(3)
Cette application est un morphisme de groupes. En effet, si q,q0 ∈S3, en composant x7→qxq avec x7→q0xq0 on obtient
x7→qq0xq0q, c’est-à-dire x7→(qq0)x(qq0).
De plus, π est continue et comme S3 est connexe l’application det◦π, qui est à valeurs dans l’espace discret {−1,+1} est
constante. Comme π(1) est l’identité de R3 (dont le déterminant est 1), on voit que cette constante est 1 et donc que π
est à valeurs dans SO(3).
Onvamaintenantdéterminerlenoyaudeπ.Supposonsqueq ∈S3 soittelqueπ(q)(x)=xpourtoutxdeR3,c’est-à-dire
pour tout quaternion pur x. En décomposant q en q =a+u (partie réelle + partie imaginaire), on a (a+u)x(a−u)=x,
c’est-à-dire
(a2−kuk2)x+2a(u∧x)+2(u.x)u=x
En supposant u non nul et en prenant pour x un vecteur non nul orthogonal à u, on obtient une base (u,x,u∧x) de R3
et l’égalité :
(a2−kuk2−1)x+2a(u∧x)=0
qui est impossible car elle implique a=0 puis kuk2 =−1. On a donc u=0, donc a=±1 car q est de norme 1. Comme
1 et −1 sont clairement dans le noyau de π, on voit que Ker(π)={−1,+1}.
On sait que tout élément de SO(3) (sauf l’identité) est une rotation autour d’un axe. Ainsi, π envoie tout quaternion de
norme 1 sur une rotation, et comme π est un morphisme de groupes, la multiplication des quaternions correspond à la
composition des rotations. Il est donc intéressant de déterminer l’axe et l’angle de la rotation π(q) en fonction de q. Ceci
nous permettra de déterminer l’axe et l’angle d’un produit de deux rotations autour d’axes distincts (faisant entre eux
un certain angle).
Toujoursavecq =a+u,notonsque(a+u)u(a−u)=(a2−kuk2+2u.u)u=(a2+kuk2)u=u(cara2+kuk2 =kqk2 =1).
Ainsi, u est fixe par la rotation π(q), dont l’axe est donc dirigé par le vecteur u (pour q =a+u6=±1 pour que u ne soit
pas nul). Il reste à trouver l’angle de π(q).
Prenons un quaternion pur x orthogonal à u. Sous l’effet de la rotation π(q), x est transformé en y = π(q)(x) et l’angle
α entre x et y est l’angle cherché. On a (toujours avec q =a+u) :
x.y =x.((a2−kuk2)x+2a(u∧x)+2(u.x)u)=(a2−kuk2)kxk2
car u.x = 0 et x.(u∧x) = u.(x∧x) = 0. En imposant la condition kxk = 1, on a kyk = 1, donc x.y = cos(α) et
cos(α) = a2−kuk2. Comme a2+kuk2 = kqk = 1, on peut poser a = cos(θ) et kuk = sin(θ), et θ est bien défini modulo
3. Enfait,un«groupedeLie».
2π, mais de toute façon compris entre 0 et π puisque sin(θ) ≥ 0. Bien sûr, θ est l’angle que fait u avec 1 dans R4. On a
donc :
cos(α)=cos(θ)2−sin(θ)2 =cos(2θ)
et donc α = ±2θ (modulo 2π), ce qui fait que α est bien déterminé modulo 2π, mais seulement au signe près. Noter
que l’expression «au signe près» n’a de sens ici que si on impose à α d’être dans l’intervalle [−π,π], ce que nous ferons
désormais.(4) Avant d’aller plus loin, remarquons que si q = ±1, la rotation π(q) est l’identité (donc d’angle 0), et si q
est pur, alors θ =π/2, donc α=π =−π (modulo 2π). Dans ces deux cas, on n’a pas de problème de signe à résoudre.
Lesigneàattribueràαrésultedetoutemanièred’uneconventiond’orientation.Ils’agiteneffetd’utiliserupourorienter
le plan (vectoriel) orthogonal à u (dans R3). Ceci peut se faire car R3 lui-même est orienté. Mais il est inutile de se casser
la tête avec ces problèmes. Il suffit en effet de prendre n’importe quelle convention ayant la propriété que si on remplace
u par −u, α sera remplacé par −α (pour la même rotation). Il y a un moyen très simple de parvenir à ce résultat qui est
de convenir que α est du signe de a. En effet, la même rotation est définie par q et par −q (π(q)=π(−q)). Or remplacer
q par −q est équivalent à remplacer simultanément a par −a et u par −u.
En résumé, la rotation π(q) (pour tout quaternion q de norme 1 non réel, c’est-à-dire tel que q = a+u avec u 6= 0) a
pour axe (orienté) la droite orientée par u, le cosinus de son angle est a2−kuk2 et le signe de l’angle (qui est bien défini
modulo 2π) est le signe de a. Dans le cas où q est pur (c’est-à-dire a = 0), l’angle de rotation est π = −π (modulo 2π),
et le signe de l’angle n’est donc plus pertinent (et la rotation est une demi-tour (symétrie) autour de l’axe défini par u).
Dans le cas où q est réel (donc q =±1), la rotation est l’identité de R3.
Intéressons-nous maintenant au produit de deux rotations d’axes non colinéaires. Soient U et V deux vecteurs non
colinéaires de R3 (on peut supposer que kUk = kVk = 1), et soient r et s les rotations d’angles α et β autour des axes
orientés définis par U et V (α et β sont donc bien définis (avec le signe) modulo 2π). Soit θ l’angle entre U et V, qui est
bien défini si on lui impose d’être dans l’intervalle [0,π]. Il s’agit de déterminer l’axe (orienté) et l’angle (avec le signe et
modulo 2π) de la rotation composée s◦r.
Déterminons d’abord un quaternion q de norme 1 tel que π(q) = r. On aura q = a+u et u = kukU. On aura aussi
r r
1+cos(α) 1−cos(α)
a2−kuk2 = cos(α) et a2+kuk2 = 1, donc a = ± = ±cos(α/2) et kuk = = |sin(α/2)|. Pour
2 2
simplifier, et quitte à remplacer U par −U, on peut supposer α positif. On a alors 0≤α/2≤π/2, et donc
q =cos(α/2)+sin(α/2)U
De la même manière, en posant q0 =b+v, on aura π(q0)=s pourvu que q0 =cos(β/2)+sin(β/2)V.
Pour calculer s◦r, il suffit donc de calculer q0q. On a :
q0q = cos(α/2)cos(β/2)−sin(α/2)sin(β/2)cos(θ)
+cos(β/2)sin(α/2)U +cos(α/2)sin(β/2)V +sin(α/2)sin(β/2)U ∧V
On voit donc que l’axe de la rotation s◦r est défini par le vecteur :
cos(β/2)sin(α/2)U +cos(α/2)sin(β/2)V +sin(α/2)sin(β/2)sin(θ)W
où W est le vecteur de norme 1 colinéaire à U ∧V et de même sens que lui. Quant au cosinus de la moitié de l’angle de
la rotation s◦r, il vaut :
cos(α/2)cos(β/2)−sin(α/2)sin(β/2)cos(θ)
On voit donc que pour θ = 0, c’est-à-dire quand U = V, l’angle de la rotation s◦r est α+β, comme on pouvait s’y
attendre,etquepourθ =π,c’est-à-direV =−U,l’angledelarotationestα−β (relativementàl’orientationdonnéepar
U), comme on pouvait également s’y attendre. Comme la fonction θ 7→cos(θ) prend toutes les valeurs comprises entre 1
et −1 quand θ varie de 0 à π, on voit qu’en composant deux rotations d’angles α et β on peut obtenir une rotation d’un
angle quelconque compris entre α+β et α−β (théorème des valeurs intermédiaires). Il suffit de faire varier l’angle entre
les axes des rotations.
On déduit facilement de ce qui précède que SO(3) est un groupe simple, c’est-à-dire qu’il ne contient pas d’autre sous-
groupe distingué que le sous-groupe réduit à 1 et lui-même. En effet, si G est un sous-groupe distingué de SO(3), non
réduit à 1, il contient une rotation non triviale c’est-à-dire d’un angle α appartenant à ]0,π]. Comme il est stable par
conjugaison, il contient toutes les rotations d’angle α. D’après ce qui précède, il contient donc toutes les rotations d’angle
compris entre 2α=α+α et 0=α−α. Par composition il contient donc toutes les rotations pour tous les angles.
4. Enfait,lavraiequestioniciestcelledusignedesin(α)etnonpasdeα.Toutefois,cesdeuxnotionsseconfondentpourαentre−πetπ.