Rotation vectorielle

Un article de Wikipédia, l'encyclopédie libre.
Cet article est une ébauche concernant la géométrie.
Vous pouvez partager vos connaissances en l’améliorant. (Comment ?).
Pour les articles homonymes, voir Rotation.

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

Sommaire

[modifier] Rotation vectorielle plane

[modifier] Écriture matricielle

Dans le plan, une rotation vectorielle est simplement définie par son angle \varphi\,. Sa matrice dans une base orthonormée directe est :

\left[\begin{matrix}\cos\varphi&-\sin\varphi \\ \sin\varphi&\cos\varphi\end{matrix}\right]


Autrement dit, un vecteur V de composantes \left[x ; y\right] a pour transformé le vecteur V' de composantes \left[x' ; y'\right] que l'on peut calculer avec l'égalité matricielle :

\left[\begin{matrix}x'\\y'\end{matrix}\right] = \left[\begin{matrix}\cos\varphi&-\sin\varphi \\ \sin\varphi&\cos\varphi\end{matrix}\right]\left[\begin{matrix}x\\y\end{matrix}\right]

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

x' = x \cos \varphi - y \sin \varphi\, et
y' = x \sin \varphi + y \cos \varphi\,    

[modifier] Écriture complexe

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

x'+ i\ y' = (x + i\ y)(\cos \varphi + i \sin \varphi)

ou encore :

z' = x'+ i\ y' = (x + i\ y)\cdot e^{\ i\varphi} = z\cdot e^{\ i\varphi}\,

[modifier] Sens de rotation

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

[modifier] Composition

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 (\mathbb R/2\pi\mathbb Z,+).

[modifier] Rotations et angles

Dans la construction axiomatique de la géométrie, c'est la définition des rotations planes qui permet de définir la notion d'angle (voir l'article Angle).

[modifier] Rotation vectorielle dans l'espace de dimension 3

[modifier] Écriture matricielle

Dans l'espace de dimension 3, une rotation vectorielle est définie par son axe \vec N orienté dont les vecteurs sont invariants par cette rotation vectorielle et par son angle \varphi\,, celui de la rotation vectorielle plane qui concerne le plan orthogonal à l'axe.

Nous supposerons que le vecteur \vec N, de coordonnées \left[n_x ; n_y ; n_z\right] dans une base orthonormée directe, est normé : \|\vec N\| = 1.

Soit \vec U un vecteur quelconque. Notons \vec V sa transformée par la rotation \left[\varphi\,\,; \vec N\right].

[modifier] Cas particulier simple

Commençons par l'étude du cas particulier où la base orthonormée directe (\vec i, \vec j, \vec k)\, est telle que \vec N = \vec k
Soient \mathbf\Pi \, le plan vectoriel orthogonal à \vec N. Compte tenu du cas particulier, le plan \mathbf\Pi \, est le plan engendré par les vecteurs \vec i et \vec j. Le vecteur \vec U se décompose en un vecteur z\vec k colinéaire à \vec N qui est invariant par la rotation, et un vecteur x\vec i + y\vec j qui subit une rotation d'angle \varphi dans le plan \mathbf\Pi, et l'on peut appliquer à x\vec i + y\vec j les formules établies dans le cas des rotations vectorielles planes. On peut donc écrire :

z' = z\,       et aussi    \left[\begin{matrix}x'\\y'\end{matrix}\right] = \left[\begin{matrix}\cos\varphi&-\sin\varphi \\ \sin\varphi&\cos\varphi\end{matrix}\right]\left[\begin{matrix}x\\y\end{matrix}\right]comme ci-dessus

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

\left[\begin{matrix}x'\\y'\\z'\end{matrix}\right] = \left[\begin{matrix}\cos\varphi&-\sin\varphi&0\\ \sin\varphi&\cos\varphi&0\\0&0&1\end{matrix}\right]\left[\begin{matrix}x\\y\\z\end{matrix}\right]

[modifier] Cas général

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

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

Finalement, l'image de \vec U par la rotation vaut :

\vec V = (\vec U \cdot \vec N) \vec N + \cos(\varphi)\vec W + \sin(\varphi)\vec N \wedge \vec W

et si on remplace \vec W par sa valeur \vec U - (\vec U \cdot \vec N) \vec N, on obtient :

\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 d'Olinde Rodrigues :

\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 du transformé \vec V d'un vecteur \vec U quelconque, dans la rotation \left[\varphi\,\,; \vec N\right] d'angle \varphi\, et d'axe \vec N normé (n^2_x+n^2_y+n^2_z=1).

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

\left[\begin{matrix}x'\\y'\\z'\end{matrix}\right] = \left[\begin{matrix}M\end{matrix}\right]\left[\begin{matrix}x\\y\\z\end{matrix}\right] avec :


M = cos\,\varphi\begin{bmatrix}1&0&0\\0&1&0\\0&0&1\end{bmatrix} + (1-\cos\varphi)\begin{bmatrix}n^2_x&n_x n_y&n_x n_z\\n_x n_y&n^2_y&n_y n_z\\n_x n_z&n_y n_z&n^2_z\end{bmatrix} +\ sin\,\varphi\begin{bmatrix}0&-n_z&n_y\\n_z&0&-n_x\\-n_y&n_x&0\end{bmatrix}



[modifier] Remarques

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\varphi\,. Par ailleurs, on remarque que :

M - {}^t M = 2\sin(\varphi)\begin{bmatrix}0&-n_z&n_y\\n_z&0&-n_x\\-n_y&n_x&0\end{bmatrix}

ce qui permet de retrouver rapidement l'axe et le sinus associés à la rotation. Géométriquement, M \vec U et {}^t M \vec U forment les deux côtés d'un losange dont le vecteur (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 de Olinde Rodrigues.

[modifier] Composition de deux rotations vectorielles

La composée de deux rotations vectorielles \left[\varphi_1\,\,; \vec N_1\right] et \left[\varphi_2\,\,; \vec N_2\right] de l'espace de dimension 3 est une rotation vectorielle. Les caractéristiques \left[\varphi_3\,\,; \vec N_3\right] de celle-ci se déterminent rapidement à partir de M3 − tM3, où M3 est le produit M2M1 des matrices de rotation initiales.

On peut également faire appel à la notion de quaternions. En effet, on peut calculer le transformé \vec V\, du vecteur \vec U\, en utilisant le produit de quaternions sous la forme suivante :

(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)

Composer des rotations revient alors à multiplier des quaternions.



[modifier] Articles connexes