Géométrie Différentielle et Application au Contrôle

Géométrique
Frédéric Jean

To cite this version:

Frédéric Jean. Géométrie Différentielle et Application au Contrôle Géométrique. Licence. France.
2022. �hal-03243924v2�

HAL Id: hal-03243924

[Link]
Submitted on 8 Mar 2022

HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est

archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents
entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non,
lished or not. The documents may come from émanant des établissements d’enseignement et de
teaching and research institutions in France or recherche français ou étrangers, des laboratoires
abroad, or from public or private research centers. publics ou privés.
 AOT 13 - Cours de l’ENSTA Paris

Géométrie Différentielle
et Application au Contrôle Géométrique

Frédéric Jean
UMA, ENSTA Paris

Notes de cours
12 octobre 2021
ii
Table des matières

1 Variétés différentiables 1
1.1 Variétés différentiables . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2 Applications différentiables . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.3 Sous-variétés . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8

2 Espaces tangent et cotangent 13

2.1 Vecteurs tangents . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.2 Dérivations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.3 Différentielle d’une application . . . . . . . . . . . . . . . . . . . . . . . . 18
2.4 Coordonnées sur l’espace tangent . . . . . . . . . . . . . . . . . . . . . . 19
2.5 Cas des sous-variétés de Rn . . . . . . . . . . . . . . . . . . . . . . . . . 22
2.6 Espace cotangent* . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23

3 Champs de vecteurs 27
3.1 Fibrés tangent et cotangent . . . . . . . . . . . . . . . . . . . . . . . . . 27
3.2 Champs de vecteurs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
3.3 Équations différentielles . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
3.4 Flots et groupes de difféomorphismes . . . . . . . . . . . . . . . . . . . . 34

4 Familles de champs de vecteurs 39

4.1 Systèmes commandés . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
4.2 Crochets et algèbre de Lie . . . . . . . . . . . . . . . . . . . . . . . . . . 42
4.3 Orbite d’une famille de champs de vecteurs . . . . . . . . . . . . . . . . . 49
4.4 Espace tangent à une orbite . . . . . . . . . . . . . . . . . . . . . . . . . 55
4.5 Distributions et mécanique non-holonome* . . . . . . . . . . . . . . . . . 58

5 Systèmes commandés affines* 63

5.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
5.2 Ensembles atteignables . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
5.3 Systèmes sans dérive . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
5.4 Systèmes affines commandables . . . . . . . . . . . . . . . . . . . . . . . 72
iv TABLE DES MATIÈRES

A Exercices 79
A.1 Champs de coordonnées . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
A.2 Véhicule articulé . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
A.3 Sphère roulant sur une sphère . . . . . . . . . . . . . . . . . . . . . . . . 82
A.4 Transfert de population pour systèmes quantiques . . . . . . . . . . . . . 86
A.5 Application du lemme de Gauss . . . . . . . . . . . . . . . . . . . . . . . 89
A.6 La sphère imaginaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89

Bibliographie 93

Index 94
Chapitre 1

Variétés différentiables

1.1 Variétés différentiables

On considère un espace topologique M (c’est-à-dire un espace de points muni d’une
topologie). On suppose que cet espace est
— à base dénombrable : la topologie de M a une base dénombrable d’ouverts. Cette
propriété est équivalente à l’existence d’un sous-ensemble dénombrable dense (par
exemple Qn pour Rn ) ;
— séparé : deux points distincts ont des voisinages distincts.
Définition. Une carte de dimension n sur M est un couple (U, φ) formé de
— un ouvert U ⊂ M ;
— un homéomorphisme φ : U → φ(U ) ⊂ Rn (un homéomorphisme est une applica-
tion continue et inversible dont l’inverse est continue).
L’ouvert U est le domaine de la carte. Pour p ∈ U , φ(p) = (x1 (p), . . . , xn (p)) ∈ Rn : φ
est ce que l’on appelle une fonction coordonnées.
Un point de M peut appartenir à deux domaines différents correspondant à deux
cartes (U, φ) et (V, ψ).
Définition. Deux cartes (U, φ) et (V, ψ) sur M sont compatibles si U ∩ V = ∅ ou si
φ ◦ ψ −1 est un difféomorphisme entre les ouverts de Rn que sont ψ(U ∩ V ) et φ(U ∩ V ).
(Rappels : Une application f : U ⊂ E → F , où E et F sont des espaces vectoriels normés
et U un ouvert de E, est un difféomorphisme de U dans f (U ) si elle est de classe C ∞ ,
inversible et si son inverse est de classe C ∞ . Une condition nécessaire et suffisante est
que f soit injective et que Df (x) soit un isomorphisme ∀x ∈ U .)
Remarques.
— A priori les dimensions des cartes n’ont pas été fixées. On pourrait donc avoir
φ(U ) ⊂ Rn et ψ(V ) ⊂ Rm avec m ̸= n. Cependant, si U ∩ V ̸= ∅, le fait que
φ ◦ ψ −1 et ψ ◦ φ−1 soient des difféomorphismes impose que les deux cartes soient
de même dimension.
2 Chap. 1 – Variétés différentiables

— Signification en coordonnées. Une carte (U, φ) donne un système local de coor-

données. Sur U ∩ V , on a donc deux systèmes de coordonnées : φ = (x1 , . . . , xn )
et ψ = (y 1 , . . . , y n ). Comme ce sont des homéomorphismes, l’application φ ◦ ψ −1
est une bijection et son inverse est ψ ◦ φ−1 . Ces deux applications s’écrivent

φ ◦ ψ −1 : y = (y 1 , . . . , y n ) 7→ (x1 = f 1 (y), . . . , xn = f n (y)),

ψ ◦ φ−1 : x = (x1 , . . . , xn ) 7→ (y 1 = g 1 (x), . . . , y n = g n (x)).

La compatibilité signifie que les fonctions f i et g i sont de classe C ∞ .

Définition. Un atlas de dimension n de M est un ensemble A = {(Uα , φα )} de cartes

de dimension n tel que :
— les ouverts Uα recouvrent M ;
— toutes les cartes de A sont compatibles deux à deux.

Un atlas permet donc de définir des coordonnées locales partout sur M . On dit que deux
atlas sont équivalents si leur union est encore un atlas, c’est-à-dire que A = {(Uα , φα )} et
A′ = {(Vβ , ψβ )} sont équivalents si toutes les cartes (Uα , φα ) et (Vβ , ψβ ) sont compatibles
deux à deux.

Définition. Une structure différentiable de dimension n sur M est une classe d’équivalence
d’atlas de dimension n de M .

En pratique on définit une structure différentiable en donnant un atlas représentant

la classe.

Définition. Une variété différentiable de dimension n est un espace topologique M séparé

et à base dénombrable muni d’une structure différentiable de dimension n.

Remarques.
— Un même espace topologique peut être muni de plusieurs structures différentiables
différentes. Par exemple, pour la topologie canonique, la sphère unité S 7 de dimen-
sion 7 peut être munie de 28 structures différentiables différentes, S 8 de 2, S 31 de
plus de 16 millions. En revanche Rn , S 1 , S 2 , S 4 , S 5 et S 6 ont une unique structure
différentiable. Pour S 3 , le problème reste ouvert.
— On peut définir de même la C k -compatibilité en demandant que φ ◦ ψ −1 soit un
difféomorphisme de classe C k . On obtient alors des variétés différentiables de classe
Ck.
1.1 Variétés différentiables 3

Commentaires sur les propriétés topologiques de M .

Espace séparé : sur un espace non séparé il n’existe pas de métrique, puisque tout
espace muni d’une distance est séparé.
De même un sous-espace compact n’est pas forcément fermé et l’image d’un com-
pact par une application continue n’est pas toujours compact. C’est pour avoir
ce type de propriété que l’on impose à une variété d’être un espace séparé. En
revanche, tout sous-espace d’un espace topologique séparé est séparé. Comme on
s’intéressera généralement à des sous-ensemble d’un RN (pour sa topologie cano-
nique), cette condition sera automatiquement satisfaite.
Base dénombrable : la classe d’équivalence d’un atlas A peut être représentée par
son atlas maximal qui est l’ensemble de toutes les cartes compatibles avec celles
de A. On veut que la topologie définie par les domaines de ces cartes ait une base
dénombrable.
Cette hypothèse est importante pour éviter un certain nombre de pathologies des
espaces topologiques : sans elle, il est par exemple possible de munir Rn d’une
topologie qui le rende homéomorphe à un Rk muni de la topologie canonique, pour
k < n quelconque. Elle n’est en revanche pas indispensable : définissons une variété
généralisée de la même façon qu’une variété mais sans exiger que sa topologie soit
à base dénombrable (beaucoup d’auteurs définissent ainsi les variétés) ; tous les
résultats de ce polycopié sont valables pour les variétés généralisées (à l’exception
notable du théorème 1.8 du plongement de Whitney).
Notons enfin qu’une variété généralisée connexe M est une variété (i.e. a une
base dénombrable) si et seulement si elle est métrisable, c’est-à-dire si elle peut
être équipée d’une distance telle que la topologie induite sur M par la distance
coı̈ncide avec la topologie originale de M (voir par exemple [5]). Nous utiliserons
ce point dans la preuve du théorème de l’orbite au § 4.3).

Exemples de variétés différentiables.

1. Rn est une variété différentiable de dimension n pour l’atlas à une seule carte
(Rn , id).
2. Tout R-espace vectoriel E de dimension n est une variété de même dimension :
tout isomorphisme φ : E → Rn définit un atlas (E, φ). De même tout ouvert
U ⊂ E de l’espace vectoriel est également une variété, l’atlas étant (U, φ).
3. L’espace euclidien En est une variété de dimension n : il est en bijection avec Rn
via le choix d’un système de coordonnées x. L’atlas à une carte (En , x) y définit
donc un structure différentiable.
Tous ces exemples sont triviaux puisqu’il s’agit d’espaces homéomorphes à Rn .
4. Le cercle S 1 ⊂ R2 , muni de la topologie induite, est une variété de dimension 1 :
cependant il n’est pas homéomorphe à R (puisque S 1 est compact). Une seule carte
4 Chap. 1 – Variétés différentiables

ne sera donc pas suffisante pour créer un atlas. On définit deux cartes (U1 , φ1 ) et
(U2 , φ2 ) :

U1 = S 1 \ {(1, 0)} U2 = S 1 \ {(−1, 0)}

φ1 : U1 → ]0, 2π[ φ2 : U2 → ] − π, π[
(cos θ, sin θ) 7→ θ (cos θ, sin θ) 7→ θ

Les domaines de ces cartes recouvrent clairement le cercle : U1 ∪ U2 = S 1 . De

plus φ1 ◦ φ−12 est un difféomorphisme, ce qui montre que les deux cartes sont
compatibles.
Ainsi {(U1 , φ1 ), (U2 , φ2 )} est un atlas et définit une structure différentiable sur S 1 .
5. La sphère S 2 ⊂ R3 est une variété de dimension 2 : on peut construire un atlas en
utilisant la projection stéréographique .
Les points S et N désignant respectivement les pôles sud et nord, on considère les
ouverts UN = S 2 \ N et US = S 2 \ S et les applications

φN : UN → R2 φS : US → R2
x y x y
p = (x, y, z) 7→ ( 1−z , 1−z ) (x, y, z) 7→ ( 1+z , 1+z )

Exercice 1.1. Montrer que les cartes (UN , φN ) et (US , φS ) forment un atlas de S 2 .

6. La sphère S n ⊂ Rn+1 est une variété de dimension n : pour y définir une structure
différentiable, on peut utiliser soit la projection stéréographique (2 cartes), soit les
projections sur les hyperplans {xi = 0} (2n + 2 cartes).

Opérations Tout sous-ensemble ouvert Ω d’une variété différentiable M est lui-même

une variété différentiable. Sa structure différentiable est définie par la restriction à Ω d’un
atlas de M , c’est-à-dire par l’atlas AΩ = {(Uα ∩ Ω, φα |Uα ∩Ω )}, où A = {(Uα , φα )} est un
atlas de M . On dit parfois que Ω est une sous-variété ouverte de M .
Exemple. L’ensemble des matrices inversibles GLn (R) est une variété en tant que sous-
variété ouverte de Mn (R) (qui est un espace vectoriel donc une variété).

Lemme 1.1. Soient M et N des variétés différentiables respectivement de dimension n

et k et d’atlas {(Uα , φα )} et {(Vβ , ψβ )}.
Alors l’espace produit M × N est une variété de dimension n + k dont la structure
différentiable est définie par l’atlas formé de toutes les cartes de la forme {(Uα × Vβ , φα ×
ψβ )}, où (φα × ψβ )(p, q) = (φα (p), ψβ (q)) ∈ Rn+k .

La preuve est laissée en exercice.

Exemples.
1.2 Applications différentiables 5

— Le tore T2 = S 1 × S 1 est une variété, de même que le tore plat de dimension n,

Tn = S 1 × · · · × S 1 .
Ce tore intervient en robotique en tant qu’espace des configurations d’un robot
bras à n articulations planes.
— Le cylindre R × S 1 et le cylindre de dimension n + 1, R × S n sont des variétés.
— L’espace des configurations d’une voiture R2 ×S 1 ×S 1 est une variété (R2 paramètre
la position planaire du centre de gravité de la voiture, un S 1 paramètre la direction
de la voiture dans le plan, et l’autre S 1 l’angle des roues).

1.2 Applications différentiables

On connaı̂t les notions de différentiabilité et de difféomorphismes pour les applications
entre espaces vectoriels normés. On va définir ces notions pour les applications entre
variétés. Le principe est toujours le même : on dira qu’une application entre variétés est
différentiable (ou est un difféomorphisme) si, lue dans une carte, elle l’est. Formalisons
cette définition.
Soient M et N des variétés différentiables de dimension n et k et F : M → N une
application. Si (U, φ) est une carte de M contenant p et (V, ψ) une carte de N contenant
F (p), avec F (U ) ⊂ V , on dit que
F φψ = ψ ◦ F ◦ φ−1 : φ(U ) ⊂ Rn → ψ(V ) ⊂ Rk
est l’application F lue dans les cartes (U, φ) et (V, ψ). Dans le cas particulier où M ou
N est égal à Rn , l’application de carte correspondante est l’identité et on note
g φ = g φid = g ◦ φ−1 : φ(U ) ⊂ Rn → V ⊂ Rk
si g : M → Rk et
hψ = hidψ = ψ ◦ h : U ⊂ Rn → ψ(V ) ⊂ Rk
si h : Rn → N .
Définition. L’application F est différentiable (ou de classe C ∞ ) en p ∈ M si il existe une
carte (U, φ) de M contenant p et une carte (V, ψ) de N contenant F (p), avec F (U ) ⊂ V ,
telles que F φψ est de classe C ∞ .
On dit que F est une application différentiable de M dans N si elle est différentiable
en tout point p ∈ M .
Cette définition est correcte car la notion de différentiabilité ne dépend pas des cartes
choisies dans les variétés. En effet, si on choisit deux systèmes de coordonnées locales
différents φ1 , φ2 (resp. ψ1 , ψ2 ) sur M (resp. N ), on a
ψ2 ◦ F ◦ φ−1 −1 −1 −1
2 = ψ2 ◦ ψ1 ◦ (ψ1 ◦ F ◦ φ1 ) ◦ φ1 ◦ φ2 ,

et les applications entre espaces vectoriels normés ψ2 ◦ ψ1−1 et φ1 ◦ φ−1 ∞

2 sont de classe C .

Donnons quelques propriétés des applications différentiables.

6 Chap. 1 – Variétés différentiables

— Toute application différentiable est continue.

sur tout ouvert de carte U , on a F |U = ψ −1 ◦ (ψ ◦ F ◦ φ−1 ) ◦ φ.)
(Car S
— Soit i∈I Ui un recouvrement ouvert de M . Alors F est différentiable si et seule-
ment si chaque restriction F |Ui , i ∈ Ui , l’est.
— La composition d’applications différentiables est différentiable.
En effet, soient F : M → M ′ , G : M ′ → M ′′ et G ◦ F : M → M ′′ . Alors, dans des
cartes de M , M ′ et M ′′ ,

φ′′ ◦ (G ◦ F ) ◦ φ−1 = (φ′′ ◦ G ◦ φ′−1 ) ◦ (φ′ ◦ F ◦ φ−1 ).

Définition. Une application F : M → N est un difféomorphisme de M sur N si F est une

bijection et si F et F −1 sont différentiables. On a alors nécessairement dim M = dim N .

Notons que (U, φ), avec U ouvert de M , est une carte de la variété si et seulement si
φ est un difféomorphisme de U sur φ(U ) ⊂ Rn .

Rang d’une application.

Rappel : soit f : Rn → Rk une application dérivable en x ∈ Rn . Le rang de f en x est
défini comme le rang de l’application linéaire Df (x) (c’est-à-dire dim ImDf (x) ou encore
n − dim ker Df (x)).
Soient F : M → N une application différentiable et p ∈ M .

Proposition 1.2 (et définition). Le rang de ψ ◦ F ◦ φ−1 en φ(p) ne dépend pas des cartes
(U, φ) de M et (V, ψ) de N telles que p ∈ U et F (p) ∈ V . Cette quantité est appelée le
rang de F en p et est notée rgp F .

Preuve. Soient (U ′ , φ′ ) et (V ′ , ψ ′ ) d’autres cartes. Sur l’intersection des domaines :

ψ ′ ◦ F ◦ φ′−1 = (ψ ′ ◦ ψ −1 ) ◦ (ψ ◦ F ◦ φ−1 ) ◦ (φ′ ◦ φ−1 )

et donc

D ψ ′ ◦ F ◦ φ′−1 (x) = D ψ ′ ◦ ψ −1 ◦ D ψ ◦ F ◦ φ−1 ◦ D φ′ ◦ φ−1 (x).

Comme D (ψ ′ ◦ ψ −1 ) et D (φ′ ◦ φ−1 ) sont des isomorphismes, on obtient

rgD ψ ′ ◦ F ◦ φ′−1 (φ′ (p)) = rgD ψ ◦ F ◦ φ−1 (φ(p)).

Lue en coordonnées locales, l’application F devient

F φψ (x1 , . . . , xn ) = (F 1 (x1 , . . . , xn ), . . . , F k (x1 , . . . , xn ))

1.2 Applications différentiables 7

et le rang de F en p est celui de la matrice jacobienne

 ∂F 1 ∂F 1

∂x1
· · · ∂xn
 .. .. 
 . . 
∂F k ∂F k
∂x1
··· ∂xn x(p)

Le rang permet de caractériser les difféomorphismes.

Lemme 1.3. Une application différentiable de M dans N est un difféomorphisme si et

seulement si elle est bijective et de rang n = dim M = dim N en tout point de M .

Preuve. La preuve est laissée en exercice (indication : appliquer le théorème d’inversion

locale à F φψ ).

Il y a d’autres classes importantes d’applications différentiables caractérisées par leur

rang.

Définition. — F : M → N est une immersion si F est différentiable et rgF = dim M

en tout point de M (autrement dit, en coordonnées locales DF φψ (x) est injective
pour tout x). Dans ce cas, on a nécessairement dim M ≤ dim N .
— F : M → N est une submersion si F est différentiable et rgF = dim N en tout
point de M (autrement dit, en coordonnées locales, DF φψ (x) est surjective pour
tout x). On a alors dim M ≥ dim N .

Une immersion n’est pas forcément injective.

Exemple. F : R → R2 , F (t) = (cos 2πt, sin 2πt), est une immersion mais n’est pas injective
puisque F (t + k) = F (t) pour tout entier k.
De même pour F (t) = 2 cos(t − π2 ), sin(2t − π) .

Considérons maintenant une immersion injective F : M → N . C’est une bijection de

M dans M f = F (M ). Si on utilise F pour munir M f d’une topologie et d’une structure
différentiable, F devient un difféomorphisme entre les variétés M et M f. Cependant la
structure différentiable et la topologie sur M ne dépendent que de M et de F . Il n’y a
f
donc aucune raison en général pour que M f soit un sous-espace de la variété N (pour la
topologie induite).
Exemple. L’application F : R → R2 , F (t) = 2 cos(2 cot t + π2 ), sin(4 cot t + π) n’est pas

un homéomorphisme de R dans F (R) muni de la topologie induite de R2 .

On introduit donc une nouvelle définition qui permet d’éviter ces problèmes.

Définition. On dit que F : M → N est un plongement si F est une immersion injective

et un homéomorphisme de M dans F (M ) pour la topologie induite.
8 Chap. 1 – Variétés différentiables

Remarquons qu’une immersion injective est déjà une bijection de M dans F (M ) et

est continue, car différentiable. Pour qu’elle soit de plus un homéomorphisme, il suffit
donc que F −1 soit continue sur F (M ).


Exemple. Les applications
F : R → R3 , F (t) = cos t, sin t, t ou G :]1, +∞[→ R2 ,
G(t) = 1t cos 2πt, 1t sin 2πt sont des plongements.

1.3 Sous-variétés
Voyons maintenant comment obtenir simplement des variétés.

Sous-variétés de Rn .

Définition. Un sous-ensemble N ⊂ Rn est une sous-variété de Rn de dimension k ≤ n

si, pour tout point x de N , il existe un ouvert Ux ⊂ Rn contenant x et un difféomorphisme
φ : Ux → φ(Ux ) ⊂ Rn tel que

φ(Ux ∩ N ) = φ(Ux ) ∩ Rk .

Autrement dit, une sous-variété de Rn est un sous-ensemble que l’on peut localement
redresser en un sous-espace vectoriel Rk .
Il est clair qu’une sous-variété de Rn est une variété, les (Ux ∩ N, φ|Ux ∩N ) formant un
atlas pour N . La structure différentiable et la topologie sont induites par celles de Rn .
Remarquons alors (le montrer) que l’inclusion de la variété N dans Rn est un plongement
(le terme sous-variété plongée est d’ailleurs souvent employé à la place de sous-variété).
De façon générale, l’image d’un plongement est une sous-variété.

Lemme 1.4. Soient U ⊂ Rk un ouvert et f : U → Rn un plongement. Alors N = f (U )

est une sous-variété de Rn de dimension k.

Preuve. Soient x0 ∈ U et f (x0 ) ∈ N = f (U ). Comme f est une immersion, sa

différentielle Df (x0 ) : Rk → Rn est de rang k. Quitte à modifier l’ordre des coordonnées
(x1 , . . . , xn ), on suppose que Df (x0 ) s’écrit par blocs
 
A
Df (x0 ) =
B

où A ∈ Mk (R) est une matrice inversible (c’est-à-dire que les dérivées partielles de f
par rapport aux k premières coordonnées sont linéairement indépendantes). On définit
l’application
φ : Rk × Rn−k → Rn
(x, z) 7→ f (x) + (0, z)
1.3 Sous-variétés 9

Cette application vérifie φ(x0 , 0) = f (x0 ) = y0 et sa différentielle en (x0 , 0) est inversible

puisqu’elle s’écrit  
A 0
Dφ(x0 , 0) =
B In−k
Il résulte alors du Théorème d’Inversion Locale qu’il existe un voisinage U(x0 ,0) ⊂ Rn de
(x0 , 0) et un voisinage Vy0 ⊂ Rn de y0 tels que φ restreinte à U(x0 ,0) est un difféomorphisme
sur Vy0 .
On voudrait alors que le difféomorphisme φ−1 permette de définir N comme une
sous-variété. Il faut pour cela qu’il vérifie

φ−1 (Vy0 ∩ N ) = φ−1 (Vy0 ) ∩ Rk = U(x0 ,0) ∩ Rk ,

c’est-à-dire que Vy0 ∩ f (U ) = φ U(x0 ,0) ∩ Rk .
Or ces deux ensembles ont pour définition

Vy0 ∩ f (U ) = {y ∈ Vy0 t.q. ∃x ∈ U, f (x) = y} ,

φ U(x0 ,0) ∩ Rk = y ∈ Vy0 t.q. ∃x ∈ U(x0 ,0) ∩ U, f (x) = y .



Utilisons maintenant que f est un plongement : f est injective (pas de croisement) et

f −1 est continue sur f (U ).
Il existe donc un voisinage Vy′0 ⊂ Vy0 de y0 tel que l’image de Vy′0 ∩ f (U ) par f −1 est
contenue dans U(x0 ,0) . En considérant maintenant le difféomorphisme φ−1 restreint à Vy′0 ,
il vérifie
φ−1 Vy′0 ∩ f (U ) = φ−1 (Vy′0 ) ∩ Rk ,

ce qui montre que N = f (U ) est une sous-variété de Rn .

Exercice 1.2. Soit U =] − π2 , π2 [×]0, 2π[⊂ R2 et f : U → R3 définie par

f (θ, ϕ) = cos θ cos ϕ, cos θ sin ϕ, sin θ .

Montrer que f (U ) (la sphère S 2 privée d’un demi-grand cercle) est une sous-variété de
R3 de dimension 2.

Enfin la façon la plus pratique d’obtenir des sous-variétés de Rn est de considérer des
pré-images de submersions.
Lemme 1.5. Soient F : Rn → Rn−k une application différentiable et y ∈ F (Rn ) ⊂
Rn−k . Si F est une submersion sur N = F −1 (y), alors N est une sous-variété de Rn de
dimension k.
Preuve. Soit x0 un point de N = F −1 (y). Quitte à changer l’ordre des coordonnées dans
Rk , on peut supposer que la différentielle de F en x0 s’écrit par blocs


DF (x0 ) = A B
10 Chap. 1 – Variétés différentiables

où A ∈ Mn−k (R) est une matrice inversible. Considérons alors l’application

φ : Rn → Rn

x 7→ F (x) − y, xn−k+1 , . . . , xn

La différentielle de φ en x0 s’écrit par blocs

 
A B
Dφ(x0 ) =
0 Ik

et est donc inversible. Il résulte alors du Théorème d’Inversion Locale qu’il existe un
voisinage Ux0 ⊂ Rn de x0 tel que l’application φ restreinte à Ux0 est un difféomorphisme.
De plus

φ(Ux0 ∩ N ) = φ {x ∈ Ux0 t.q. F (x) = y}

= {(0, . . . , 0, xn−k+1 , . . . , xn ) t.q. x ∈ Ux0 }
= φ(Ux0 ) ∩ Rk .

Ce lemme est très utilisé dans le cas suivant. Soient f 1 (x), . . . , f p (x) des fonctions de
classe C ∞ de Rn dans R et S la “surface” définie par les zéros de ces fonctions :

S = {x ∈ Rn t.q. f 1 (x) = · · · = f p (x) = 0}.

∂f i


Si, pour tout point x de S, la jacobienne ∂x j i,j est de rang p, alors S est une sous-variété

de Rn de dimension n − p.
Exemples.
— Soit f (x) = (x1 )2 + · · · + (xn )2 , x ∈ Rn . La sphère f −1 (r2 ) de rayon r est une
sous-variété de dimension n − 1 car la différentielle de f ne s’y annule pas.
2 2 2
— Soit f (x, y, z) = xa2 ± yb2 ± zc2 . La surface f −1 (1), qui peut être soit un ellipsoı̈de
soit un hyperboloı̈de à une ou deux nappes, est une sous-variété de dimension 2
de R3 .
— Soit f (x, y, z) = x2 + y 2 − z 2 . Le cône de sommet origine, f −1 (0), n’est pas une
sous-variété de R3 car la différentielle de f à l’origine est nulle. En revanche le
cône privé de l’origine est une sous-variété.

Sous-variétés de variétés.
Définition. Une partie N d’une variété M de dimension n est une sous-variété de M de
dimension k ≤ n si, pour tout point q de N , il existe une carte (U, φ) de M contenant x
telle que
φ(U ∩ N ) = φ(U ) ∩ Rk .
1.3 Sous-variétés 11

La carte (U, φ) est dite adaptée à N . Elle vérifie

φ(U ∩ N ) = (x1 , . . . , xn ) ∈ φ(U ) t.q. xk+1 = · · · = xn = 0 .



Les deux résultats suivants se montrent presque de la même façon que les lemmes 1.5
et 1.4 (il faut en plus passer par des cartes des variétés M et N ).

Lemme 1.6. Soient M , N des variétés de dimension n, k et F : N → M un plongement.

Alors W = F (N ) est une sous-variété de M de dimension k.

Lemme 1.7. Soient M , N des variétés de dimension n, k, F : M → N une submersion

et y ∈ F (M ). Alors W = F −1 (y) est une sous-variété de M de dimension n − k.

Exercice 1.3. Montrer que le groupe spécial linéaire SLn (R) (l’ensemble des matrices
(n × n) de déterminant égal à un) est une variété.

Exercice 1.4. Montrer que le groupe orthogonal

O(n) = M ∈ Mn (R) t.q. M T M = I



est une variété. Quelle est sa dimension ?

On aura besoin ultérieurement (§ 4) d’une notion moins forte que celle de sous-variété
pour caractériser les sous-ensembles d’une variété. On vient de voir dans le lemme 1.4
qu’une sous-variété peut être vue comme l’image d’un plongement. On introduit alors la
définition suivante, suggérée dans la discussion page 7.

Définition. Un sous-ensemble W d’une variété M est une sous-variété immergée de M

de dimension k ≤ n si il existe une immersion injective f : N → M , où N est une variété
de dimension k, dont l’image f (N ) est égale à W .

Remarques.
— Une sous-variété immergée peut aussi être définie comme une variété contenue
dans M telle que l’inclusion i : W → M est une immersion (alors que c’est un
plongement pour les vraies sous-variétés).
— Des ouverts suffisamment petits de W sont des “vraies” sous-variétés de M (nous
laissons la preuve de cette affirmation en exercice). En revanche W lui-même n’en
est pas forcément une, sa topologie étant en général plus forte que celle induite
par M (voir page 7).
12 Chap. 1 – Variétés différentiables

Des variétés aux sous-variétés. Soit M une variété et f : M → RN un plongement.

L’ensemble f (M ) est alors une sous-variété de RN . Les variétés M et f (M ) sont alors
difféomorphes, c’est-à-dire indistinguables du point de vue de la géométrie différentielle.
Inversement, est-il possible de plonger toute variété dans un espace RN ? De façon
équivalente, est-il possible de considérer n’importe quelle variété comme une sous-variété
d’un RN ? La réponse est oui, et on peut de plus préciser N .

Théorème 1.8 (Plongement de Whitney). Toute variété de dimension n admet un plon-

gement sur une sous-variété fermée de R2n+1 .

Admis.
Chapitre 2

Espaces tangent et cotangent

2.1 Vecteurs tangents

Le but de ce chapitre est de formaliser dans une variété la notion de direction de
déplacement ou encore de mouvements possibles à partir d’un point donné.

Dans Rn et En Regardons d’abord comment cette notion apparaı̂t dans un espace

vectoriel, Rn par exemple. Dans Rn , les directions possibles à partir d’un point x sont
tous les vecteurs de Rn . L’ensemble des déplacements est donc Rn × Rn , l’ensemble des
couples (x, v) formés d’un point de départ x et d’une direction v. Cette notation est
cependant trompeuse : dans le produit Rn × Rn , les deux Rn jouent des rôles différents
puisque l’un est un espace de points et l’autre un espace de directions. La distinction
entre ces deux rôles est également mise en évidence dans l’expression de la différentielle
d’une application f : Rn → Rp :

Df : Rn × Rn → Rp
(x, v) 7→ Df (x)[v]

Le cas de l’espace euclidien En est plus intéressant. Un mouvement est considéré ici
comme une courbe M (t) dans En . La vitesse en M0 = M (0) est un vecteur Ṁ (0) : ce
n’est donc pas un élément de En , qui est un espace de points. En fait, si on regarde tous
les mouvements possibles M (t) à partir de M0 , leurs tangentes Ṁ (0) forment un espace
−−→
vectoriel de dimension n. Il s’agit de l’ensemble des segments orientés M0 P .
Notons TM0 En cet espace vectoriel et appelons-le l’espace tangent à En en M0 . À
chaque point de En correspond un espace tangent différent, qui sont isomorphes deux à
deux puisque de même dimension. De plus, et c’est propre à En et aux espaces affines, il
y a un isomorphisme naturel de TM En dans TM ′ En induit par la translation de M à M ′ .
Cette propriété ne sera plus vraie pour les variétés.
14 Chap. 2 – Espaces tangent et cotangent

Dans une variété Considérons maintenant une variété différentiable M et un point p

de M . On s’intéresse aux courbes dans M qui sont différentiables et qui passent par p

c : ] − ε, ε[ → M , c(0) = p.
t 7→ c(t)

Définition. Deux courbes c1 et c2 sont tangentes au point p si c1 (0) = c2 (0) = p et si il

existe une carte locale (U, φ) telle que p ∈ U et
d
d

φ ◦ c1 (0) = φ ◦ c2 (0).
dt dt
La définition est indépendante de la carte choisie. En effet si (V, ψ) est une autre carte
autour de p, on a
d d
(ψ ◦ φ−1 ) ◦ (φ ◦ c1 ) (0)


ψ ◦ c1 (0) =
dt dt
d
= D(ψ ◦ φ−1 ) ◦


φ ◦ c1 (0)
dt
−1 d
d

= D(ψ ◦ φ ) ◦ φ ◦ c2 (0) = ψ ◦ c2 (0)
dt dt
On définit ainsi une relation d’équivalence (c’est-à-dire une relation qui est transitive,
symétrique et réflexive) sur l’ensemble des courbes passant par p : c1 ∼ c2 si elles sont
tangentes en p.
Définition. Un vecteur tangent à M en p est une classe d’équivalence de courbes tan-
gentes en p.
L’espace tangent à M en p, noté Tp M , est l’ensemble des vecteurs tangents à M en p.

Exemple. Dans Rn , il est clair que deux courbes c1 , c2 sont tangentes au point x dès que
ċ1 (0) = ċ2 (0). Il y a donc un isomorphisme canonique entre l’ensemble des classes de
courbes tangentes Tx Rn et l’ensemble des directions ċ(0). Ce qui est propre à Rn c’est
que cet isomorphisme ne dépend pas du point x.

On peut montrer que Tp M est un espace vectoriel en utilisant une carte. La structure
vectorielle n’apparaı̂t cependant pas clairement. De plus la définition de Tp M fait inter-
venir un espace très gros, l’ensemble des courbes passant par p, qui n’est pas aisé à mani-
puler. Nous allons voir maintenant qu’on peut donner une autre définition – équivalente
– des vecteurs tangents qui résoudra ces difficultés.

2.2 Dérivations
Considérons l’ensemble des fonctions à valeurs réelles, de classe C ∞ , définies sur un
ouvert de M contenant un voisinage de p, dans lequel on identifie les fonctions qui sont
2.2 Dérivations 15

égales sur un voisinage de p (on obtient ainsi des germes de fonction). On note C ∞ (p)
cet ensemble. Notons que c’est une algèbre, c’est-à-dire un espace vectoriel muni d’une
opération interne (la multiplication).
Sur cet ensemble de fonctions on définit des opérateurs.
Définition. Une dérivation en p est une application linéaire Dp : C ∞ (p) → R qui vérifie
la règle de Leibniz. Autrement dit, Dp est une dérivation si, pour tous réels α, β et toutes
fonctions f , g dans C ∞ (p),
(i) Dp · (αf + βg) = αDp · f + βDp · g (linéarité),
(ii) Dp · (f g) = g(p)Dp · f + f (p)Dp · g (Leibniz).
L’ensemble D(p) des dérivations en p forme un espace vectoriel pour les opérations :
(Dp + Dp′ ) · f = Dp · f + Dp′ · f


Dp · (αf ) = αDp · f

Remarque. Toute dérivation vérifie Dp · cte = 0 (car Dp · 1 = Dp · (1 × 1) = 0).

Nous allons montrer que l’espace vectoriel tangent Tp M s’identifie à D(p). La première
étape est de déterminer la dimension de l’espace vectoriel D(p). Nous avons besoin pour
cela du résultat suivant.
Lemme 2.1 (Lemme d’Hadamard). Soit (U, φ), φ = (x1 , . . . , xn ), une carte de M centrée
en p. Pour toute fonction g ∈ C ∞ (p), il existe χ1 , . . . , χn ∈ C ∞ (p) telles que
n
X
g = g(p) + xi χ i
i=1
Pn
(autrement dit g(q) = g(p) + i=1 xi (q)χi (q) pour tout q ∈ U ).
Remarque. Dire qu’une carte (U, φ) est centrée en p signifie simplement que φ(p) = 0.

Preuve. Quitte à réduire U , on suppose φ(U ) convexe (par exemple une boule). On
travaille sur la fonction g lue dans la carte φ, c’est-à-dire g φ = g ◦ φ−1 . C’est une fonction
de classe C ∞ de φ(U ) ⊂ Rn dans R.
On peut calculer g φ en utilisant la formule de Taylor avec reste intégral :
Z 1
φ 1 n φ d φ 1
g (tx , . . . , txn ) dt

g (x , . . . , x ) − g (0) =
0 dt
n
X Z 1 ∂g φ
= xi i
(tx)dt.
i=1 0 ∂x

∂g φ
R1
Posons alors χφi (x) = 0 ∂xi
(tx)dt. On a obtenu
n
X
g φ (x) = g φ (0) + xi χφi (x).
i=1
16 Chap. 2 – Espaces tangent et cotangent

Comme g = g φ ◦ φ et φ(p) = 0, il suffit de composer la formule ci-dessus par φ pour

achever la démonstration.

Utilisons ce lemme pour caractériser les éléments de D(p). Fixons une carte (U, φ)
centrée en p. Une dérivation s’écrit alors
n

X
χi (p)Dp · xi + xi (p)Dp · χi


Dp · g = Dp · g(p) +
i=1
n
X
= χi (p)Dp · xi . (2.1)
i=1

Ainsi la donnée de Dp est équivalente à la donnée des réels Dp · xi , i = 1, . . . , n.

Lemme 2.2. dim D(p) = n = dim M .

Remarque. Ce lemme montre que, sur Rn , toute dérivation est une dérivée directionnelle.
En effet, à toute direction v ∈ Rn est associée une dérivation en x

1

g 7→ lim g(x + tv) − g(x) , (2.2)
t→0 t

qui est la dérivée directionnelle en x dans la direction v. L’ensemble des dérivées direc-
tionnelles en x est ainsi un sous-espace vectoriel de D(x) de dimension n et est donc égal
à D(x).

Faisons maintenant le lien entre les dérivations et les vecteurs tangents.

Proposition 2.3. Soient g ∈ C ∞ (p) et Xp un vecteur tangent en p. Alors la dérivée

d
dt
(g ◦ c)(0) est la même pour toutes les courbes c(s) passant par p et appartenant à la
classe d’équivalence Xp .

Preuve. Choisissons des coordonnées locales φ et écrivons g ◦ c comme la composée de

g φ = g ◦ φ−1 avec cφ = φ ◦ c. On obtient alors

d
d
(g ◦ c)(0) = D g ◦ φ−1 ◦ (φ ◦ c)(0).
dt dt
d
La proposition résulte alors du fait que dt
(φ◦c)(0) ne dépend que de la classe d’équivalence
Xp .

On note Xp · g la valeur de cette dérivée.

Proposition 2.4. L’application g 7→ Xp · g est une dérivation.

2.2 Dérivations 17

Preuve. Il suffit de choisir des coordonnées locales φ et un représentant c(t) de la classe

Xp . La linéarité et Leibniz se déduisent aisément de l’expression
d φ φ

Xp · g = g ◦ c (0).
dt

Cette dérivation est en fait une généralisation des dérivées directionnelles. En effet
on a vu que dans Rn un vecteur tangent vx est associé canoniquement à une direction
v = ċ(0). La dérivation g 7→ vp ·g est alors clairement égale à la dérivée directionnelle (2.2).
D’ailleurs on appellera parfois g 7→ Xp · g dérivée directionnelle de g dans la direction Xp .
Théorème 2.5. L’ensemble des vecteurs tangents Tp M s’identifie à l’espace vectoriel
D(p) de dimension n des dérivations en p.
Cette identification permet de définir une structure vectorielle sur l’espace tangent
Tp M (appelé également, en conséquence, espace vectoriel tangent). Notons que cette
structure vectorielle coı̈ncide avec celle que l’on peut obtenir à partir de la lecture dans
une carte.
Preuve. Soit Ψ : Tp M → D(p) qui à une vecteur tangent Xp fait correspondre la
dérivation définie par Xp · g. Ainsi, si c est un représentant de la classe d’équivalence
Xp = [c]p , Ψ(Xp ) est la dérivation Xp · g = dtd (g ◦ c)(0).
Montrons d’abord que Ψ est injective. Soient Xp = [c]p et Xp′ = [c′ ]p des vecteurs de
Tp M tels que Ψ(Xp ) = Ψ(Xp′ ), c’est-à-dire que, pour toute fonction g ∈ C ∞ , on a
d d
g ◦ c′ (0).



g ◦ c (0) =
dt dt
Fixons des coordonnées locales φ = (x1 , . . . , xn ). On a, pour chaque i,
d i
d i ′

x ◦ c (0) = x ◦ c (0),
dt dt
ce qui implique ċφ (0) = ċ′φ (0) et donc Xp = Xp′ . Ainsi l’application Ψ est injective.
Montrons qu’elle est aussi surjective. Soit Dp une dérivation. On a vu que, dans une
carte φ donnée, Dp est déterminée par les réels di = Dp · xi , i = 1, . . . , n. Considérons
alors la courbe c(t) = φ−1 ◦ (td1 , . . . , tdn ) et notons Xp son vecteur tangent en p.
On a
d i

X p · xi = x ◦ c (0) = di = Dp · xi ,
dt
ce qui montre que Dp = Ψ(Xp ) et donc que Ψ est surjective.

Dans la suite nous identifierons systématiquement Tp M et D(p). Ainsi le terme vec-

teur tangent désignera indifféremment la classe d’équivalence de courbes tangentes ou la
dérivation associée alors que l’espace tangent Tp M sera employé à la place de D(p).
18 Chap. 2 – Espaces tangent et cotangent

2.3 Différentielle d’une application

Soient M et N des variétés différentiables de dimension n et k et F : M → N une
application différentiable. Si g : N → R est une fonction sur N , F permet de lui faire
correspondre une fonction sur M , F ∗ g = g ◦ F , appelée image réciproque de g par F . On
définit ainsi l’application :

F ∗ : C ∞ F (p) → C ∞ (p)

g 7→ F ∗ g = g ◦ F

Remarquons que le sens de F ∗ est l’inverse de celui de F .

Proposition 2.6 (et définition). L’application dFp : Tp M → TF (p) N définie par

dFp (Xp ) · g = Xp · (F ∗ g), ∀g ∈ C ∞ F (p) ,

est linéaire.
On l’appelle la différentielle de F en p (ou encore application linéaire tangente à F en p).

Remarque. L’élément dFp (Xp ) de TF (p) N est donc un vecteur tangent à N . On peut
également le caractériser de la façon suivante : soit la courbe c(s) ⊂ M passant par p de
vecteur tangent Xp en p. Alors dFp (Xp ) est le vecteur tangent en F (p) de la courbe F ◦ c
incluse dans N .

Preuve. La linéarité de dFp découle immédiatement de celle de Xp . En effet

dFp λXp + µYp · g = λXp + µYp · (F ∗ g)

= λXp · (F ∗ g) + µYp · (F ∗ g)



= λdFp Xp · g + µdFp Yp · g.

La définition de la différentielle ne fait appel qu’à des propriétés locales de la variété

et de l’application. Comme localement une variété est difféomorphe à un ouvert de Rn ,
toute les propriétés des applications différentiables dans les espaces vectoriels normés sont
vraies localement pour les applications différentiables sur les variétés (c’est le principe de
base pour obtenir des résultats locaux dans les variétés). Citons les plus importantes de
ces propriétés.

Théorème 2.7 (Théorème de composition). Soient F : M → N une application différentiable

en p ∈ M et G : N → W une application différentiable en F (p) ∈ N . Alors G ◦ F est
différentiable en p et
d(G ◦ F )p = dGF (p) ◦ dFp
2.4 Coordonnées sur l’espace tangent 19

Corollaire 2.8. Si F : M → N est un difféomorphisme, alors, pour tout p ∈ M , dFp est

un isomorphisme.
La réciproque à ce corollaire n’est vraie que localement. Elle nécessite d’introduire la
notion de difféomorphisme local.
Définition. Une application F : M → N est un difféomorphisme local en p s’il existe un
voisinage U ⊂ M de p et un voisinage V ⊂ N de F (p) tels que l’application F |U : U → V
est un difféomorphisme.

Théorème 2.9 (Théorème d’inversion locale). Soit F : M → N une application différentiable

en p ∈ M telle que dFp : Tp M → TF (p) N est un isomorphisme.
Alors F est un difféomorphisme local en p. De plus la réciproque du difféomorphisme
F |U a pour différentielle

−1
d F |−1


U F (p) = dF p .
Ce théorème a une conséquence importante pour la détermination de coordonnées lo-
cales. En effet, un système de coordonnées locales n’est rien d’autre qu’un difféomorphisme
local de M dans Rn . Ainsi une application différentiable φ : M → Rn définit des coor-
données locales en p si et seulement si dφp est un isomorphisme.
Notons enfin que l’on pourrait également énoncer le théorème des fonctions implicites
(mais nous n’en avons pas besoin dans ce cours).

2.4 Coordonnées sur l’espace tangent

En utilisant les différentielles des cartes, nous allons étendre le calcul en coordonnées
locales à l’espace tangent. On commence par traiter le cas de Rn .

Description de Tx Rn On a vu plusieurs propriétés de Tx Rn : il est canoniquement

isomorphe à Rn et il peut être identifié à l’ensemble des dérivées partielles en x. Il est
temps maintenant d’en donner une description facile à utiliser.
Soit x ∈ Rn . Considérons les dérivées partielles en x, c’est-à-dire les dérivations sur
n
R
∂ ∂g
i x
: g 7→ (x), i = 1, . . . , n.
∂x ∂xi
Ces dérivées partielles forment une base de l’ensemble des dérivées directionnelles en x et
donc une base de Tx Rn dite base naturelle. Ainsi tout vecteur tangent vx ∈ Tx Rn s’écrit
∂ n ∂
vx = v 1 + · · · + v .
∂x1 x ∂xn x
Ce vecteur est aussi la classe d’équivalence des courbes c(t) passant par x telles que
ċ(0) = (v 1 , . . . , v n ).
L’identification canonique Tx Rn ≃ Rn est définie alors comme l’isomorphisme vx 7→
(v , . . . , v n ).
1
20 Chap. 2 – Espaces tangent et cotangent

Coordonnées sur Tp M Soient p ∈ M et (U, φ), φ = (x1 , . . . , xn ), une carte de M dont

le domaine contient p. Comme φ est un difféomorphisme de U sur φ(U ), sa différentielle

−1
dφp : Tp M → Tφ(p) Rn est inversible (corollaire 2.8) et dφp = d(φ−1 )φ(p) est un
n
isomorphisme de Tφ(p) R sur Tp M .
Nous allons nous servir de cet isomorphisme pour définir des coordonnées sur Tp M .
Soient x = φ(p) ∈ Rn et ( ∂x∂ 1 x , . . . , ∂x∂n x ) la base canonique de Tx Rn . L’image de cette
base définit des vecteurs tangents à M en p, que l’on note

∂ −1

 ∂ 
=d φ φ(p) ∂xi x
.
∂xi p

Cestangents ( ∂x∂ 1 p , . . . , ∂x∂n p ) forment une base de Tp M , appelée base naturelle associée
aux coordonnées locales φ.
Remarquons que, si g ∈ C ∞ (p), on a

∂  ∂  ∂ ∂g φ
· g = d φ−1 · (g ◦ φ−1 ) =



φ(p) ∂xi x
· g = x
φ(p) ,
∂xi p ∂xi ∂xi

où g φ = g ◦ φ−1 est la fonction g lue dans les coordonnées xi .

Lemme 2.10. Dans la base naturelle associée aux coordonnées locales φ, un vecteur
tangent Xp ∈ Tp M s’écrit

∂ ∂
Xp = X 1 + · · · + Xn , avec X i = Xp · xi .
∂x1 p ∂xn p
Remarque. Comparer avec la formule 2.1, page 16.

∂xi i
Pn ∂x
Preuve. Il suffit de remarquer que ∂xj
= δij et que Xp · xi = j=1 X j ∂x j.

Exercice 2.1. Considérons dans R2 la courbe c(t) = x0 cos t − y0 sin t, x0 sin t + y0 cos t .
Calculer son vecteur tangent. L’écrire également en coordonnés polaires.

Changement de coordonnées. Le lemme ci-dessus permet de trouver assez simple-

ment les formules de changement de coordonnées dans l’espace tangent. Soient φ =
(x1 , . . . , xn ) et ψ = (y 1 , . . . , y n ) des coordonnées locales autour de p et Xp un vec-
teur tangent
Pn en p de coordonnées X 1 , . . . , X n dans la base naturelle associée à φ :
i ∂
Xp = i=1 X ∂xi p . On note

ψ ◦ φ−1 = y 1 (x), . . . , y n (x)

l’application ψ lue dans la carte φ.

2.4 Coordonnées sur l’espace tangent 21

Dans la base naturelle associée à ψ, Xp a pour coordonnées Y j = Xp · y j , j = 1, . . . , n.

En écrivant tout dans les coordonnées φ, on obtient
n j φ

n j
i ∂ (y ) i ∂y (x)
X
X
j
Y = X φ(p) = X
i=1
∂xi i=1
∂xi x=φ(p)

ou encore, en écriture matricielle,

   
Y1  i  X1
 ..  ∂y (x)  .. 
 . = ∂xj  . . (2.3)
n φ(p)
Y 1≤i,j≤n Xn

On peut enfin calculer la différentielle d’une application F en coordonnées locales.

Soit F : M → N une application différentiable en p ∈ M . Choisissons φ = (x1 , . . . , xn )
des coordonnées locales sur M centrées en Ppn et ψi = (y 1 , . . . , y k ) des coordonnées locales
sur N centrées en F (p). Enfin soit Xp = i=1 X ∂x∂ i p un élément de Tp M .
Pour tout g ∈ C ∞ F (p) , on écrit

n φ


X ∂ (g ◦ F )
dFp (Xp ) · g = Xp · (g ◦ F ) = Xi (0). (2.4)
i=1
∂xi

Écrivons (g ◦ F )φ = g ψ ◦ F φψ , où F φψ = ψ ◦ F ◦ φ−1 et g ψ = g ◦
ψ −1 sont les applications

F et g lues dans les cartes. Notons F φψ (x) = F 1 (x), . . . , F k (x) . Pour i = 1, . . . , n, on a

k k
∂ (g ◦ F )φ X ∂g ψ ∂F j X ∂F j ∂
(0) = (0) (0) = (0) · g.
∂xi j=1
∂y j ∂xi j=1
∂xi ∂y j F (p)

En insérant cette valeur dans (2.4), qui est valable pour tout g ∈ C ∞ F (p) , on obtient

l’expression suivante pour la différentielle de F en p :

k n
!
j
X X ∂F ∂
dFp (Xp ) = i
(0)X i .
j=1 i=1
∂x ∂y j F (p)

On peut réécrire cette expression de façon matricielle en faisant intervenir la matrice

jacobienne de F φψ .
Lemme 2.11. Si Xp = ni=1 X i ∂x∂ i p dans la base naturelle de Tp M associée aux co-
P
Pk j ∂
ordonnées locales φ, alors dFp (Xp ) = j=1 Y ∂y j p dans la base naturelle de TF (p) N
associée à ψ, avec    
Y1 X1
 ..  φψ  . 
 .  = JF (0)  ..  .
Yk Xn
22 Chap. 2 – Espaces tangent et cotangent

Cette formulation en coordonnées permet de calculer le rang de F (tel que défini par
la proposition 1.2).

Corollaire 2.12. Le rang de F en p est égal à la dimension de dFp (Tp M ).

Ceci nous donne d’ailleurs une définition intrinsèque du rang d’une application.

2.5 Cas des sous-variétés de Rn

L’espace tangent à une sous-variété est un sous-espace vectoriel de l’espace tangent
à la variété ambiante : si N est une sous-variété de M et p ∈ N , alors Tp N ⊂ Tp M .
L’espace tangent à une sous-variété de Rn est donc un sous-espace vectoriel de Tx Rn . Il
sera en fait plus pratique de le voir comme un sous-espace vectoriel de Rn en utilisant
l’identification canonique Tx Rn ≃ Rn . Sa caractérisation dépend de la façon dont a été
définie la sous-variété.

Cas général. Une sous-variété N de Rn est définie au voisinage de chacun de ses points
x par un difféomorphisme φ : Ux → φ(Ux ) ⊂ Rn tel que φ(Ux ∩ N ) = φ(Ux ) ∩ Rk .
En tant que sous-espace vectoriel de Rn , l’espace tangent à N est alors

Tx N = Dφ−1 φ(x) Rk × {0}Rn−k .

 

Plongement. Soit N une sous-variété de Rn définie comme N = F (U ), où F : U ⊂

Rk → Rn est un plongement. Alors, en tant que sous-espace vectoriel de Rn , son espace
tangent est
Tx N = DF F −1 (x) [Rk ].

Exemple. Considérons une courbe c(t) = (x1 (t), . . . , xn (t)) dans Rn , de classe C ∞ , qui
est un plongement sur un intervalle ouvert I ⊂ R. L’ensemble de la courbe C = c(I) est
alors une sous-variété. En t0 ∈ I, la différentielle Dc(t0 ) a pour matrice
 
ẋ1 (t0 )
ċ(t0 ) =  ... 
 
ẋn (t0 )

et l’espace tangent à C est Tc(t0 ) C = Rċ(t0 ).

Ainsi l’espace tangent correspond à l’espace vectoriel engendré par la tangente à la
courbe, comme on s’y attendait.
On peut également retrouver ce résultat en construisant une base de Tc(t0 ) C. Soit dtd |t0
la base de Tt0 R. Pour i = 1, . . . , n, on a
d d i

· xi = = ẋi (t0 ),

Dc(t0 ) t0
x (t) t0
dt dt
2.6 Espace cotangent* 23

et l’image de la base de Tt0 R est

n
d X ∂
ẋi (t0 )

Dc(t0 ) t0
= .
dt i=1
∂xi

C’est bien le vecteur vitesse (ou vecteur tangent) de la courbe.

Submersions Soit N ⊂ Rn une sous-variété définie comme N = F −1 (y), où F : U ⊂

Rn → Rn−k est une submersion.

Proposition 2.13. En tant que sous-espace vectoriel de Rn , l’espace tangent à la sous-

variété N est Tx N = ker DF (x).

Preuve. Par définition, Tx N est l’ensemble des vecteurs tangents ċ(0) pour toutes les
courbes différentiables c(t) tracées sur N passant par x. Une telle courbe vérifie en par-
ticulier F ◦ c(t) = y = cte, c’est-à-dire

d
 
F ◦ c (0) = DF (x) ċ(0) = 0.
dt
Ainsi Tx N est inclus dans ker DF (x). Comme ces espaces vectoriels ont même dimension,
ils sont égaux.

C’est la façon la plus pratique de déterminer les espaces tangents.

Exemple. La sphère de dimension n : S n = F −1 (1), où F : Rn+1 → R est la forme
quadratique F (x) = ⟨x, x⟩ = (x1 )2 + · · · + (xn+1 )2 , qui est une submersion. Comme
DF (x)[h] = 2⟨x, h⟩, l’espace tangent à la sphère est

Tx S n = ker DF (x) = (Rx)⊥ ,

c’est-à-dire l’ensemble des vecteurs orthogonaux à x.

Exercice 2.2. Déterminer les espaces tangents à O(n) et SO(n).

Exercice 2.3. Déterminer l’espace tangent à SLn (R).

2.6 Espace cotangent*

On définit maintenant le dual de l’espace tangent à une variété M .
24 Chap. 2 – Espaces tangent et cotangent

Définition. — Une 1-forme (ou covecteur) en p ∈ M est une forme linéaire sur Tp M ,
c’est-à-dire une application linéaire

ωp : Tp M → R
Xp 7→ ωp (Xp )

On note ωp (Xp ) = ⟨ωp , Xp ⟩, le crochet étant ici le crochet de dualité.

— L’espace cotangent à M en p, noté Tp∗ M , est l’espace vectoriel des 1-formes en p.
C’est l’espace vectoriel dual de Tp M (c’est-à-dire Tp∗ M = (Tp M )∗ ).

Exemple. Soit g : M → R une fonction différentiable sur M . Alors, en identifiant Tt R

avec R, la différentielle de g en p, dgp : Tp M → R, peut être vue comme une 1-forme. Cet
exemple justifie les notations ci-dessous en coordonnées.

Rappelons d’abord que, si e1 , . . . , en est la base d’un espace vectoriel V , il existe une
unique base duale e∗1 , . . . , e∗n du dual V ∗ telle que e∗i (ej ) = δij .
Considérons maintenant des coordonnées locales φ = (x1 , . . . , xn ) en p et dxip : Tp M →
R la différentielle de la i-ème coordonnée (on identifie de nouveau Tt R avec R). Par
définition, dxip (Xp ) = Xp · xi . En particulier, pour tout couple i, j,

∂ ∂
∂xi
⟨dxip , ⟩ = dx i
p = = δij .
∂xj p ∂xj p ∂xj p

Ainsi dx1p , . . . , dxnp est une base de Tp∗ M , duale de la base ∂

∂x1 p
, . . . , ∂x∂n p de Tp M .
Dans cette base, toute 1-forme de Tp∗ M s’écrit
n
X ∂
ωp = ωi dxip , où ωi = ωp ( ),
i=1
∂xi p
Pn
puisque ωp ( ∂x∂ i p ) = j=1 ωj ⟨dxjp , ∂x∂ i p ⟩). De même, tout vecteur de Tp M s’écrit
n
X ∂
Xp = Xi où X i = dxip (Xp ).
i=1
∂xi p

Les formules de changement de coordonnées pour les 1-formes s’obtiennent comme

les formules (2.3) pour les vecteurs tangents : on choisit des coordonnées locales φ =
(x1 , . . . , xn ) et ψ = (y 1 , . . . , y n ) et on note φ ◦ ψ −1 = x1 (y), . . . , xn (y) l’application φ

lue dans la carte ψ. Soit ωp une 1-forme, ω1 , . . . , ωn ses coordonnées dans la base de Tp∗ M
associée à φ et µ1 , . . . , µn ses coordonnées dans la base associée à ψ. En écrivant que
µj = ωp ( ∂y∂ j p ), on obtient, en écriture matricielle,

∂xi (y)
 



µ1 · · · µn = ω1 · · · ωn .
∂y j ψ(p) 1≤i,j≤n
2.6 Espace cotangent* 25

Remarque. Cette formule est l’inverse de celle pour les vecteurs tangents. En effet, en
notant A la matrice des dérivées partielles des y i par rapport à x (c’est-à-dire la matrice
qui intervient dans (2.3)), on obtient

Y = AX et µ = ωA−1 ,

X, Y étant des vecteurs colonnes et µ, ω des vecteurs lignes.

On a vu dans la section 2.3 qu’à une application différentiable correspond une application
linéaire entre espaces tangents, la différentielle. De la même façon, il lui correspond également
une application linéaire entre espaces cotangents. Soit F : M → N une application différentiable
et p ∈ M .

Définition. L’application réciproque de F en p est l’application linéaire F ∗ : TF∗ (p) N → Tp∗ M

telle que
F ∗ ωF (p) (Xp ) = ωF (p) dFp (Xp ) .

C’est l’application duale de dFp :

⟨ F ∗ ω p , Xp ⟩ = ⟨ωF (p) , dFp (Xp )⟩,

le membre de droite utilisant le crochet de dualité dans Tp M alors que celui de gauche utilise
celui dans TF (p) M .
26 Chap. 2 – Espaces tangent et cotangent
Chapitre 3

Champs de vecteurs

Jusqu’à maintenant nous avons vu les notions d’espace et de vecteur tangents en

un point, qui correspondent aux mouvements infinitésimaux à partir d’une configuration
donnée. En considérant toutes les configurations possibles, nous sommes maintenant en
mesure de définir des mouvements sur toute la variété : ces mouvements vont être intro-
duits comme des systèmes dynamiques, c’est-à-dire des équations différentielles. Du point
de vue des espaces de configurations, cela signifie que l’on caractérise les mouvements par
la donnée en tout point de la vitesse.
L’outil principal sera ainsi la notion de champ de vitesse, ou plutôt de champ de
vecteurs : il s’agit d’une application assignant à chaque point p de la variété un vecteur
Xp de l’espace tangent. Avant d’aller plus loin dans la définition, il faut préciser l’espace
auquel appartient l’image de cette application.

3.1 Fibrés tangent et cotangent

On s’intéresse à l’ensemble de tous les vecteurs tangents en tous les points de la variété.
Définition. L’ensemble T M = {(p, Xp ), p ∈ M, Xp ∈ Tp M } est appelé le fibré tangent
de la variété M .
De même, l’ensemble T ∗ M = {(p, ωp ), p ∈ M, ωp ∈ Tp∗ M } est appelé le fibré cotangent
de M .
Le fibré tangent est l’union des espaces tangents
[ [
TM = {p} × Tp M ou simplement Tp M,
p∈M p∈M

mais il faut bien préciser que cette union est disjointe : on ne peut pas additionner des
éléments Xp et Yp′ appartenant à des espaces tangents différents.
On appelle projection canonique sur T M la projection
π: TM → M
(p, Xp ) 7→ p
28 Chap. 3 – Champs de vecteurs

et la fibre au-dessus de p la pré-image π −1 (p) = {p} × Tp M d’un point p.

Théorème 3.1. Le fibré tangent (resp. cotangent) a un structure naturelle de variété
différentiable de dimension 2n.
Preuve. Soit (U, φ) une carte de M et π −1 (U ) l’ensemble des fibres au-dessus d’un point
de U , i.e. π −1 (U ) = {p} × Tp M, p ∈ U . L’application
Φ : π −1 (U ) ⊂ T M → φ(U ) × Rn ⊂ R2n
(p, Xp ) 7→ (φ(p), dφp (Xp ))
est une bijection d’après le corollaire 2.8 (notons que l’on a utilisé ici l’identification
Tx Rn ≃ Rn puisque dφp (Xp ) ∈ Tφ(p) Rn est identifié à ses coordonnées dans la base
naturelle).
Nous allons montrer que les couples (π −1 (U ), Φ) forment un atlas de T M et définissent
donc une structure différentiable. La preuve se fait en deux étapes.
(i) Il faut d’abord munir T M d’une topologie : on la définit en posant que les appli-
cations Φ (pour toutes les cartes de M ) sont des homéomorphismes (c’est-à-dire
que les ouverts de T M sont les parties W ⊂ T M telles que Φ(W ∩ π −1 (U )) est
ouvert dans R2n ). Pour que cela soit possible, il suffit de montrer que, si Φ et Ψ
sont les applications correspondant à des cartes (U, φ) et (V, ψ), alors
Ψ ◦ Φ−1 : φ(U ∩ V ) × Rn → ψ(U ∩ V ) × Rn
est un homéomorphisme. Or, sur son ensemble de définition,
Ψ ◦ Φ−1 = ψ ◦ φ−1 , D(ψ ◦ φ−1 ) ,

et est donc un difféomorphisme (pour la structure différentiable de R2n ). C’est

donc a fortiori un homéomorphisme.
(ii) L’expression ci-dessus montre de plus que les couples (π −1 (U ), Φ) sont des cartes
compatibles entre elles. Comme elles recouvrent T M elles forment donc un atlas
et définissent une structure différentiable.
On montre de même que T ∗ M est une variété en utilisant les applications
Φ′ : (p, ωp ) 7→ (φ(p), (φ−1 )∗ ωp ).

La projection canonique π : T M → M apparaı̂t maintenant comme une application

différentiable. C’est de plus une submersion surjective.
Remarquons également que l’on peut définir un difféomorphisme h : π −1 (U ) ⊂ T M →
U × Rn en posant h(p, Xp ) = (p, dφp (Xp )), qui est de plus linéaire sur les fibres. Ainsi
le fibré tangent est difféomorphe localement à un produit : π −1 (U ) ≃ U × Rn (h est ce
qu’on appelle une trivialisation locale).
En revanche, en général T M n’est pas globalement trivial, c’est-à-dire qu’il n’existe
pas de difféomorphisme de T M dans M × Rn linéaire le long des fibres.
3.1 Fibrés tangent et cotangent 29

Exemples.
— Le fibré tangent à Rn admet une trivialisation globale T Rn ≃ Rn × Rn via l’iden-
tification canonique Tx Rn ≃ Rn .
— Le fibré tangent au cercle S 1 admet une trivialisation globale car il est difféomorphe
au cylindre : T S 1 ≃ S 1 × R. En revanche le fibré tangent T S 2 n’admet pas de
trivialisation globale.
— Est-ce que O(n) (ou SO(n)) admet une trivialisation globale ?

Remarquons que l’application Φ permet de définir des coordonnées locales sur T M en

utilisant les coordonnées de l’espace tangent construites dans la section 2.4. En effet, si
(U, φ) est une carte de M , avec φ = (x1 , . . . , xn ), le difféomorphisme Φ sur π −1 (U ) s’écrit

Φ(p, Xp ) = (x1 , . . . , xn , X 1 , . . . , X n ) ∈ R2n ,

avec (x1 , . . . , xn ) = φ(p) et X i = Xp · xi . Ce sont bien des coordonnées locales sur T M .

En s’appuyant sur la construction de Φ, on peut prolonger une application entre deux
variétés en un application entre leurs fibrés tangents.

Définition. Soit F : M → N une application différentiable. On définit le prolongement

(ou différentielle) de F comme l’application différentiable :

dF : TM → TN
(p, Xp ) 7→ (F (p), dFp (Xp ))

Il est clair que Φ est le prolongement de l’application de carte φ : Φ = dφ. Le prolon-

gement vérifie de plus les propriétés suivantes
— Le diagramme ci-dessous est commutatif :

dF
TM −→ TN
π ↓ ↓π
F
M −→ N

c’est-à-dire F ◦ π = π ◦ dF .
— La restriction de dF aux fibres est linéaire car dF Tp M = dFp .
— Si F : M → W et G : W → N sont des applications différentiables, alors d(G◦F ) =
dG ◦ dF .
— Si F : M → N est un difféomorphisme, alors le prolongement dF : T M → T N est
un difféomorphisme également et (dF )−1 = d(F −1 ).
La preuve des deux dernières propriétés est laissée en exercice : il s’agit essentiellement
de passer dans des cartes et d’utiliser la preuve du théorème 3.1.
30 Chap. 3 – Champs de vecteurs

3.2 Champs de vecteurs

Considérons une variété M de dimension n.
Définition. Un champ de vecteurs différentiable (ou champ de vecteurs) sur M est une
application différentiable X : M → T M qui, à un point p ∈ M , associe un couple formé
de p et d’un vecteur tangent à M en p : X(p) = (p, Xp ). Autrement dit, π ◦ X = idM .

On notera X (M ) l’ensemble de tous les champs de vecteurs sur M .

De même qu’un vecteur tangent en p définit une dérivation sur l’ensemble des germes
C (p), un champ de vecteur définit une dérivation sur l’ensemble C ∞ (M ) des fonctions
∞

de M dans R de classe C ∞ . En effet l’application

C ∞ (M ) → C ∞ (M )
g 7→ X · g, t.q. X · g(p) = Xp · g

définie par un champ de vecteurs X est linéaire et vérifie la règle de Leibniz. L’ensemble
X (M ) s’identifiera donc avec l’espace vectoriel de dimension infinie des dérivations sur
C ∞ (M ).

Exemple (Champs de vecteurs sur Rn ). En utilisant la base naturelle de Tx Rn , tout champ

de vecteurs sur Rn s’écrit comme
n
 X
i ∂ 
X(x) = x, X (x) i x , (3.1)
i=1
∂x

où X i est une fonction différentiable de Rn dans R. En particulier, pour i = 1, . . . , n,

l’application
∂  ∂ 
: x →
7 x,
∂xi ∂xi x
est un champ de vecteur sur Rn . Bien entendu, en tant que dérivation, ce champ est tout
∂g
simplement la dérivée partielle par rapport à xi , c’est-à-dire l’application g 7→ ∂xi . On

pourra donc réécrire le champ de vecteur X donné en (3.1) comme la dérivation

n
∂ 1 n ∂
X ∂g
X=X 1
+ · · · + X n
: g →
7 Xi i .
∂x ∂x i=1
∂x

Enfin un champ de vecteurs sur Rn peut aussi être considéré comme une application
de Rn → Rn en l’identifiant à x 7→ (X 1 (x), . . . , X n (x)).

Considérons maintenant une carte (U, φ) de M , avec φ = (x1 , . . . , xn ). Sur le domaine

U de la carte, tout champ de vecteurs X sur M s’écrit
n
 X ∂ 
X : p ∈ U 7→ p, X i (p) i p , (3.2)
i=1
∂x
3.2 Champs de vecteurs 31

où X i = X · xi ∈ C ∞ (M ).
En particulier, comme U est lui-même une variété, l’application de U sur π −1 (U ) :

∂ M −1 ∂  ∂ 
= dφ ◦ ◦ φ : p →
7 p, i p ∈ Tp M = Tp U
∂xi ∂xi ∂x
définit un champ de vecteurs sur U . On l’appelle champ de coordonnées.
En utilisant les champs de coordonnées, le champ X donné par (3.2) s’écrit également,
en tant que dérivation sur C ∞ (M ),
n
∂ M1 n ∂
M X
i ∂
M
X=X + ··· + X : g 7→ X · g.
∂x1 ∂xn i=1
∂xi

Définissons le champ X lu dans la carte φ comme le champ de vecteur sur φ(U ) ⊂ Rn

∂ ∂
X φ (x) = dφ ◦ X ◦ φ−1 (x) = X 1 (x) 1
+ · · · + X n (x) n .
∂x ∂x

Considérons alors d’autres coordonnées locales ψ = (y 1 , . . . , y n ), dans lesquelles le champ

X se lit
∂ ∂
X ψ (y) = Y 1 (y) 1 + · · · + Y n (y) n .
∂y ∂y
Alors, en notant y(x) = ψ ◦ φ−1 (x), on obtient la formule suivante pour le changement
de coordonnées :
Y j (y) = X φ · yj φ ◦ ψ −1 (y).

Exemple. Champ de gravité d’un objet de masse m : la variété est ici M = R3 \{0} et le
champ de vecteurs est

m 1 ∂ 2 ∂ 3 ∂

X(x) = 3 x +x +x
r ∂x1 ∂x2 ∂x3
p
où r = (x1 )2 + (x2 )2 + (x3 )2 .

Comment trouver d’autres champs de vecteurs ? Une première possibilité est de res-
treindre les champs de vecteurs sur Rn aux sous-variétés de Rn . Le résultat suivant indique
sous quelle condition une telle restriction est possible.

Lemme 3.2. Soient N une sous-variété de M et X un champ de vecteurs sur M tel que,

∀p ∈ N, X(p) ∈ Tp N.

Alors la restriction X|N de X à N est un champ de vecteurs sur N .

32 Chap. 3 – Champs de vecteurs

Preuve. La restriction X|N : p ∈ N 7→ X(p) ∈ Tp N est une application de N dans T N .

Il suffit donc de montrer qu’elle est différentiable.
Comme N est une sous-variété de M , on peut choisir des cartes (U, φ) sur M (dim M =
n) et (V, ψ) sur N (dim N = k) telles que V = U ∩ N , ψ(V ) = φ(U ) ∩ Rk et ψ = ϖ ◦ φ|V ,
où ϖ est la projection canonique de Rn sur Rk .
On obtient alors une carte (π −1 (U ), dφ) sur T M et une carte (π −1 (V ), dψ) sur T N ,
telles que π −1 (V ) ⊂ π −1 (U ) et dψ = dϖ ◦ dφ|π−1 (V ) . Notons dφ = (x1 , . . . , xn , α1 , . . . , αn ).
Un couple (p, Yp ) appartient à π −1 (V ) si et seulement si
dφ(p, Yp ) = (x1 , . . . , xk , 0, . . . , 0, α1 , . . . , αk , 0, . . . , 0),
et dans ce cas dψ(p, Yp ) = (x1 , . . . , xk , α1 , . . . , αk ).
Ainsi, si dφ(X) = (x, α1 (x), . . . , αn (x)) sont les composantes du champ de vecteurs
X, on a αk+1 (x) = · · · = αn (x) = 0 pour tout x ∈ Rk × {0}, puisque p = φ−1 (x) ∈ N et
X(p) ∈ Tp N . Ceci implique que X|N a pour coordonnées
dψ(X|N ) = (x1 , . . . , xk , α1 (x), . . . , αk (x))
et est différentiable puisque les αi le sont.

Exercice 3.1. On définit des champs de vecteurs sur R4 :

∂ 1 ∂ 4 ∂ 3 ∂
X = −x2 + x + x − x
∂x1 ∂x2 ∂x3 ∂x4
∂ ∂ ∂ ∂
Y = −x3 1 − x4 2 + x1 3 + x2 4
∂x ∂x ∂x ∂x
4 ∂ 3 ∂ 2 ∂ 1 ∂
Z = −x +x −x +x .
∂x1 ∂x2 ∂x3 ∂x4
Montrer que leurs restrictions à S 3 sont des champs de vecteurs sur S 3 dont les valeurs
en tout point forment une base orthonormée de l’espace tangent.

Remarque. Il est rare en général qu’il existe une famille de champs de vecteurs sur une
variété dont les valeurs en tout point forment une base de l’espace tangent. Une variété
possédant cette propriété est dite parallélisable. Les seules sphères parallélisables sont S 1 ,
S 3 et S 7 . Aucune sphère de dimension paire ne peut l’être puisque tout champ de vecteurs
s’y annule au moins une fois (“on ne peut pas peigner la sphère”). En revanche tout ouvert
de carte est parallélisable (avec les champs de coordonnées). Notons enfin qu’une variété
est parallélisable si et seulement si son fibré tangent admet une trivialisation.
On peut enfin se demander si il est possible de “transporter” un champ de vecteurs :
étant donnée une application différentiable F : M → N et un champ de vecteurs X sur
M , existe-t-il un champ de vecteurs sur N dont les valeurs soient les images de celles de
X par dF ?
En général la réponse est non (prendre par exemple M = R2 , N = R, F = la projection
sur la première coordonnée et X(x) = (x2 , 0)). Ce n’est possible que dans certains cas,
en particulier quand F est un difféomorphisme.
3.3 Équations différentielles 33

Définition. Soient F : M → N un difféomorphisme et X un champ de vecteurs sur M .

Alors le transport de X par F , noté F∗ X, est le champ de vecteurs sur N défini par

F∗ X = dF (X) ◦ F −1 ou F∗ X(q) = dFF −1 (q) X(F −1 (q)) .

Exemple. Les champs de coordonnées s’écrivent

∂ M ∂
i
= φ−1 ∗ i .
∂x ∂x

3.3 Équations différentielles

Nous allons maintenant étudier le mouvement sur une variété.

Définition. On appelle équation différentielle sur la variété M une équation de la forme

q̇ = X(q), q ∈ M,

où X est un champ de vecteurs sur M .

La donnée d’une équation différentielle est ainsi équivalente à celle du champ de

vecteurs X. Une solution de l’équation différentielle, c’est-à-dire une courbe intégrale du
champ X, est une courbe c(t) ∈ M , définie sur un intervalle J ⊂ R, telle que

dc

(t) = X c(t) , ∀t ∈ J.
dt
Comme pour les équations différentielles dans les espaces vectoriels normés (voir le
cours AO102), nous allons avoir des résultats d’existence et d’unicité pour les solutions
et une dépendance différentiable par rapport aux conditions initiales.

Théorème 3.3 (Existence et unicité des solutions). Soit q̇ = X(q), q ∈ M , une équation
différentielle sur M . Pour tout point p ∈ M , si η > 0 est suffisamment petit il existe
une unique courbe intégrale cp (t) de X, définie pour t ∈ ] − η, η[, satisfaisant la condition
initiale cp (0) = p.

Preuve. Le résultat est local : il suffit donc de le prouver dans une carte, c’est-à-dire
dans un ouvert de Rn . C’est alors un résultat classique (voir le cours AO102 [3]).

Dans la suite, nous parlerons de la courbe intégrale de X issue de p : il s’agit de la

solution cp (t) pour laquelle l’intervalle de définition ]a, b[ ⊃ ] − η, η[ est maximal.
Mentionnons ici une conséquence du théorème d’existence et d’unicité qui sera utile
pour les sous-variétés.
34 Chap. 3 – Champs de vecteurs

Corollaire 3.4. Soient q̇ = X(q) une équation différentielle sur M et N ⊂ M une

sous-variété de M tels que, pour tout q ∈ N , X(q) ∈ Tq N . Alors la courbe intégrale de
X issue d’un point p ∈ N est incluse dans N pour t suffisamment petit.
Preuve. D’après le lemme 3.2, la restriction X|N est un champ de vecteurs sur N .
L’équation différentielle q̇ = X|N (q), q ∈ N , a donc une solution cp (t) ⊂ N pour p ∈ N
et t suffisamment petit.
Or cette courbe est aussi solution de l’équation différentielle q̇ = X(q) sur M , pour la
condition initiale p. Le théorème d’unicité montre que c’est bien une courbe intégrale de
X.

Théorème 3.5 (Dépendance par rapport aux conditions initiales). La solution d’une
équation différentielle q̇ = X(q) sur M dépend de façon différentiable de la condition
initiale : pour tout p ∈ M , il existe un voisinage Up ⊂ M de p et un intervalle I = {|t| <
εp } tels que l’application

ϕ : I × Up → M
(t, q) 7→ ϕ(t, q) = cq (t)
est différentiable, ainsi que p 7→ εp .
Preuve. Comme pour le théorème 3.3, il suffit de se restreindre à une carte et d’utiliser
les résultats sur Rn .

3.4 Flots et groupes de difféomorphismes

Étant donné un point p ∈ M , on peut définir pour chaque t ∈ I une application

ϕt : Up → M
q 7→ ϕt (q) = ϕ(t, q).
Il résulte des deux précédents théorèmes que pour chaque t l’application ϕt est un
difféomorphisme sur Up , autrement dit un difféomorphisme local sur M . La famille de ces
difféomorphismes locaux ϕt est appelée le flot de X.
Remarque. Il est parfois pratique de se représenter X comme le champ de vitesse d’un
fluide se déplaçant sur M . Le flot transporte alors une particule de liquide de la position
q au temps 0 à la position ϕt (q) au temps t.
Par définition de ϕ, et donc de ϕt , le flot peut également être caractérisé comme la
famille de difféomorphismes locaux solution de
∂ϕt
= X ◦ ϕt , ϕ0 = id.
∂t
Cette caractérisation permet de montrer que les flots se transportent par conjugaison.
3.4 Flots et groupes de difféomorphismes 35

Lemme 3.6. Soient F : M → N un difféomorphisme et X un champ de vecteurs sur M ,

de flot ϕt . Le flot du transport F∗ X de X par F est F ◦ ϕt ◦ F −1 , c’est-à-dire le conjugué
du flot de X.

Preuve. Posons ϕ′t = F ◦ ϕt ◦ F −1 . Comme ϕ′0 = idN , il suffit de vérifier que

∂ϕ′t
= F∗ X ◦ ϕ′t .
∂t

Avec les notations du théorème 3.5, ϕ′t (q) = F ◦ ϕ(t, F −1 (q)), et la dérivée partielle par
rapport à t s’écrit

∂ϕ′t ∂ϕ
(q) = dF ◦ (t, F −1 (q))
∂t ∂t
= dF ◦ X ◦ ϕt ◦ F −1 (q))
= dF ◦ X ◦ F −1 ◦ F ◦ ϕt ◦ F −1 (q)) = F∗ X ◦ ϕ′t (q).

Le flot forme un groupe local à un paramètre de difféomorphismes de M . Un tel groupe

est défini par les propriétés suivantes :
(i) ϕt est un difféomorphisme local ∀t ∈ I ;
(ii) t 7→ ϕt est différentiable ;
(iii) ϕ0 = id ;
(iv) ϕt ◦ ϕs = ϕs ◦ ϕt = ϕt+s , ∀t, s, t + s ∈ I.
Pour le flot, les trois premières propriétés viennent des théorèmes 3.3 et 3.5. Pour
la quatrième, on remarque que c(t) = ϕt+s (q) et c′ (t) = ϕt ϕs (q) sont solutions de la

même équation différentielle avec la condition initiale c(0) = c′ (0) = ϕs (q), ce qui montre
l’égalité. Résumons ceci dans une proposition.

Proposition 3.7. Tout champ de vecteurs sur M engendre un groupe local à un paramètre
de difféomorphismes sur M .

Champs de vecteurs complets.

Définition. Un champ de vecteurs X sur M est complet si, pour tout p ∈ M , la courbe
intégrale de X issue de p est définie pour tout t ∈ R.

Remarque. Même dans R tous les champs ne sont pas complets. Par exemple X(x) = x,
x0
x ∈ R est complet mais X(x) = x2 ne l’est pas (la solution est de la forme x(t) = 1−tx0
).
36 Chap. 3 – Champs de vecteurs

Un champ de vecteurs complet engendre un flot constitué de difféomorphismes (glo-

baux) de M . Ce flot forme un groupe à un paramètre de difféomorphismes sur M (défini
comme un groupe local par les propriétés (i)–(iv) dans lesquelles on remplace local par
global et I par R).
Dans la suite de ce cours nous ferons souvent l’hypothèse de complétude pour simplifier
les résultats et les preuves. Il est généralement aisé de retrouver le résultat sans cette
hypothèse. Il faut remarquer toutefois qu’en mécanique la plupart des champs de vecteurs
hamiltoniens ne sont pas complets. On peut également noter qu’il existe une condition
suffisante simple : si M est une variété compacte, tout champ sur M est complet.
On se pose maintenant la question de la réciproque de le proposition 3.7 : étant donné
un groupe à un paramètre de difféomorphismes sur M , est-ce le flot d’un champ de
vecteurs ?
Lemme* 3.8. Tout groupe à un paramètre de difféomorphismes ϕt sur M est le flot d’un
champ de vecteurs, appelé générateur infinitésimal du groupe et obtenu comme
∂ϕt
X= .
∂t t=0

Preuve. Notons ϕt les difféomorphismes de ce groupe. En tout point p ∈ M , on définit

un vecteur tangent à M :
∂ϕt
X(p) = (p) ∈ Tp M.
∂t t=0
Comme la famille ϕt dépend de façon différentiable de t et de p (propriétés (i) et (ii) des
groupes à un paramètre), p 7→ X(p) définit bien un champ de vecteurs X sur M .
Montrons maintenant que le flot de X est ϕt . Comme ϕ0 = id (propriété (iii)), il suffit
de montrer que ∂ϕ∂t
t
= X ◦ ϕt . Or, d’après la propriété (iv),
∂ϕt+s ∂ϕt+s ∂

(p) = (p) = ϕs ϕt (p) .
∂t ∂s ∂s
En choisissant s = 0 dans cette relation, on obtient
∂ϕt ∂


(p) = ϕs ϕt (p) = X ϕt (p) ,
∂t ∂s s=0

ce qui prouve le lemme.

Le tableau ci-dessous résume les relations entre les champs de vecteurs et les groupe
de difféomorphismes.

Champ de vecteurs =⇒ Groupe local à 1 paramètre

̸
=⇒ Groupe global à 1 paramètre
Champ de vecteurs complet ⇐⇒ Groupe global à 1 paramètre
3.4 Flots et groupes de difféomorphismes 37

Exemples.
— Dans M = R3 : soient a ∈ R3 et le groupe de translations à 1 paramètre ϕt (x) =
x + ta.
Le générateur infinitésimal de ce groupe est X(x) = a, ou plus exactement

∂ ∂ ∂
X(x) = a1 1
+ a2 2 + a3 3 .
∂x ∂x ∂x
Inversement, considérons l’équation différentielle dans R3 : ẋ = X(x) = a. La
courbe intégrale issue d’un point x0 est ϕt (x0 ) = x0 + ta.
— Dans M = GL2 (R), on considère les applications
 
1 t
ϕt (A) = A, t ∈ R.
0 1

L’ensemble de ces applications forme un groupe à un paramètre de difféomorphismes

sur GL2 (R) (le montrer). Le générateur infinitésimal de ce groupe est
 
d 0 1
X(A) = ϕt (A) t=0 = A.
dt 0 0

Inversement, l’équation différentielle sur GL2 (R)

 
0 1
Ȧ = A
0 0

est linéaire et a donc pour solution issue de A0 :

h 0 1  i 
1 t

A(t) = exp t A0 = A0 .
0 0 0 1

— Plus généralement, si M = Rn (ou un ouvert de Rn ) et si le champ de vecteurs X

est linéaire, c’est-à-dire X(x) = Ax, alors X est complet et son flot est ϕt (x) =
exp(tA)x (voir le cours AO102 [3]).
38 Chap. 3 – Champs de vecteurs
Chapitre 4

Familles de champs de vecteurs

Nous avons vu que, pour un système dynamique, c’est-à-dire un système dont la loi
d’évolution ou dynamique est donnée par une équation différentielle, le futur ϕ(t, p), t > 0,
est complètement déterminé par l’état présent p = ϕ(0, p).
Nous allons maintenant étudier des systèmes dont la dynamique est certes contrainte,
mais possède un certain nombre de degrés de libertés. Comme nous l’avons fait jusqu’à
maintenant, nous considérons que l’espace des configurations est une variété M et nous
supposons qu’à partir de chaque configuration, le système peut suivre un certain nombre
de directions de l’espace tangent. Ces directions peuvent être vues comme les valeurs d’une
famille de champs de vecteurs (en supposant bien entendu qu’elles soient suffisamment
régulières) ou, en introduisant une paramétrisation de cette famille, comme un système
commandé.

4.1 Systèmes commandés

Soit M une variété différentiable de dimension n.

Définition. Un système commandé est une famille d’équations différentielles

q̇ = Xu (q), q ∈ M, u ∈ U, (4.1)

où
— l’ensemble de commande U est un sous-ensemble (quelconque) de Rm ; m est donc
la taille de la commande u, ou encore le nombre de commandes scalaires ui ;
— les Xu forment une famille de champs de vecteurs sur M paramétrés par u ; au-
trement dit, on a une application X : M × U → T M telle que, pour u ∈ U fixé,
Xu (q) = X(q, u) est un champ de vecteurs sur M .

On peut également représenter un système commandé sans faire apparaı̂tre explicite-

ment les commandes, en insistant au contraire sur la famille de champs de vecteurs. Pour
40 Chap. 4 – Familles de champs de vecteurs

cela, on définit la famille F de champs de vecteurs associée au système commandé

F = {Xu , u ∈ U }.

Le système commandé s’écrit alors q̇ ∈ F. Il est important de noter que cette notation
ici est strictement équivalente à l’équation (4.1). Elle sera cependant plus pratique pour
énoncer un certain nombre de résultats et elle peut de plus prendre un sens plus général
(voir § 4.3).
Exemple (La voiture (modèle simplifié)). Nous prenons comme exemple type tout au long
de ce cours un modèle simplifié de la voiture. Dans ce modèle, la voiture est constituée d’un
seul essieu se déplaçant dans le plan. Une configuration est donnée par les coordonnées
planaires (x1 , x2 ) du milieu de l’essieu et par l’angle θ que fait la voiture avec l’axe 0x1 .
L’espace des configurations est donc la variété tridimensionnelle M = R2 × S 1 .
Les deux commandes sont la vitesse tangentielle u1 et la vitesse angulaire u2 . Nous
supposerons pour l’instant que l’ensemble de commande est R2 tout entier (nous verrons
au chapitre 5 d’autres possibilités pour U ).
Le système commandé prend alors la forme

q̇ = Xu (q) = u1 X1 + u2 X2 , q = (x, y, θ) ∈ M, u = (u1 , u2 ) ∈ R2 ,

où X1 et X2 , qui correspondent respectivement aux déplacements rectilignes et aux rota-

tions sur place, sont donnés en coordonnées par
∂ ∂ ∂
X1 (q) = cos θ + sin θ , X2 (q) = .
∂x1 ∂x2 ∂θ
La famille de champs de vecteurs associée à ce système commandé est l’ensemble des com-
binaisons linéaires à coefficients constants de X1 et X2 , c’est-à-dire F = Vect{X1 , X2 }.

On appellera loi de commande une fonction u(t) : J ⊂ R → U . Dans ce cours nous

nous limiterons aux lois de commande constantes par morceaux : nous obtiendrons ainsi
une bonne formulation géométrique tout en évitant des considérations techniques délicates
sur les mesures.
Une trajectoire du système commandé est une solution de la famille d’équations
différentielles (4.1), c’est-à-dire une courbe c(t) dans M , définie sur un intervalle J ⊂ R,
pour laquelle il existe une loi de commande u(t), t ∈ J, constante par morceaux et telle
que
dc

(t) = Xu(t) c(t) ,
dt
pour tout t ∈ J excepté aux points de discontinuité de u(t).
On peut écrire une trajectoire du système en fonction des flots des champs de vecteurs
Xu . En effet, soient p un point de M et u(t), t ∈ [0, T ], une loi de commande constante
par morceaux, caractérisée par ses points de discontinuité 0 < t1 < · · · < ts < T et sa
4.1 Systèmes commandés 41

valeur ui ∈ U sur chacun des intervalles ]ti , ti+1 [, i = 0, . . . , s (on pose t0 = 0 et ts+1 = T ).
Quand elle existe (problème de complétude), la trajectoire issue de p et correspondant à
la loi de commande u(t) est donnée par
u
c(t) = ϕut−t
i
i
◦ ϕtii−1 u0
−ti−1 ◦ · · · ◦ ϕt1 −t0 (p), pour t ∈ [ti , ti+1 ],

où ϕut i désigne le flot du champ de vecteurs Xui .

Dans ce cours nous nous intéresserons essentiellement au problème de l’atteignabilité :
quelle est la structure de l’ensemble des configurations atteignables à partir d’une confi-
guration donnée ? L’écriture des trajectoire en fonction des flots permet de bien définir
cet ensemble.

Définition. Pour un système commandé q̇ = Xu (q), u ∈ U , sur M , l’ensemble atteignable

à partir d’un point p ∈ M , noté Ap (ou Ap (F)), est l’ensemble des points pouvant être
atteints par une trajectoire issue de p, c’est-à-dire

Ap = ϕutkk ◦ · · · ◦ ϕut11 (p) ti ≥ 0, ui ∈ U, k ∈ N .



Essayons de comprendre quelle est la nature de cet ensemble.

Considérons d’abord le cas où l’ensemble U ne contient qu’une commande : le système
se réduit alors à une équation différentielle et l’ensemble atteignable à partir de p est la
demi-trajectoire {ϕt (p), t ≥ 0}. Ainsi Ap est contenu dans une variété de dimension 1
(voir § 4.3).
Prenons ensuite le cas où la commande ne peut prendre que deux valeurs, U = {1, 2}.
On a alors n o
Ap = ϕ2t2k ◦ ϕ1t2k−1 ◦ · · · ◦ ϕ2t2 ◦ ϕ1t1 (p) ti ≥ 0, k ∈ N .

La commutativité des flots ϕ1t et ϕ2s va alors jouer un rôle essentiel.

Supposons en effet que ces flots commutent, ce qui s’écrit ϕ1t ◦ ϕ2s = ϕ2s ◦ ϕ1t pour tous
t, s. Dans l’expression de Ap , on peut alors regrouper tous les flots ϕ1 d’un coté et tous
les flots ϕ2 de l’autre, ce qui donne

Ap = ϕ2s ◦ ϕ1t (p) t, s ≥ 0 .



Dans ce cas l’ensemble atteignable est donc un morceau de variété de dimension 2.

Que se passe-t-il quand les flots ne commutent pas (ce qui est généralement le cas) ?
A-t-on encore un résultat du type “k valeurs pour la commande ⇒ Ap inclus dans une
variété de dimension k ” ? Évidemment non, comme le montre l’exemple suivant.
Exemple. Reprenons l’exemple de la voiture. Le système commandé est donné ici par
deux champs de vecteurs sur la variété tridimensionnelle M = R2 × S 1 . Il est cependant
clair que l’on peut joindre toute paire de points de M en alternant les flots de X1 et de
X2 (on peut s’en convaincre sur un dessin). En particulier Ap = M pour tout point p.
42 Chap. 4 – Familles de champs de vecteurs

Ceci est possible parce que les flots de X1 et de X2 ne commutent pas. On peut en
effet montrer que ϕ1t ◦ ϕ2t ̸= ϕ2t ◦ ϕ1t en calculant directement la composée ϕ2−t ◦ ϕ1−t ◦ ϕ2t ◦ ϕ1t ,
par exemple à partir de 0 :
ϕ1
t tϕ2 ϕ1−t
0 −→ (t, 0, 0) −→ (t, 0, t) −→ (t(1 − cos t), −t sin t, t)
ϕ2−t
−→ (t(1 − cos t), −t sin t, 0) = γ(t).

Cette courbe γ(t) est bien différente de 0, ce qui montre la non commutativité des flots.
On constate de plus que, pour t petit,

γ(t) = (0, −t2 , 0) + o(t2 ).

√
Ainsi le vecteur tangent ċ(0) à la courbe c(t) = γ( t), t ≥ 0, est la direction transverse
au plan engendré par X1 (0), X2 (0).

Pour pouvoir continuer l’étude des ensembles atteignables, il est nécessaire de ca-
ractériser la commutation des flots. C’est ce que nous allons faire dans la section suivante.
L’objectif est de faire cette caractérisation à partir des champs de vecteurs, sans avoir à
calculer les flots. C’est un des principes essentiels de l’étude des équations différentielles et
des systèmes commandés (voir le cours AO102) : il faut caractériser le comportement des
solutions à partir des équations elles-mêmes, c’est-à-dire à partir des champs de vecteurs
qui les constituent, sans calculer les solutions (ce que l’on ne sait généralement pas faire).

4.2 Crochets et algèbre de Lie

Rappelons (§ 3.2) qu’un champ de vecteurs X ∈ X (M ) peut être considéré comme
une dérivation sur l’ensemble C ∞ (M ) des fonctions à valeurs réelles sur M :

X : C ∞ (M ) → C ∞ (M )
g 7→ X · g.
On définit alors le produit de deux champs de vecteurs X et Y comme l’opérateur

XY : C ∞ (M ) → C ∞ (M )
g 7→ (XY ) · g = X · (Y · g).
Cet opérateur n’est pas une dérivation, XY n’est donc pas un champ de vecteurs. En
effet, il ne vérifie pas la règle de Leibniz :

XY · (f g) = X · (f Y · g + gY · f )
= f XY · g + gXY · f + (X · f )(Y · g) + (X · g)(Y · f ) . (4.2)
| {z }
termes en trop

Il s’agit en fait d’un opérateur du second ordre.

4.2 Crochets et algèbre de Lie 43

∂ ∂ ∂2g
Exemple. Dans Rn , avec X = ∂x1
et Y = ∂x2
, on a XY · g = ∂x1 ∂x2
.

Il est cependant possible de définir un champ de vecteurs à partir du produit.

Définition. Le crochet de Lie de deux champs de vecteurs X, Y ∈ X (M ) est le champ
de vecteurs
[X, Y ] = XY − Y X.

Pour que cette définition soit correcte il faut montrer que XY − Y X est bien un
champ de vecteurs, c’est-à-dire une dérivation sur C ∞ (M ). Or la linéarité est évidente et
la règle de Leibniz se déduit du calcul (4.2) :

[X, Y ] · (f g) = XY · (f g) − Y X · (f g)
= f XY · g + gXY · f − f Y X · g − gY X · f
= f [X, Y ] · g + g[X, Y ] · f.

PnEn coordonnées
i ∂
locales, un champ de vecteurs étant donné sous la forme X(x) =
i=1 X (x) ∂xi , les composantes du crochet de Lie [X, Y ](x) sont
n 
i
X ∂Y i ∂X i 
[X, Y ] (x) = X (x) j (x) − Y (x) j (x) = X · Y i (x) − Y · X i (x).
j j

j=1
∂x ∂x
i

Matriciellement, en notant JX(x) = ∂X ∂xj
(x) la matrice jacobienne du vecteur des com-
posantes de X, on obtient
     
[X, Y ]1 (x) X 1 (x) Y 1 (x)
..  .   . 
 = JY (x)  ..  − JX(x)  ..  .
 
 .
[X, Y ]n (x) X n (x) Y n (x)

Exemples.
— Si X = ∂x∂ i et Y = ∂x∂ j , alors [X, Y ] = 0. Plus généralement si X et Y sont
constants dans un système de coordonnées, c’est-à-dire X 1 (x), . . ., X n (x) et Y 1 (x),
. . ., Y n (x) tous constants, alors [X, Y ] = 0 sur le domaine des coordonnées.
— Dans Rn , considérons des champs de vecteurs linéaires, c’est-à-dire que, dans la
base des ∂x∂ i , ils s’écrivent matriciellement comme

X(x) = Ax et Y (x) = Bx, x ∈ Rn

A et B étant des matrices (n × n) à coefficients constants. Le crochet de Lie de

ces champs s’écrit alors

[X, Y ](x) = JY (x)X(x) − JX(x)Y (x) = (BA − AB)x.

Autrement dit le crochet de Lie s’écrit en fonction du commutateur BA − AB des

matrices A et B.
44 Chap. 4 – Familles de champs de vecteurs

Exercice 4.1. Calculer [X1 , X2 ] pour les champs de vecteurs X1 et X2 du système com-
mandé correspondant√à la voiture. Comparer la valeur [X1 , X2 ](0) avec le vecteur tangent
à la courbe c(t) = γ( t) définie page 42.

Comme nous venons de le voir pour la voiture, le commutateur des flots peut être
relié au crochet de Lie des champs de vecteurs correspondants.

Proposition 4.1. Soient X1 , X2 des champs de vecteurs sur M de flots respectifs ϕ1t , ϕ2s
et p un point de M . On considère la courbe

γ(t) = ϕ2−t ◦ ϕ1−t ◦ ϕ2t ◦ ϕ1t (p)

pour t suffisamment petit.

Alors, dans des coordonnées locales φ = (x1 , . . . , xn ) centrées en p,

γ φ (t) = t2 [X1 , X2 ]φ (p) + o(t2 ),

où l’on a identifié le champ [X1 , X2 ]φ sur Rn avec le vecteur de ses coordonnées.

Preuve. Pour t suffisamment petit, on peut se placer dans un système local de coor-
données centré en p (c’est-à-dire que p a pour coordonnées 0), et on identifie vecteurs
tangents et champs de vecteurs aux vecteurs colonnes de leurs composantes. On va faire
un développement limité de γ(t) et montrer que

γ(t) = t2 [X1 , X2 ](p) + o(t2 ),

ce qui prouvera la proposition.

Nous commençons par calculer

ϕ1t ◦ ϕ2t γ(t) = ϕ2t ◦ ϕ1t (p) − ϕ1t ◦ ϕ2t (p).

Remarquons tout d’abord que, si ϕt désigne le flot d’un champ de vecteurs X,

dϕt t2 d2 ϕt
ϕt (x) = x + t (x) + (x) + o(t2 )
dt t=0 2 dt2 t=0
2
t
= x + tX(x) + JX(x)X(x) + o(t2 ).
2
En choisissant ϕ = ϕ1 et x = 0 puis ϕ = ϕ2 et x = ϕ1t (p), on obtient

t2
ϕ2t ϕ1t (p)


◦ = tX1 (p) + JX1 (p)X1 (p) + t X2 (p) + tJX2 (p)X1 (p) +
2
t2
+ JX2 (p)X2 (p) + o(t2 ).
2
4.2 Crochets et algèbre de Lie 45

On peut faire le même calcul en échangeant les rôles de X1 et de X2 , ce qui nous conduit
à

ϕ2t ◦ ϕ1t (p) − ϕ1t ◦ ϕ2t (p) = t2 JX2 (p)X1 (p) − JX1 (p)X2 (p) + o(t2 )

= t2 [X1 , X2 ](p) + o(t2 ),

et finalement

γ(t) = ϕ2−t ◦ ϕ1−t ϕ2t ◦ ϕ1t (p) − ϕ1t ◦ ϕ2t (p)

= t2 [X1 , X2 ](p) + o(t2 ).

Cette proposition met en évidence le lien entre le crochet et la non-commutation des

flots (ce qui est le but de cette section !). Nous avons en fait la caractérisation suivante,
dont la preuve est donnée un peu plus loin page 47.

Lemme 4.2. Deux champs de vecteurs X1 , X2 sur M satisfont [X1 , X2 ] ≡ 0 si et seule-

ment si leurs flots ϕ1t , ϕ2s commutent pour tous t, s suffisamment petits.

Une autre propriété importante du crochet de Lie est qu’il est préservé par transport
par un difféomorphisme.

Lemme 4.3. Soient F : M → N un difféomorphismes entre variétés et X, Y des champs

de vecteurs sur M . Alors
F∗ [X, Y ] = [F∗ X, F∗ Y ].

Preuve. Rappelons que, en tant que dérivation, le champ F∗ X s’écrit

F∗ X : C ∞ (N ) → C ∞ (N )
g 7→ X · (g ◦ F ) ◦ F −1 .

Le transport par F du crochet [X, Y ] vérifie donc

F∗ [X, Y ] · g = [X, Y ] · (g ◦ F ) ◦ F −1 , ∀g ∈ C ∞ (N )


= X · Y · (g ◦ F ) − Y · X · (g ◦ F ) ◦ F −1 , ∀g ∈ C ∞ (N ).

Or Y · (g ◦ F ) est égal à F∗ Y ◦ F (de même avec X). L’égalité ci-dessus devient donc,
pour toute fonction g ∈ C ∞ (N ),

F∗ [X, Y ] · g = F∗ X · (F∗ Y · g) − F∗ Y · (F∗ X · g) = [F∗ X, F∗ Y ] · g,

ce qui montre que les champs de vecteurs F∗ [X, Y ] et [F∗ X, F∗ Y ] sont égaux.
46 Chap. 4 – Familles de champs de vecteurs

Dérivées de Lie.* Le crochet de Lie est un cas particulier d’une notion plus générale,
la dérivée de Lie. Cette notion prolonge celle de dérivée directionnelle.

Définition. Soient g ∈ C ∞ (M ) une fonction et X ∈ X (M ) un champ de vecteurs sur

M , de flot ϕt . La dérivée de Lie de g par rapport à X, notée LX g, est la fonction sur M
dont la valeur en un point p est

1


LX g(p) = lim g ϕt (p) − g(p) .
t→0 t

Cette dérivée de Lie LX g n’est en fait rien d’autre que l’action de l’opérateur de
dérivation X sur g,
LX g = X · g.
En effet, par définition,
d

LX g(p) = g ◦ ϕt (p) t=0 .
dt
Comme le vecteur tangent en 0 de la courbe c(t) = ϕt (p) est X(p) = Xp , il résulte des
propositions 2.3 et 2.4 que LX g(p) = Xp · g.
On peut de la même façon définir la dérivée de Lie d’un champ de vecteurs Y . L’idée
ici est de comparer les valeurs de Y le long du flot ϕt de X avec la valeur en p. Cette
comparaison ne pouvant se faire qu’entre des objets appartenant au même espace vecto-
riel, on utilise la différentielle du flot pour ramener les valeurs Y (ϕt (p)) dans Tp M . D’où
la définition suivante.

Définition. Soient X, Y ∈ X (M ) des champs de vecteurs sur M et ϕt le flot de X. La

dérivée de Lie de Y par rapport à X, notée LX Y , est le champ de vecteurs sur M dont
la valeur en un point p est

1
 d

LX Y (p) = lim dϕ−t Y (ϕt (p)) − Y (p) = ϕ−t ∗ Y (p) t=0
.
t→0 t dt

Remarque. On peut aussi écrire LX Y (p) comme :

1
 d

LX Y (p) = lim Y (p) − dϕt Y (ϕ−t (p)) = − ϕt ∗ Y (p) t=0
.
t→0 t dt

La dérivée de Lie coı̈ncide avec le crochet de Lie de X et Y .

Lemme 4.4. Pour tous X, Y ∈ X (M ), LX Y = [X, Y ].

Preuve. Faisons au préalable la remarque suivante, qui généralise le lemme d’Hadamard

(lemme 2.1) :
4.2 Crochets et algèbre de Lie 47

si g : ] − ε, ε[ ×M → R est une fonction différentiable satisfaisant g(0, p) = 0 pour tout

p ∈ M , il existe une fonction différentiable h : ] − ε, ε[ ×M → R telle que
∂g
g(t, p) = th(t, p) et (0, p) = h(0, p).
∂t
Il suffit en effet de poser Z 1
∂g
h(t, p) = (st, p)ds.
0 ∂s
Venons-en à la preuve du lemme. Considérons g ∈ C ∞ (M ). D’après la remarque
précédente, il existe une fonction h(t, x) telle que
∂(g ◦ ϕ−t )
g ◦ ϕ−t (q) = g(q) + th(t, q) et h(0, q) = (q) = −X · g(q).
∂t t=0

Ainsi,




dϕ−t Y (ϕt (p)) · g = Y · (g ◦ ϕ−t ) ϕt (p) = Y · g ϕt (p) + tY · h(t, ϕt (p)),
et
1

LX Y (p) · g = lim dϕ−t Y (ϕt (p)) − Y (p) · g
t→0 t
1

= lim Y · g ϕt (p) − Y · g(p) + lim Y · h(t, ϕt (p))
t→0 t t→0
= LX (Y · g)(p) − Y · (X · g)(p) = [X, Y ](p) · g.

Remarque. Comme nous l’avons vu dans la preuve précédente, la valeur en un point de

la dérivée de Lie d’une fonction peut être représentée, selon le contexte, par l’une des
notations suivantes : LX g(p), X · g(p), X(p) · g ou Xp · g. Ces notations sont bien entendu
toutes équivalentes mais mettent l’accent sur des propriétés différentes de l’expression.
L’interprétation du crochet de Lie comme une dérivée de Lie permet de prouver le
lemme de la page 45 sur la commutation des flots.
Preuve (du lemme 4.2). Soient X1 , X2 ∈ X (M ) et ϕ1t , ϕ2s leurs flots respectifs. Ces deux
flots commutent si et seulement si, pour tous t, s suffisamment petits,
ϕ1−t ◦ ϕ2s ◦ ϕ1t = ϕ2s .
Posons θs = ϕ1−t ◦ ϕ2s ◦ ϕ1t . C’est un groupe local à un paramètre de difféomorphismes. Il
existe donc localement un champ Zt générateur infinitésimal de ce groupe :
∂θs ∂ 1
ϕ−t ◦ ϕ2s ◦ ϕ1t


Zt = =
∂s s=0 ∂s s=0
= dϕ1−t ◦ X2 ◦ ϕ1t = ϕ1−t ∗ X2 .
48 Chap. 4 – Familles de champs de vecteurs

Ainsi les flots commutent si et seulement si ce champ Zt est indépendant de t, c’est-à-dire

Zt = Z0 = X2 . Or
d d 1 d
Zt = (ϕ−t ∗ X2 ) = (ϕ1−t−s ∗ X2 )
dt dt ds s=0
d
= dϕ1−t ∗ (ϕ1−s ∗ X2 ) = dϕ1−t ∗ LX1 X2 ,
ds s=0

ce qui implique que Zt est indépendant de t si et seulement si LX1 X2 = [X1 , X2 ] = 0.

Algèbres de Lie. On a déjà vu que l’ensemble X (M ) des champs de vecteurs sur M

est un espace vectoriel pour l’addition point par point :

(αX + βY )(q) = αX(q) + βY (q), ∀q ∈ M, ∀α, β ∈ R.

En fait X (M ), muni de cette addition et de la loi interne (ou produit) défini par le crochet
de Lie est une algèbre 1 sur R.
De plus le crochet de Lie satisfait les deux propriétés suivantes :
(a) l’anti-symétrie : [X, Y ] = −[Y, X],
(b) l’identité de Jacobi : [X, [Y, Z]] + [Z, [X, Y ]] + [Y, [Z, X]] = 0.
Toute algèbre sur R dont la loi interne satisfait ces deux propriétés est appelée une algèbre
de Lie. Ainsi X (M ) est une algèbre de Lie.
Exercice 4.2. Vérifier que Mn (R) muni de la loi interne [A, B] = BA−AB est une algèbre
de Lie.

Une sous-algèbre de Lie de X (M ) est un sous-espace vectoriel de X (M ) qui est stable

par le crochet de Lie. C’est donc également une algèbre de Lie. En particulier, si F une
famille de champs de vecteurs sur M , on appelle (sous-)algèbre de Lie engendrée par F,
notée Lie(F), le plus petit sous-espace vectoriel S de X (M ) qui contient F et qui est
stable par le crochet de Lie, c’est-à-dire [X, Y ] ∈ S pour tout X, Y ∈ S.
Prenons par exemple une famille à deux éléments, F = {X1 , X2 }. Les champs X1 et
X2 appartiennent à Lie(X1 , X2 ), donc leur crochet [X1 , X2 ] aussi, ainsi que leurs crochets
avec le crochet [X1 , [X1 , X2 ]], [X2 , [X1 , X2 ]], etc. . . et finalement

Vect [Xi1 , [Xi2 , . . . , [Xik−1 , Xik ] . . . ]], k ≥ 1, ij = 1 ou 2 ⊂ Lie(X1 , X2 )

(la relation de Jacobi permet de montrer que l’on a en fait égalité entre ces deux espaces).
Il est clair qu’en général Lie(F) est un sous-espace vectoriel de X (M ) de dimension
infinie, même si F ne contient que deux champs. En revanche si F ne contient qu’un
élément, F = {X}, alors Lie(X) = Vect{X} car [X, X] = 0.
1. Rappelons qu’une algèbre sur R est un R-espace vectoriel E muni d’une application bilinéaire de
E × E → E.
4.3 Orbite d’une famille de champs de vecteurs 49

En un point p ∈ M , on note

Liep (F) = {X(p), X ∈ Lie(F)}.

C’est un sous-espace vectoriel de l’espace tangent Tp M , c’est donc un espace vectoriel de

dimension finie.
Exemples.
— Soit F = {X1 , X2 } sur M = R2 , avec X1 = ∂x∂ 1 , X2 = ∂x∂ 1 + x1 ∂x∂ 2 . Leur crochet
de Lie est [X1 , X2 ] = ∂x∂ 2 et tous les autres crochets sont nuls. Ainsi
 ∂ ∂ 1 ∂
Lie(X1 , X2 ) = Vect , , x ,
∂x1 ∂x2 ∂x2
est de dimension 3 et, pour tout x ∈ R2 , Liex (X1 , X2 ) = Tx R2 est de dimension 2.
— Prenons la même famille que dans l’exemple précédent avec X1 = ∂x∂ 1 , X2 =
2
∂
∂x1
+ f (x1 ) ∂x∂ 2 , où f (t) = e−1/t . Cette fonction f est de classe C ∞ et toutes ses
dérivées sont nulles en 0. En calculant les crochets successifs, on obtient
 ∂ ∂
Lie(X1 , X2 ) = Vect 1
, f (k) (x1 ) 2 , k ≥ 0 ,
∂x ∂x
qui est de dimension infinie. En revanche, si x1 ̸= 0, Liex (X1 , X2 ) = Tx R2 est de
dimension 2 alors que, si x1 = 0, Liex (F) = R ∂x∂ 1 est de dimension 1.

Exercice 4.3. Calculer Lie(X1 , X2 ) dans les cas suivants :

— X1 = ∂x∂ 1 et X2 = ∂x∂ 2 + x1 ∂x∂ 3 dans R3 ;
— X1 = ∂x∂ 1 et X2 = ∂x∂ 2 + (x1 )2 ∂x∂ 3 dans R3 ;
— X1 et X2 champs de vecteurs du système commandé de la voiture.

4.3 Orbite d’une famille de champs de vecteurs

Nous allons maintenant généraliser les notions de flot et de courbe intégrale d’un
champ de vecteurs et les définir pour des familles de champs de vecteurs. Ces notions
sont primordiales pour l’étude de l’ensemble atteignable d’un système commandé, comme
nous le verrons au chapitre 5.
Considérons donc une famille F de champs de vecteurs sur M . C’est un sous-ensemble,
a priori quelconque, de X (M ).
Remarque. Les champs de vecteurs de F, ainsi que tous les champs que nous rencon-
trerons dans cette section, seront supposés complets. Cette hypothèse ne modifie pas
fondamentalement les résultats (ni les preuves) mais permet en revanche de simplifier
leur énoncé. Pour retrouver les résultats généraux, il faudra systématiquement remplacer
les définitions et résultats globaux utilisant les flots par les équivalents locaux et vérifier
que les preuves restent valables (ce sera toujours le cas).
50 Chap. 4 – Familles de champs de vecteurs

Définition. Le groupe de difféomorphismes engendré par F, noté G(F), est le groupe

engendré par les flots des éléments de F :
n o
G(F) = ϕX tk
k
◦ · · · ◦ ϕX1
t1 ti ∈ R, X i ∈ F, k ∈ N ,

où ϕX
t désigne le flot du champ de vecteurs X.

Définition. Soit p ∈ M . L’orbite de F issue de p est l’ensemble

n o
Orbp (F) = ϕX tk
k
◦ · · · ◦ ϕX1
t1 (p) ti ∈ R, X i ∈ F, k ∈ N ,

ou encore Orbp (F) = {φ(p), φ ∈ G(F)}.

Quand F = {Xu , u ∈ U } est la famille de champs associée à un système com-

mandé (4.1), l’orbite Orbp (F) ressemble beaucoup à l’ensemble atteignable Ap (F), avec
une différence importante : dans l’orbite, il est permis de se déplacer en avant et en
arrière dans chaque direction Xu ∈ F alors que dans Ap (F), seul le mouvement en avant
(c’est-à-dire pour un temps positif) est autorisé. On a donc Ap (F) ⊂ Orbp (F).
Bien entendu les deux notions coı̈ncident si les champs −Xu appartiennent aussi à F.
Proposition 4.5. Soit F la famille associée à un système commandé. Si F est une
famille symétrique, c’est-à-dire si

X ∈ F ⇒ −X ∈ F,

l’orbite de F et l’ensemble atteignable par le système commandé coı̈ncident pour tout

point p ∈ M : Orbp (F) = Ap (F).

Nous nous intéressons maintenant à la nature des orbites. Est-ce que ce sont des
sous-variétés de M ? ou au moins des variétés ? Commençons par le cas le plus simple,
celui d’une famille à un élément X.
Dans ce cas, le groupe engendré par F est simplement le flot de X et l’orbite issue
d’un point p est la courbe intégrale de X issue de ce point p : Orbp (F) = {ϕt (p), t ∈ R}.
Cette orbite peut être :
— soit réduite à un point : si X(p) = 0, Orbp (F) = {p} ;
— soit une sous-variété de M difféomorphe à S 1 (si X(p) ̸= 0 et t 7→ ϕt (p) n’est pas
injective) ;
— soit une sous-variété immergée de dimension 1 (si X(p) ̸= 0 et t 7→ ϕt (p) est
injective) : en effet, dans ce cas Orbp (F) est l’image de l’application t 7→ ϕt (p),
qui est une immersion injective (car ∂ϕ ∂t
t
(p) = X(ϕt (p)) ̸= 0)
(la preuve de ces résultats est laissée en exercice). Dans les trois cas, Orbp (F) est une sous-
variété immergée. Est-ce une sous-variété ? En général non, comme le montre l’exemple
ci-dessous (rappelons qu’une sous-variété est l’image d’un plongement alors qu’une sous-
variété immergée est l’image d’une immersion injective, voir § 1.3).
4.3 Orbite d’une famille de champs de vecteurs 51

Exemple (Enroulements denses sur le tore). Sur le tore T2 , représenté ici comme R2 /Z2 ,
on considère le champ de vecteurs (donné dans les coordonnées locales induites par R2 )
∂ ∂
X(x) = α + β
∂x1 ∂x2
où α et β sont des réels dont le rapport β/α est irrationnel. L’orbite de X issue d’un
point quelconque x0 s’écrit en coordonnées locales x1 (t) = x10 + αt, x2 (t) = x20 + βt ; cette
orbite est dense dans le tore. Il n’est donc pas possible que ce soit une sous-variété de T2 .
Ainsi, pour une famille F quelconque, le meilleur résultat qu’on puisse espérer pour
les orbites est qu’elles soient des sous-variétés immergées (c’est déjà beaucoup, puisque
ça implique une structure de variété). C’est ce que nous allons montrer dans le reste de
cette section.
Théorème 4.6 (Théorème de l’orbite). Les orbites d’une famille F de champs de vecteurs
sur M sont des sous-variétés immergées connexes de M .
L’espace tangent à une orbite est

Tq Orbp (F) = Vect φ∗ X(q) X ∈ F, φ ∈ G(F) , pour q ∈ Orbp (F).
Preuve*. La preuve 2 de ce théorème est longue et, parfois, difficile... On y développe
cependant des outils qui seront utiles par ailleurs, tout particulièrement la construction de
coordonnées locales. Pour plus de clarté, nous allons présenter cette preuve en plusieurs
étapes.
Fixons un point p ∈ M et son orbite Orbp (F). Pour q ∈ M , on désigne par Πq le
sous-espace vectoriel de Tq M :

Πq = Vect φ∗ X(q) X ∈ F, φ ∈ G(F) .
Voici le plan de la démonstration.
(i) Immersion locale Pour chaque point q ∈ Orbp (F), on construit une application
Gq : Rm → M , où m = dim Πq , à valeurs dans Orbp (F), qui est une immersion
sur un voisinage Uq ⊂ Rm de 0.
(ii) Topologie On munit Orbp (F) d’une topologie T dont la base est formée des
Gq (U ), U ⊂ Uq ouvert de Rm .
(iii) Structure différentiable L’ensemble des cartes (Gq (U ), G−1 q ) forme un atlas
et définit donc une structure différentiable sur Orbp (F).
(iv) Variété connexe Orbp (F) est un espace topologique connexe, séparé et à base
dénombrable. Muni de la structure différentiable du (iii), il forme donc une variété
connexe, dont l’espace tangent en un point q est Πq .
(v) Sous-variété immergée L’inclusion i : Orbp (F) → M est une immersion, ce
qui fait de la variété Orbp (F) une sous-variété immergée de M .
Nous allons maintenant détailler ces différents points.
2. Il y a plusieurs façons différentes de prouver le théorème de l’orbite. La preuve que nous donnons
ici est, à peu de choses près, celle établie par A. Agrachev et Y. Sachkov [1].
52 Chap. 4 – Familles de champs de vecteurs

(i) Immersion locale Commençons par montrer que Πq , le candidat à l’espace tangent,
est de dimension constante sur une orbite.
Lemme 4.7. Pour tout q ∈ Orbp (F), dim Πq = dim Πp .
Preuve. Soit q un point de Orbp (F). Il existe donc un difféomorphisme ψ ∈ G(F) tel
que q = ψ(p). Considérons d’autre part un élément φ∗ X(p) de l’espace vectoriel Πp .
Comme ψ est un difféomorphisme, dψp : Tp M → Tq M est un isomorphisme et

dψp (φ∗ X(p)) = dψ ◦ dφ(X) ◦ φ−1 ◦ ψ −1 (q) = (ψ ◦ φ)∗ X(q) ∈ Πq .

Ainsi dψp (Πp ) ⊂ Πq , ce qui implique que dim Πp ≤ dim Πq .

Comme Orbp (F) = Orbp (F), on peut échanger les rôles de p et de q dans ce raison-
nement et donc dim Πp = dim Πq .

Notons m = dim Πp . Étant donné un point q ∈ Orbp (F), il est possible de choisir des
champs de vecteurs Y1 , . . . , Ym , tous de la forme ψ∗ X, dont les valeurs en q engendrent
Πq , c’est-à-dire
Vect⟨Y1 (q), . . . , Ym (q)⟩ = Πq .
On considère alors l’application

Gq : Rm → M
(t1 , . . . , tm ) 7→ ϕm 1
tm ◦ · · · ◦ ϕt1 (q)

où ϕit est le flot du champ de vecteurs Yi .

Montrons d’abord que Gq est une immersion locale autour de 0, c’est-à-dire que sa
différentielle en 0 est de rang m. Pour cela, on calcule l’image de la base canonique de
T Rm par la différentielle de Gq :
∂
∂Gq
dt Gq = (t) = dϕm i+1 i−1 1
tm ◦ · · · ◦ dϕti+1 (Yi ) ◦ ϕti−1 ◦ · · · ◦ ϕt1 (q).
∂ti ∂ti

Comme chaque flot vérifie ϕj0 = id et dϕj0 = Id, on obtient

∂Gq
(t = 0) = Yi (q).
∂ti
Les vecteurs tangents Y1 (q), . . . , Ym (q) étant linéairement indépendants, l’application
d0 Gq : T0 Rm → Tq M est de rang m, donc injective. Il existe ainsi un voisinage Uq ⊂ Rm
de 0 sur lequel dGq est injective, ce qui implique que Gq Uq est une immersion. Autrement
dit, Gq (Uq ) est une sous-variété de M de dimension m.
De plus cette sous-variété Gq (Uq ) est incluse dans Orbp (F). En effet, chaque champ
−1
Yi est de la forme ψ∗ X, avec X ∈ F et ψ ∈ G(F) et son flot ϕit = ψ ◦ ϕX t ◦ψ appartient
au groupe G(F). La composition de tels flots appartient également au groupe G(F) et,
pour tout t, Gq (t) ∈ Orbq (F) = Orbp (F).
4.3 Orbite d’une famille de champs de vecteurs 53

On peut enfin déterminer l’espace tangent à la sous-variété Gq (Uq ) de M . Nous allons

montrer pour cela que dt Gq (Tt Rm ) = ΠGq (t) pour tout t ∈ Uq .
Remarquons d’abord que, comme Gq (t) ∈ Orbp (F), il résulte du lemme 4.7 que l’es-
pace vectoriel ΠGq (t) est de dimension m, qui
est également le rang de dt Gq . Il suffit donc
de montrer que, pour i = 1, . . . , m, dt Gq ∂t∂ i ∈ ΠGq (t) . Or

∂

dt Gq = dϕm i+1 i−1 1
tm ◦ · · · ◦ dϕti+1 (Yi ) ◦ ϕti−1 ◦ · · · ◦ ϕt1 (q)
∂ti
= dϕm i+1 i+1 m m 1
tm ◦ · · · ◦ dϕti+1 (Yi ) ◦ ϕ−ti+1 ◦ ϕ−tm ◦ ϕtm ◦ · · · ◦ ϕt1 (q)


= ϕ∗ Yi Gq (t) ∈ ΠGq (t)
i+1
puisque ϕ = ϕm
tm ◦ · · · ◦ ϕti+1 appartient au groupe G(F). Ainsi,

Tq′ Gq (Uq ) = Πq′ ∀q ′ ∈ Gq (Uq ).

En résumé, nous avons montré que Gq (Uq ) est une sous-variété de M qui est incluse
dans Orbp (F) et dont l’espace tangent en un point q ′ est Πq′ .

(ii) Topologie Rappelons d’abord quelques définitions. Une topologie sur un espace E
est un ensemble de parties de E, les ouverts, qui satisfont les trois propriétés suivantes :
— toute réunion d’ouverts est un ouvert ;
— toute intersection finie d’ouverts est un ouvert ;
— E et ∅ sont des ouverts.
Une base d’une topologie est une famille d’ouverts telle que tout ouvert de la topologie
est réunion d’éléments de la base. Une base est donc caractérisée par les deux propriétés
suivantes :
— l’intersection finie d’éléments de la base est égale à la réunion d’éléments de la
base ;
— E et ∅ sont inclus dans la réunion des éléments de la base.
Il nous faut donc montrer que les parties Gq (U ), U ⊂ Uq ouvert de Rm , forment la
base d’une topologie de Orbp (F).
S La deuxième propriété d’une base est trivialement vérifiée : ∅ = Gq (∅) et Orbp (F) =
q∈Orbp (F ) Gq (Uq ). Pour établir la première propriété, il suffit de montrer que, pour tout
point q̄ appartenant à un élément Gq (U ), il existe un autre élément Gq̄ (U ) inclus dans
Gq (U ). En effet on pourra alors écrire toute intersection finie comme réunion d’éléments
de la base :
k
\ [
Gqj (Uj ) = Gq̄ (U ).
j=1 q̄∈∩

Fixons q̄ ∈ Gq (U ). On choisit des champs Y 1 , . . . , Y m de la forme φ∗ X dont les valeurs

en q̄ forment une base de Πq̄ et on définit l’application Gq̄ : t 7→ ϕ̄m 1 i
tm ◦ · · · ◦ ϕ̄t1 (q̄), où ϕ̄ti
est le flot de Y i .
54 Chap. 4 – Familles de champs de vecteurs

Pour tout point q ′ ∈ Gq (U ), on a Y 1 (q ′ ) ∈ Πq′ = Tq′ Gq (U ). Il résulte alors du

corollaire 3.4 que ϕ̄1t1 (q̄) appartient à G
q (U ) pour t1 suffisamment petit. En répétant
cet argument, on montre que ϕ̄2t2 ϕ̄1t1 (q̄) appartient à Gq (U ) pour t1 , t2 suffisamment
petits, . . . et ainsi que Gq̄ (t) ∈ Gq (U ) pour t suffisamment petit. En notant U le voisinage
de 0 dans Rm correspondant à “t suffisamment petit”, on a bien Gq̄ (U ) ⊂ Gq (U ).
Ainsi les Gq (U ), q ∈ Orbp (F), U ⊂ Uq voisinage de 0, définissent une topologie T sur
Orbp (F). On supposera dorénavant que l’orbite Orbp (F) est munie de cette topologie.

(iii) Structure différentiable Considérons l’ensemble des couples (Gq (U ), G−1 q ), pour
q ∈ Orbp (F) et U ⊂ Uq . Chacun de ces couples est une carte sur Orbp (F), les applications
G−1
q étant clairement des homéomorphismes. Ces cartes sont de plus compatibles entre
elles puisque G−1 ′
q ◦ Gq ′ est un difféomorphisme sur U ∩ U . Ces couples forment donc un
atlas et définissent une structure différentiable sur Orbp (F).

(iv) Variété connexe Étudions maintenant les propriétés de Orbp (F) en tant qu’es-
pace topologique (rappelons que Orbp (F) est munie de la topologie T définie au point
(ii)).
Tout d’abord, Orbp (F) est connexe, et même connexe par arcs. En effet, si q1 et q2
sont des points de Orbp (F), on peut écrire q2 sous la forme q2 = ϕX X1
tk ◦ · · · ◦ ϕt1 (q1 ), ce
k

qui implique que q1 et q2 peuvent être joints par le chemin constitué de la trajectoire de
X1 issue de q1 suivie pendant un temps t1 , puis de la trajectoire de X2 pendant un temps
t2 , etc. De plus, chacun des flots t 7→ ϕX
t , X ∈ F, étant une application continue de R
dans Orbp (F), le chemin ainsi construit est continu, ce qui prouve la connexité par arcs.
D’autre part, il est clair que Orbp (F) est un espace topologique séparé. Muni de la
structure différentiable du (iii), il forme donc une variété généralisée (voir le commentaire
page 3), c’est-à-dire une variété différentiable dont la topologie n’est pas forcément à base
dénombrable. Or nous venons de voir que Orbp (F) est connexe. Pour montrer que c’est
une variété il suffit donc de montrer que c’est un espace métrisable (voir encore page 3).
Rappelons d’abord que la variété M elle-même est métrisable : en effet, toute variété
différentiable peut être munie d’une métrique riemannienne. On peut donc choisir une
distance dM sur M , les boules ouvertes BM (q, ε) associées à cette distance étant des
ouverts de M . Cette distance induit une notion de longueur : si γ(t), t ∈ [a, b], est un
chemin continu de M , la longueur P de γ, notée ℓ(γ), est définie comme le supremum de
toutes les sommes de la forme N i=0 dM (γ(ti ), γ(ti+1 )) où a = t0 ≤ t1 ≤ · · · ≤ tN +1 = b
est une partition finie de [a, b] (de façon équivalente, si dM est une distance riemannienne,
ℓ(γ) est la longueur riemannienne de γ).
Définissons alors une distance d sur Orbp (F) en posant, pour q1 , q2 ∈ Orbp (F),

d(q1 , q2 ) = inf ℓ(γ),

l’infimum étant pris sur tous chemins γ reliant q1 à q2 et contenus dans Orbp (F). L’orbite
4.4 Espace tangent à une orbite 55

étant connexe par arcs, d est à valeurs finies et on vérifie aisément que c’est une distance 3 .
Il reste à montrer que d induit la topologie T , ce qui revient à montrer que toute boule
(définie par d) suffisamment petite contient un ouvert de T et que tout ouvert suffisam-
ment petit de T contient une boule. Il suffit de montrer ces propriétés au voisinage d’un
point quelconque de l’orbite.
Soit donc q ∈ Orbp (F) et Gq (Uq ) un voisinage de q dans Orbp (F) qui est une sous-
variété de M (nous avons vu au point (i) qu’il existait de tels voisinages).
Considérons une boule B(q, ε), avec ε > 0. Remarquons que tout point de Gq (Uq )
s’écrit ϕm 1
tm ◦ · · · ◦ ϕt1 (q) et peut donc être vu comme le point final d’un chemin dont
le paramètre évolue de 0 à |t1 | + · · · + |tm |. Or il existe une constante C telle que la
longueur de tout chemin de ce type soit inférieure à C∥t∥, où t = (t1 , . . . , tm ) (c’est une
conséquence triviale du lemme de Gronwall). On peut donc choisir un réel η > 0 tel que
∥t∥ < η implique d(q, Gq (t)) < ε. Autrement dit, en posant U = {∥t∥ < η} ⊂ Rm , la
boule B(q, ε) contient l’ouvert Gq (U ) de la topologie T .
Inversement, puisque Gq (Uq ) est une sous-variété (plongée) de M , les voisinages ou-
verts de q pour la topologie T suffisamment petits sont engendrés par les intersections
BM (q, ε) ∩ Gq (Uq ), avec ε suffisamment petit. Il suffit donc de montrer qu’un tel ouvert
BM (q, ε) ∩ Gq (Uq ) contient une boule pour la distance d. Quitte à réduire ε, on suppose
de plus d(q, ∂Gq (Uq )) > ε. Tout chemin sur Orbp (F) de longueur inférieure à ε issu de q
est donc contenu dans Gq (Uq ). D’autre part, puisque d ≥ dM , la boule B(q, ε) est incluse
dans BM (q, ε). Par conséquent l’ouvert BM (q, ε) ∩ Gq (Uq ) contient la boule B(q, ε).
En conclusion, nous avons montré que Orbp (F) était métrisable. C’est donc une variété
différentiable. Elle est également connexe et son espace tangent en un point q est Πq .

(v) Sous-variété immergée L’inclusion i : Orbp (F) → M est une immersion, car,
sur chaque ouvert de carte Gq (Uq ) de Orbp (F), la restriction de i s’écrit
i|Gq (Uq ) = Gq ◦ G−1
q |Gq (Uq ) ,

c’est-à-dire comme composition d’une immersion et d’un difféomorphisme. Finalement la

variété Orbp (F) est bien une sous-variété immergée de M .

Remarque. Il est important de noter que les orbites ne sont pas toutes de même dimension.
C’était déjà clair pour les orbites d’un champ de vecteurs qui peuvent être de dimension
0 ou 1.

4.4 Espace tangent à une orbite

Quand la famille F est réduite à un unique champ de vecteurs X, l’espace tangent
à une orbite Orbp (F) est simplement la droite vectorielle RX(p) ⊂ Tp M . En revanche
3. Rappelons qu’une distance sur un ensemble N est une fonction d : N ×N → R+ telle que d(x, y) = 0
si et seulement si x = y, qui est symétrique et qui vérifie l’inégalité triangulaire.
56 Chap. 4 – Familles de champs de vecteurs

quand F est constituée de plusieurs champs de vecteurs, il devient difficile de déterminer

la structure de l’espace tangent Tp Orbp (F) en terme de champs de vecteurs de F : la
principale raison est que cet espace n’est pas, en général, déterminé par les propriétés
locales des éléments de F (voir ci-dessous l’exemple page 57).
Le but de cette section est d’étudier Tp Orbp (F) et de voir dans quels cas – et comment
– il dépend des champs de vecteurs de F.
Il est déjà clair que tout champ de vecteurs de F est tangent aux orbites de F. En
général il y en a d’autres : les combinaisons linéaires à coefficients constants de ces champs
bien sûr, mais également leurs crochets de Lie.

Lemme 4.8. Si X et Y sont des champs de vecteurs tangents à une sous-variété N de

M , leur crochet de Lie [X, Y ] est aussi tangent à N .

Preuve. D’après le lemme 3.2, les restrictions X|N et Y |N de X et Y à N sont des

champs de vecteurs sur N . De plus ϕX s (p) est contenu dans N pour tout p ∈ N et est égal
X|N
à ϕs (p) (corollaire 3.4). Écrivons alors le crochet de Lie comme une dérivée de Lie :

1 X
dϕ−t ◦ Y (ϕX


[X, Y ](p) = LX Y (p) = lim t (p)) − Y (p) .
t→0 t

X|
Dans cette expression, ϕXt (p) peut être remplacé par ϕs
N
(p) et Y par Y |N , car Y n’est
évalué qu’en des points ϕX
t (p) ∈ N . Ainsi, [X, Y ](p) = [X|N , Y |N ](p) est un vecteur
tangent à N .

Corollaire 4.9. Pour tout q ∈ Orbp (F), Lieq F ⊂ Tq Orbp (F).

Cette inclusion donne déjà un premier résultat sur les orbites dans un cas particulier
qui sera très utile pour l’étude des systèmes commandés (voir théorème 5.5).

Lemme 4.10. Soient M une variété connexe et F ⊂ X (M ) une famille de champs

de vecteurs telle que, pour tout q ∈ M , Lieq F = Tq M . Alors F a une unique orbite
Orbp (F) = M .

Preuve. Comme dim Lieq F = dim M , chaque orbite est une variété de même dimension
que M , et est donc ouverte dans M . Or M est l’union des orbites disjointes et M est
connexe : il y a donc une unique orbite égale à M tout entière.

Plus généralement, à quelle condition les espaces tangents aux orbites sont donnés par
les valeurs de l’algèbre de Lie engendrée par les champs de F ? Le lemme suivant donne
une condition nécessaire et suffisante.

Lemme* 4.11. Tp Orbp (F) = Liep (F) pour tout p ∈ M si et seulement si ϕX t ∗ Y (p) ∈
X
Liep (F) pour tout flot ϕt d’un élément X de F, tout t ∈ R et tout élément Y de Lie(F).
4.4 Espace tangent à une orbite 57

Preuve. Rappelons que Liep (F) ⊂ Tp Orbp (F), l’égalité est donc équivalente à l’inclusion
dans l’autre sens. Nous allons montrer d’abord que la condition énoncée dans le lemme
est nécessaire puis qu’elle est suffisante.
=⇒ : Supposons que, ∀p ∈ M , Tp Orbp (F) = Liep (F) et considérons Y ∈ Lie(F),
X ∈ F et p ∈ M . Comme Y est tangent aux orbites, la courbe σt (s) = ϕX Y X
t ◦ ϕs ◦ ϕ−t (p)
est incluse dans l’orbite Orbp (F) de p. Son vecteur tangent en s = 0 appartient donc à
l’espace tangent Tp Orbp (F), égal à Liep (F) par hypothèse, c’est-à-dire
dσt
(s = 0) = ϕX
t ∗ Y (p) ∈ Tp Orbp (F) = Liep (F).
ds
⇐= : Supposons maintenant que ϕX t ∗ Y (q) ∈ Lieq (F) pour tout q ∈ M , X ∈ F,
t ∈ R et Y ∈ Lie(F) et considérons ϕ∗ X(p) un élément de Tp Orbp (F), où X ∈ F et
ϕ = ϕktk ◦ · · · ◦ ϕ1t1 appartient à G(F) (c’est-à-dire que chaque ϕit est le flot d’un champ
Xi ∈ F).
Comme X ∈ Lie(F), il résulte de l’hypothèse que ϕ1t1 ∗ X(q) ∈ Lieq (F) pour tout q. Par


conséquent, ϕ2t2 ∗ ϕ1t1 ∗ X (q) appartient aussi à Lieq (F), et, en itérant, ϕ∗ X(q) également.
L’espace tangent à l’orbite Tp Orbp (F) est donc inclus dans Liep (F), ce qui conclut la
preuve.

Exemple. Donnons un exemple de famille pour laquelle Liep (F) ̸= Tp Orbp (F). Prenons
∂ ∂
la famille F = {X1 , X2 } dans R2 avec X1 = ∂x et X2 = f (x) ∂y où

0 si x ≤ 0
f (x) = −1/x2 .
e si x > 0

Il est clair (faire un dessin) que F a une unique orbite égale à R2 tout entier. En revanche,
∂
en q = (x, y) avec x ≤ 0, l’algèbre de Lie de F est Lieq (F) = R ∂x , et est de dimension 1.

La condition énoncée dans le lemme 4.11 permet d’obtenir un théorème qui a des
implications importantes, tant du point de vue théorique (voir le théorème de Frobenius
dans la section suivante) que pour les applications aux systèmes commandés.
Théorème 4.12. Soit F une famille de champs de vecteurs telle que la dimension de
Liep (F) est une constante k indépendante du point p ∈ M .
Alors toutes les orbites de F sont des sous-variétés immergées de dimension k et leur
espace tangent en tout point est Tp Orbp (F) = Liep (F).
Preuve*. La seule chose à montrer est que, sous les hypothèses du théorème, la condition
énoncée dans le lemme 4.11 est vérifiée. On considère donc un champ X ∈ F, de flot ϕ,
un point p ∈ M et on va montrer que, ∀Y ∈ Lie(F) et ∀t, on a ϕt ∗ Y (p) ∈ Liep (F).
Choisissons k champs de vecteurs Y1 , . . . , Yk dans Lie(F) dont les valeurs en p forment
une base de Liep (F). Comme la dimension de Lieq (F) est constante, Y1 (q), . . . , Yk (q)
forment une base de Lieq (F) pour tout point q dans un voisinage U de p.
58 Chap. 4 – Familles de champs de vecteurs

Remarquons alors que, pour t suffisamment petit, le point qt = ϕ−t (p) appartient à
U et que le vecteur tangent Y (qt ) est combinaison linéaire des vecteurs Yi (qt ). Comme
ϕt ∗ Y (p) = dϕt (Y (qt )), ceci implique que ϕt ∗ Y (p) est combinaison linéaire des ϕt ∗ Yi (p).
Ainsi, pour t suffisamment petit, il suffit de montrer que ϕt ∗ Yi (p) ∈ Liep (F) pour i =
1, . . . , k.
Considérons la courbe Wi (t) = ϕt ∗ Yi (p) dans Tp M . Nous allons montrer qu’elle est
solution d’une équation différentielle linéaire dans Tp M . En effet :
dWi d

(t) = ϕt ∗ Yi (p) = ϕt ∗ LX Yi (p) = ϕt ∗ [X, Yi ](p).
dt dt
Comme [X, Yi ] est un élément de Lie(F) et que ϕ−t (p) appartient à U , il existe des
fonctions différentiables fij sur M telles que
k
X
[X, Yi ](ϕ−t (p)) = fij (ϕ−t (p))Yj (ϕ−t (p)).
j=1

Par conséquent
k
dWi X j
(t) = fi (ϕ−t (p))ϕt ∗ Yj (p),
dt j=1

ou encore Ẇi (t) = kj=1 aji (t)Wj (t).

P

Notons A(t) la matrice k × k de composantes aji (t) et M (t) = mji (t) la matrice k × k

solution de l’équation différentielle linéaire matricielle


Ṁ (t) = A(t)M (t),
M (0) = Ik
(autrement dit M (t) = RA (t, 0), où RA (t, 0) est la P résolvante de l’équation x′ (t) =
k j P j
A(t)x(t) dans R ). On a alors, pour tout i, Wi (t) = j mi (t)Wj (0) = j mi (t)Yj (p), ce
qui implique que tous les Wi (t) appartiennent à Liep (F) pour t suffisamment petit.
On a donc montré que, pour tous X ∈ F, Y ∈ Lie(F) et p ∈ M , ϕX t ∗ Y (p) ∈ Liep (F)
pour t suffisamment petit.
Pour t quelconque, la conclusion s’obtient en utilisant la compacité de la courbe ϕX s (p),
X
0 ≤ s ≤ t, et en écrivant ϕt ∗ Y (p) sous la forme
ϕX X X



t ∗ Y (p) = ϕtl ∗ · · · ϕt1 ∗ Y (p).

4.5 Distributions et mécanique non-holonome*

Définition. Une distribution ∆ ⊂ T M sur M est une famille de sous-espaces vectoriels
∆q ⊂ Tq M dépendant de façon différentiable de q ∈ M .
4.5 Distributions et mécanique non-holonome* 59

Autrement dit, on peut toujours trouver localement des champs de vecteurs formant
une base de ∆. Ce qui implique que la dimension des sous-espaces ∆q est constante sur
M . Cette dimension dim ∆q est appelée le rang de la distribution. On parle parfois d’une
distribution comme d’un champ de sous-espaces vectoriels tangents à M .
Exemples.
— Le fibré tangent T M est l’unique distribution de rang n.
— Un champ de vecteurs sur M définit une distribution de rang 1 si il ne s’annule
jamais sur M .
— Plus généralement, une famille F de champs de vecteurs peut engendrer une distri-
bution si les sous-espaces Fq = Vect{X(q), X ∈ F} sont de dimension constante
sur M .
Inversement, on peut associer à toute distribution ∆ la famille des champs de vecteurs
qui y sont tangents

F∆ = X ∈ X (M ), X(q) ∈ ∆q ∀q ∈ M

Étant donnée une distribution, il est naturel de se demander si elle n’est pas l’espace
tangent à une variété.
Définition. Une distribution ∆ ⊂ T M est intégrable si, pour tout point p ∈ M , il existe
une sous-variété immergée Np ⊂ M contenant p et telle que Tq Np = ∆q pour tout point
q ∈ Np . On appelle Np une variété intégrale de ∆ en p.

Si ∆ est une distribution intégrable, tout champ de F∆ est tangent aux variétés
intégrales Np . Ainsi, localement, l’orbite de F∆ par p est incluse dans Np . Vu le corol-
laire 4.9, ceci implique que Liep (F∆ ) = ∆p . Cette propriété est l’involutivité.
Définition. Une distribution ∆ est involutive si, pour tous champs de vecteurs X et Y
dans F∆ , on a [X, Y ](p) ∈ ∆p pour tout p ∈ M .

Cette propriété est en fait équivalente à l’intégrabilité, comme l’indique le théorème

suivant.
Théorème 4.13 (Théorème de Frobenius). Une distribution est intégrable si et seulement
si elle est involutive.
Preuve. Nous venons de montrer que l’intégrabilité implique l’involutivité. Inversement,
supposons que ∆ est une distribution involutive. Alors l’algèbre de Lie Lieq (F∆ ) = ∆q
est de rang constant sur M . Vu le théorème 4.12, ceci implique que toutes les orbites de
F∆ sont des variétés intégrales de ∆, ce qui achève la preuve.

Remarques.
— Notons que les théorèmes de Frobenius et 4.12 sont équivalents puisque, sous les
hypothèses de ce dernier, l’algèbre de Lie Lie(F) est une distribution involutive.
60 Chap. 4 – Familles de champs de vecteurs

— Le théorème de Frobenius apparaı̂t ici comme une conséquence (ou plutôt un cas
particulier) du théorème de l’orbite. Historiquement, ces deux théorèmes sont ap-
parus dans l’ordre inverse, puisque le théorème de l’orbite date des années 70 alors
que le théorème de Frobenius est un résultat classique de la géométrie différentielle
et date du début du 20ème siècle. La preuve classique de ce dernier théorème est
nettement plus simple que celle du théorème de l’orbite et a son intérêt propre : elle
fait intervenir les propriétés des familles de champs de vecteurs formant localement
des champs de coordonnées (voir l’exercice A.1).

Application à la mécanique non-holonome La théorie des distributions a des ap-

plications importantes en mécanique lagrangienne (voir les cours AOT11 et AOT12).
Dans notre contexte, un système lagrangien est un couple (M, L) composé de :
— un espace des configurations, noté M : c’est une variété différentiable et son fibré
tangent T M est appelé l’espace des phases ;
— un lagrangien L ∈ C ∞ (T M ), c’est-à-dire une fonction différentiable de T M dans
R.
Rappelons que des coordonnées locales q 7→ (x1 , . . . , xn ) sur M induisent une base
( ∂x∂ 1 |q , . . . , ∂x∂n |q ) sur Tq M et donc des coordonnées locales sur T M

(q, v) 7→ (x1 , . . . , xn , ẋ1 , . . . , ẋn ).

On appelle mouvement du système lagrangien une courbe q(t) ∈ M qui satisfait le

principe de d’Alembert-Lagrange, que l’on écrit localement : pour tout t et pour tout
ξ = (ξ 1 , . . . , ξ n ) ∈ Tq(t) M ,
n     
X d ∂L dq ∂L dq
i
q(t), (t) − i q(t), (t) ξ i = 0.
i=1
dt ∂ ẋ dt ∂x dt

Comme ξ est arbitraire ici, le principe de d’Alembert-Lagrange est équivalent à l’équation

d’Euler-Lagrange :
   
d ∂L dq ∂L dq
q(t), (t) − i q(t), (t) = 0 pouri = 1, . . . , n.
dt ∂ ẋi dt ∂x dt

Considérons maintenant un système lagrangien sous contraintes : la vitesse des mou-

vements du système lagrangien appartient à un sous-ensemble ∆q de Tq M en tout point.
Pour simplifier, on suppose que ∆ forme une distribution de rang m sur M . On distingue
alors deux cas.
— Si ∆ est intégrable (c’est-à-dire involutive par le théorème de Frobenius), on dit
que les contraintes sont holonomes. Dans ce cas, tout point q ∈ M est inclus dans
une unique sous-variété intégrale Nq de dimension m avec T Nq = ∆|Nq et chaque
mouvement du système lagrangien reste dans une variété intégrale.
4.5 Distributions et mécanique non-holonome* 61

— Si ∆ n’est pas intégrable, on dit que les contraintes sont non-holonomes. Chaque
mouvement du système reste alors dans une orbite de la distribution (c’est-à-dire
une orbite de la famille F∆ ), qui est une variété de dimension k > m.
Supposons de plus que Lie(F∆ ) soit une distribution (c’est-à-dire qu’elle est de rang
constant), ce qui est automatiquement vérifiée pour les contraintes holonomes. Nous
pouvons alors écrire les contraintes en coordonnées locales.
Soit k la dimension des orbites de F∆ (cette dimension ne dépend pas de l’orbite
d’après le théorème 4.12) : k = m si les contraintes sont holonomes, m < k ≤ n sinon.
Il existe alors (voir exercice A.1) un système local de coordonnées (x1 , . . . , xn ) et des
fonctions fm+1 , . . . , fk ∈ C ∞ (Rn+m ) tels que :
— le plan tangent à l’orbite Tx Orbx (F∆ ) a pour équation

ẋk+1 = · · · = ẋn = 0;

— la distribution ∆ a localement pour équation :



 ẋm+1 = fm+1 (x1 , . . . , xn , ẋ1 , . . . , ẋm )
 ..

.

 ẋk = fk (x1 , . . . , xn , ẋ1 , . . . , ẋm )
 ẋk+1 = · · · = ẋn = 0


Autrement dit, ∆x est paramétré par ẋ1 , . . . , ẋm .

Dans ces coordonnées, le principe de d’Alembert-Lagrange s’écrit : pour tout t et pour
tout ξ = (ξ 1 , . . . , ξ m ),
m     
X d ∂L dq ∂L dq
i
q(t), (t) − i q(t), (t) ξ i +
i=1
dt ∂ ẋ dt ∂x dt
k     
X d ∂L dq ∂L dq
i
q(t), (t) − i q(t), (t) fi (x, ξ) = 0.
i=m+1
dt ∂ ẋ dt ∂x dt

Quand les contraintes sont holonomes (m = k), ceci implique que les mouvements du
systèmes satisfont
   
 d ∂L q(t), dq (t) − ∂L q(t), dq (t) = 0 pour i = 1, . . . , m,

dt ∂ ẋi dt ∂xi dt
m+1 n

q(t) ∈ Nq(0) = {x = · · · = x = 0}.

Le principe de d’Alembert-Lagrange reste alors équivalent à l’équation d’Euler-Lagrange

sur chaque variété intégrale de la distribution des contraintes.
En revanche, pour les contraintes non-holonomes, l’équation d’Euler-Lagrange n’est
plus valable.
62 Chap. 4 – Familles de champs de vecteurs
Chapitre 5

Systèmes commandés affines*

5.1 Introduction
Nous allons utiliser les outils de géométrie différentielle que nous avons introduits pour
analyser les systèmes dynamiques commandés, et plus particulièrement ceux provenant
de la robotique.
Nous nous limiterons dans ce cours aux systèmes commandés dits affines.
Définition. Un système commandé affine, ou simplement système affine, est un système
commandé de la forme

q̇ = X0 (q) + u1 X1 (q) + · · · + um Xm (q), q ∈ M, u ∈ U, (5.1)

où X0 , . . . , Xm sont des champs de vecteurs sur la variété M et U un sous-ensemble de

Rm (a priori quelconque). Le champ X0 est appelé la dérive du système.

Ces systèmes représentent une classe très générale de systèmes commandés. En effet :
— on peut considérer un système de commandé comme une contrainte sur la vitesse
du type q̇ ∈ Vq , où Vq ⊂ Tq M . Si les sous-ensembles Vq forment une distribution
(voir § 4.5), ou plus généralement un champ de sous-espaces affines, choisir des
champs de vecteurs engendrant Vq donne une caractérisation du type (5.1).
— Du point de vue de la mécanique classique, qui est généralement celui utilisé en
robotique, toute commande résulte de l’action d’une force ou d’un couple. Elle agit
donc linéairement sur l’accélération du système. Le système commandé doit alors
s’écrire dans l’espace des phases et prend la forme (5.1). Par exemple, le système
ẍ = u dans R se réécrit dans T R ≃ R2 sous la forme :
     
ẋ y 0
q̇ = = +u = X0 (q) + uX1 (q).
ẍ 0 1

— Un autre point de vue, typique de la robotique, est de considérer que l’on dispose de
m actionneurs, l’action de chacun d’eux sur le système étant décrite par q̇ = Xi (q).
64 Chap. 5 – Systèmes commandés affines*

En les combinant et en tenant compte de la vitesse X0 (q) du système, on obtient

à nouveau un système affine.
— Enfin, prenons un système commandé quelconque (dans Rn ou en coordonnée
locales) ẋ = f (x, u) utilisant des lois de commande dérivables. En posant v = u̇ et
y = (x, u), on retrouve un système affine ẏ = X0 (y) + vX1 (y).
Dans tout le chapitre nous nous appuierons sur les conventions et notations intro-
duites dans la section 4.1. Ainsi, nous considérerons des lois de commande constantes par
morceaux. L’ensemble des points atteignables par le système à partir de p à l’aide de ces
lois de commande est noté Ap . Notons cependant que cette hypothèse sur la classe des
lois de commande n’est pas restrictive : pour un système affine, tout point atteignable
à partir de p au moyen de commandes mesurables l’est aussi au moyen de commandes
constantes par morceaux (nous ne montrerons pas ce résultat dans ce cours, sa preuve
sort du cadre de la géométrie différentielle étudié ici).
Définition. Un système commandé est dit commandable si Ap = M quel que soit le
point p ∈ M .

Le but de ce chapitre est d’étudier les ensembles atteignables des systèmes affines et
d’en déduire des conditions de commandabilité.

5.2 Ensembles atteignables

Considérons donc un système affine de la forme (5.1). Ce système correspond à une
famille de champs de vecteurs
n m
X o
F(U ) = X0 + u i Xi , u∈U .
i=1

L’ensemble atteignable Ap est alors un sous-ensemble de l’orbite de F(U ) passant par p.

Dans ce chapitre nous
supposerons toujours que F(U ) est Lie-déterminée, c’est-à-dire
que l’algèbre Liep F(U ) a une dimension constante sur M . D’après le théorème 4.12,
cette hypothèse assure que toutes les orbites de la famille F(U ) sont
des sous-variétés
immergées de même dimension dont l’espace tangent est Liep F(U ) .

la famille F(U )|Orbp (F ) ,
Quitte à se restreindre à une orbite Orbp (F) et à considérer
nous supposerons dorénavant que la dimension de Liep F(U ) est égale à la dimension
de M en tout point (l’orbite de tout point est donc la variété M tout entière).
La propriété fondamentale du sous-ensemble Ap de M est la suivante.


Théorème 5.1. Soit un système affine tel que Liep F(U ) = Tp M en tout point p ∈ M .
Alors, ∀p ∈ M ,
Ap ⊂ intAp .
Autrement dit l’ensemble atteignable a un intérieur non vide. Il peut être :
5.2 Ensembles atteignables 65

— un ensemble ouvert, donc une variété de dimension n ;

— une variété à bord : intAp est une variété de dimension n et ∂Ap une variété de
dimension n − 1 ;
— une variété de dimension n avec un bord ayant des singularités (coins, cusps).
En revanche Ap ne peut pas être :
— un sous-ensemble d’une variété de plus petite dimension ;
— un ensemble dont les points du bord sont isolés des points intérieurs.
Exercice 5.1. Quelle est la structure de l’ensemble atteignable A0 pour les systèmes affines
dans R2 définis par
∂ ∂ ∂
  
X0 = y ∂x X0 = ∂x X0 = ∂x
∂ , ∂ , ∂ , avec U = [−1, 1].
X1 = ∂y X1 = ∂y X1 = x ∂y

Preuve (du théorème 5.1). Fixons p ∈ M et notons F = F(U ). Prenons un point q ′ ∈ Ap

et montrons qu’il appartient aussi à intAp .
— Si dim M = dim Lieq′ (F) ̸= 0, il existe Y1 ∈ F tel que Y1 (q ′ ) ̸= 0. Soit ϕ1s1 le flot
de Y1 . La courbe s1 7→ ϕ1s1 (q ′ ), s1 ∈]0, ε[, est une variété de dimension 1 pour ε > 0
suffisamment petit et est incluse dans Ap .
Si dim M = 1, alors ϕ1s1 (q ′ ) ∈ intAp pour s1 > 0 suffisamment petit, d’où le
théorème.
— Si dim M = dim Lieq′ (F) > 1, alors, arbitrairement près de q ′ , on peut trouver un
point q1 sur la courbe ϕ1s1 (q ′ ) et un champ Y2 ∈ F tels que Y2 (q1 ) n’est pas tangent
à la courbe ϕ1s1 (q ′ ) :
q1 = ϕ1t1 (q ′ ), t11 suffisamment petit,

1
Y1 (q1 ) et Y2 (q1 ) linéairement indépendants
(en effet, sinon tous les champs de F sont tangents à la courbe pour s1 petit, donc
tous leurs crochets aussi, ce qui entraı̂nerait dim Lie(F)|courbe = 1).
Il existe donc un voisinage V (t11 , 0) de (t11 , 0) dans R2 tel que l’application
R2+ ∩ V (t11 , 0) → Ap
(s1 , s2 ) 7→ ϕ2s2 ◦ ϕ1s1 (q ′ )
est une immersion. Son image est donc une surface (une sous-variété de dimension
2) de M . Si dim M = 2, la preuve est terminée.
— Si dim M > 2, par le même raisonnement que précédemment, on peut trouver un
point q2 = ϕ2t2 ◦ ϕ1t1 (q ′ ) arbitrairement proche de q ′ sur la surface et Y3 ∈ F tels
2 1
que Y3 (q2 ) ne soit pas tangent à la surface. L’application
R3+ ∩ V (t11 , t22 , 0) → Ap
(s1 , s2 , s3 ) 7→ ϕ3s3 ◦ ϕ2s2 ◦ ϕ1s1 (q ′ )
est une immersion et son image une sous-variété de dimension 3.
66 Chap. 5 – Systèmes commandés affines*

Ainsi, par récurrence, on construit une sous-variété de M de dimension n – c’est-à-

dire un ouvert – incluse dans Ap et arbitrairement proche de q ′ . Ceci montre bien que
q ′ ∈ intAp .

Remarque. Cette preuve montre également que, pour tout ouvert Ω ⊂ M et tout point
p ∈ Ω, Ap ∩ Ω contient un ouvert.

Corollaire 5.2. Sous les hypothèses du théorème précédent, si Ap = M pour un point

p ∈ M , alors Ap = M .

Preuve. Fixons q0 ∈ M . Nous allons montrer qu’il appartient à Ap . Considérons pour

cela la famille de champs de vecteurs

−F = {−X, X ∈ F}.

Comme −F satisfait les mêmes propriétés que F, son ensemble atteignable vérifie, pour
tout q ∈ M ,
Aq (−F) ⊂ intAq (−F).
Prenons alors un point quelconque q ∈ intAq0 (−F) et un voisinage de ce point V (q) ⊂
Aq0 (−F). Comme Ap est dense dans M , l’intersection de Ap et de V (q) est non vide, ce
qui implique que celle de Ap et de Aq0 (−F) est non vide également. Il existe donc un
point q ′ dans cette intersection, qui peut s’écrire de deux façons différentes :

q ′ = ϕYtkk ◦ · · · ◦ ϕYt11 (p), Yi ∈ F, ti > 0,

= ϕ−W
sl
l
◦ · · · ◦ ϕ−W
s1 (q0 ),
1
Wi ∈ F, si > 0.

−1
Comme ϕ−W
si
i
= ϕW
si
i
, cette égalité s’écrit
Yk Y1
q0 = ϕW Wl
s1 ◦ · · · ◦ ϕsl ◦ ϕtk ◦ · · · ◦ ϕt1 (p)
1

ce qui implique q0 ∈ Ap .

Ce corollaire est important pour l’étude de la commandabilité, puisqu’il permet de

remplacer Ap par son adhérence.
Enfin, nous allons expliciter la dépendance de Ap par rapport à la famille F. Énonçons
d’abord un résultat sur l’approximation d’un flot.

Théorème 5.3. Soient X et Y deux champs de vecteurs sur une variété M , p un point
de M , λ et µ des réels positifs et T > 0 tel que le flot ϕt (p) du champ λX + µY soit défini
sur l’intervalle [0, T ].
Alors il existe une suite γ N (t) de trajectoires de la famille {X, Y }, avec γ N (0) = p,
qui converge uniformément vers ϕt (p) sur [0, T ] quand N → +∞.
5.3 Systèmes sans dérive 67

Nous admettrons ce résultat ici, sa preuve sort du cadre de ce cours (on peut la trouver
dans [1, 4]).
Ce théorème implique que l’on peut ajouter à F toute combinaison linéaire positive
de ses éléments sans changer la adhérence de Ap .
Lemme 5.4. Pour un système affine Lie-déterminé,


Ap (F) = Ap cone(F) ,
où cone(F) est le cône convexe positif engendré par F :
k
nX o
cone(F) = λi Xi , λi ∈ C ∞ (M ), λi ≥ 0, Xi ∈ F, k ∈ N ⊂ X (M ).
i=1

5.3 Systèmes sans dérive

Un cas particulier important est celui des systèmes affines sans dérive :
q̇ = u1 X1 (q) + · · · + um Xm (q), q ∈ M, u ∈ U.
La famille associée à ce système est F(U ) = { m
P
i=1 ui Xi , u ∈ U }.

5.3.1 Conditions suffisantes de commandabilité

Pour les systèmes affines sans dérive, la commandabilité se déduit essentiellement du
théorème suivant, qui est une conséquence directe du lemme 4.10.
Théorème 5.5 (Théorème de Chow). Soit F = {Xu , u ∈ U } une famille symétrique
telle que
dim Lieq F = dim M, pour tout q ∈ M. (5.2)
Alors, pour le système commandé correspondant, l’ensemble atteignable à partir d’un
point p est la composante connexe de son espace de configuration contenant p.
Remarque. Ce théorème a été établi en 1938 indépendamment par Chow et Rashevsky.
Nous avons choisi de le présenter ici comme une conséquence du théorème de l’orbite (qui
date de 1973) mais il en serait plutôt le précurseur. La preuve originale du théorème de
l’orbite s’appuie d’ailleurs sur la preuve du théorème de Chow.
La condition (5.2) est appelée condition du rang (ou condition de Chow). Notons
qu’une famille qui satisfait la condition du rang est Lie-déterminée. Inversement, on a
vu au § 5.2 que, quitte à se restreindre à une orbite, on peut toujours supposer qu’une
famille Lie-déterminée satisfait la condition du rang.
Le théorème de Chow donne directement une condition de commandabilité pour les
systèmes sans dérive quand U est Rm tout entier.
68 Chap. 5 – Systèmes commandés affines*

Proposition 5.6. Soit un système affine sans dérive sur une variété M connexe et dont
l’ensemble de commande est U = Rm . Alors, si dim Liep (X1 , . . . , Xm ) = dim M pour tout
p ∈ M , le système est commandable.

Preuve. La famille associée au système est ici F(Rm ) = Vect(X1 , . . . , Xm ) et son algèbre
de Lie est égale à l’algèbre de Lie engendrée par X1 , . . . , Xm . Ainsi F(Rm ) est symétrique
et satisfait la condition du rang, et on obtient la conclusion par le théorème de Chow.

Ce résultat reste valable avec des hypothèses beaucoup moins restrictives sur l’en-
semble de commande.

Théorème 5.7. Soit un système affine sans dérive sur une variété connexe M . Supposons
que les deux conditions suivantes sont vérifiées
— l’enveloppe convexe de l’ensemble de commande U contient un voisinage de l’ori-
gine dans Rm ,
— dim Liep (X1 , . . . , Xm ) = dim M pour tout p ∈ M .
Alors le système est commandable.

Preuve. Quitte à restreindre U , on peut supposer que U est une boule B(0, ε) dans Rm .
La famille
F(U ) est alors symétrique et satisfait la condition du rang puisque les algèbres
Lie F(U ) et Lie(X1 , . . . , Xm ) sont égales. La commandabilité résulte alors du théorème
de Chow.

Exemple. Reprenons l’exemple de la voiture (voir page 40), en autorisant maintenant

seulement 4 commandes, (±1, 0) et (0, ±1). L’enveloppe convexe de cet ensemble de
commande étant [−1, 1]2 et Lie(X1 , X2 ) étant toujours de dimension 3, le système est
commandable.

5.3.2 Problème de la sphère en roulement

Intéressons-nous à un problème de manipulation d’objet (c’est une classe de problèmes
classique en robotique). On pose une balle de tennis sur une table et on la fait rouler
entre la table et la paume de la main ouverte (on peut aussi imaginer une bille posée sur
la paume de la main que l’on fait rouler avec le pouce). On suppose que la balle roule
sans glisser ni pivoter. Peut-on amener la balle n’importe où sur la table avec n’importe
quelle orientation ?
Commençons par formaliser cette question comme un problème de commandabilité.
On considère la balle comme la sphère unité et la table comme un plan horizontal dans
l’espace euclidien E3 . Soit O un point fixe quelconque du plan et C le centre de la sphère.
On choisit deux repères dans E3 :
— (O, e1 , e2 , e3 ) est un repère orthonormé fixe, où (e1 , e2 ) sont horizontaux et e3 est
dirigé vers le haut ;
5.3 Systèmes sans dérive 69

— (C, a1 , a2 , a3 ) est un repère orthonormé mobile lié à la sphère.

Nous noterons en minuscule les coordonnées d’un point dans la base fixe (p ∈ R3
représente donc les coordonnées du point P ∈ E3 ).
Une configuration (position et orientation) de la sphère est caractérisée par :
— la matrice orthogonale R de transition entre la base fixe (e1 , e2 , e3 ) et la base mobile
(a1 , a2 , a3 ), c’est-à-dire R = (a1 a2 a3 ) dans la base (e1 , e2 , e3 ) ;
— le centre C de la sphère, appartenant à un plan horizontal E2 .
L’espace des configurations est donc la variété de dimension 5

M = E2 × SO(3).

Écrivons maintenant les équations du mouvement. Soit P le point de contact de la

sphère et du plan à l’instant t. La condition du roulement sans glissement se traduit en
−→
coordonnées fixes par ṗ(t) = 0. Comme p = c + CP , on obtient
−→
dCP
ċ(t) + (t) = 0.
dt
−→
D’autre part, P étant un point de la sphère, CP (t) = R(t)p0 , et l’équation du mouvement
s’écrit
ċ(t) = −Ṙ(t)p0 .
Avant de continuer le calcul, détaillons la structure de T M . Comme M est un produit
de variété, son espace tangent est le produit des espaces tangents : T(C,R) M = TC E2 ×
TR SO(3), et, d’après l’exercice 2.2,

TR SO(3) = {ΩR, Ω matrice anti-symétrique 3 × 3}.

Puisque Ṙ ∈ TR SO(3), il existe une courbe Ω(t) de matrices antisymétriques telle que
Ṙ(t) = Ω(t)R(t).
L’équation du mouvement devient donc
−→
ċ(t) = −Ω(t)R(t)p0 = −Ω(t)CP (t).
−→
Enfin, CP (t) = −e3 puisque P est le point de contact à l’instant t, ce qui donne

ċ(t) = Ω(t)e3 .

Pour mettre cette équation sous la forme d’un système commandé, choisissons une
base (A1 , A2 , A3 ) de l’ensemble des matrices anti-symétriques :
     
0 0 0 0 0 1 0 −1 0
A1 = 0 0 −1 , A2 =  0 0 0 , A3 = 1 0 0
0 1 0 −1 0 0 0 0 0
70 Chap. 5 – Systèmes commandés affines*

Les champs de vecteurs A1 R, A2 R, A3 R sur SO(3) forment une base de TR SO(3) en tout
point. Définissons maintenant des champs de vecteurs sur M :

X1 (c, R) = A1 e3 ⊕ A1 R = −e2 ⊕ A1 R
X2 (c, R) = A2 e3 ⊕ A2 R = e1 ⊕ A2 R
X3 (c, R) = A3 e3 ⊕ A3 R = 0 ⊕ A3 R

En notant q = (c, R) les points de M , on en déduit que tout mouvement de la sphère

satisfait la contrainte 
q̇ ∈ Vect X1 (q), X2 (q), X3 (q)
ou, de façon équivalente,

q̇ = u1 X1 (q) + u2 X2 (q) + u3 X3 (q), u ∈ R3 .

On obtient ainsi un système commandé affine sans dérive.

Nous n’avons tenu compte pour le moment que de la contrainte de roulement sans
glissement. Il faut rajouter celle de non pivotement qui s’écrit tout simplement u3 = 0.
Le système commandé décrivant le roulement de la sphère sur un plan est donc :

q̇ = u1 X1 (q) + u2 X2 (q), u ∈ R2 , (5.3)

et la question initiale devient : est-ce que ce système est commandable ?

D’après la proposition 5.6, si l’algèbre de Lie engendrée par X1 et X2 est partout de
dimension 5, le système (5.3) sera commandable. C’est ce que nous allons montrer.
Commençons par calculer le crochet de Lie de X1 et X2 :

[X1 , X2 ]M = −[e2 , e1 ]E2 ⊕ [A1 R, A2 R]SO(3) .

Comme les champs e1 et e2 sont constants, leur crochet est nul : [e2 , e1 ] = 0. D’autre
part, pour le crochet dans SO(3) des champs ξi = Ai R,
d ξ1

[ξ1 , ξ2 ] = Lξ1 ξ2 = ϕ ξ2 .
dt −t ∗ t=0

Comme ξ1 = A1 R est linéaire, son flot ϕξs1 = exp(sA1 )R est également linéaire. Sa
différentielle est dϕξs1 R · Ω = exp(sA1 )Ω et donc

ϕξ−t1 ∗ ξ2 = exp(−tA1 )A2 exp(tA1 )R,

ce qui implique que

[A1 R, A2 R] = (−A1 A2 + A2 A1 )R = −[A1 , A2 ]R

où [A, B] = AB − BA désigne le commutateur des matrices (c’est en fait un résultat très
général sur les groupes de Lie).
5.3 Systèmes sans dérive 71

Or les commutateurs des matrices Ai sont :

[A1 , A2 ] = A3 , [A3 , A1 ] = A2 , [A2 , A3 ] = A1 ,

ce qui donne finalement

[X1 , X2 ] = 0 ⊖ A3 R = −X3 .
Le même calcul montre immédiatement que

[X1 , [X1 , X2 ]] = 0 ⊕ A2 R
[X2 , [X1 , X2 ]] = 0 ⊖ A1 R.

Il est clair que les champs X1 , X2 , [X1 , X2 ], [X1 , [X1 , X2 ]] et [X2 , [X1 , X2 ]] engendrent
toujours un espace de dimension 5, ce qui montre que dim Lieq (X1 , X2 ) = 5 pour tout
q ∈ M.
Ainsi le système (5.3) est commandable : il est toujours possible d’amener la balle
dans une position donnée avec une orientation donnée.

5.3.3 Autres notions de commandabilité

Les systèmes sans dérive ont des propriétés de commandabilité plus fortes que la
propriété habituelle Ap = M .
Étant donné un système commandé sur M , un point p ∈ M et un réel T > 0, on note
Ap,T l’ensemble des points atteignables à partir de p en temps inférieur ou égal à T . Le
système est dit fortement commandable si, pour tous p ∈ M et T > 0, on a Ap,T = M .
Alors un système affine sans dérive
— est fortement commandable si l’enveloppe convexe de U est égale à Rm tout entier ;
— n’est pas fortement commandable quand U est borné.
(La preuve des ces deux résultats est laissée en exercice).
Une autre propriété des systèmes affines sans dérive est la commandabilité locale : on
peut joindre deux points d’un ouvert sans sortir de l’ouvert.
Lemme 5.8. Soit q̇ = m
P
i=1 ui Xi , q ∈ M , u ∈ U , un système affine sans dérive satisfai-
sant les hypothèses du théorème 5.7.
Alors, pour tout ouvert connexe Ω ⊂ M , deux points quelconques de Ω peuvent être
joints par une trajectoire du système qui reste dans Ω.
Preuve. L’ouvert Ω est une sous-variété ouverte et connexe de M . Les restrictions Xi |Ω
sont donc des champs de vecteurs sur Ω qui satisfont, pour tout p ∈ Ω,

Liep (X1 |Ω , . . . , Xm |Ω ) = Liep (X1 , . . . , Xm ).

Donc le système affine sans dérive q̇ = m

P
i=1 ui Xi |Ω , q ∈ Ω, u ∈ U vérifie les hypothèses
du théorème 5.7 : il est donc commandable, ce qui montre le résultat.
72 Chap. 5 – Systèmes commandés affines*

Application à la planification des mouvements Imaginons un robot – par exemple

une voiture – évoluant au milieu d’obstacles, ou plus généralement un système localement
commandable placé dans un milieu encombré. On se demande à quelle condition sur
les obstacles le système reste commandable, et, le cas échéant, comment trouver une
trajectoire reliant deux configurations données.
On commence d’abord par représenter les obstacles comme des zones interdites de
l’espace des configurations du robot, c’est-à-dire un fermé O ⊂ M (fermé car on n’autorise
pas les configurations au contact des obstacles).
Puisque notre système est localement commandable, il suffit que M \O soit connexe
pour que deux configurations q0 et q1 ∈ M \O quelconques soient joignables sans entrer
en collision avec un obstacle. Trouver une trajectoire reliant q0 et q1 se fait alors en deux
temps :
— trouver un arc dans M \O reliant q0 à q1 (ce n’est en général pas une trajectoire
du système) ;
— approcher cet arc par une trajectoire, c’est-à-dire trouver une trajectoire restant
dans un voisinage tubulaire de l’arc et joignant q0 et q1 .
Cette méthode est toujours réalisable. En effet, la première étape est possible dès que
M \O est connexe (rappelons qu’un ouvert connexe est connexe par arcs). Pour la seconde
étape, comme M \O est ouvert, il existe un voisinage Ω de l’arc contenu dans M \O.
Comme le système est supposé localement commandable et que Ω peut être supposé
connexe (il contient l’arc), q0 et q1 sont joignables par une trajectoire restant dans Ω.
Remarque. Comparer avec le cas de la voiture de Dubins, page 75.

5.4 Systèmes affines commandables

Revenons maintenant au cas général d’un système affine avec dérive :

q̇ = X0 (q) + u1 X1 (q) + · · · + um Xm (q), q ∈ M, u ∈ U.

Nous allons voir différentes conditions de commandabilité pour ces systèmes.

5.4.1 Commandes non bornées

Lemme 5.9. Soit un système affine sur une variété connexe M . Supposons que les deux
conditions suivantes sont vérifiées
— U = Rm (pas de restriction sur les commandes),
— dim Liep (X1 , . . . , Xm ) = dim M pour tout p ∈ M .
Alors le système est commandable (et même fortement commandable).
Preuve. Commençons par remarquer que, si F désigne l’adhérence d’une famille F de
champs de vecteurs, on a :
Ap (F) ⊂ Ap (F).
5.4 Systèmes affines commandables 73

Ce résultat est en fait la conséquence d’une version de la dépendance continue par rapport
aux paramètres pour les équations différentielles : soit {X N } une suite de champs de
vecteurs sur M qui converge vers X quand N → +∞ et ϕX t (p) le flot de X à partir de p
XN
défini pour t ∈ [0, T ] ; alors les flots ϕt (p) sont définis et convergent uniformément vers
ϕXt (p) sur [0, T ].
= F(U ) = {X0 + i ui Xi , u ∈ Rm }. Pour tout λ > 0,
P
Revenons maintenant à FP
les champs
P de vecteurs λX0 + i λui Xi appartiennent au cône cone(F), donc les champs
λX0 + i vi Xi aussi : il suffit de prendre ui = vi /λ, ce qui est toujours possibleP puisque
les ui ne sont pas bornés. En faisant tendre λ vers 0, on obtient que les champs i vi Xi ,
v ∈ Rm , appartiennent à cone(F).
Ainsi cone(F) est la famille associée à un système sans dérive dont l’algèbre de Lie
est de dimension n = dim M partout. Ceci implique Ap (cone(F)) = M (proposition 5.6).
Nous avons remarqué au début de cette preuve que Ap (cone(F)) est inclus dans
Ap (cone(F)), qui est lui-même égal à Ap (F) d’après le lemme 5.4. On a donc Ap (F) = M
ce qui implique la commandabilité par le corollaire 5.2.
La preuve de la commandabilité forte est laissée en exercice.

Remarque. Le résultat n’est bien entendu plus valable si il y a des restrictions sur les
commandes. Considérons par exemple dans R le système commandé

ẋ = x + u, u ≥ −1.

Dès que x > 1, on a ẋ > 0 et donc x ne fait qu’augmenter. Ce système n’est pas
commandable.

Exemple (Sphère sur un plateau en rotation). Reprenons l’exemple de la sphère mais sup-
posons maintenant qu’elle est posée sur un plateau en mouvement. Prenons par exemple
le cas où le plateau est en rotation autour d’un axe vertical avec une vitesse angulaire ω.
En reprenant les équations du mouvement du § 5.3.2, on obtient

ṗ(t) = ωA3 p(t).

Or A3 p(t) = A3 c(t) puisque p = c − e3 et A3 e3 = 0. Ainsi,

ċ(t) = ωA3 c(t) + Ω(t)e3 .

On pose X0 (c, R) = ωA3 c ⊕ 0 et on obtient le système commandé suivant, décrivant une

sphère roulant sans glisser ni pivoter sur un plateau en rotation :

q̇ = X0 (q) + u1 X1 (q) + u2 X2 (q), q ∈ E2 × SO(3), u ∈ R2 .

Comme on a déjà montré que X1 et X2 engendrent une algèbre de Lie de dimension 5

partout, il résulte du lemme 5.9 que le système est commandable, et ce indépendamment
de la vitesse de rotation ω du plateau.
74 Chap. 5 – Systèmes commandés affines*

5.4.2 Champs récurrents

D’après la section précédente, la commandabilité dépend des tailles comparées de la
dérive et des champs commandés (penser à un bateau avançant dans le courant). Il arrive
cependant que, malgré une “forte dérive”, on arrive à commander le système.
L’idée est la suivante : pour obtenir la commandabilité il faut être capable de “remon-
ter le flot” de la dérive, c’est-à-dire d’atteindre les points situés sur le flot de −X0 . Cela
peut se faire localement quand les commandes ne sont pas bornées (section précédente),
mais également par un détour. L’exemple le plus simple est le cas où le flot de la dérive
est périodique, car on a alors ϕ−X t
0
(p) = ϕX
T −t (p).
0

Pour mettre en œuvre cette idée, nous allons introduire quelques définitions.

Définition. Soit X un champ de vecteurs complet sur M . Un point p ∈ M est dit

récurrent pour X si, pour tout t > 0 et tout voisinage V (p) de p, il existe q ∈ V (p) et
t′ > t tels que ϕX
t′ (q) ∈ Vp .

Autrement dit, toutes les trajectoires de X ne peuvent quitter le voisinage d’un

point récurrent, certaines doivent revenir dans le voisinage pour des temps arbitraire-
ment grands. Bien entendu, si ϕX t (p) est périodique, p est récurrent pour X.

Définition. On dit qu’un champ de vecteurs complet X est récurrent si tous les points
de M sont récurrents pour X.

Bien que cette propriété paraisse restrictive, il y a en fait beaucoup de champs de

vecteurs récurrents dans les applications, en particulier en mécanique 1 .

Théorème 5.10. Soit un système commandé affine sur une variété connexe M . Suppo-
sons que les conditions suivantes sont vérifiées
1. l’enveloppe convexe de l’ensemble de commande U contient un voisinage de l’ori-
gine dans Rm ;
2. dim Liep (X0 , . . . , Xm ) = dim M pour tout p ∈ M ;
3. la dérive X0 est un champ récurrent.
Alors le système est commandable.

Remarque. Quitte à faire une translation de l’ensemble de commande, on peut remplacer

l’hypothèse que X0 est récurrent par la propriété suivante : il existe une commande ū ∈ U
appartenant
P à l’intérieur de l’enveloppe convexe de U telle que le champ correspondant
X0 + i ūi Xi est récurrent.

1. D’après un théorème de Poincaré, sur une variété orientée de volume borné, tous les champs de
vecteurs dont le flot préserve la forme volume sont récurrents (voir par exemple [1]).
5.4 Systèmes affines commandables 75

P
Preuve. Notons F = F(U ) = {X0 + i ui Xi , u ∈ U }. Nous allons d’abord montrer que

Ap (F ∪ −X0 ) = Ap (F) (5.4)

Pour cela, il suffit de montrer que, ∀p ∈ M et t > 0, ϕ−X t

0
(p) est la limite d’une suite de
points atteignables à partir de p.
Choisissons un ouvert V ⊂ intAp arbitrairement proche de p. L’ouvert ϕt−X0 (V ) est
donc arbitrairement proche de ϕ−X t
0
(p). Comme X0 est récurrent, il existe t′ > t et
−X0
q ′ ∈ ϕ−X X0 ′
t (V ) tels que ϕt′ (q ) ∈ ϕt (V ).
−X
Or q = ϕt (q), où q ∈ V ⊂ intAp , donc ϕX
′ 0 0 ′ X0
t′ (q ) = ϕt′ −t (q) appartient à Ap (rappe-
lons que t′ − t > 0). Ainsi, dans tout voisinage de ϕ−X t
0
(p) il y a un point atteignable,
d’où la relation (5.4).
En appliquant maintenant le lemme 5.4 à F ∪ −X0 , (5.4) devient


Ap (F) = Ap cone(F ∪ −X0 ) .

Or le cône cone(F∪−X0 ) contient toutes les combinaisons linéaires des champs X0 , . . . , Xm .

Comme Vect{X0 , . . . , Xm } est une famille symétrique qui satisfait la condition du rang,
on a


Ap cone(F ∪ −X0 ) = M
et on conclut grâce au corollaire 5.2.

Exemple (Voiture de Dubins). La voiture de Dubins est une voiture ne pouvant aller qu’en
marche avant, à une vitesse fixée v, et dont le rayon de braquage est borné par ρ > 0.
Ses mouvements sont donc gouvernés par le système affine

q̇ = X1 + uX2 , q ∈ R2 × S 1 , |u| ≤ v/ρ,

où
∂ ∂ ∂
X1 (q) = v cos θ 1
+ v sin θ 2 , X2 (q) = .
∂x ∂x ∂θ
On a déjà vu au Chapitre 4 que Lieq (X1 , X2 ) est de dimension 3 en tout point. De plus
l’ensemble de commande U = [−v/ρ, v/ρ] contient un voisinage de l’origine. Enfin, pour
tout ū ∈]0, v/ρ[, le champ de vecteurs

Xū = X1 + ūX2 ∈ F(U )

est récurrent car toutes ses trajectoires sont des cercles de rayon v/ū > ρ.
En appliquant la remarque suivant le théorème 5.10, on obtient que la voiture de
Dubins est commandable.
76 Chap. 5 – Systèmes commandés affines*

5.4.3 Satellite à un rotor

On considère ici un satellite dont l’orbite est fixée et nous cherchons à commander
son orientation. Le système dynamique associé est donc celui d’un solide dont le centre
de gravité est fixé.
L’état du système est déterminé par sa position, une matrice de SO(3), et par sa
vitesse, un élément de T SO(3). Seule la vitesse nous intéresse ici.
Rappelons que, si on note so(3) l’ensemble des matrices symétriques 3 × 3, on a

TR SO(3) = {RΩ, Ω ∈ so(3)} et TI SO(3) = so(3).

D’autre part, à toute matrice Ω ∈ so(3) on peut associer le vecteur de R3 composé de ses
coordonnées dans la base (A1 , A2 , A3 ) :
   1
0 −ω 3 ω 2 ω
Ω =  ω3 0 −ω 1  ↔ ω = ω 2  .
−ω 2 ω 1 0 ω3

L’action de l’opérateur Ω sur un vecteur de x ∈ R3 s’écrit alors

Ωx = ω ∧ x.

Le vecteur ω associé à la vitesse Ω = R−1 Ṙ du satellite est appelé vitesse angulaire (le
mouvement correspond alors à une rotation autour de l’axe ω à une vitesse angulaire
∥ω∥).
En terme de vitesse angulaire, les équations du satellite sont données par les équations
d’Euler : 
 I1 ω̇ 1 = (I2 − I3 )ω 2 ω 3
I2 ω̇ 2 = (I3 − I1 )ω 3 ω 1
I3 ω̇ 3 = (I1 − I2 )ω 1 ω 2


où I1 , I2 , I3 sont les moments principaux d’inertie du satellite. Cette équation s’écrit de
façon plus compacte en fonction du moment angulaire µ = (µ1 , µ2 , µ3 ), avec µi = Ii ω i :

µ̇ = µ ∧ ω.

Rajoutons maintenant une action extérieure, exercée au moyen d’un moment Q.

L’équation sur le moment angulaire devient :

µ̇ = µ ∧ ω + Q.

Dans le cas du satellite, l’action se fait par le biais de 3 paires de réacteurs situés à
l’extérieur du satellite. Chaque paire exerce un moment d’amplitude constante que l’on
peut diriger dans un sens ou dans l’autre ou bien annuler (commandes du type bang-bang).
Autrement dit, pour i = 1, 2, 3,

Qi = ui bi , bi ∈ R3 , ui = −1, 0 ou 1.
5.4 Systèmes affines commandables 77

Le système commandé prend finalement la forme

µ̇ = X0 (µ) + u1 X1 (µ) + u2 X2 (µ) + u3 X3 (µ), µ ∈ R3 , u ∈ U,
où  
a1 µ 2 µ 3
X0 (µ) = a2 µ1 µ3  , Xi (µ) = bi , U = {−1, 0, 1}3
a3 µ 1 µ 2
I2 −I3 I3 −I1 I1 −I2
où on a posé a1 = I2 I3
, a2 = I3 I1
, a3 = I1 I2
.
Cas trivial : si le satellite est une sphère de densité uniforme, les constantes ai sont
toutes nulles. Le système est alors affine sans dérive, et tous les champs de vecteurs sont
constants. Il est commandable si et seulement si b1 , b2 , b3 sont linéairement indépendants.
Il faut donc au minimum 3 paires de réacteurs pour commander le satellite.
Dans le cas général, on cherche quel est le nombre minimal de paires de réacteurs
nécessaires pour rendre le système commandable.
Remarquons tout d’abord que la dérive est un champ de vecteurs récurrent de R3 . En
effet, X0 a deux intégrales de mouvement :
— la norme du moment angulaire :
(µ1 )2 + (µ2 )2 + (µ3 )2 = cte;
— l’énergie cinétique
1  (µ1 )2 (µ2 )2 (µ3 )2 
+ + = cte.
2 I1 I2 I3
Une trajectoire µ(t) de X0 est donc contenue dans l’intersection de la sphère ∥µ(t)∥2 = cte
et de l’ellipsoı̈de d’énergie constante. En traçant les trajectoires pour une énergie donnée,
on constate que par tout ouvert passe une trajectoire périodique, ce qui implique que X0
est un champ récurrent.
Pour montrer la commandabilité, il reste donc à montrer que l’algèbre de Lie engendrée
par X0 et les champs commandés est de dimension 3.
Nous allons montrer qu’une seule paire de réacteurs suffit. On s’intéresse donc au
système :
µ̇ = X0 (µ) + uX1 , X1 = b.
La preuve se déroule en deux étapes.


1er point Si le plan Vect b, X0 (b) est de dimension 2 et non invariant par X0 , alors
Lieµ (X0 , X1 ) est de dimension 3 pour tout µ.
En effet, puisque le plan est non invariant, il existe un vecteur c = αb + βX0 (b)
tel que X0 (c), b et X0 (b) sont linéairement indépendants. Or b = X1 et X0 (b) =
1
2
[X1 , [X1 , X0 ]]. Le vecteur c est donc un champ de vecteurs (constant) X2 appar-
tenant à Lie(X0 , X1 ), et X0 (c) = 12 [X2 , [X2 , X0 ]] aussi. Les 3 champs de vecteurs
constants b, X0 (b) et X0 (c) sont ainsi linéairement indépendants et appartiennent
à Lie(X0 , X1 ).
78 Chap. 5 – Systèmes commandés affines*

2ème point Si au moins 2 coordonnées de a = (a1 , a2 , a3 ) sont non nulles, il existe

b tel que Lie(X0 , X1 ) soit de dimension 3 partout.
En effet, si c = αb + βX0 (b), on a X0 (c) = α2 X0 (b) + β 2 (a1 a2 a3 b1 b2 b3 )b + αβd où
d est le vecteur  
a1 (a2 b23 b1 + a3 b1 b22 )
d =  a2 (a3 b21 b2 + a1 b2 b23 )  .
a3 (a1 b22 b3 + a2 b21 b3 )
La dimension de Lie(X0 , X1 ) est égale à 3 partout si b, X0 (b) et d sont des vecteurs
linéairement indépendants, c’est-à-dire si


∆ = det b, X0 (b), d ̸= 0.

Il suffit alors de calculer ∆ et de vérifier que si au moins deux ai sont non nuls, il
existe b tel que ∆ ̸= 0.
Finalement, si le satellite n’est pas une sphère, les moments d’inertie ne sont pas tous
égaux, donc au moins deux ai sont non nuls. Les deux points précédents montrent alors
que la vitesse du satellite est commandable avec une seule paire de réacteurs.
Annexe A

Exercices

A.1 Champs de coordonnées

Soient M une variété de dimension n, p un point de M et X1 , . . . , Xn des champs de
vecteurs sur M dont les valeurs en p sont linéairement indépendantes. On note ϕit le flot
de Xi , pour i = 1, . . . , n.
Sur un voisinage V0 de 0 dans Rn , on définit l’application

ψ: V0 ⊂ Rn → M
x = (x1 , . . . , xn ) 7→ ϕ1x1 ◦ · · · ◦ ϕnxn (p)

Le but de ce problème est d’étudier les relations entre deux familles de champs de
vecteurs sur M : les champs Xi et les champs de coordonnées ∂x∂ i .
1. Rappeler pourquoi ψ −1 définit des coordonnées (c’est-à-dire une carte) sur un
voisinage U ⊂ M de p. On notera V = ψ −1 (U ) dans la suite.
∂
2. Rappeler pourquoi les champ de vecteurs ∂xi
définis sur U vérifient :

∂ ∂ψ −1

(q) = ψ (q) , ∀q ∈ U.
∂xi ∂xi

3. Montrer que, pour tout x ∈ V ,

∂ψ

1
(x) = X1 ψ(x) .
∂x
En déduire le résultat suivant :
Redressement d’un champ de vecteurs. Soit X un champ de vecteurs sur
M . Si X(p) ̸= 0, il existe des coordonnées locales (x1 , . . . , xn ) autour de p telles
que
∂
X= .
∂x1
80 Chap. A – Exercices

4. Montrer que, si [X1 , X2 ] = 0, alors, pour tout x ∈ V ,

∂ψ

(x) = X 2 ψ(x) .
∂x2

5. En utilisant la question 4, montrer le résultat suivant :

Redressement d’une famille de champs de vecteurs.
Soient X1 , . . . , Xk des champs de vecteurs sur la variété M dont les valeurs en p
sont linéairement indépendantes (k ≤ n). Alors les deux propriétés suivantes sont
équivalentes :
— [Xi , Xj ] = 0 pour tout i, j ∈ {1, . . . , k} sur un voisinage de p ;
— il existe des coordonnées locales (x1 , . . . , xn ) autour de p pour lesquelles les
champs X1 , . . . , Xk sont des champs de coordonnées (c’est-à-dire Xi = ∂x∂ i
pour i = 1, . . . , k).

On considère maintenant la famille de champs de vecteurs F = {X1 , . . . , Xm }, m ≤ n,

et on note Orb(q) l’orbite de F passant par un point q ∈ M . On suppose également que,
∀q ∈ M , Lieq (F) est de dimension constante k ≥ m et que

Lieq (F) = Vect{X1 (q), . . . , Xk (q)}.

6. Montrer que, pour tout x ∈ V ,

∂ψ

(x) = ϕ∗ X i ψ(x) ,
∂xi
où ϕ est un difféomorphisme de M .
7. Montrer que, pour i = 1, . . . , k et x ∈ V ,

∂ψ 


(x) ∈ Vect X 1 ψ(x) , . . . , X k ψ(x) .
∂xi

8. Montrer que, pour q ∈ U ,

 ∂ ∂ 
Vect (q), . . . , (q) = Vect X1 (q), . . . , Xk (q) .
∂x1 ∂xk

9. Montrer que, dans les coordonnées locales (x1 , . . . , xn ), l’orbite Orb(p) a pour
équation
{xk+1 = · · · = xn = 0}.

10. Soit ∆ la distribution sur M définie par ∆q = Vect{X1 (q), . . . , Xm (q)}. Quelle est
l’équation de cette distribution en coordonnées locales ?
A.2 Véhicule articulé 81

A.2 Véhicule articulé

On considère un véhicule formé de 2 essieux reliés l’un à l’autre par une barre articulée.
On note d (resp. ℓ) la distance entre l’articulation et le milieu de l’essieu avant (resp.
arrière). Le véhicule est une “traction-avant”, c’est-à-dire que l’on commande la vitesse
tangentielle v0 et la vitesse angulaire ω0 de l’essieu avant. La vitesse et le rayon de
braquage de l’essieu avant sont bornés, ce qui s’écrit
ρ|ω0 | ≤ |v0 | ≤ vmax .
Le véhicule se déplace dans un plan. Ses configurations sont donc déterminées par :
— (x, y) les coordonnées planaires du centre de l’essieu avant ;
— θ l’angle entre l’axe des x et la direction de l’essieu avant ;
— φ l’angle entre les directions des essieux.
L’équation du mouvement pour ce véhicule est donnée par


 ẋ = v0 cos θ
 ẏ = v sin θ
0

 θ̇ = ω0
φ̇ = −v0 sin φ − ω0 (1 + dℓ cos φ)


Le but du problème est de déterminer l’ensemble des configurations atteignables à

partir d’une configuration donnée en fonction des paramètres d et ℓ.
1. Déterminer l’espace des configurations M du problème et montrer que c’est une
variété.
∂ ∂ ∂ ∂
2. Donner un atlas de M . Montrer que les champs de vecteurs ∂x , ∂y , ∂θ et ∂φ sont
définis dans toutes les cartes de l’atlas.
3. Écrire l’équation du mouvement comme un système commandé dans M . Préciser
l’ensemble U de commande. Expliquer pourquoi les champs de vecteurs X1 et X2
apparaissant dans ce système sont définis globalement sur M .
4. Montrer, que pour toute configuration q ∈ M , l’ensemble Aq des configurations
atteignables à partir de q est une variété. À quelle condition le système est-il
commandable ?
5. Calculez les crochets de Lie de X1 et X2 de longueur inférieure à 3 : [X1 , X2 ],
[X1 , [X1 , X2 ]] et [X2 , [X1 , X2 ]].
6. En distinguant les cas d = ℓ et d ̸= ℓ, montrer que, en tout point q ∈ M ,
Lieq (X1 , X2 ) est égal à
Vect {X1 (q), X2 (q), [X1 , X2 ](q), [X1 , [X1 , X2 ]](q), [X2 , [X1 , X2 ]](q)} .

7. Déterminer les différents ensembles atteignables Aq selon les valeurs de d et ℓ.

8. Le système est-il commandable ? Donnez une explication “mécanique” des résultats.
82 Chap. A – Exercices

A.3 Sphère roulant sur une sphère

L’objet de ce problème est d’étudier l’ensemble des configurations atteignables par
une sphère roulant sur une autre sphère. Après avoir établi un résultat théorique en
préliminaire (Partie I), nous étudierons le cas du roulement sans glissement dans la partie
II puis celui du roulement sans glissement ni pivotement dans la partie III.

Partie I : Préliminaires
Soient N et Ne deux variétés différentiables de dimensions respectives n et ñ, n ≤ ñ, et
F :Ne → N une application différentiable et surjective. On dit que des champs de vecteurs
X e ) et X ∈ X (N ) sont reliés par F si, pour tout q ∈ N et tout p ∈ F −1 (q) ⊂ N
e ∈ X (N e,
1
on a dF ◦ X(p) = X(q).
e
Soient Xe ∈ X (Ne ) et X ∈ X (N ) deux champs reliés par F , de flots respectifs ϕet et ϕt .
1º) Montrer que la courbe t 7→ F ϕet (p) est solution de l’équation différentielle dans

N q̇ = X(q). En déduire ϕt ◦ F = F ◦ ϕet .

Considérons maintenant les systèmes commandés (Σ) et (Σ),
e

(Σ) : q̇ = Xu (q), q ∈ N, u ∈ U et (Σ) eu (p), p ∈ N

e : ṗ = X e , u ∈ U,

eu ∈ X (N
où X e ) et Xu ∈ X (N ) sont des champs paramétrés par u ∈ U . On suppose que,
pour tout u fixé, X eu et Xu sont reliés par F .
2º) Soit c̃(t), t ∈ [0, T ], une trajectoire de (Σ)
e correspondant à une loi de commande


u(t) (constante par morceaux). Montrer que F c̃(t) est une trajectoire de (Σ)
correspondant à la même loi de commande u(t).
3º) Inversement, montrer que si c(t) est une trajectoire de (Σ), il existe une trajectoire


c̃(t) de (Σ)
e correspondant à la même loi de commande et telle que c(t) = F c̃(t) .
4º) Notons Aq ⊂ N (resp. Aep ⊂ N e ) les ensembles atteignables pour le système (Σ)
e Soit q ∈ N . Montrer que, quel que soit p ∈ F −1 (q), Aq = F (Aep ). En
(resp. Σ).
déduire que, si le système (Σ)
e est commandable, alors (Σ) l’est aussi.
5º) Montrer que la réciproque de cette dernière implication est fausse (on donnera un
contre-exemple avec N e = R3 et N = R2 ). Donner l’exemple d’une propriété des
ensembles Aep qui impliquerait la non commandabilité de (Σ).

Partie II : Sphère en roulement sans glissement sur une autre

sphère
Dans l’espace euclidien E3 , on considère une sphère fixe de rayon ρ > 0 et de centre
O. Le choix d’un repère orthonormé fixe (O, e1 , e2 , e3 ) définit des coordonnées sur E3 : on
notera m (en minuscules) les coordonnées d’un point M .
1. Si ñ = n et si F est un difféomorphisme, ceci implique X = F∗ X.
e
A.3 Sphère roulant sur une sphère 83

Une sphère mobile de rayon 1 et de centre C roule sur la sphère fixe. On choisit un
repère orthonormé mobile (C, a1 , a2 , a3 ) attaché à cette sphère.
Une configuration du système est un couple (x, R) où
1
— x ∈ S 2 , x = ρ+1 c, où c désigne les coordonnées du centre C de la sphère mobile ;
— R ∈ SO(3) est la matrice de transition entre le repère fixe et le repère mobile.
On note M = S 2 ×SO(3) l’espace des configurations et on identifie Tq M , le sous-espace
tangent à M en q, avec un sous-espace vectoriel de R3 × M3 (R).
Enfin, so(3) désigne l’ensemble des matrices antisymétriques 3 × 3 et (A1 , A2 , A3 ) est
la base de so(3) définie par
     
0 0 0 0 0 1 0 −1 0
A1 = 0 0 −1 , A2 =  0 0 0 , A3 = 1 0 0 .
0 1 0 −1 0 0 0 0 0
On notera Ωu la matrice de so(3) ayant pour coordonnées u = (u1 , u2 , u3 ), c’est-à-dire
Ωu = u1 A1 + u2 A2 + u3 A3 .
6º) Montrer que les mouvements de la sphère mobile roulant sans glisser sur la sphère
fixe sont solutions du système commandé
(ẋ, Ṙ) = Xu (x, R), u ∈ R3 , (x, R) ∈ M, (A.1)
1


avec Xu (x, R) = αΩu x, Ωu R et α = ρ+1 ∈ ]0, 1[.
f = SO(3) × SO(3), on considère le système commandé
Sur la variété M

(Ṫ , Ṙ) = Xeu (T, R), u ∈ R3 , (T, R) ∈ M

f, (A.2)


avec X
eu (T, R) = αΩu T, Ωu R .
7º) Montrer que, si le système (A.2) est commandable, alors le système (A.1) l’est
aussi.
Indication : utiliser l’application T 7→ T e3 de SO(3) dans S 2 et appliquer les
résultats de la partie I, question 4.

Intermède : Champs de vecteurs sur M

f
L’étude du système (A.2) requiert l’analyse de champs de vecteurs sur M
f dits “inva-
riants à droite”.
Notons e = (I, I) l’identité de M f. L’espace tangent à M
f en l’identité est Te M
f =
so(3) × so(3). C’est une algèbre de Lie pour le crochet de Lie
 

f = AC − CA, BD − DB .
(A, B), (C, D) Te M

Pour q = (P, Q) ∈ M f → M
f, on définit la multiplication à droite rq : M f (r pour
“right”) par
rq (T, R) = (T P, RQ).
84 Chap. A – Exercices

8º) Montrer que rq est un difféomorphisme. Expliciter rq−1 .

On dit qu’un champ de vecteurs X sur M f est invariant à droite si, pour tout q ∈ Mf,
rq ∗ X = X.
9º) Montrer que si deux champs X, Y ∈ X (M f) sont invariants à droite, leur crochet
de Lie [X, Y ] l’est aussi.
10º) Montrer qu’un champ invariant à droite est déterminé de façon unique par sa
valeur en e. En déduire que l’application φ : X 7→ X(e) est un isomorphisme entre
l’ensemble des champs de vecteurs invariants à droite et Te Mf.
11º) Montrer à l’aide de la question précédente que, si X est un champ invariant à
droite, il existe (A, B) ∈ Te M
f tel que

X(T, R) = (AT, BR).

12º) Calculer le crochet de Lie [X, Y ] de deux champs X et Y invariants à droite et

montrer que
φ([X, Y ]) = −[φ(X), φ(Y )]Te M
f.

Considérons maintenant une famille F de champs invariants à droite.

13º) Montrer que Lie(F) est isomorphe à Lie(φ(F)), la sous-algèbre de Lie de Te M
f
engendrée par les éléments φ(X), X ∈ F.
14º) En conclure que Orbq (F) = M f si et seulement si dim Lie(φ(F)) = 6.

Retour au problème de la sphère (suite de la question 7).

15º) Montrer à l’aide de la question 14 que le système (A.2) est commandable si et
seulement si l’algèbre de Lie
Lie(ξ1 , ξ2 , ξ3 ) ⊂ Te M
f (A.3)

engendrée par les éléments ξi = (αAi , Ai ) de Te Mf est de dimension 6.

16º) Calculer les crochets de Lie [ξi , ξj ]Te M
f et donner une base de l’algèbre de Lie (A.3).
Conclure pour le problème de la sphère roulant sans glisser sur une autre sphère.

Partie III : Roulement sans glissement ni pivotement

Supposons maintenant que la sphère mobile roule sans glisser ni pivoter sur la sphère
fixe. La contrainte de non pivotement se traduit dans le système commandé (A.1) par
une condition sur la commande : u⊥x. Les mouvements de la sphère mobile roulant sans
glissement ni pivotement sur la sphère fixe sont donc solution du système commandé :
u ∈ (Rx)⊥ ⊂ R3 , (x, R) ∈ M.
(ẋ, Ṙ) = Xu (x, R), (A.4)
1


Rappelons que Xu (x, R) = αΩu x, Ωu R et que α = ρ+1 ∈ ]0, 1[.
A.3 Sphère roulant sur une sphère 85

17º) Soit F 1 : SO(3) → S 2 l’application définie par F 1 (T ) = T e3 . Montrer que, pour

tout u ∈ (Rx)⊥ , il existe v = v(u) ∈ (Re3 )⊥ tel que le champ Yev (T ) = αT Ωv sur
SO(3) et le champ Yu (x) = αΩu x sont reliés par F 1 (voir la partie A.3 pour la
définition de champs reliés).
Indication : on utilisera l’égalité ΩT v = T Ωv T −1 que l’on démontrera.
18º) Soient F : Mf → M l’application définie par F (T, S) = (T e3 , T S −1 ) et X
ev le champ
de vecteurs sur M f défini par


ev (T, S) = αT Ωv , (α − 1)SΩv .
X

Montrer que, pour tout u ∈ (Rx)⊥ , il existe v ∈ (Re3 )⊥ tel que les champs X
ev (T, S)
et Xu sont reliés par F (v = v(u) est le même que dans la question précédente).
19º) Montrer que, v étant fixé, X
ev est un champ invariant à gauche sur M f, c’est-à-dire
que, pour tout q ∈ Mf, on a lq X
∗ v = Xv , où lq désigne la multiplication à gauche
e e
(l pour “left”) :
l(P,Q) (T, R) = (P T, QR).

Nous admettrons que le résultat de la question 13 est valable également pour les champs
invariants à gauche (la preuve est la réplique exacte de celle faite pour les champs inva-
riants à droite) : si F est une famille de champs invariants à gauche, Lie(F) est isomorphe
à Lie(φ(F)), la sous-algèbre de Lie de Te M f engendrée par les éléments φ(X), X ∈ F.

20º) En déduire que toutes les orbites de la famille Fe = {Xev : v ∈ (Re3 )⊥ } sont des
variétés de même dimension k, où k est la dimension de l’algèbre de Lie

L = Lie(ζ1 , ζ2 ) ⊂ Te M
f

engendrée par les éléments ζi = αAi , (α − 1)Ai de Te M
f.

Cas ρ ̸= 1 : les deux sphères ont un rayon différent.

21º) Montrer que k = 6. En déduire que le système (A.4) est commandable.

Cas ρ = 1 : les deux sphères ont même rayon.

22º) Montrer que k = 3. En déduire que le système (A.4) n’est pas commandable
(utiliser les questions 4 et 5).
23º) Déterminer Orbe (F)
e et montrer qu’elle est difféomorphe à SO(3).
24º) Donner la forme de l’ensemble atteignable A(x,R) du système (A.4). Interprétez.
86 Chap. A – Exercices

A.4 Transfert de population par laser

pour des systèmes quantiques
Cet exercice est tiré de l’article “Optimal Control in laser-induced population
transfer for two- and three-level quantum systems”, de Ugo Boscain, Gre-
goire Charlot, Jean-Paul Gauthier, Stephane Guerin et Hans–Rudolf Jaus-
lin, Journal of Mathematical Physics, n. 43, pp. 2107–2132 (2002). Cet ar-
ticle est disponible sur la page personnelle de l’un des auteurs, à l’adresse
http ://[Link]/ boscain/[Link].

On s’intéresse à un système quantique à 3 niveaux d’énergie E1 , E2 , E3 . La dyna-

mique du système est gouvernée par l’équation de Schrödinger suivante (dans un système
d’unités où ℏ = 1) :
dψ
i (t) = Hψ(t), t ∈ R, ψ ∈ C3 ,
dt
où H est le hamiltonien du système. Le vecteur d’état ψ(t) = (ψ1 (t), ψ2 (t), ψ3 (t)) ∈ C3
satisfait de plus |ψ1 (t)|2 + |ψ2 (t)|2 + |ψ3 (t)|2 = 1, |ψi (t)|2 étant la probabilité de mesurer
l’énergie Ei à l’instant t.
Ce système quantique est commandé au moyen de deux lasers, c’est-à-dire au moyen de
deux champs monochromatiques de fréquence ω1 = E2 − E1 et ω2 = E3 − E2 commençant
à l’instant t0 et finissant à l’instant t1 . Le hamiltonien du système est alors
 
E1 u1 (t)eiω1 t 0
H = u1 (t)e−iω1 t E2 u2 (t)eiω2 t  ,
0 u2 (t)e−iω2 t E3

où les lois de commande u1 (·) : R → R et u2 (·) : R → R sont des fonctions continues
par morceaux (u1 et u2 dépendent en fait des enveloppes des champs monochromatiques
et des coefficients de couplage entre les niveaux d’énergie dans le système quantique).
Les commandes u1 (t) et u2 (t) sont nulles pour t < t0 et t > t1 . En particulier, pour ces
valeurs de t, les probabilités |ψi (t)|2 sont constantes.
Le problème du transfert de population est d’amener le système d’un niveau d’énergie
donné à un autre. Nous le formulons de la façon suivante :

Supposons |ψ1 (t)|2 = 1 pour t < t0 .

Trouver u1 (·) et u2 (·) tels que |ψ3 (t)|2 = 1 pour t > t1 .

Préliminaire : élimination de la dérive

La dynamique que l’on vient de décrire apparaı̂t comme un système commandé affine
dépendant du temps. Nous allons commencer par le rendre autonome et sans dérive.
A.4 Transfert de population pour systèmes quantiques 87

1. Vérifier que la fonction φ(t) = (eiE1 t ψ1 (t), eiE2 t ψ2 (t), eiE3 t ψ3 (t)) est solution de

dφ
i (t) = H̄φ(t), t ∈ R,
dt
où  
0 u1 (t) 0
H̄ = u1 (t) 0 u2 (t) .
0 u2 (t) 0
Remarquons que |φi (t)|2 = |ψi (t)|2 est aussi la probabilité de mesurer l’énergie Ei
à l’instant t.
2. En posant φ1 = x1 + ix2 , φ2 = x3 + ix4 , φ3 = x5 + ix6 , réécrire le système sous la
forme
ẋ = u1 X1 (x) + u2 X2 (x), x ∈ S 5 , u ∈ R2 . (A.5)

Problème réduit dans S 5

Le problème du transfert de population prend maintenant la forme suivante :

Partant d’une condition initiale quelconque dans le cercle

1
Sin = {x ∈ S 5 : x21 + x22 = 1},

trouver une trajectoire du système (A.5) atteignant le cercle

1
Sbut = {x ∈ S 5 : x25 + x26 = 1}.

On s’intéresse surtout à l’existence d’une telle trajectoire. On va donc chercher à montrer

1 1
que l’intersection Ax ∩ Sbut est non-vide pour tout x ∈ Sin (Ax est l’ensemble atteignable
à partir de x pour le système (A.5)). Pour cela, il faut commencer par déterminer les
ensembles Ax .

3. Déterminer Lie(X1 , X2 ).
4. Montrer que la dimension de Liex (X1 , X2 ) est donnée par
(
dim Liex (X1 , X2 ) = 3 si x ∈ S 5 \Q
,
dim Liex (X1 , X2 ) = 2 si x ∈ Q

où Q = {x ∈ S 5 : F (x) = 0}, avec

F (x) = (x4 x6 + x3 x5 , x1 x6 − x2 x5 , x1 x3 + x2 x4 )

(ne pas détailler les calculs).

88 Chap. A – Exercices

1 1
Remarquons que Sin et Sbut sont contenus dans Q. Nous allons donc maintenant nous
intéresser plus particulièrement à cet ensemble.

5. L’application F : S 5 → R3 est-elle une submersion ? Peut-on en conclure quelque

chose sur Q ?
6. Soit U1 = {x ∈ S 5 : x1 ̸= 0}. Montrer qu’il existe une submersion g1 : U1 → R2
telle que Q ∩ U1 = g1−1 (0).
7. Montrer, à l’aide de la question précédente, que Q est une sous-variété de dimension
3 de S 5 .
8. Montrer que les restrictions X1 |Q et X2 |Q de X1 et X2 à Q sont des champs de
vecteurs sur Q.
9. Montrer que, pour tout x ∈ S 5 , l’ensemble atteignable Ax est une sous-variété
immergée de S 5 et préciser sa dimension selon que x appartient ou non à Q.
Peut-on conclure sur la faisabilité du transfert de population ?

Considérons maintenant la sphère

S02 = {x ∈ S 5 : x21 + x24 + x26 = 1}.

10. Montrer que les restrictions X1 |S02 et X2 |S02 sont des champs de vecteurs sur S02 .
11. En déduire que Ax = S02 pour tout x ∈ S02 , et que le problème du transfert de
population a une solution si la condition initiale est dans S02 ∩ Sin
1
. Quels sont les
vecteurs d’états que l’on peut atteindre avec de telles conditions initiales ?

Pour α ∈ R, on définit l’application ϕα : S 5 → S 5 par

 
x1
  x2 

Rα 0 0  
 x3 
 
cos α sin α
ϕα (x) =  0 Rα 0    où Rα = .
 x4  − sin α cos α
0 0 Rα   
x5 
x6

12. Montrer que {ϕα , α ∈ R} est un groupe de difféomorphismes à un paramètre dont

on précisera le générateur infinitésimal.
13. Montrer que ϕα commute avec les flots de X1 et de X2 .
14. Pour α fixé, on pose xα = (cos α, sin α, 0, 0, 0, 0). Déduire de ce qui précède que

Axα = ϕα (Ax0 ).
1 1
15. Déterminer Ax ∩ Sbut pour x ∈ Sin . En conclure que le problème de transfert de
population a toujours une solution (en supposant |ψ1 (t)|2 = 1 pour t < t0 ).
A.5 Application du lemme de Gauss 89

Questions complémentaires

16. Supposons maintenant que, pour t < t0 , la probabilité de mesurer l’énergie E1 est
1 − ε2 , et celle de mesurer E2 est ε2 . Est-il toujours possible de transférer toute
la population sur le niveau d’énergie E3 ? Si non, à quelle condition sur φ(t0 ) un
vecteur d’état particulier peut-il être transféré sur le niveau d’énergie E3 ?
17. Proposer une loi de commande constante par morceaux permettant de réaliser le
transfert de population (dans le cas où |ψ1 (t)|2 = 1 pour t < t0 ).
18. Montrer que Q est localement difféomorphe à S 1 × S 2 (Q est ce que l’on appelle
un fibré en sphères de S 1 ).

A.5 Application du lemme de Gauss

Soit (M, g) une variété riemannienne, d la distance riemannienne associée et p un point
de M . Pour λ > 0, on note Sλ = S(p, λ) la sphère de rayon λ centrée en p, c’est-à-dire
Sλ = {q ∈ M : d(p, q) = λ}.
1. Montrez que, pour λ suffisamment petit, Sλ est une sous-variété de dimension n−1
de M (utiliser l’application expp ).
2. Soient q ∈ Sλ et γ la géodésique issue de p telle que γ(1) = q (λ est supposé
arbitrairement petit). Montrez que γ̇(1) est orthogonal à Sλ , c’est-à-dire que
∀wq ∈ Tq Sλ , gq (γ̇(1), wq ) = 0.

A.6 La sphère imaginaire

Dans Rn+1 on considère la forme bilinéaire (appelée métrique de Minkowski )
h(x, y) = x1 y 1 + · · · + xn y n − xn+1 y n+1
et l’espace
Hn = {x ∈ Rn+1 : h(x, x) = −1},
appelé espace hyperbolique ou sphère imaginaire (c’est en quelque sorte la sphère de rayon
i).
1. Montrez que Hn est une sous-variété de Rn+1 de dimension n et montrer que Tx Hn
est l’ensemble des vecteurs v tels que h(x, v) = 0.
2. Pour x ∈ Hn et u, v ∈ Tx Hn , on pose
gx (u, v) = h(u, v).
Montrez que g est une métrique riemannienne sur Hn .
90 Chap. A – Exercices

On se restreint maintenant au cas n = 2 et on ne considère que la composante connexe

supérieure de H2 , c’est-à-dire

H = {x ∈ R3 : h(x, x) = −1 et x3 > 0}.

Le but du reste de l’exercice est de trouver les géodésiques de (H, g). On va procéder exac-
tement comme on l’a fait en cours pour la sphère S n : montrer d’abord que la dérivée
covariante est une sorte de projection orthogonale de la dérivée usuelle sur l’espace tan-
gent, puis en déduire une équation des géodésiques facile à résoudre.
Fixons un point x̄ ∈ H. On peut alors écrire R3 = Tx̄ H⊕Rx̄. Dans cette décomposition,
tout w ∈ R3 s’écrit
w = v − h(x̄, w)x̄, où v ∈ Tx̄ H.
Enfin, on choisit une base orthonormée (v1 , v2 ) de Tx̄ H (muni du produit scalaire gx̄ ) et
on appelle y les coordonnées de R3 dans la base (v1 , v2 , x̄), i.e.

w = y 1 v1 + y 2 v2 + y 3 x̄.

3. Montrez que l’application

π: H → R2
w = (y 1 , y 2 , y 3 ) 7→ (y 1 , y 2 )

est un difféomorphisme local en x̄.

L’application π définit donc des coordonnées (y 1 , y 2 ) sur H.
∂ ∂
4. Calculez les champs de coordonnées ∂y 1
et ∂y 2
sur H.
5. Montrez que Γkij (x̄) = 0 ∀i, j, k, où les Γkij sont les symboles de Christoffel associés
aux coordonnées (y 1 , y 2 ).
6. Soit c : I → H une courbe dans H et Y un champ le long de c. Montrez à l’aide
des questions précédentes que, pour tout t ∈ I,
dY dY

∇ċ Y (t) = (t) + h c(t), (t) c(t).
dt dt

7. Écrivez l’équation des géodésiques de H (sans coordonnées, comme on l’a fait pour
S n ).
8. Soit (x, v) ∈ T H. Montrez que la courbe
v
γv (t) = cosh(∥v∥t)x + sinh(∥v∥t)
∥v∥

est la géodésique issue de (x, v) (on suppose v ̸= 0).

Commentaires.
A.6 La sphère imaginaire 91

— Les géodésiques sont donc les intersections de H avec des plans passant par l’ori-
gine, exactement comme le sont les grands cercles pour S 2 .
— La projection stéréographique

π : H → D2
x1 x2
x 7→ ( 1+x 3 , 1+x3 )

4
est une isométrie de (H, g) dans le disque hyperbolique (D2 , (1−∥x∥2 )2 ⟨·, ·⟩), qui est

lui-même isométrique au demi-plan de Lobatchevski. Montrez-le si vous avez le

temps et donnez l’équation des géodésiques dans D2 .
92 Chap. A – Exercices
Bibliographie

[1] A. A. Agrachev and Y. L. Sachkov. Control Theory from the Geometric Viewpoint.
Springer-Verlag, 2004.
[2] W. Boothby. An Introduction to Differentiable Manifolds and Riemannian Geometry.
Academic Press, 1986.
[3] F. Jean. Systèmes dynamiques. Stabilité et commande. Cours AO102, ENSTA.
[4] V. Jurdjevic. Geometric Control Theory. Cambridge University Press, 1997.
[5] M. Spivak. A Comprehensive Introduction to Differential Geometry. Publish or Perish,
Inc., second edition, 1979.
[6] Y. Talpaert. Geo Diff. Cépaduès-Editions, 1993.
Index

Notations Xu . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
Ap . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41 X (M ) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
C ∞ (p) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 [X, Y ] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
C ∞ (M ) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
A
cone(F) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
algèbre de Lie . . . . . . . . . . . . . . . . . 48, 56, 57
cφ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
∂ application
| . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
∂xi p
∂ différentiable . . . . . . . . . . . . . . . . . . . . . . . 5
∂xi
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
linéaire tangente . . . voir différentielle
dF . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
lue dans des cartes . . . . . . . . . . 5, 21, 24
dFp . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18 réciproque . . . . . . . . . . . . . . . . . . . . . . . . 25
D(p) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 atlas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
F . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40 de T M . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
F∆ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59 maximal . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
F ∗ ωF (p) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
F ∗ g . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18 B
F∗ X . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33 base
F φψ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 naturelle . . . . . . . . . . . . . . . . . . . . . . . . . . 20
G(F) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50 base dénombrable voir espace topologique
gφ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
C
intAp . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64 cône engendré par F . . . . . . . . . . . . . . . . . . 67
Lie(F) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48 carte . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1, 6
Liep (F) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49 adaptée . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
LX g . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46 centrée . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
LX Y . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46 compatibles . . . . . . . . . . . . . . . . . . . . . . . . 1
Orbp (F) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50 domaine de . . . . . . . . . . . . . . . . . . . . . . . . 1
ϕt (q) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34 sur T M . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
ϕ(t, q) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34 champ
rgp F . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6 d’espaces vectoriels . . . . . . . . . . . . . . . . 59
T ∗ M . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27 de coordonnées . . . . . . . . . . . . . . . . . . . . 79
Tp∗ M . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24 de coordonnées . . . . . . . . . . . . . . . . . . . . 31
T M . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27 de vecteurs . . . . . . . . . . . . . . . . . . . . 30, 79
Tp M . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14 de vecteurs complet . . . . . . . . . . . . . . . 35
Xp · g . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16 récurrent . . . . . . . . . . . . . . . . . . . . . . . . . . 74
X φ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31 commandable . . . . . . . . . . . . . . . . . . . . . . . . . 64
INDEX 95

fortement . . . . . . . . . . . . . . . . . . . . . . . . . 71 différentielle . . . . . . . . . . . . . . . . . . . . . . . 33
localement . . . . . . . . . . . . . . . . . . . . . . . . 71 espace
commande . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39 cotangent . . . . . . . . . . . . . . . . . . . . . . . . . 24
bang-bang . . . . . . . . . . . . . . . . . . . . . . . . 76 des configurations . . . . . . . . . . . . . . . . . 60
ensemble de . . . . . . . . . . . . . . . . . . . . . . . 39 des phases . . . . . . . . . . . . . . . . . . . . . . . . 60
loi de . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40 euclidien . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
condition du rang . . . . . . . . . . . . . . . . . . . . . 67 tangent . . . . . . . . . . . . . . . . . . . . . . . . 14, 17
conjugué . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35 topologique . . . . . . . . . . . . . . . . . . 1, 3, 54
coordonnées . . . . . . . . . . . . . . . . . . . . . . . . 1, 79
coordonnées . . . . . . . . . . . . . . . . . . . . . . . . . . 19 F
changement de . . . . . . . . . . . . . 20, 24, 31 famille
de dFp . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21 associée . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
pour une distribution . . . . . . . . . . . . . 61 symétrique . . . . . . . . . . . . . . . . . . . . 50, 67
sur Tp M . . . . . . . . . . . . . . . . . . . . . . . . . . 20 fibré
sur T M . . . . . . . . . . . . . . . . . . . . . . . 29, 60 cotangent . . . . . . . . . . . . . . . . . . . . . . . . . 27
courbes tangent . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
intégrales . . . . . . . . . . . . . . . . . . . . . . . . . 33 fibre . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
tangentes . . . . . . . . . . . . . . . . . . . . . . . . . 14 flot . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
covecteur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24 approximation . . . . . . . . . . . . . . . . . . . . 66
crochet de Lie . . . . . . . . . . . . . . . . . . . . . . . . . 43 commutation . . . . . . . . . . . . . . . . . . . . . . 45
1-forme . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
D
dérivée G
de Lie . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46 générateur infinitésimal . . . . . . . . . . . . . . . 36
directionnelle . . . . . . . . . . . . . . . . . . 16, 46 germe . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
dérivation . . . . . . . . . . . . . . . . . . . . . . . . . 30, 42 groupe
en un point . . . . . . . . . . . . . . . . . . . . . . . 15 de difféomorphismes . . . . . . . . . . . . . voir
ensemble des . . . . . . . . . . . . . . . . . . . . . . 15 difféomorphismes
dérive . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63 orthogonal . . . . . . . . . . . . . . . . . voir O(n)
difféomorphisme . . . . . . . . . . . . . . . . . . 1, 6, 7 spécial linéaire . . . . . . . . . . voir SLn (R)
difféomorphisme . . . . . . . . . . . . . . . . . . . . . . 19
groupe de . . . . . . . . . . . . . . . . . . . . . . . . . 35 H
local . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19 Hadamard
différentielle . . . . . . . . . . . . . . . . . . . . . . . 18, 29 lemme . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
distribution . . . . . . . . . . . . . . . . . . . . . . . 59, 80 holonome . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
homéomorphisme . . . . . . . . . . . . . . . . . . . . . . 1
E
ensemble I
atteignable . . . . . . . . . . . . . . . . . . . . 41, 64 identification canonique . . . . . . . . . . . . . . . 19
de commande . . . . . . . . voir commande immersion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
équation locale . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
d’Euler . . . . . . . . . . . . . . . . . . . . . . . . . . . 76 intégrable . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
d’Euler-Lagrange . . . . . . . . . . . . . . . . . . 60 involutive . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
96 INDEX

J so(3) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
jacobienne . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7 sous-variété
de Rn . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
L de variété . . . . . . . . . . . . . . . . . . . . . . . . . 10
lagrangien . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60 immergée . . . . . . . . . . . . . . . . . . . . . . . . . 11
Leibniz ouverte . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
règle de . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 plongée . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
Lie-déterminée . . . . . . . . . . . . . . . . . . . . . . . . 64 sous-variété
espace tangent . . . . . . . . . . . . . . . . . . . . 22
M
immergée . . . . . . . . . . . . . . . . . . . . . . . . . 50
moment angulaire . . . . . . . . . . . . . . . . . . . . . 76
sphère en roulement . . . . . . . . . . . . . . . 68, 73
N structure différentiable . . . . . . . . . . . . . . . . . 2
non-holonome . . . . . . . . . . . . . . . . . . . . . . . . . 61 structure différentiable
d’une orbite . . . . . . . . . . . . . . . . . . . . . . . 54
O de T M . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
O(n) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11 submersion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
orbite . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50 pré-image de . . . . . . . . . . . . . . . . . . . 9, 11
espace tangent . . . . . . . . . . . . . . . . 51, 56 système
affine . . . . . . . . . . . . . . . . . . . . . . . . . . 63, 72
P
commandé . . . . . . . . . . . . . . . . . . . . . . . . 39
planification des mouvements . . . . . . . . . 72
de coordonnées . . . . . voir coordonnées
plongement . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
lagrangien . . . . . . . . . . . . . . . . . . . . . . . . 60
de Whitney . . . . . . . . . . . . . . . . . . . . . . . 12
sans dérive . . . . . . . . . . . . . . . . . . . . . . . . 67
image d’un . . . . . . . . . . . . . . . . . . . . . 8, 11
point récurrent . . . . . . . . . . . . . . . . . . . . . . . . 74 T
principe de d’Alembert-Lagrange . . . . . . 60 théorème
projection d’existence et d’unicité . . . . . . . . . . . . 33
canonique . . . . . . . . . . . . . . . . . . . . . . . . . 27 d’inversion locale . . . . . . . . . . . . . . . . . . 19
stéréographique . . . . . . . . . . . . . . . . . . . . 4 de Chow . . . . . . . . . . . . . . . . . . . . . . . . . . 67
prolongement . . . . . . . . . . . . . . . . . . . . . . . . . 29 de composition . . . . . . . . . . . . . . . . . . . . 18
R de dépendance etc . . . . . . . . . . . . . . . . 34
rang . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6, 22 de Frobenius . . . . . . . . . . . . . . . . . . . . . . 59
d’une application linéaire . . . . . . . . . . . 6 de l’orbite . . . . . . . . . . . . . . . . . . . . . . . . . 51
d’une distribution . . . . . . . . . . . . . . . . . 59 topologie .............................. 53
base d’une . . . . . . . . . . . . . . . . . . . . . . 1, 53
S d’une orbite . . . . . . . . . . . . . . . . . . . . . . . 51
S n . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4, 23 de T M . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
1 induite . . . . . . . . . . . . . . . . . . . . . . . . . 8, 11
S . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3, 29
2 n
S . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4, 29 T R . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
séparé . . . . . . . . . . . voir espace topologique Tx Rn . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14, 19
SLn (R) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11 trajectoire . . . . . . . . . . . . . . . . . . . . . . . . . 40, 72
SO(3) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69 transport . . . . . . . . . . . . . . . . . . . . . . . . . . 33, 45
INDEX 97

flot du . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
trivialisation . . . . . . . . . . . . . . . . . . . . . . . . . . 28

V
variété . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2, 12
de classe C k . . . . . . . . . . . . . . . . . . . . . . . . 2
variété
à bord . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
intégrale . . . . . . . . . . . . . . . . . . . . . . . . . . 59
parallélisable . . . . . . . . . . . . . . . . . . . . . . 32
vecteur tangent . . . . . . . . . . . . . . . . 14, 17, 27
en coordonnées . . . . . . . . . . . . . . . . . . . . 19
vitesse angulaire . . . . . . . . . . . . . . . . . . . . . . 76
voiture . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
de Dubins . . . . . . . . . . . . . . . . . . . . . . . . . 75
modèle simplifié . . . . . . . . . . . 40, 41, 68