Rotation vectorielle

Soit E un espace vectoriel euclidien. Une rotation vectorielle de E est un élément du groupe spécial orthogonal SO(E). Si on choisit une base orthonormée de E, sa matrice dans cette base est orthogonale directe.

## Rotation vectorielle plane

[modifier | modifier le code]

### Écriture matricielle

[modifier | modifier le code]

Article détaillé : Matrice de rotation.

Dans le plan vectoriel euclidien orienté, une rotation vectorielle est simplement définie par son angle  φ {\displaystyle \varphi \,} . Sa matrice dans une base orthonormée directe est :

( cos ⁡ φ − sin ⁡ φ sin ⁡ φ cos ⁡ φ ) {\displaystyle {\begin{pmatrix}\cos \varphi &-\sin \varphi \\\\\sin \varphi &\cos \varphi \end{pmatrix}}} [1].

Autrement dit, un vecteur  U → {\displaystyle {\vec {U}}} de composantes  ( x , y ) {\displaystyle (x,y)} a pour image le vecteur  V → {\displaystyle {\vec {V}}} de composantes  ( x ′ , y ′ ) {\displaystyle (x',y')} que l'on peut calculer avec l'égalité matricielle :

( x ′ y ′ ) = ( cos ⁡ φ − sin ⁡ φ sin ⁡ φ cos ⁡ φ ) ( x y ) {\displaystyle {\begin{pmatrix}x'\\\y'\end{pmatrix}}={\begin{pmatrix}\cos \varphi &-\sin \varphi \\\\\sin \varphi &\cos \varphi \end{pmatrix}}{\begin{pmatrix}x\\\y\end{pmatrix}}} ,

c'est-à-dire que l'on a :

x ′ = x cos ⁡ φ − y sin ⁡ φ {\displaystyle x'=x\cos \varphi -y\sin \varphi \,}

y ′ = x sin ⁡ φ + y cos ⁡ φ {\displaystyle y'=x\sin \varphi +y\cos \varphi \,} .

### Exemple

[modifier | modifier le code]

Si par exemple  cos ⁡ φ = 0 , 8 {\displaystyle \cos \varphi =0{,}8} et  sin ⁡ φ = 0 , 6 {\displaystyle \sin \varphi =0{,}6} ,  φ {\displaystyle \varphi } désigne un des angles du triangle rectangle de côtés 3, 4 et 5. On peut multiplier les exemples fournissant des matrices à coefficients rationnels en utilisant à chaque fois un triplet pythagoricien.

### Écriture complexe

[modifier | modifier le code]

Ceci peut être rapproché de la formule suivante, écrite avec des nombres complexes :

x ′ + i   y ′ = ( x + i   y ) ( cos ⁡ φ + i sin ⁡ φ ) {\displaystyle x'+i\ y'=(x+i\ y)(\cos \varphi +i\sin \varphi )}

ou encore :

z ′ = x ′ + i   y ′ = ( x + i   y ) ⋅ e   i φ = z ⋅ e   i φ {\displaystyle z'=x'+i\ y'=(x+i\ y)\cdot e^{\ i\varphi }=z\cdot e^{\ i\varphi }\,} .

### Sens de rotation

[modifier | modifier le code]

Lorsque  φ {\displaystyle \varphi } est compris entre  0 {\displaystyle 0} et  π {\displaystyle \pi } et si le plan est orienté de façon usuelle, la rotation se fait dans le sens trigonométrique (ou « sens inverse des aiguilles d'une montre » ). On dit que la rotation est sénestre. Si  φ {\displaystyle \varphi } est compris entre  − π {\displaystyle -\pi } et  0 {\displaystyle 0} , la rotation se fait dans le sens des aiguilles d'une montre. Elle est dite dextre.

### Composition

[modifier | modifier le code]

La composée de deux rotations vectorielles est une rotation vectorielle dont l'angle est la somme des angles des deux rotations, ce qu'on traduit en disant que le groupe des rotations vectorielles est isomorphe au groupe  ( R / 2 π Z , + ) {\displaystyle (\mathbb {R} /2\pi \mathbb {Z} ,+)} .

### Rotations et angles

[modifier | modifier le code]

Dans la construction axiomatique de la géométrie reposant sur l'algèbre linéaire, c'est la définition des rotations planes qui permet de définir la notion d'angle [1] (voir aussi l'article Angle).

## Rotation vectorielle dans l'espace de dimension 3

[modifier | modifier le code]

### Écriture matricielle

[modifier | modifier le code]

Article détaillé : Matrice de rotation.

Dans l'espace euclidien orienté de dimension 3, une rotation vectorielle est définie par :

* un vecteur unitaire N → {\displaystyle {\vec {N}}} , qui détermine son axe : la droite des vecteurs invariants par cette rotation vectorielle est engendrée et orientée par ce vecteur ;
* son angle  φ {\displaystyle \varphi \,} , celui de la rotation vectorielle plane associée, restriction de cette rotation au plan  Π {\displaystyle \mathbf {\Pi } \,} orthogonal à l'axe.

L'orientation de ce plan est déterminée par le choix de l'orientation de l'axe. Les couples  ( N → , φ ) {\displaystyle ({\vec {N}},\varphi )} et  ( − N → , − φ ) {\displaystyle (-{\vec {N}},-\varphi )} représentent donc la même rotation de l'espace.

Nous noterons  ( n x , n y , n z ) {\displaystyle \left(n_{x},n_{y},n_{z}\right)} les coordonnées du vecteur unitaire  N → {\displaystyle {\vec {N}}} dans une base orthonormée directe  ( i → , j → , k → ) {\displaystyle ({\vec {i}},{\vec {j}},{\vec {k}})\,} fixée :

n x 2 + n y 2 + n z 2 = ‖ N → ‖ 2 = 1 {\displaystyle n_{x}^{2}+n_{y}^{2}+n_{z}^{2}=\|{\vec {N}}\|^{2}=1} .

Soit  U → {\displaystyle {\vec {U}}} un vecteur quelconque. Notons  V → {\displaystyle {\vec {V}}} son image par la rotation  ( N → , φ ) {\displaystyle ({\vec {N}},\varphi )} .

#### Cas particulier simple

[modifier | modifier le code]

Commençons par l'étude du cas particulier  N → = k → {\displaystyle {\vec {N}}={\vec {k}}} .

Le plan  Π {\displaystyle \mathbf {\Pi } \,} est alors le plan engendré par les vecteurs  i → {\displaystyle {\vec {i}}} et  j → {\displaystyle {\vec {j}}} . Le vecteur  U → {\displaystyle {\vec {U}}} se décompose en un vecteur  z k → {\displaystyle z{\vec {k}}} colinéaire à  N → {\displaystyle {\vec {N}}} qui est invariant par la rotation, et un vecteur  x i → + y j → {\displaystyle x{\vec {i}}+y{\vec {j}}} qui subit une rotation d'angle  φ {\displaystyle \varphi } dans le plan  Π {\displaystyle \mathbf {\Pi } } , et l'on peut appliquer à  x i → + y j → {\displaystyle x{\vec {i}}+y{\vec {j}}} les formules établies dans le cas des rotations vectorielles planes. On peut donc écrire :

z ′ = z {\displaystyle z'=z\,}        et     ( x ′ y ′ ) = ( cos ⁡ φ − sin ⁡ φ sin ⁡ φ cos ⁡ φ ) ( x y ) {\displaystyle {\begin{pmatrix}x'\\\y'\end{pmatrix}}={\begin{pmatrix}\cos \varphi &-\sin \varphi \\\\\sin \varphi &\cos \varphi \end{pmatrix}}{\begin{pmatrix}x\\\y\end{pmatrix}}}    comme ci-dessus,

ce qui peut s'écrire sous la forme synthétique :

( x ′ y ′ z ′ ) = ( cos ⁡ φ − sin ⁡ φ 0 sin ⁡ φ cos ⁡ φ 0 0 0 1 ) ( x y z ) {\displaystyle {\begin{pmatrix}x'\\\y'\\\z'\end{pmatrix}}={\begin{pmatrix}\cos \varphi &-\sin \varphi &0\\\\\sin \varphi &\cos \varphi &0\\\0&0&1\end{pmatrix}}{\begin{pmatrix}x\\\y\\\z\end{pmatrix}}}

#### Cas général

[modifier | modifier le code]

Si le vecteur unitaire  N → {\displaystyle {\vec {N}}} est quelconque par rapport à la base orthonormée directe  ( i → , j → , k → ) {\displaystyle ({\vec {i}},{\vec {j}},{\vec {k}})\,} qui sert à exprimer les composantes, le raisonnement est plus délicat.

Le vecteur  U → {\displaystyle {\vec {U}}} se décompose en la somme de  ( U → ⋅ N → ) N → {\displaystyle ({\vec {U}}\cdot {\vec {N}}){\vec {N}}} , colinéaire à  N → {\displaystyle {\vec {N}}} et invariant par la rotation, et de  W → = U → − ( U → ⋅ N → ) N → {\displaystyle {\vec {W}}={\vec {U}}-({\vec {U}}\cdot {\vec {N}}){\vec {N}}} , élément de  Π {\displaystyle \mathbf {\Pi } \,} et qui va subir une rotation dans ce plan. Le vecteur directement orthogonal à  W → {\displaystyle {\vec {W}}} dans le plan et de même norme est  N → ∧ W → {\displaystyle {\vec {N}}\wedge {\vec {W}}} , de sorte que l'image de  W → {\displaystyle {\vec {W}}} dans la rotation d'angle  φ {\displaystyle \varphi } est  ( cos ⁡ φ ) W → + ( sin ⁡ φ ) N → ∧ W → {\displaystyle (\cos \varphi ){\vec {W}}+(\sin \varphi ){\vec {N}}\wedge {\vec {W}}} .

Finalement, l'image de  U → {\displaystyle {\vec {U}}} par la rotation vaut :

V → = ( U → ⋅ N → ) N → + ( cos ⁡ φ ) W → + ( sin ⁡ φ ) N → ∧ W → {\displaystyle {\vec {V}}=({\vec {U}}\cdot {\vec {N}}){\vec {N}}+(\cos \varphi ){\vec {W}}+(\sin \varphi ){\vec {N}}\wedge {\vec {W}}}

et si on remplace  W → {\displaystyle {\vec {W}}} par sa valeur  U → − ( U → ⋅ N → ) N → {\displaystyle {\vec {U}}-({\vec {U}}\cdot {\vec {N}}){\vec {N}}} , on obtient :

V → = ( U → ⋅ N → ) N → + ( cos ⁡ φ ) ( U → − ( U → ⋅ N → ) N → ) + ( sin ⁡ φ ) N → ∧ U → {\displaystyle {\vec {V}}=({\vec {U}}\cdot {\vec {N}}){\vec {N}}+(\cos \varphi )({\vec {U}}-({\vec {U}}\cdot {\vec {N}}){\vec {N}})+(\sin \varphi ){\vec {N}}\wedge {\vec {U}}}

d'où finalement la formule de rotation de Rodrigues [2] :

V → = ( cos ⁡ φ )   U → + ( 1 − cos ⁡ φ ) ( U → ⋅ N → )   N → + ( sin ⁡ φ ) ( N → ∧ U → ) {\displaystyle {\vec {V}}=(\cos \varphi )\ {\vec {U}}+(1-\cos \varphi )({\vec {U}}\cdot {\vec {N}})\ {\vec {N}}+(\sin \varphi )\,\,\left({\vec {N}}\wedge {\vec {U}}\right)}

La formule encadrée ci-dessus donne l'expression vectorielle de l'image  V → {\displaystyle {\vec {V}}} d'un vecteur  U → {\displaystyle {\vec {U}}} quelconque, par la rotation  ( N → , φ ) {\displaystyle ({\vec {N}},\varphi )} .

On peut présenter le même résultat sous la forme matricielle équivalente suivante :

( x ′ y ′ z ′ ) = M ( x y z ) {\displaystyle {\begin{pmatrix}x'\\\y'\\\z'\end{pmatrix}}=M{\begin{pmatrix}x\\\y\\\z\end{pmatrix}}}

avec :

M = ( cos ⁡ φ ) ( 1 0 0 0 1 0 0 0 1 ) + ( 1 − cos ⁡ φ ) ( n x 2 n x n y n x n z n x n y n y 2 n y n z n x n z n y n z n z 2 ) +   ( sin ⁡ φ ) ( 0 − n z n y n z 0 − n x − n y n x 0 ) {\displaystyle M=(\cos \varphi ){\begin{pmatrix}1&0&0\\\0&1&0\\\0&0&1\end{pmatrix}}+(1-\cos \varphi ){\begin{pmatrix}n_{x}^{2}&n_{x}n_{y}&n_{x}n_{z}\\\n_{x}n_{y}&n_{y}^{2}&n_{y}n_{z}\\\n_{x}n_{z}&n_{y}n_{z}&n_{z}^{2}\end{pmatrix}}+\ (\sin \varphi ){\begin{pmatrix}0&-n_{z}&n_{y}\\\n_{z}&0&-n_{x}\\\\-n_{y}&n_{x}&0\end{pmatrix}}}

### Remarques

[modifier | modifier le code]

La matrice M est appelée matrice de rotation. C'est une matrice orthogonale directe, ce qui signifie que ses colonnes forment une base orthonormée directe, ou encore que sa matrice transposée est égale à sa matrice inverse et que son déterminant vaut 1.

Inversement, étant donné une matrice de rotation quelconque, on retrouve facilement le cosinus de l'angle de rotation. En effet, la trace de la matrice (c'est-à-dire la somme de ses éléments diagonaux) est égale à  1 + 2 cos ⁡ φ {\displaystyle 1+2\cos \varphi \,} . Par ailleurs, on remarque que :

M − t M = 2 ( sin ⁡ φ ) ( 0 − n z n y n z 0 − n x − n y n x 0 ) {\displaystyle M-{}^{t}M=2(\sin \varphi ){\begin{pmatrix}0&-n_{z}&n_{y}\\\n_{z}&0&-n_{x}\\\\-n_{y}&n_{x}&0\end{pmatrix}}}

ce qui permet de retrouver rapidement l'axe et le sinus associés à la rotation. Géométriquement,  M U → {\displaystyle M{\vec {U}}} et  t M U → {\displaystyle {}^{t}M{\vec {U}}} forment les deux côtés d'un losange dont le vecteur  ( M − t M ) U → = 2 ( sin ⁡ φ ) N → ∧ U → {\displaystyle (M-{}^{t}M){\vec {U}}=2(\sin \varphi ){\vec {N}}\wedge {\vec {U}}} est la diagonale, orthogonale à l'axe de rotation. C'est le losange d'Olinde Rodrigues.

### Utilisation des quaternions

[modifier | modifier le code]

Article détaillé : Quaternions et rotation dans l'espace.

On peut également faire appel à la notion de quaternions. En effet, on peut calculer l'image  V → {\displaystyle {\vec {V}}\,} du vecteur  U → {\displaystyle {\vec {U}}\,} en utilisant le produit de quaternions sous la forme suivante :

( 0 ,   V → ) = ( 0 ,   R ( φ , N → ) ( U → ) ) = ( cos ⁡ φ 2 ,   sin ⁡ φ 2   N → ) ⋅ ( 0 ,   U → ) ⋅ ( cos ⁡ φ 2 ,   − sin ⁡ φ 2   N → ) {\displaystyle (0,\ {\vec {V}})=\left(0,\ \mathbf {R} _{\left(\varphi ,{\vec {N}}\right)}({\vec {U}})\right)=(\cos {\frac {\varphi }{2}},\ \sin {\frac {\varphi }{2}}\ {\vec {N}})\cdot (0,\ {\vec {U}})\cdot (\cos {\frac {\varphi }{2}},\ -\sin {\frac {\varphi }{2}}\ {\vec {N}})}

### Composition de deux rotations vectorielles

[modifier | modifier le code]

La composée  R 2 ∘ R 1 {\displaystyle R_{2}\circ R_{1}} de deux rotations vectorielles  R 1 = ( N → 1 , φ 1 ) {\displaystyle R_{1}=({\vec {N}}_{1},\varphi _{1})} et  R 2 = ( N → 2 , φ 2 ) {\displaystyle R_{2}=({\vec {N}}_{2},\varphi _{2})} de l'espace de dimension 3 est une rotation vectorielle. Les caractéristiques  ( N → 3 , φ 3 ) {\displaystyle ({\vec {N}}_{3},\varphi _{3})} de celle-ci se déterminent à partir de  M 3 − t M 3 {\displaystyle M_{3}-{}^{t}M_{3}} , où  M 3 {\displaystyle M_{3}} est le produit  M 2 M 1 {\displaystyle M_{2}M_{1}} des matrices de rotation initiales, ou bien à partir du produit des quaternions définissant chacune des rotations, ou bien en composant les formules de Rodrigues relatives à chaque rotation.

On trouve que[3] :

cos ⁡ ( φ 3 2 ) = cos ⁡ ( φ 1 2 ) cos ⁡ ( φ 2 2 ) − sin ⁡ ( φ 1 2 ) sin ⁡ ( φ 2 2 ) ( N → 1 ⋅ N → 2 ) {\displaystyle \cos({\frac {\varphi _{3}}{2}})=\cos({\frac {\varphi _{1}}{2}})\cos({\frac {\varphi _{2}}{2}})-\sin({\frac {\varphi _{1}}{2}})\sin({\frac {\varphi _{2}}{2}})({\vec {N}}_{1}\cdot {\vec {N}}_{2})}
sin ⁡ ( φ 3 2 ) N → 3 = cos ⁡ ( φ 1 2 ) sin ⁡ ( φ 2 2 ) N → 2 + cos ⁡ ( φ 2 2 ) sin ⁡ ( φ 1 2 ) N → 1 + sin ⁡ ( φ 1 2 ) sin ⁡ ( φ 2 2 ) N → 2 ∧ N → 1 {\displaystyle \sin({\frac {\varphi _{3}}{2}}){\vec {N}}_{3}=\cos({\frac {\varphi _{1}}{2}})\sin({\frac {\varphi _{2}}{2}}){\vec {N}}_{2}+\cos({\frac {\varphi _{2}}{2}})\sin({\frac {\varphi _{1}}{2}}){\vec {N}}_{1}+\sin({\frac {\varphi _{1}}{2}})\sin({\frac {\varphi _{2}}{2}}){\vec {N}}_{2}\wedge {\vec {N}}_{1}}

## Rotations en dimension 4

[modifier | modifier le code]

Article détaillé : Rotation en quatre dimensions.

Les matrices du groupe orthogonal SO(4) peuvent de même se mettre sous forme canonique (après diagonalisation dans C) ; on montre qu'il existe deux plans vectoriels orthogonaux tels que dans une base orthonormale constituée de deux vecteurs de chaque plan, la matrice s'écrive

( cos ⁡ α − sin ⁡ α 0 0 sin ⁡ α cos ⁡ α 0 0 0 0 cos ⁡ β − sin ⁡ β 0 0 sin ⁡ β cos ⁡ β ) {\displaystyle {\begin{pmatrix}\cos \alpha &-\sin \alpha &0&0\\\\\sin \alpha &\cos \alpha &0&0\\\0&0&\cos \beta &-\sin \beta \\\0&0&\sin \beta &\cos \beta \end{pmatrix}}} .

On voit donc que la rotation est composée de deux rotations planes, et ne possède en particulier pas de vecteur fixe (pas d'« axe ») sauf si l'un des angles α ou β est nul (dans ce cas, on peut parler, par analogie avec le cas tridimensionnel, de rotation « autour » d'un plan). Si  α ≠ β {\displaystyle \alpha \neq \beta } , les deux plans sont uniques, et ce sont les seuls plans globalement invariants par la rotation ; dans le cas où  α = ± β {\displaystyle \alpha =\pm \beta } (rotations dites  isoclines), tous les plans engendrés par un vecteur et son image sont globalement invariants.

## Expression d'une rotation en algèbre géométrique

[modifier | modifier le code]

Dans une algèbre géométrique(algèbre de Clifford), on considère deux vecteurs n et m tels que n²=1 et m²=1, et la rotation composée de la symétrie par rapport à l'hyperplan orthogonal à n, et de la symétrie par rapport à l'hyperplan orthogonal à m. Alors[4], on a l'expression suivante de la rotation :

a → n m a m n {\displaystyle \quad a\to nmamn}

ce qui, en notant  R = n m {\displaystyle R=nm} et en notant  R † = m n {\displaystyle R^{\dagger }=mn} , se réécrit

a → R a R † {\displaystyle \quad a\to RaR^{\dagger }}

R {\displaystyle R} s'appelle un tourneur (rotor en anglais).

## Notes et références

[modifier | modifier le code]

1. ↑ a et b Jean Dieudonné, Algèbre linéaire et géométrie élémentaire, Paris, Hermann, 1964, p.113 pour l'étude mathématique et voir aussi la préface : "je pense en particulier aux invraisemblables confusions et paralogismes auxquels donne lieu une notion aussi simple que celle d' "angle" quand on la prend du point de vue traditionnel, alors, que, du point de vue de l'algèbre linéaire, ce n'est pas autre chose que l'étude du groupe des rotations dans le plan.",p. 13
2. ↑ Olinde Rodrigues, « Des lois géométriques qui régissent les déplacements d'un corps solide dans l'espace, et de la variation des coordonnées provenant de ces déplacements considérés indépendammet des causes qui peuvent les produire », Journal de mathématiques pures et appliquées, 1840, p. 380-440, plus spécialement p. 403
3. ↑ Olindes Rodrigues, op. cit., plus spécialement p. 408
4. ↑ en:Chris J. L. Doran et Anthony N. Lasenby, Geometric algebra for physicists, Cambridge University Press, 2003 (ISBN 978-0-521-48022-2), p. 44

## Voir aussi

[modifier | modifier le code]

### Articles connexes

[modifier | modifier le code]

* Rotation affine :
* plane
* dans l'espace
* Groupe orthogonal
* Angles d'Euler
* SO(4) (en), le groupe des rotations de l'espace à quatre dimensions
* William Rowan Hamilton
* Harmonique sphérique

### Lien externe

[modifier | modifier le code]

Utilisation de la DCM

* Portail de la géométrie

*[p.]: page