République Algérienne Démocratique et Populaire

Ministère de l’Enseignement Supérieur et de la Recherche Scienti…que

Faculté des Sciences

Département de Mathématiques
POUR L’OBTENTION DU GRADE DE LICENCE EN
MATHEMATIQUES
Option : Equations di¤érentielles
Sous le thème
Equations di¤érentielles ordinaires

————————

Présenté par : Gheziel Mama et Bendahma Amina

Encadré par M.MESSIRDI BACHIR
Année Universitaire
2012-2013

www.mathonec.com
Dédicaces:
Nous dédions ce modeste travail à :

Nos chers parents.

Nos frères et soeures.

Nos grandes familles.

Et à tous ceux qui nous ont aidées et encouragées pour …nir ce travail.

www.mathonec.com
 Remerciements
Nous tenons, à exprimer notre reconnaissance et notre profonde gratitude à Monsieur MES-
SIRDI BACHIR, Maître de conférence au Département de Mathématiques à la Faculté des
Sciences de Tlemcen qui nous a, tout au long de ces travaux, conseillées, soutenues, le plus
souvent dirigées sans jamais rechigner de son temps combien précieux et ce après nous avoir
proposées ce judicieux thème de recherche.
De même que nous sommes obligées de la disponibilité, de la compétence ainsi que de la dévo-
tion avec lesquelles nos professeurs YEBEDRI MOUSTAFA, DIB MOHAMED et BOUGUIMA
ont exercé leur métier, nous enrichissant sans cesse et nous permettant de valoriser nos e¤orts
durant ce dernier semestre.

www.mathonec.com
Table des Matières

I Introduction 7

1 Dé…nitions et préliminaires 12
1.0.1 Un outil indisponsable . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
1.0.2 <<Résolution>> des équations di¤érentielles . . . . . . . . . . . . 15

2 Problème de Cauchy général 16

2.0.3 Dé…nitions et notations . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.0.4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17

2.0.5 Existence locale de solutions : . . . . . . . . . . . . . . . . . . . . . . . . . 20

2.0.6 Unicité locale des solutions : . . . . . . . . . . . . . . . . . . . . . . . 24
2.0.7 Application 1 (du lemme de Gronwall) . . . . . . . . . . . . . . . . 27
2.1 Prolongement des solutions locales, solutions maximales : . . . . . . . . 30
2.1.1 Cas particulier des EDO linéaires . . . . . . . . . . . . . . . . . . . 31
2.1.2 Dépendance des solutions en fonction des conditions initiales: . . 33

2.1.3 Continuité en fonction des conditions intiales . . . . . . . . . . . . 33

2.1.4 Dépendance des solutions en fonction des conditions initiales . . 34
2.1.5 Di¤érentiabilité en fonction des conditions initiales . . . . . . . . 37

3 EQUATIONS LINEAIRES 42

3.1 Existence globale . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42

3.2 Résolvante . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45

www.mathonec.com
 3.2.1 Cas des équations linéaires autonomes . . . . . . . . . . . . . . . . 48

3.2.2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
3.2.3 Cas des équations linéaires en dimension …nie . . . . . . . . . . . . 49
3.3 Formule de Duhamel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51

3.4 Equations linéaires autonomes . . . . . . . . . . . . . . . . . . . . . . . . . 52

3.4.1 Dimension …nie . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52

3.4.2 Dimension in…nie . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
3.5 Equations di¤érentielles linéaires à coe¢ cients périodiques . . . . . . 63

4 EQUATIONS AUTONOMES 68

4.1 Champs de vecteurs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68

4.2 Flot . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
4.2.1 Orbites . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72

4.2.2 Redressement du ‡ot . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73

4.3 Portrait de phases . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76

4.4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
4.5 Ensemble ! limite . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79

4.5.1 Propriétés générales . . . . . . . . . . . . . . . . . . . . . . . . . . . 79

4.5.2 Equations dans le plan . . . . . . . . . . . . . . . . . . . . . . . . . 82

5 STABILITE DES SOLUTIONS STATIONNAIRES 85

5.1 Théorie de Lyapunov . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85

5.2 Points …xes hyperboliques . . . . . . . . . . . . . . . . . . . . . . . . . . . . 94

5.3 Variétés invariantes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97

5.4 Introduction aux bifurcations . . . . . . . . . 101

www.mathonec.com
5.4.1 Bifurcation col-noeud . . . . . . . . . . . . . . . . . . . . 102
5.4.2 Bifurcation de Hopf . . . . . . . . . . . . . . . . . . . . . . . . . . . 104

www.mathonec.com
 Partie I

Introduction

www.mathonec.com
 On peut diviser le monde des équations di¤érentielles (EDO) en deux:
le monde familier, qui correspond en gros aux équations linéaires, et le monde étrange.
Le monde familier:
La plus simple:
x0 = ax:

Plus généralement,
x0 = a(t)x + b(t);

ou bien
x0 = Ax

en dimension supérieure. La caractéristique principale : on exprime les solutions avec des

formules.
Le monde étrange:

Exemple 0.1 loi de la dynamique et loi de la graviation(Newton).ceci permet de modéliser le

système solaire par une EDO.Cette EDO est non linéaire: on peut résoudre le problème des deux
corps (ce qu’a fait Newton), mais pas au-delà. Exemples de solutions complexes (animation).
Hors de portée de ce cours...
requins et sardines (Volterra 1920). En l’absence d’ interractions

X 0 = ax

et
y0 = by;

le nombre de rencontres est proportionnelle à xy; on obtient

x0 = ax cxy

y0 = by + dxy

On ne peut pas résoudre, mais on sait néanmoins décrire le comportement qualitatif des

www.mathonec.com
solutions. Et déjà, dire qu’elles éxistent!

Exemple 0.2 Petites oscillations du pendule. On ne sait pas résoudre l’équation

y 00 = sin(y)

.on peut linéariser, et éspérer que l’équation linéarisée décrit le comportement des petites
oscillations, mais comment le justi…er?

La forme la plus générale d’une équation di¤érentielle ordinaire (en abrégé EDO) est

F (t; u; u0 ; :::; u(k) ) = 0:

où u est une fonction inconnue de la variable réelle t à valeurs dans Rn ou plus généralement
dans un espace de Banach X; u0 ; :::; u(k) désignent les dérrivées successives de u, et F est une
fonction donnée, supposée <<régulière>>(on précisera comment par la suite) sur I U U1
::: U k où I est un intèrvalle ouvert de R,U; U 1 ; :::; U k sont des ouverts connexes de X.On
ne s’intéressera dans ce cours qu’à des équations di¤érentielles résolues, pour lesquelles il éxiste
une fonction G, régulière sur I U U1 ::: Uk 1 telle que

F (t; u; u0 ; :::; u(k) ) = 0 , u(k) = G(t; u; u0 ; :::; u(k 1)

):

On observe de plus que

u(k) = G(t; u; u0 ; :::; u(k 1)

) , U 0 = G(t; U );

www.mathonec.com
 0 1
u
B C
B C
B u0 C
B C
B C
B : C
U := B
B
C;
C
B : C
B C
B C
B : C
@ A
u(k 1)

0 1 0 1
0 I ::: 0 0
B C B C
B C B C
B : C B 0 C
B C B C
B C B C
B : : : : : : : C B : C
G(t; U ) := B
B
CU + B
C B
C
C
B : C B : C
B C B C
B C B C
B I C B : C
@ A @ A
0 ::: 0 G(t; u; u0 ; :::; u(k 1) )

la fonction G étant évidemment aussi régulière que G. On supposera donc sans perte de
généralité k = 1.
Désormais, on considère une équation dite d’ordre 1, de la forme

du
= f (t; u)
dt

où u est une fonction inconnue de la variable réelle t à valeurs dans un espace de Banach
X, et f est une fonction donnée sur I U , ouvert connexe non vide de R X.
Lorsque f ne dépend pas de t, l’équation di¤érentielle est dite autonome.

Remarque 0.1 On peut toujours se ramener, par une astuce, à une équation non autonome.
En e¤ et, il su¢ t de considérer l’équation étendue

0 1 0 1
d @ u A @ f (t; u) A
=
ds t 1

Cette approche, parfois utile, est malgré tout arti…cielle. Il faut savoir étudier directement
certaines propriétés des équations dites à coe¢ cients variables, où f dépend vraiment de t: Ce

10

www.mathonec.com
sera le cas au moins dans les deux premiers chapitres.

11

www.mathonec.com
Chapitre 1

Dé…nitions et préliminaires

Dé…nition 1.1 Soit I R un intervalle ouvert de R et F : I Rd ! Rd une application . On

appelle solution de l’équation di¤ érentielle

y 0 = F (t; y) (1.1)

tout couple (J; y) où J I est un sous intervalle de I et y une fonction dérivable dé…ne sur
J telle que

8t 2 J; y 0 (t) = F (t; y(t)):

Remarque 1.1 _On peut tout à fait étendre ces dé…nitions au cas où F est dé…nie sur un
ouvert quelconque de R Rd :

Remarque 1 _L’équation 1:1) est appelée du premier ordre car elle ne fait intervenir que les
dérivées premières de la fonction inconnue.
_La forme la plus générale d’une équation di¤ érentielle est :

F (t; y; y 0 ) = 0;

mais nous ne traiterons pas ce cas ici. Notons que le théorème des fonctions implicites,
permet de se ramener au cas résolu dans un certain nombre de cas.

12

www.mathonec.com
Dé…nition 1.2 (solution maximale) [1]
^ y^) véri…ant
On dit que (J; y) est une solution maximale de 1:1) s’il n’éxiste pas de solution (J;
J J^ et y^jJ = y.

Dé…nition 1.3 Toute solution (J; y) de 1:1) dé…nie sur l’intervalle J = I tout entier est dite
globale.

1.0.1 Un outil indisponsable

la théorie des équations di¤érentielles utilise abondamment le lemme (ou l’inégalité) de Gron-
wall. Ce chapitre introductif est l’occasion de le rappeler.Parfois, on appelle (à tort) lemme de
Gronwall le fait suivant : si une fonction u à valeurs dans R:de classe C 1 , satisfait une inégalité
di¤érentielle du type :
u0 (t) b(t) + a(t)u(t) , pour tout t 2 [0; T ]
alors
Rt
a(s)ds
Rt Rt
a(s)ds
u(t) e 0 u(0) + 0 b( )e d ; pour tout t 2 [0; T ] :
En e¤et, l’inégalité di¤érentielle implique
Rt Rt
d a(s)ds a(s)ds
dt (e u(t)) e b(t);
0 0

et donc par intégration entre 0 et t on obtient immédiatement l’inégalité annoncée. Le

lemme de Gronwall est un peu plus subtil, puisqu’il suppose une inégalité intégrale et non
une inégalité di¤érentielle (la seconde impliquant la première mais l’iverse est faux). Or les
éstimations à priori que l’on obtient en général sont plutôt du type intégral, d’où l’intérêt de
ce lemme, dont la preuve est néanmoins élémetaire.

Lemme 1.1 (Gronwall)

si u 2 C([0; T ]; R+ ) est telle qu’il existe a et b 2 C([0; T ]; R+ ) avec

Z t
u(t) b(t) + a( )u( )d ; pourtoutt 2 [0; T ];
0

13

www.mathonec.com
 alors
Z t Rt
a(s)ds
u(t) b(t) + b( )a( )e d
0

;pour tout t 2 [0; T ]:

Démonstration:
La seule astuce consiste à majorer l’intégrale du second membre

Z t
v(t) = a( )u( )d
0

par la méthode décrite précédemment.Comme

v 0 (t) = a(t)u(t) a(t)(b(t) + v(t))

par hypothèse, on a donc

d Rt Rt
a(s)ds a(s)ds
(e 0 v(t)) e 0 a(t)b(t);
dt

d’où après intégration (notons que v(0) = 0 par dé…nition):

Z t Rt
a(s)ds
v(t) a( )b( )e d ;
0

pour tout t 2 [0; T ]:

En majorant de cette façon v(t) dans l’inégalité de départ, on obtient le résultat.
Cette version du lemme de Gronwall est donnée surtout pour mettre en évidence la méthode
de calcul. On pourrait donner une autre version de l’inégalité obtenue, en integrant par parties
si b est dérivable :

Rt Z t Rt
a(s)ds
u(t) b(0)e 0 + b0 ( )e a(s)ds
d
0

, pour tout t 2 [0; T ] :

En particulier, si b est constante (b(t) = k; 8t 2 [0; T ]), on obtient simplement:

Rt
a(s)ds
u(t) ke 0 :

14

www.mathonec.com
 En pratique, il est conseillé de refaire rapidement le calcul pour éviter les erreurs.

1.0.2 <<Résolution>> des équations di¤érentielles

Dans le chapitre II, on va s’intéresser à l’éxistence, l’unicité, et la dépendance des solutions par
rapport aux <<conditions initiales>> u(t0 ) = u0 : on (re) verra des résultats essentiellement
théoriques, car il y a trés peu d’équations di¤erentielles dont on connait explicitement les so-
lutions, en dehors des équations linéaires à coe¢ cients constants (dont le calcul des solutions
se ramène à un calcul de primitive). Le chapitre III sera consacré aux propriétés spéci…ques
des équations linéaires, en insistant sur le cas à coe¢ cients variables.A partir du chapitre IV,
on s’attaquera aux propriétés qualitatives des équations di¤érentielles :sans prétendre résoudre
ces équations, on peut en e¤et obtenir beaucoup d’informations sur le comportement de leurs
solutions.(Cette idée générale remonte à Poincré).On étudiera notamment l’éxistance et la sta-
bilité de solutions particulières, comme les solutions stationnaires et les solutions périodiques
(qui jouent un grand rôle dans les applications).

15

www.mathonec.com
Chapitre 2

Problème de Cauchy général

2.0.3 Dé…nitions et notations

On note E = RN ; N 1 et x = (x1 ; :::; xN ) un élément de E. Sa norme jjxjjE sera l’une des

normes usuelles sur RN (on rappelle que toutes les normes sont équivalentes sur E ).
Soit D un ouvert de R E et f : D ! E une fonction continue. Pour tout (t; x)2 D; on
notra
f (t; x) = (f1 (t; x); :::; fN (t; x))

où chaque fonction fi est continue de D dans R: la notation (a; b) recouvre tous les intervalles
de R de la forme [a; b]; ]a; b]; [a; b[; ]a; b[. de même , on utilisera par exemple la notation (a; b]
pour ne pas préciser la nature de l ’intervalle en a:
Dans ce chapitre, on s’intéresse aux équations di¤erentielles ordinaires (notée en abrégé
EDO)du premier ordre, sous forme normale (ou résolue):

x0 (t) = f (t; x) (2.1)

Commençons par préciser la notion de solution pour ce type d’équation:

Dé…nition 2.1 Une solution de 2:1est un couple ('; J) où J est un intervalle de R et ' =
('1 ; :::; 'N ) est une fonction dérivable sur J à valeurs dans E telle que (t; '(t)) 2 D pour tout
t 2 J et

16

www.mathonec.com
 '0i (t) = fi (t; '(t)); 8t 2 J; 8i = 1; :::; N:

On remarque tout de suite que, f et ' étant deux fonctions continues sur J et ' est de
classe C 1 surJ:(' 2 C 1 (J)):

2.0.4

Exemple 2.1 pour 0 1

x1
N = 2; x = @ A 2 R2
x2

et 0 10 1
a(t) b(t) x
f (t; x) = M (t)x = @ A@ 1 A
c(t) d(t) x2
0 10 1
a(t) b(t) x1
f (t; x) = M (t)x = @ A@ A
c(t) d(t) x2

où a(t); b(t); c(t) et d(t) sont des fonctions réelles continues, l’équation x0 (t) = f (t; x(t)) est
appelée équation linéaire du premier ordre.

Exemple 2.2 Pour

N = 1; f (t; x) = a(t)x + b(t)x

où a(t); b(t) sont des fonctions continues et 2 Rnf0; 1g;l’EDO 2:1 est une équation de
Bernoulli.

Exemple 2.3 L’équation

x0 (t) = a(t)x2 (t) + b(t)x(t) + c(t)

pour laquelle
N = 1; f (t; x) = a(t)x2 + b(t)x + c(t)

17

www.mathonec.com
où a(t); b(t) et c(t) sont trois fonctions continues, l’EDO 2:1 est une équation , de Riccati.Le
problème 2:1 peut avoir de nombreuses solutions sur un intervalle donné. Par exemple, pourN =
1; D = R R et f (t; x) 1, l’EDO

x0 (t) = 1;

admet '(t) = t + c comme solution sur R pour tout c 2 R:On introduit la notion de problème
de Cauchy:

Dé…nition 2.2 Soit (t0 ; x0 ) 2 D:Résoudre le problème de Cauchy:

x0 (t) = f (t; x);

x(t0 ) = x0

consiste à déterminer un couple ('; J) où J est un intrvalle de R contenant t0 et ' une fonction
dérivable (en fait C 1 ) de J dans E telle que (t; '(t)) 2 D pour tout t 2 J;

'0 (t) = f (t; '(t) (2.2)

pour tout t 2 J et
'(t0 ) = x0

En intégrant l’EDO du problème de Cauchy 2.2 entre t0 et t et compte tenu de la condition

x(t0 ) = x0 ; on obtient :

Z t
'(t) = x0 + f (s; '(s))ds (2.3)
t0

où

Z t Z t Z t
f (s; '(s))ds = ( f1 (s; '(s)ds; :::; fN (s; '(s))ds):
t0 t0 t0

18

www.mathonec.com
Réciproquement, toute fonction ' véri…ant( 2.3) est bien une solution C 1 de ( 2.2)
Nous utiliserons souvent l’équivalence entre les deux formulations( 2.2) et ( 2.3) dans la
suite du chapitre.
En reprenant l’exemple précédent du problème de cauchy, on voit que si l’on ajoute à
l’équation( 1.1) la condition x(0) = 0, on doit …xer la constante c et l’unique solution dé…nie
sur R est (t) = t: Nous verrons plus tard cependant que ce n’est pas toujours aussi simple et
que sans hypothèses supplémentaires, notamment sur la fonction f , l’unicité de la solution n’est
pas assurée en général pour le problème de cauchy.
Les formulations( 2.1) et( 2.2) bien que ne faisant intervenir que la dérivée première de
x(t); recouvrent en fait une large classe de problèmes. En e¤ et, il est souvent possible de mettre
sous la forme 2.1 des EDO dans lesquelles apparaissent des dérivées à un ordre quelconque.
Considérons pour simpli…er que les fonctions x(t) sont à valeurs dans R (ie:N = 1). soit D un
ouvert deR Rp avec p 1 et f : D ! R une fonction continue. En notant x(k) (t) la derivée
k–ème de x(t); toute équation di¤ érentielle ordinaire du p–ème ordre associée à f qui s’écrit :

x(p) (t) = f (t; x(t); x0 (t); :::; x(p 1)

(t)); (2.4)

peut se mettre sous la forme( 2.1) .

Notons en e¤ et
x1 (t) = x(t) et xi+1 (t) = x(i) (t)

pour i = 1; :::; p 1 et introduisons

0 1 0 1
x1 (t) x2 (t)
B C B C
B C B C
B x2 (t) C B : C
B C B C
B C B C
B : C B : C
X(t) = B
B
C 2 Rp et F (t; X) = B
C B
C
C
B : C B : C
B C B C
B C B C
B : C B xp (t) C
@ A @ A
xp (t) f (t; x1 (t); :::; xp(t) )

L’EDO 2.4 devient alors :

X 0 (t) = F (t; X(t))

19

www.mathonec.com
Le problème de Cauchy correspondant à l’équation ci-dessus s’obtient en ajoutant une condition
du type X(t0 ) = X0 où t0 2 R et X0 2 Rp : Ceci correspond à la donnée de x(t0 ); x0 (t0 ); :::; x(p 1) (t
0 ):

Nous allons introduire maintenant la notion de solution approchée. Ces solutions ne seront
en général pas aussi régulières que les solutions exactes dont nous venons de parler.

Dé…nition 2.3 On dira qu’une fonction à valeus dans E est C 1 par morceaux sur un inter-
valle J de R si:
1. est continue sur R
2.Il existe un ensemble …ni S = ft1 ; : : : ; tp g de points de J tels que soit C 1 sur J n S et
0 0
limt!t+ (t) et limt!t (t) existent mais ne coïncident pas forcément.
i i

Nous avons alors :

Dé…nition 2.4 soit " > 0 et J un intervalle de R. On dira que 2 C(J) est une solution
" approchée de 2.1 si :
1. (t; (t)) 2 D; 8t 2 J:
0
2. est C 1 par morceaux sur J (on note S les points où n’est pas dé…nie).
3. jj 0 (t) f (t; (t)jjE "; 8t 2 JnS:

2.0.5 Existence locale de solutions :

Historiquement, il a été d’abord démontré que le problème de Cauchy ( 2.2) admettait lo-
calement une solution unique. Ce résultat, que nous le démontrerons plus loin, nécessite une
hypothèse plus forte que la simple continuité sur la fonction f . En introduisant une suite de
solutions approchées, le mathématicien italien Peano en s’appuyant sur un résultat d’Ascoli a
démontré :

Théorème 2.1 soit (t0 ; x0 ) 2 D et soit a > 0 et b > 0 tels que le cylindre C = fjt t0 j a;
jjx x0 jj bg soit inclus dans D. On note

M = sup jjf (t; x)jjE

(t;x)2C

20

www.mathonec.com
et
b
= min(a; )
M

alors pour tout " > 0 il existe une solution " approchée du problème de cauchy sur ( 2.2)
l’intervalle [t0 ; t0 + ] :

Dé…nition 2.5 Un ensemble de fonctions F dé…nies sur un intervalle réel I et à valeurs dans
E;est dit uniformément équicontinu si pour tout " > 0, il existe " > 0 tel que:

jjf (t) f (T )jjE " 8f 2 F

et
8t; T 2 I; jt Tj ":

Nous aurons besoin également du:

Lemme 2.1 (Ascoli)

soit I un intervalle borné de R et C(I) l’espace vectoriel des fonctions continues sur I à
valeurs dans E, muni de la norme

jjf jj1 = sup jjf (t)jjE :

t2I

Alors tout sous ensemble F de C(I); borné et uniformément équicontinu est relativement com-
pact.

Démonstration :
Si F est un ensemble …ni, c’est évident. On peut donc supposer que F contient une suite
de fonctions deux à deux distinctes (fn )n . Soit " > 0 …xé.
L’uniforme équicontinuité de F nous assure l’existence de " > 0 tel que

"
jjf (t) f (T )jjE ; sijt Tj ":
3

21

www.mathonec.com
 L’intervalle I étant borné, on peut le recouvrir par un nombre …ni d’intervalles ouverts
I1 ; :::; Ip de longueur inférieure à " :Dans chacun de ces intervalles, on choisit des réel t1 ; :::; tp :
La suite (fn (t1 ))n est bornée dans E: Elle admet donc une sous suite convergente (fn1 (t1 ))n1 :
De la même façon, la suite ( fn1 (t2 ))n1 admet également une sous suite convergente notée
(fn2 (t2 ))n2 :En poursuivant ce procédé d’extraction, on aboutit à la construction de p suites,
toutes convergentes (fnp (ti ))np ; i = 1; :::; p: Notons fen = fnp : Pour le " choisi plus haut, il existe
donc un rang N" tel que

"
jjfem (ti ) fen (ti )jjE ; 8i = 1; :::; p; 8n; m N" :
3

Montrons que (fen )n est uniformément convergente sur I. Pour tout t 2 I, il existe ti tel que
jt ti j ": Pour n; m N" ; on a alors:

jjfem (t) fen (t)jjE jjfem (t) fen (ti )jjE + jjfem (ti ) fen (ti )jjE + jjfen (ti ) fen (t)jjE

" " "

+ + = ";
3 3 3

Ce qui prouve que (fen )n est de Cauchy uniforme.

Théorème 2.2 (Ascoli-Peano):

Soient (t0 ; x0 ) 2 D , a > 0 et b > 0 tels que le cylindre C = fjt t0 j a; jjx x0 jj bg soit
inclus dans D: On note
M = sup jjf (t; x)jjE
(t;x)2C

et
b
= min(a; );
M

alors il existe (au moins) une solution ' du problème de Cauchy 2.2 sur l’intervalle [t0 ; t0 + ]:

Démonstration:

22

www.mathonec.com
 Soit ("n )n une suite décroissante de réels positifs tendant vers 0. Pour chaque "n , et d’après
le théorème précédent, il existe une solution "n –approchée du problème 2.2, notée 'n et dé…nie
sur [t0 ; t0 + ]: Toujours d’après la démonstration du théorème précédent , chacune de ces
solutions véri…e :

jj'n (t) 'n (T )jjE M jt T j; 8t; T 2 [t0 ; t0 + ]:

Cette inégalité entraîne d’une part que l’ensemble des termes de la suite ('n )n est unifor-
mément équicontinu de C([t0 ; t0 + ]): D’autre part, en choisissant T = t0 ; on montre
que
jj'n (t)jjE jjx0 jjE + b

pour tout t 2 [t0 ; t0 + ]: C’est à dire que la suite ('n )n est uniformément bornée. Les hy-
pothèses du Lemme d’Ascoli sont donc véri…ées et il existe une suite extraite ('nk )n convergeant
vers une fonction continue sur [t0 ; t0 + ] et notée ':
Posons 8
< n (t) = '0n (t) f (t; 'n (t)) si '0n (t) existe
:
n (t) =0 sinon

Par dé…nittion des solutions " approchées,

jj n (t)jjE "n

pour tout n et pour tout t 2 [t0 ; t0 + ]: En réecrivant maintenant les solutions sous forme
intégrale, on obtient que :

Z t
'n (t) = x0 + f (s; 'n (s)) + n (s)ds:
t0

Comme 'nk converge uniformément vers ' et que f est uniformément continue sur le com-
pact C; f (t; 'nk (t)) converge uniformément vers f (t; '(t)) . Si on ajoute la convergence uniforme
de nk vers 0, on peut passser à la limite dans l’égalité ci-dessus pour obtenir :

Z t
'(t) = x0 + f (s; '(s))ds;
t0

23

www.mathonec.com
 ce qui prouve que ' est de classe C 1 sur [t0 ; t0 + ] et est solution du problème de
Cauchy 2.2.
On utilisera en général le théorème sous la forme simpli…ée suivante :

Corollaire 2.1 Soit f 2 C(D): Alors pour tout (t0 ; x0 ) 2 D; il existe un voisinage de t0 dans R
sur lequel le problème de Cauchy 2.2 admet une solution. A
de nous intéresser au problème de l’unicité des solutions, voici un exemple simple qui montre
qu’un problème de Cauchy peut avoir une in…nité de solutions même sur un voisinage arbi-
trairement petit autour de la condition de Cauchy.

Exemple 2.4 Le problème de Cauchy

p
x0 (t) = jx(t)j; x(0) = 0;

où
p
f (x; t) = jxj

est dé…nie et continue sur D = R R et (0; x(0)) 2 D; admet sur [0; 1] et pour tout
0 c 1 la solution 'c dé…nie par :

1
'c (t) = (t c)2 ; c t 1; 'c (t) = 0; 0 t < c:
4

2.0.6 Unicité locale des solutions :

On commence par un lemme technique mais qui sera très utile dans la suite :

Lemme 2.2 (Gronwall)

Soit une fonction continue dé…nie sur un intevalle [a; b] R et à valeurs dans R+ :On
suppose qu’il existe t0 2 [a; b] et trois constantes positives A; B et C telles que :

Z t
(t) A + Bj (s)dsj + Cjt t0 j; 8t 2 [a; b]:
t0

Alors, pour tout t 2 [a; b] on a l’estimation :

24

www.mathonec.com
 C Bjt
(t) AeBjt t0 j
+ (e t0 j
1):
B

Démonstration :
Si t = t0 ; la conclusion est évidente. Intéressons nous dans un premier temps aux valeurs
de t telles que t0 < t < b et posons :

Z r
(r) = (s)ds; 8r 2 [t0 ; b]:
t0

La continuité de entraîne que 2 C 1 ([t0 ; b]) et

Z r
Br 0 0 Br Br
( (r)e ) =( (r) B (r))e = ( (r) B (s)ds)e :
t0

On en déduit, d’après l’inégalité véri…ée par , que :

Br 0 Br
( (r)e ) e (A + C(r t0 )); 8r 2 [t0 1):

puis, en intégrant cette relation entre t0 et t, on obtient :

Z t
Bt A Bt0 Bt Br
(t)e (e e )+C e (r t0 )dr
B t0

:
Une intégration par parties nous permet de calculer le dernier terme :

Z t
Br 1 B(t t0 ) 1
e (r t0 )dr = e (t t0 );
t0 B B

ce qui nous permet d’estimer la fonction comme suit :

Z t
A B(t t0 ) C 1 B(t t0 ) 1
j (s)dsj = (t) (e 1) + ( e (t t0 )):
t0 B B B B

En combinant cette relation avec l’inégalité donnée dans les hypothèses, on obtient la con-
clusion du lemme dans le cas t0 < t b: Pour le cas a t < t0 , on procède de la même façon
mais en posant

25

www.mathonec.com
 Z t0
(r) = (s)ds:
r

La bonne propriété qui va assurer l’unicité pour le problème de Cauchy 2.2 est le caractère
lipschitzien de la fonction f . Précisons cette notion :

Dé…nition 2.6 On dira que f est lipschitzienne en x (uniformément par rapport à t), et on
notera f 2 Lip(D), s’il existe k > 0 tel que :

jjf (t; x1 ) f (t; x2 )jjE kjjx1 x2 jjE ; 8(t; x1 ); (t; x2 ) 2 D:

Remarquons que cette notion n’entaîne pas que f est continue sur D comme le prouve
l’exemple suivant :
D = R2 ;
f (t; x) = 1

si t > 0 et
f (t; x) = 0

si t 0:
En revanche, si f est lipschitzienne au sens classique, c’est à dire s’il existe k > 0 tel que :

jjf (t1 ; x1 ) f (t2 ; x2 )jjE k(jt1 t2 j + jjx1 x2 jjE ); 8(t1 ; x1 )2; (t1 ; x2 ) 2 D;

alors f est en particulier lipschitzienne en x uniformément en t.

L’exercice suivant fournit un exemple simple de fonction lipschitzienne.

Exercice 1 :Montrer que si f 2 C 1 (D) ou si Dx f est continue en (t; x), alors f est lipschitzi-
enne en x uniformément par rapport à t sur tout compact convexe K inclus dans D.
Une application importante du lemme de Gronwall concerne les solutions " approchées :

26

www.mathonec.com
2.0.7 Application 1 (du lemme de Gronwall)

Soit f 2 Lip(D) \ C(D) avec pour constante de Lipschitz k > 0: soient '1 ; '2 deux solutions
respectivement "1 et "2 approchées de 2.1 sur un même intervalle (a; b)
et telles que, pour un certain a < t0 < b on ait :

"1 + "2 kjt

jj'1 (t) '2 (t)jjE ekjt t0 j
+ (e t0 j
1): (2.5)
k

Démonstration :
En reprenant la démonstration du théorème 2.2, on a l’écriture des solutions " approchées
sous forme intégrale :

Z t
'i (t) = 'i (t0 ) + f ('i (s); s) + i (s)ds; i = 1; 2;
t0

avec
jj i (s)jjE "i ; i = 1; 2:

On pose
(t) = jj'1 (t) '2 (t)jjE

et on a

Z t
(t) jj'1 (t0 ) '2 (t0 )jjE + jjf ('1 (s); s) f ('2 (s); s)jjE ds + ("1 + "2 )jt t0 j:
t0

Or, f étant Lipschitzienne en x uniformément en t :

jjf ('1 (s); s) f ('2 (s); s)jjE kjj'1 (s) '2 (s)jjE = k (s):

On applique alors l’inégalité du lemme de Gronwall pour obtenir 2.5:

Exercice 2 Soient ' et deux fonctions continues de

Exemple 2.5 Exercice 3 On suppose qu’il existe des constantes positives A; B telles que [a; b]
dans R+ et t0 2 [a; b]:

27

www.mathonec.com
 Z t
'(t) A + Bj (s)'(s)dsj; 8t 2 [a; b]:
t0

Montrer qu’alors:

Z t
'(t) A exp( (s)dsj):
t0

Théorème 2.3 (Cauchy-Lipschitz) :

Soit f 2 C(D) \ Lip(D) et avec les mêmes notations que pour le théorème 2.2, il existe une
unique solution du problème de Cauchy 2.2 sur l’intervalle [t0 ; t0 + ]:

Démonstration :
Le théorème 2.2 nous assure l’existence d’au moins une solution '1 sur l’intervalle considéré.
On note '2 une éventuelle autre solution et on introduit la fonction continue

(t) = jj'1 (t) '2 (t)jjE ;

dé…nie sur [t0 ; t0 + ]: Les fonctions '1 et '2 étant des solutions du problème de Cauchy,
elles s’écrivent, sous la forme intégrale 2.3 :

Z t
'i (t) = x0 + f (s; 'i (s))ds; 8t 2 [t0 ; t0 + ]; i = 1; 2
t0

:
On obtient en particulier que :

Z t Z t
(t) = jj f (s; '1 (s)) f (s; '2 (s))dsjjE jjf (s; '1 (s)) f (s; '2 (s))jjE ds;
t0 t0

puis, f étant lipschitzienne en x uniformément en t sur D; il existe k > 0 telle que

jjf (s; '(s)) f (s; '2 (s))jjE kjj'1 (s) '2 (s)jjE

28

www.mathonec.com
pour tout s 2 [t0 ; t0 + ]; ce qui nous donne :

Z t
(t) kj (s)jds:
t0

On conclut ensuite en appliquant le lemme de Gronwall avec A = 0; B = k et C = 0:

Introduisons une nouvelle dé…nition:

Dé…nition 2.7 On dira que f est localement lipschitzienne sur D si pour tout (t; x) 2 D il
existe une boule B = f(t0 ; x0 ) 2 D; jjx x0 jjE < "; jt t0 j < "g D et une constante k > 0 telles
que f soit lipschitzienne sur B: On note alors f 2 Liploc (D):

Remarque 2.1 Dans cette dé…nition la constante de Lipschitz n’est valable que localement.

Exercice 4 Montrer que si f 2 C 1 (D) alors f 2 C(D) \ Liploc (D):Du théorème précédent, on
déduit que l’unicité de la solution est en fait globale, ce qui peut se traduire par : lorsque la
fonction f 2 Liploc (D) alors les graphes de deux solutions distincts de 2.1 ne peuvent se croiser.

Théorème 2.4 Soit f 2 C(D) \ Liploc (D) et soient ('1 ; J1 ) et ('2 ; J2 ) deux solutions de 2.1
telles que J1 \ J2 6= ?: S’il existe un point t0 de J1 \ J2 tel que '1 (t0 ) = '2 (t2 ) alors '1 '2
sur J1 \ J2 :

Démonstration:
Les ensembles J1 et J2 sont des intervalles. Il en est donc de même de J = J1 \ J2 qui est en
particulier connexe et non vide puisqu’il contient t0 : On note I = ft 2 J tels que '1 (t) = '2 (t)g:
Montrons que I est ouvert et fermé dans J ce qui entraînera que I = J. Les fonctions '1 et
'2 étant continues, on en déduit que I est fermé. Soit t1 2 I; notons x1 = '1 (t) = '2 (t): Alors
(t1 ; x1 ) 2 D et selon le théorème 1.3, le problème de Cauchy :

x0 (t) = f (t; x(t)); x(t1 ) = x1 ;

29

www.mathonec.com
 admet une unique solution sur [t0 ; t0 + ]: On en déduit que '1 '2 sur cet intervalle
et que ]t0 ; t0 + [ I et donc que I est ouvert.
On considèrera à partir de maintenant que l’on a toujours f 2 Liploc (D) \ C(D) ou plus
simplement f 2 Lip(D) \ C(D); c’est à dire qu’il existe toujours une solution unique pour le
problème de Cauchy ( 2.2)

2.1 Prolongement des solutions locales, solutions maximales :

commençons par poser quelques dé…nitions :

Dé…nition 2.8 Soit ('1 ; J1 ) et ('2 ; J2 ) deux solutions de 2.1. On dit que ('2 ; J2 ) prolonge
('1 ; J1 ) si J1 J2 et '1 '2 sur J1 :

Dé…nition 2 Dé…nition 2.9 Une solution ('; J ) de 2.1 est dite maximale si elle n’admet
aucun prolongement.

Théorème 2.5 (Existence d’une solution maximale)

Soit f 2 C(D) \ Liploc (D): Alors par tout point (t0 ; x0 ) 2 D il passe une unique solution
maximale au problème de Cauchy (2.2).

Démonstration:
Considérons l’ensemble S de tous les couples ('; J) de solutions au problème de Cauchy
2.2. Si ('1 ; J1 ) et ('2 ; J2 ) sont deux tels couples alors J1 \ J2 n’est pas vide car il contient t0
et '1 '2 sur J1 \ J2 d’après le théorème 1.4.
Soit I la réunion de tous les intervalles J. Sur I on peut donc dé…nir la fonction par
' sur J pour tout ('; J) 2 S. Cette fonction est la solution maximale cherchée.
La question à laquelle nous allons répondre maintenant est :pourquoi une solution maximale,
dé…nie sur un intervalle borné, ne peut-elle être prolongée sur un intervalle plus grand?

Théorème 2.6 On suppose que est un ouvert de RN et que D =]a; b[ : Soient f 2 C(D) \
Liploc (D) et (t0 ; x0 ) 2 D: Si ('; (T ; T+ )) est une solution maximale du problème de Cauchy
(PC), alors on a l’alternative suivante :

30

www.mathonec.com
 –ou bien T+ = b
–ou bien T < b et pour tout compact K de il existe t < T+ tel que '(t) 2
= K:
–Enoncé analogue pour T :

Démonstration:
Supposons que T+ < b et qu’il existe un compact K tel que '(t) 2 K pour tout t 2 (t0 ; T+ ):
Alors, comme f 2 C(D); il existe M > 0 tel que jjf (t; x)jjE M pour tout (t; x) 2 [t0 ; T+ ] K:
Soit (tn )n une suite croissante tendant vers T+ et telle que t0 < tn < T+ pour tout n: En
écrivant la solution ' sous forme intégrale, on obtient que :

Z tm
jj'(tm ) '(tn )jjE jjf (s; '(s))jjE ds M jtm tn j; 8m > n:
tn

La suite (tn )n étant de Cauchy, il en est de même pour ('(tn ))n qui est donc convergente.
Notons x1 = limn!1 '(tn ): Alors x1 2 K et on a donc (T+ ; x1 ) 2 D: La solution du
problème de Cauchy

x0 (t) = f (t; x(t)); x(T+ ) = x1 ;

admet selon le théorème 1.2 une solution locale qui prolonge ' au dela de T+ : Ceci contredit
la maximalité de ':On procède de façon analogue pour T :
Si est borné, ce théorème se traduit par : Le point de R E de coordonées (t; '(t)) tend
vers un point de la frontière du cylindre ]a; b[ quand t ! T+ et t ! T :

2.1.1 Cas particulier des EDO linéaires

Un cas particulier intéressant est celui où f est une application linéaire en x: Considérons I un
intervalle de R et N N fonctions réelles continues :

mij : t 2 I ! mij (t) 2 R:

Notons alors M (t) la matrice carrée N N dont les coe¢ cients sont les coe¢ cients sont les
fonctions mij (t): On notera $(E) l’espace vectoriel des matrices carrées N N dont la norme
naturelle est :

31

www.mathonec.com
 jjM jj$(E) = max jjM xjjE :
jjxjjE =1

On dira qu’une EDO 2.1 est linéaire si elle s’écrit :

x0 (t) = M (t)x(t) (2.6)

Théorème 2.7 Si mij 2 C(I) pour tout i; j = 1; :::; N; alors pour tout t0 2 I et x0 2 E;
le système (2.6) admet une unique solution maximale dé…nie sur I tout entier et telle que
x(t0 ) = x0 :

Démonstrtion:
Il est clair que
f (t; x) = M (t)x
O
est une fonction continue sur I E. Soit I(t0 ) un voisinage compact de t0 dans I (si t0 2 I ; on
peut choisir un intervalle de la forme [t0 ; t0 + ]; > 0; sinon, t0 est une extremité de I et
on peut choisir [t0 ; t0 + ] par exemple). Les fonctions aij étant continues sur I(t0 ), la fonction

k(t) = jjM (t)jj$(F )

est elle aussi continue sur I(t0 ). On peut alors considérer

k = max k(t)
t2I(t0 )

et f est uniformément lipschitzienne en x, de constante de lipschitz k sur I(t0 ) E:Selon le

théorème 1.3, il existe une unique solution locale ' au problème de Cauchy considéré. D’autre
part, en appliquant l’inégalité de Gronwall 2.5 avec '1 = ' et '2 = 0 ; on obtient l’estimation :

jj'(t)jjE jjx0 jjE ekjt t0 j

:

32

www.mathonec.com
 La solution reste donc bornée sur tout intervalle borné et suivant le théorème 1.6, elle peut
donc être prolongée sur I tout entier.
Nous reviendrons largement sur les EDO linéaires dans le chapitre suivant qui leur sera
dédié.

2.1.2 Dépendance des solutions en fonction des conditions initiales:

Lorsqu’elle est unique, la solution d’un problème de Cauchy:

x0 (t) = f (t; x(t));

x( ) = ;

peut être vue comme une fonction de la variable t dépendant du paramètre ( ; ) 2 R E:

Pour mettre en exergue cette dépendance, on écrira cette solution '(t; ; ): Noter que l’unicité
est fondamentale : si le problème de Cauchy ci-dessus admettait plusieurs solutions, la valeur
de '(t; ; ) ne serait pas dé…nie de façon univoque et on ne pourrait pas parler de fonction!

En considérons alors la fonction ' dé…nie sur un ouvert de R R E; on peut se poser la

question de la régularité de la fonction par rapport à l’ensemble de ses variables.

2.1.3 Continuité en fonction des conditions intiales

Si est une fonction continue dé…nie sur un intervalle I R et à valeurs dans E; on notera :

0 0
T ( ; ) = f(t; x) 2 I E; jjx (t)jjE < ; 8t 2 Ig;

le tube ouvert centré en (t; (t)): On peut alors énoncer le théorème.

Théorème 2.8 soit f 2 Lip(D) \ C(D) et soit une solution de 2.1 dé…nie sur un intervalle
I = [a; b]: Alors il existe > 0 tel que T ( ; ) D et pour tout ( ; ) 2 T ( ; ) il existe une
unique solution ' à 2.1 dé…nie sur I et véri…ant '( ; ; ) = : En outre, ' est continue sur
V =]a; b[ T ( ; ):

33

www.mathonec.com
 Démonstration:
Par dé…nition d’une solution de 2.1, l’ensemble f(t; (t)); t 2 Ig D: Cet ensemble étant
compact et le domaine D étant ouvert, il est possible de trouver 1 > 0 tel que T ( ; 1) D.
Choisissons alors <e k(b a) où k désigne la constante de Lipschitz de f sur D. Pour
1

tout ( ; ) 2 T ( ; );il existe d’après le théorème 1.3 une unique solution locale à l’ EDO 2.1
véri…ant '( ) = . L’estimation de Gronwall 2.5 fournit, sur l’intervallenous d’exercice de
'(:; ; ) :

jj'(t; ; ) (t)jjE ekjt j

jj ( )jjE < 1:

Ceci prouve que le point (t; '(t; ; )) reste à l’intérieur de T ( ; 1) D et donc, d’après le
théorème 1.6, '(:; ; ) peut être prolongée sur tout l’intervalle [a; b]:
On peut donc a¢ rmer que toute solution qui passe par un point de T ( ; ) est entièrement
contenue dans T ( ; 1 ):

2.1.4 Dépendance des solutions en fonction des conditions initiales

Pour montrer la continuité de la fonction ', nous allons montrer qu’elle est la limite uniforme
d’une suite de fonctions continues sur V . Introduisons pour cela la suite ('n )n dé…nie de la
façon suivante :

'0 (t; ; ) = (t) + ( );

Z t
'n+1 (t; ; ) = + f (s; 'n (s; ; ))ds; 8(t; ; ) 2 V; 8n > 0;

et montrons par récurrence qu’elle a les bonnes propriétés. Il est clair que '0 est continue
sur V et que pour tout n; la continuité sur V . Montrons maintenant que, pour tout n 1:
(t; 'm (t; ; )) 2 T ( ; 1 );

jt jm
jj'm (t; ; ) 'm 1 (t; ; )jjE km jj ( )jjE ; 81 m n; (2.7)
m!

8(t; ; ) 2 V:

34

www.mathonec.com
 D’une part :

jj'0 (t; ; ) (t)jjE = jj ( )jjE < 1;

Ce qui prouve que (t; '0 (t; ; )) 2 T ( ; 1) D pour tout (t; ; ) 2 V . D’autre part,
comme :

Z t
(t) = ( ) + f (s; (s))ds; 8t 2 [a; b];

on obtient que

Z t
jj'1 (t; ; ) '0 (t; ; )jjE j jjf (s; '0 (s; ; ))_f (s; (s))jjE dsj

Z t
j kjj'0 (s; ; ) (s)jjE dsj

= kjt jjj ( )jjE ;

ce qui est la relation 2.7 pour n = 1: On obtient également :

jj'1 (t; ; ) (t)jjE jj'1 (t; ; ) '0 (t; ; )jjE + jj'0 (t; ; ) (t)jjE

(1 + kjt j)jj ( )jjE ekjt j

jj ( )jjE 1;

ce quio prouve que (t; '1 (t; ; )) 2 T ( ; 1) pour tout (t; ; ) 2 V . Supposons maintenant
que 2.7 soit véri…ée pour un certain rang n. Alors

Z t
jj'n+1 (t; ; ) 'n (t; ; )jjE j jjf (s; 'n (s; ; )) f (s; 'n 1 (t; ; ))jjE dsj:

D’après 2.7, (t; 'n (t; ; )) et (t; 'n 1 (s; ; )) sont dans T ( ; 1) D, on peut donc écrire
que :

35

www.mathonec.com
 Z t Z t
j jjf (s; 'n (s; ; )) f (s; 'n 1 (s; ; ))jjE dsj j kjj'n (s; ; ) 'n 1 (s; ; )jjE dsj:

On utilise l’inégalité de 2.7, pour obtenir :

Z t Z t
js jn
j kjj'n (s; ; ) 'n 1 (s; ; )jjE dsj j k n+1 jj ( )jjE dsj
n!

jt jn+1
= k n+1 jj ( )jjE :
(n + 1)!

L’inigalité triangulaire fournie :

n
X
jj'n (t; ; ) (t)jjE jj'p (t; ; ) 'p 1 (t; ; )jjE + jj'0 (t; ; ) (t)jjE
p=1

n
X jt jp
= kp jj ( )jjE ekjt j
jj ( )jjE 1;
p!
p=0

et la propriété 2.7 est démontrée pour tout n 1; On en déduit que, pour tout n 0 et
tout m 1:

n+m
X jt jp
jj'n+m (t; ; ) 'n (t; ; )jjE kp jj ( )jjE
p!
p=n+1

n+m
X jb ajp
kp :
p!
p=n+1

La suite ('n )n est donc de Cauchy uniforme sur V . Sa limite, notée 'e est continue sur V .
En passant à la limite quand n ! 1 dans la relation :

Z t
'n+1 (t; ; ) = + f (s; 'n (t; ; ))ds;

On obtient que

36

www.mathonec.com
 Z t
e (t; ; ) = +
' e (t; ; ))ds;
f (s; '

e
ce qui prouve que ' ':
Remarquer que dans la démonstration de ce théorème, on a prouvé une nouvelle fois
l’éxistence (locale) d’une solution de 2.2.

2.1.5 Di¤érentiabilité en fonction des conditions initiales

Nous allons maintenant étudier la di¤érentiabilité de la solution ' en fonction des conditions
initiales ( ; ). Nous noterons Dx f la matrice (quand elle existe) :

0 1
@f1 @f1
B @x1 (t; x)::: @xN (t; x) C
B C
B : : C
B C
B C
Dx f = B : : C ; i = 1; :::; N ;
B C
B C
B : : C
@ A
@fN @fN
@x1 (t; x)::: @xN (t; x)

Théorème 2.9 En reprenant les mêmes notations et sous les mêmes hypothèses que dans le
théorème 1.8 précédent et en supposant de plus que Dx f existe et que Dx f 2 C(D), alors
' 2 C 1 (v):

Démonstration :
Prouver que ' est de classe C 1 est équivalent à prouver que toutes ces dérivées partielles
#' =#t; #' =# ; #' =# 1 ; :::; #' =# N
existent et sont continues sur V . La relation

Z t
'(t; ; ) = + f (s; '(t; ; ))ds;

et la continuité de ' assurée par le théorème précédent, nous donne immédiatement l’existence
et la contnuité de #' =#t : Montrons que #' =# 1
existe et est continue sur V . Pour cela, …xons

1 et comme dans le théorème précédent et considérons ( ; ) un point de T ( ; ): Comme

T ( ; ) est un ouvert, si l’on pose h = (h1 ; 0; :::; 0)T 2 E; alors, pour tout h1 proche de 0;

h = +h=( 1 + h1 ; 2 ; :::; N ) 2 T ( ; ):

37

www.mathonec.com
 En…n considérons :

'(t; ; h) '(t; ; )
h (t; ; )= :
h1

On peut alors énoncer :

Lemme 2.3 Pour tout " > 0; il existe " > 0 tel que si jh1 j " alors pour tout ( ; ) 2 T ( ; );
la fonction h (:; ; ) est une solution " approchée du problème de Cauchy linéaire :

y 0 (t) = Dx f (t; '(t; ; ))y(t);

y( ) = e1 ;

où e1 = (1; 0; :::; 0)T 2 E:

Démonstration du Lemme:
Posons:

h (t; ; ) = '(t; ; h) '(t; ; ):

En appliquant l’inégalité de Gronwall 2.5 à '(:; ; h) et '(:; ; ); on obtient que :

jj h (t; ; )jjE jj h ( ; ; )jjE ejt j

jh1 je(b a)
: (2.8)

Ainsi, h tend uniformement vers 0 sur V quand h ! 0: D’autre part, ' étand une solution
de 2.1, on déduit que :

0
h (t; ; ) = f (t; '(t; ; h )) f (t; ; '(t; ; )):

Le Théorème des accroisements …nis appliqué à la fonction C 1 :

Ft; ; ;h : s 2 [0; 1] 7! f (t; (1 s)'(t; ; h) + s'(t; ; h ));

38

www.mathonec.com
 nous assure de l’existence de

s1 = s1 (t; ; ; h) 2 [0; 1]

tel que :

0
h (t; ; ) = Ft; ; ;h (1) Ft; ; ;h (0) = Ft;0 ; ;h (s1 )

= Dx f (t; (1 s1 )'(t; ; ) + s1 '(t; ; )) h (t; ; ):

Cette relation devient :

0
h (t; ; ) = (Dx f (t; '(t; ; )) + h (t; ; )) h (t; ; )): (2.9)

si l’on pose :

h (t; ; ) = Dx f (t; (1 s1 )'(t; ; ) + s1 '(t; ; h )) Dx f (t; '(t; ; h )):

Les vecteurs '(t; ; ) et '(t; ; h) étant dans la boule de centre (t) et de rayon 1 par
convexité il en est de même de(1 s1 )'(t; ; ) + '(t; ; h )et donc (t; (1 s1 )'(t; ; ) +
s1 '(t; ; h )) 2 T( ; 1 ): Or T ( ; 1) est compacte, inclus dans D et Dx f est continue sur
D. Donc Dx f est uniformément continue sur T ( ; 1 ):

Pour tout " > 0; il existe e" > 0 tel que

jjDx f (t; x) e)jjF

Dx f (t; x "e (b a)
e 2 E; jjx
; 8x; x ejj
x e":

Un calcul simple conduit à l’estimation :

jj(1 s1 )'(t; ; ) + s1 '(t; ; h) '(t; ; h )jjE = js1 jjj h (t; ; )jjE

jh1 je(b a) :

En posant " = e" e (b a) ; on obtient donc que :

39

www.mathonec.com
 (b a)
jj h (t; ; )jjF "e ; 8(t; ; ) 2 V; 8jh1 j ": (2.10)

Comme

h = h =h1 ;

en rassemblant les estimations 2.8; 2.9 et 2.10, on obtient :

0
jj h (t; ; ) Dx f (t; '(t; ; )) h (t; ; )jjE "; 8(t; ; ) 2 V; 8jh1 j E:

Il su¢ t de véri…er (ce qui est évident) que

(t; ; ) = e1

pour avoir la conclusion du lemme.

Considérons maintenant la solution (t; ; ) du problème de Cauchy linéaire :

y 0 (t) = Dx f (t; '(t; ; ))y(t);

y( ) = e1 :

D’après le théorème 1.7, cette solution existe pour tout t 2 [a; b] et d’après le lemme, pour
tout " > 0, il existe E > 0 tel que h (t; ; ) soit une solution " approchée si jh1 j E: En
utilisant l’estimation de Gronwall 2.5, on en déduit que :

" k(b a)
jj h (t; ; ) (t; ; )jjE (e 1);
k

pour tout (t; ; ) 2 V: En d’autres termes, h (t; ; ) ! (t; ; ) uniformément sur V . On

en conclut que
@' =@ 1
=

40

www.mathonec.com
existe et est une fonction continue sur V car c’est la limite uniforme des fonctions h qui sont
continues sur V:

41

www.mathonec.com
Chapitre 3

EQUATIONS LINEAIRES

Dans ce chapitre, on s’intéresse exclusivement aux équations di¤érentielles du type

du
= A(t)u + b(t); (3.1)
dt

avec b 2 C(I; X) et A 2 C(I; $(X)) (où $(X) désigne l’espace des applications linéaires
continues dans l’espace de Banach X). On dit que sont des équations <<linéaires>> (car
A(t)est linéaire ) avec <<terme source>> (b). On quali…e parfois 3.1 simplement d’équation
linéaire (par abus de langage, le terme correct étant plutôt a¢ ne).

3.1 Existence globale

D’après le théorème II.3, les solutions de 3.1 sont globales.

Théorème 3.1 Pour b 2 C(I; X) et A 2 C(I; $(X))); quel que soit (t0 ; u0 ) 2 I X; il existe
une unique fonction u 2 C 1 (I; X) solution de 3.1 telle que

u(t0 ) = uO :

On peut aussi donner une démonstration de ce théorème en faisant appel au théorème de

point …xe unique.

42

www.mathonec.com
Théorème 3.2 Soit T un opérateur sur un espace de Banach dont une puissance est contrac-
tante. Alors T admet un point …xe unique.

Démonstration:
Ce resultat est une conséquence facile du théorème de point …xe de Banach-picard.
Par hypothèse, il existe un entier k tel que T k soit contractant et donc admette un point
…xe unique u. Alors :

T k (T (u)) = T (T k (u)) = T (u):

Donc T (u) est aussi un point …xe de T k : à cause de l ’unicité de ce point …xe on a ainsi

T (u) = u:

De plus, tout point …xe de T étant évidemment un point …xe de T k : u est l unique point …xe
de T:
Démonstration du théorème III.1 :
Pour résoudre le problème de cauchy de donnée ‹‹initiale›› u(t0 ) = u0 ; on cherche u tel que

Z t
u(t) = u0 + (A(s)u(s) + b(s))ds
t0

c’est-à-dire de façon équivalente

u(t) = v(t) + u(t);

avec
Z t
v(t) = (A(s)u(s) + b(s))ds:
t0

Pour simpli…e, on suppose (sans perte de généralité, il su¢ t de translater toutes les fonctions
en jeu) t0 = 0: soient et tels que <0< et [ ; ] I: Dé…nissons alors l’opérateur :

T : C([ ; ]; X) ! C([ ; ]; X)

43

www.mathonec.com
 v 7! T v

Z t0
(T v)(t) = (A(s)(v(s) + u0 ) + b(s))ds:
0

Quel que soient v et w dans C(j ; ]; X) on a

Z t
(T v)(t) (T w)(t) = A(s)(v(s) w(s))ds
0

et donc

jj(T v)(t) (T w)(t)jj Ct max jjv wjj;

s2[0;t]

avec
C := max jjA(s)jj:
s2[ ; ]

Montrons par récurrence que

jj(T n v)(t) (T n w)(t)jj C n max jjv wjj

s2[0;t]

pour tout t 2 [ ; ]: C’est vrai à l’ordre 1 comme on vient de le voir. Supposons l’inégalité
vraie à l’ordre n. Alors

jj(T n+1 v)(t) (T n+1 w)(t)jj = jjT ((T n v)(t)) T ((T n w)(t))jj

Z t
C jj(T n v)(s) (T n w)(s)jjds
0

Z t
sn
C max jjv wjj Cn ds
s2[0;t] 0 n!

d’aprés l’hypothèse de récurrence. En évaluant l’intégrale on obtient immédiatement l’inégalité

à l’ordre n + 1 :

44

www.mathonec.com
 tn+1
jj(T n+1 )(t) (T n+1 !)(t)jj Cn+1 jj !jj:
(n + 1)!

Par conséquent, on a

( )n
jjT n T n !jj Cn jj !jj
n!

quels que soient et !. Donc T n est contractant pour n assez grand et d’aaprès le théorème
III.2, T admet un point …xe unique.

3.2 Résolvante

Considérons l’équation di¤érentielle linéaire dans $(X) :

dM
= A(t) M (3.2)
dt

D’après le théorème III.1, ses solutions maximales sont globales.

Dé…nition 3.1 On appelle résolvante de l’équation di¤ érentielle

du
= A(t)u
dt

l’application
R:I I!$

(t; t0 ) ! R(t; t0 )(X)

où
t ! R(t; t0 )

est la solution du problème de Cauchy pour (3.2) avec la condition initiale

R(t0 ; t0 ) = IX :

45

www.mathonec.com
 On observe que pour tout u0 2 X, l’application

u : t ! R(t; t0 )u0

est la solution du problème de Cauchy

pour (3.1) et la condition initiale u(t0 ) = u0 . Autrement dit, le ‡ot de (3.1) est donné par

t0
(t; u0 ) = R(t; t0 )u0 :

Par construction, R est continûment di¤ érentiable par rapport à t. De plus, le lemme
II.3(application simple du lemme de Gronwall) montre que R est lipschitzienne par rapport
à t0 . Donc en particulier, R est une fonction continue de (t; t0 ). Comme l’ensemble des
isomorphismes de $(X) est un ouvert et R(t; t0 ) = IX est trivialement un isomorphisme,
R(t; t0 ) est un isomorphisme pour tout t voisin de t0 . En fait, ceci est vrai pour tout t 2 I
comme conséquence de la :

proposition 3.1 La résolvante est telle que

R(t; s) R(s; t0 ) = R(t; t0 )

quels que soient t; s; t0 2 I:

Démonstration :
Soit u0 2 X et
u(t) = R(t; t0 )u0 :

Alors la solution du problème de Cauchy

d
dt = A(t)
(s) = u(s)

n’est autre que u, mais elle s’exprime aussi comme

46

www.mathonec.com
 : t ! (t) = R(t; s)u(s) = R(t; s) R(s; t0 )u0 :

D’où

R(t; t0 )u0 = R(t; s) R(s; t0 )u0 :

Ceci étant vrai quel que soit u0 , on en déduit la formule voulue.

Corollaire 3.1 Pour tout (t; t0 ) 2 I I; R(t; t0 ) est un isomorphisme, d’inverse R(t0 ; t). De
plus, R estcontinûment di¤ érentiable comme fonction de deux variables, et ses dérivées partielles
satisfons

@R @R
(t; s) = A(t) R(t; s)et (t; s) = R(t; s) A(s):
@t @s

Démonstration :
La première partie est conséquence immédiate de la proposition III.1:

R(t0 ; t) R(t; t0 ) = R(t0 ; t0 ) = IX etR(t; t0 ) R(t0 ; t) = R(t; t) = IX :

L’expression de la dérivée partielle de R par rapport à sa première variable découle di-

rectement de la dé…nition de R : elle est continue comme composée de fonctions continues
(on rappelle que R est continue comme fonction des deux variables). Pour montrer que R est
dérivable par rapport à sa seconde variable, on utilise la première partie : comme

1
R(t; s) = R(s; t) ;

R est dérivable par rapport à s par composition de s ! R(s; t) et de l’application

Isom(X) ! Isom(X)

47

www.mathonec.com
 1
M ! M ;

dont on sait qu’elle est (continûment) di¤érentiable. Pour exprimer e¤ectivement la dérivée
partielle de R par rapport à sa seconde variable, on dérive simplement l’identité

R(t; s) R(s; t) = IX

en utilisant l’expression de la dérivée partielle de R par rapport à sa première variable, ce

qui donne :

@R
(t; s) R(s; t) + R(t; s) A(s) R(s; t) = 0:
@s

Comme R(s; t) est un isomorphisme, on en déduit l’expression annoncée.

Attention ! En général, on ne connaît pas explicitement la résolvante.

3.2.1 Cas des équations linéaires autonomes

Si A est indépendant de t, la résolvante de 3.1 s’exprime à l’aide de l’exponentielle, dont on

rappelle la

Dé…nition 3 Pour tout A 2 $(X);

+1
X
A 1 n
e = A :
n!
n=0

(La série ci-dessus est normalement convergente dans $(X)).

Lemme 3.1 Pour tout A 2 $(X); l’application t ! etA est continûment dérivable (et même
de classe C 1 ), et
d tA
e = AetA = etA A
dt

48

www.mathonec.com
 La démonstration est laissée en exercice. Elle utilise de façon cruciale la propriété de
l’exponentielle :

e(t+s)A = etA esA :

Le lemme III.3 montre que la résolvante de l’équation di¤érentielle linéaire autonome

du
= Au
dt

est donnée par

R(t; s) = e(t s)A

:

3.2.2

3.2.3 Cas des équations linéaires en dimension …nie

Dimension 1.

Les équations di¤érentielles (d’ordre 1) linéaires scalaires, pour lesquelles X = R, se résolvent

explicitement : la résolvante de

du
= a(t)u
dt

où l’inconnue u est à valeurs réelles et a : I ! R est continue est donnée par

Rt
a( )d
R(t; s) = e s :

Dimension 2 et plus

Il est tentant de généraliser la formule précédente. Cependant, si les matrices A(t) (ou les
endomorphismes) ne commutent pas deux à deux, une telle formule ( que l’on se garde volon-

49

www.mathonec.com
tairement d’écrire) est fausse : ceci est lié au fait que eA+B n’est pas le produit de eA et eB
lorsque les matrices A et B ne commutent pas .
Tout ce que l’on peut dire est que la résolvante est liée à une base de solutions.

proposition 3.2 Soit A : I ! Mn (R) continue, et R : I ! Mn (R) la résolvante de

du
= A(t)u; (3.3)
dt

c’est-à-dire que pour tout (t; s) 2 I I;

@R
(t; s) = A(t)R(t; s)etR(s; s) = In :
@t

Si (e1 ; :::; en ) est une base de Rn ; alors les applications t ! R(t; t0 )ei (pour i 2 f1; :::; ng)
forment une famille indépendante de solutions de 3.3, soit B : t 2 I ! B(t) 2 GLn (R) dont les
vecteurs colonnes sont
u1 (t); :::; un (t): Alors la résolvante de 3.3 est donnée par

1
R(t; s) = B(t)B(s) :

Démonstration :
La première partie découle de la dé…nition et du corollaire III.1, qui montre que R(t; s) 2
GLn (R). La preuve de la seconde partie consiste simplement à remarquer que l’application
t ! B(t)B(s) 1 est solution du même problème de Cauchy que t ! R(t; s):

Remarque 3.1 On peut aussi montrer que R(t; s) est inversible en véri…ant que son détermi-
nant satisfait une équation di¤ érentielle scalaire linéaire (d’ordre 1). C’est l’objet du résultat
suivant, parfois appelé formule de Liouville.

proposition 3.3 Soit

(t; s) := det(R(t; s)):

50

www.mathonec.com
Alors

Rt
tr(A( ))d
(t; s) = e s :

Démonstration :
C’est un petit calcul qui utilise la di¤érentielle du déterminant, dont on rappelle qu’elle
donnée par :

d(det)(M ) H = tr((comM )0 H);

Par conséquent, la règle de dérivation des fonctions composées montre que :

@ @R
(t; s) = tr((comR(t; s))t (t; s)):
@t @t

D’où

@
(t; s) = tr((comR(t; s))t A(t)R(t; s)) = tr((comR(t; s))t R(t; s)A(t)) = (t; s)tr(A(t))
@t

d’après la formule
(comR)t R = (detR)I:

3.3 Formule de Duhamel

Ce paragraphe est consacré à un résultat très facile à démontrer et néaumoins fondamental.

Il exprime le principe de superposition bien connu selon lequel : la solution générale d’une
équation linéaire avec terme source est donnée pr la solution générale de l’équation homogène
plus une solution particulière de l’équation avec terme source.

proposition 3.4 (Formule de Duhamel)

51

www.mathonec.com
 Soit A 2 C(I; $(X)) et R la résolvante de l’équation di¤ érentielle homogène :

du
= A(t)u
dt
t0
Soit b 2 C(I; X) et le ‡ot en t0 2 I de l’équation

du
= A(t)u + b(t):
dt

Alors pour tout (t; ) 2 I X;

Z t
t0
(t; ) = R(t; t0 ) + R(t; s) b(s)ds:
t0

Démonstration :
On véri…e par un calcul facile que l’application

Z t
t! R(t; s) b(s)ds
t0

est solution de l’équation avec le terme source b. Comme t ! R(t; t0 ) est solution de
l’équation homogène, la somme des deux est, par linéarité, solution de l’équation avec terme
source. De plus elle vaut à t = t0 : Donc c’est l’unique solution cherchée.

3.4 Equations linéaires autonomes

L’objectif ici est d’introduire certains outils utiles dans la suite du chapitre. On commence par
le cadre plus élémentaire des équations en dimension …nie.

3.4.1 Dimension …nie

On considère dans ce paragraphe A 2 Mn (R) et l’équation di¤érentielle linéaire homogène

associée

52

www.mathonec.com
 du
= Au: (3.4)
dt

On sait que les solutions sont dé…nies sur R tout entier. Que peut-on dire de leur com-
portement à l’in…ni? Cela dépend à l’évidence des propriétés de la matrice A, que l’on va pour
l’occasion identi…er à un endomorphisme de Cn : L’analyse de 3.4 est intimement liée à l’algèbre
linéaire, dont voici quelques rappels (pour information le vocabulaire anglophone, lorsqu’il n’est
pas intuitif, est donné en note).

Dé…nition 3.2 Les valeurs propres de A sont les nombres complexes pour lesquels le noyeau
de (A I) est non triviale. Si est une valeur propre de A, on dit que l’espace Ker(A I)
est le sous-espace propre associé. Si k est le plus grand entier tel que

Ker(A I)k 1
Ker(A I)k

avec inclusion stricte,

l’espace Ker(A I)k est applé sous-espace caractéristique se :
- La dimension du sous-espace propre Ker(A I) est appelée multiplicité géométrrique de
:
- La dimension du sous-espace caractéristique Ker(A I)k est appelée multiplicité al-
gébrique de :
- Lorsque les multiplicités géométriques et algébriques coincident, la valeur propre est dite
semi-simple.
- lorsque les multiplicités géométriques et algébriques coincident et valent 1, la valeur propre
est dite simple.

Théorème 3.3 - Les valeurs propres de A sont les zéros du polynôme caractéristique :

P ( := det( i A):

La multiplicité algébrique d’une valeur propre est égale à son ordre comme zéro du polynôme
caractéristique.

53

www.mathonec.com
 Les sous-espaces caractéristiques sont en somme directe et leur somme est Cn :
-(Cayley-Hamilton) Les valeurs propres de A sont aussi les zéros du polynôme minimal,
dé…ni comme le générateur unitaire de l’idéal des polynômes anulateurs de A.
- Si toute les valeurs propres sont semi simples, les sous-espaces caractéristiques et les sous-
espaces propres sont égaux, et la matrice A est diagonalisable.
Pour la démonstration, se repporter à un ouvrage d’algèbre linéaire...
Attention! Le vocabulaire se télescope un peu car on appelle aussi résolvante de A la
matrice (A I) 1 lorsque n’est pas une valeur propre de A. Le terme est donc à manipuler
avec précaution.
Le comportement asymptotique des solutions de 3.4 est donné par le résultat suivant, que
l’on donne sur R+ et qui admet évidemment un analogue sur R (changer t en t et A en
A).

proposition 3.5 L’application t ! etA est bornée sur R+ si et seulement si les valeurs propres
de A sont touts de partie réelle négative ou nulle et si les valeurs propres imaginaires pures sont
semi-simple. Elle tend vers 0 en +1 si et seulement si les valeurs propres de A sont toutes de
partie réelle strictement négative.

Démonstration :
Soit f 1 ; :::; pg l’ensemble des valeurs propres de A (avec p n et i 6= j pour i 6= j). A
chaque valeur propre j on associe le sous-espace caractéristique Ej , et j la projection sur Ej
parallèlement à i6=j Ej . Par dé…nition de Ej , il existe un entier kj tel que

kj
(A j I) j = 0;

et d’après le théorème III.3, on a

p
X
I= j:
j=1

Par suite

54

www.mathonec.com
 p p p kj 1 k
X X X X t
tA tA t t(A j I) t k
e = e j = e j
e j = e j
(A j I) j;
k!
j=1 j=1 j=1 k=0

ce qui tend bien vers 0 en +1 si les j sont toutes de partieréelle strictement négative. Par
ailleurs, si les j sont de partie réelle négative ou nulle et si les kj valent tous 1, alors etA est
bornée sur R+ . Les réciproques sont faciles à véri…er :
-S’il existe de partie réelle strictement positive tel que

A =

pour un vecteur non nul, alors

etA = et

est non borné.

-Si 2 iR est une valeur propre non semi-simple, il existe des vecteurs non nuls et ! tels
que

A = ; A! = ! + ;

On montre par réccurence que

Am ! = m
!+m m 1

est de norme supérieure ou égale à tjj jj jj!jj pour t 0, ce qui est évidemment non
borné.
-En…n, si 2 iR est une valeur propre même simple,

A =

avec non nul entaîne que

etA = et

55

www.mathonec.com
 est de norme égale à jj jj, ce qui ne tend pas vers 0:

Dé…nition 3.3 Une matrice n’ayant pas de valeurs propre imaginaire pure est dite hyper-
bolique.

proposition 3.6 Soit A 2 Mn (R) une matrice hyperbolique. Alors les ensembles

E s = f! 2 Rn ; lim etA ! = 0g;

t!+1

E u = f! 2 Rn ; lim etA ! = 0g
t! 1

sont des sous-espaces vectoriels de Rn stables par A, appelés respectivement sous-espace

stable et sous-espace instable, et il existe des projecteurs u et s commutent avec A tels que
u + s = I satisfaisant les estimations suivantes. En notant (A) l’ensemble des valeurs
propres de A et

= maxfRe < 0; 2 (A)get = minfRe > 0; 2 (A)g;

pour tout " > 0; il existe des constantes b" et c" telles que

jjetA s
jj b" e ( ")t
8t 0

et
jjetA u
jj c" e( ")t
8t 0

Démonstration :
Le fait que E s et E u soient des sous-espaces vectoriels de Rn stables par A est évident
(rappelons que etA commute avec A). On va commencer par dé…nir les projecteurs u et s,

et on véri…era ensuite que E s = Im( s) et E u = Im( u ): On reprend pour cela les mêmes
notations que dans la preuve de la proposition III.5. Soient

56

www.mathonec.com
 X X
s u
:= j et := j:
j;Re j <0 j;Re j >0

Noter que ces projecteurs sont bien à valeurs dans Rn : en e¤et, les valeurs propres non
réelles de A sont deux à deux conjuguées, et j0 = j si j0 = j donc j + j0 est à valeurs
dans Rn :
Comme les s u
j; et sont bien des projecteurs commutant avec A. Et comme A n’a pas
de valeur propre imaginaire pure on a

p
X
u s
+ = j = I:
j=1

De plus, on a l’expression explicite :

kj 1 k
X X t
tA s t k s
e = e j
(A j I) ;
k!
j;Re j <0 k=0

d’où, pour t 0;

K k
X
tA s k s t t
jje jj p max jj(A j I) jjjj jje
j p;k k k!
k=0

avec
K = max kj 1:
j

Le dernier terme, polynômial en t, étant majoré par une constante fois e"t pour " > 0, on
en déduit l’estimation de jjetA s jj. L’estimation de jjetA s jj s’obtient exactement de la même
manière. On déduit directement de ces estimations les inclusions

s
Im( ) E s et Im( u
) Eu:

Réciproquement, soit par exemple ! 2 E s : Comme u commute avec A, il commute aussi

avec etA , donc

u
jj !jj = jjetA u
e tA
!jj c" e( ")t
jje tA
!jj

57

www.mathonec.com
 pour tout t 0. En prenant " < , le second membre tend vers 0 lorsque t tend vers 1
(pour le dernier terme, on utilise simplement la dé…nition de E s ). Donc u! = 0, c’est à dire
que
s s
!= ! 2 Im( ):

Ainsi E s Im( s ). On montre de façon analogue que E u est inclus dans Im( u ):

3.4.2 Dimension in…nie

Pour simpli…er la présentation, on suppose dans ce paragraphe que X est un C-espace de

Banach, et que A est un opérateur linéaire et continu sur X, ce que l’on note A 2 $(X):
On a vu au paragraphe précédent des résultats sur le comportement de etA lorsque X est
de dimension …nie, qui reposent pour une large part sur la réduction des matrices. Lorsque
X est de dimension in…nie, tout devient plus délicat. Tout d’abord, le spectre de A n’est plus
exclusivement constitué de valeurs propres.

Dé…nition 3.4 On appelle specter de A 2 $(X) l’ensemble (A) des nombres complexes
tels que (A IX ) 2
= Isom(X) (ensemble des isomorphismes de X). Le spectre penctuel est
inclus dans le spectre : c’est l’ensemble des valeurs propres, c’est à dire des nombres 2 C
pour lesquels (A IX ) est non injectif.

Remarque 3.2 - L’ensemble des isomorphismes étant un ouvert, l’ensemble résolvant

de A :

(A) := Cn (A)

est un ouvert, par continuité de l’application

C ! $(X)

(A) := Cn (A)

58

www.mathonec.com
 est un ouvert, par continuité de l’application

C ! $(X)

! (A IX ):

- Le spectre d’un opérateur A 2 $(X) est borné :

(A) f ;j j jjjAjjjg:

En e¤ et, si
j j > jjjAjjj; (A IX )

est inversible d’inverse

+1
1 1 X 1 m
(A IX ) = mA ;
m=0

cette série étant normalement convergente.

Plus précisément, on a

1
(A) f ;j j rA := lim jjjAm jjj m g (3.5)

et
(A) \ f ; j j = rA g =
6 ?:

- Ainsi, (A) est un compact non vide de C:

Attention ! Si A était un opérateur non borné (c’est à dire linéaire mais pas continu),
on pourrait encore dé…nir son spectre mais il pourrait être vide ou non borné.

Lemme 3.2 L’application

RA : (A) ! $(X)

59

www.mathonec.com
 1
! (A IX )

est holomorphe (c’est à dire dérivable comme fonction de variables complexes).

Théorème 3.4 Si (Bm ) est une famille d’éléments d’un C-espace de Banach ßtel que pour
tout ' 2ß
’,

supj'(Bm )j < 1
m

alors
supjjBm jj < 1:
m

Démonstration :
Cela revient simplement à combiner le théorème de Hahn-Banach avec le théorème de
Banach-Steihaus : le premier montre en e¤et que

jjBm jjß= jjbm jj$(ß’;C)

où
hbm ; 'i := '(Bm )8' 2 ß
0;

et le second dit que

supj hbm ; 'i j < 1 8' 2 ß

0 ) supjjbm jj$(ß’,C) < 1:
m m

Grâce à l’astuce mentionnée plus haut (utiliser des formes linéaires continues ' pour revenir
dans C), on peut faire passer les résultats classiques de la théorie des fonctions de variable
complexe aux fonctions à valeurs dans un algèbre de Banach. En particulier, on a le

60

www.mathonec.com
Théorème 3.5 Soit un ouvert connexe non vide de C et ßune algèbre de banach. Si
f : !ßholomorphe elle est analytique, c’est à dire développable en série entière (et récipro-
quement). De plus, si et sont des lacets homotopes dans , alors

Z Z
f (z)dz = f (z)dz:
0

Remarque 3.3 Dans le cas de l’algèbre ß=$(X) et de la fonction RA , la première partie

de l’énoncé est facile à véri…er. On peut même préciser l’observation faite au lemme III.2
en

proposition 3.7 L’application RA est analytique. De plus, RA n’admet pas de prolongement

analytique en dehors de (A):

Démonstration :
Pour ; 0 2 (A), on a

RA ( 0) = (IX ( 0 )RA ( 0 ))RA ( ):

Par suite, on a

+1
X
n n+1
RA ( ) = ( 0 ) (RA ( 0 )) ;
n=0

cette série étant convergente pour

j 0j < 1=jjjRA ( 0 )jjj:

Ceci montre que RA est analytique dans (A). <<Inversement>> si la série ci-dessus converge,
alors sa somme est

61

www.mathonec.com
 +1
X
n n+1 1
SA ( ) := ( 0 ) (RA ( 0 )) = (IX ( 0 )RA ( 0 )) RA ( 0 ):
n=0

Or

RA ( 0 )(A IX ) = IX ( 0 )RA ( 0 ):

Donc
SA ( )(A IX ) = IX :

On montre de la même façon que

(A IX )SA = IX :

D’où 2 (A) et
RA ( ) = SA ( ):

Par conséquent, RA n’admet pas de prolongement analytique en dehors de (A):

Pour s’en convaincre, on fait un raisonnement classique mais un peu acrobatique. Supposons
qu’une fonction analytique f prolonge RA au voisinage d’un point du bord de (A): Par
dé…nition, n’appartient pas à l’ouvert (A). Soit r le rayon de convergence du développement
en série entière de f au point . Il existe un point 0 2 (A) \ f ; j j < r=2g; et le
développement en série entière de f au point 0 a un rayon au moins égal à r j 0 j>j 0 j:
En particulier, est dans le disque de convergence. Or, f coincide avec RA au voisinage de 0

et son développement en série entière est donc

+1
X
n n+1
f( ) = ( 0 ) (RA ( 0 )) :
n=0

D’après ce qui précède, tous les points du disque de convergence devraient être dans (A),
mais ce n’est pas le cas de .
Nous sommes maintenant armés pour aborder l’extension de la proposition III.6 en dimen-
sion in…nie. Comme en dimension …nie, on dit d’un opérateur qu’il est hyperbolique s’il n’a pas
de spectre imaginaire pur.

62

www.mathonec.com
Théorème 3.6 Si A est un opérateur hyperbolique, il existe des projecteurs continus s et
u comme dans la proposition III.6.

3.5 Equations di¤érentielles linéaires à coe¢ cients péri-

odiques

Soit X un espace de Banach (surR ou C) et A 2 C(R; X) périodique de période T . On

considère l’équation di¤érentielle linéaire

du
= A(t)u;
dt

et
R:R R ! Isom(X)

sa résolvante :

@R
(t; s) = A(t) R(t; s) et R(s; s) = IX.
@t

Dé…nition 3.5 Quel que soit s 2 R, on appelle l’opérateur

C(s) := R(s + T; s)

opérateur de monodromie (au point s).

proposition 3.8 Les opérateurs de monodromie sont tous conjugués. Pour qu’il existe une
solution T -périodique non triviale de l’équation di¤ érentielle, il faut et il su¢ t que 1 soit valeur
propre de ces opérateurs.

63

www.mathonec.com
Théorème 3.7 ((Floquet-Lyapunov)
Si X est un C-espace vectoriel de dimension …nie et A 2 C(R; $(X)) est T -périodique, il
existe une application T -périodique Q : R ! Isom(X) et B 2 $(X) telle que u est solution de

du
= A(t)u
dt

si et seulement si
t ! v(t) := Q(t)u(t)

est solution de

dv
= Bv:
dt

Remarque 3.4 Dans le cadre de ce théorème, on a inversement

u(t) = P (t)v(t);

avec

1 tB
P (t) = Q(t) = R(t; 0)e :

La version réelle du théorème est le

Théorème 3.8 (Floquet-lyapunov)[?]

Si X est un R-espace vectoriel de dimension …nie et A 2 C(R; $(X)) est T -périodique, il

existe une application 2T -périodique Q : R ! Isom(X) et B 2 $(X) telle que u est solution de

du
= A(t)u
dt

si et seulement si

t ! v(t) := Q(t)u(t)

64

www.mathonec.com
est solution de

dv
= Bv:
dt

Dé…nition 3.6 Les valeurs propres des opérateurs de monodromie sont appelées multiplica-
teurs caractéristiques. Les nombres tels que e T est un multiplicateur caractéristique sont
appelés exposants de Floquet.

Théorème 3.9 Les solutions de

du
= A(t)u
dt

sont bornées sur R+ si et seulement si les multiplicateurs caractéristiques sont tous de module
inférieure ou égale à 1 sont semi-simple.
Toutes les solutions tendent vers 0 en +1 si et seulement si les multiplicateurs caractéris-
tiques sont tous de module strictement inférieur à 1.
Considérons maintenant l’équation avec terme source :

du
= A(t)u + b(t);
dt

avec b périodique de même période que A

proposition 3.9 Si 1 n’est pas dans le spectre des opérateurs de monodromie de l’équation
homogène, l’équation avec terme source admet au moins une solution périodique.

Remarque 3.5 Lorsque b 0 la solution périodique en question est identiquement nulle !

Théorème 3.10 Supposons X de dimension …nie. Pour que l’équation

65

www.mathonec.com
 du
=A(t)u + b(t)
dt

avec A et b périodique de même période, admette une solution périodique, il su¢ t qu’elle
admette une solution brnée sur R+ :

Soit P l’application a¢ ne dé…nie par

ZT
P (v) = R(T; 0)(v + R(0; s)b(s)ds);
0

où R est la résolvante de l’équation homogène. Toute solution u de l’équation avec teme

source véri…e d’après la formule de Duhamel (itérée), pour tout n 2 N,

u(nT ) = P n (u0 )

où
u0 = u(0):

Donc s’il existe une solution bornée, c’est qu’il existe u0 2 X tel que la suite P n (u0 ) soit
bornée.

On va raisonner par l’absurde et montrer que la non existence de solutions périodiques

empêche P n (u0 ) d’être bornée. S’il n’y a pas de solution périodique c’est que l’opérateur P n’a
pas de point …xe. Autrement dit, le vecteur

ZT
y := R(T; 0) R(0; s)b(s)ds
0

n’appartient pas à
E := Im(IX R(T; 0));

qui est fermé puisqu’on a supposé l’espace de dimension …nie. Donc il existe ' 2 E ? tel que
'(y) 6= 0 (sinon on aurait y 2 (E ? )? = E = E). Par construction ceci implique

66

www.mathonec.com
 '(v R(T; 0)v) = 0

pour tout v 2 X;
c’est à dire
'(v) = '(R(T; 0)v)

et donc
'(v) = '(R(T; 0)n v)

pour tout entier n. Or une récurrence facile montre que

n
X1
P n (v) = R(T; 0)n v + R(T; 0)k y:
k=0

On déduit

'(P n (v)) = '(v) + n'(y);

ce qui n’est pas borné puisque '(y) 6= 0:

Remarque 3.6 En dimension in…nie, la même démonstration fonctionne, à condition de sup-

poser

Im(IX R(T; 0)) = Im(IX R(T; O)):

67

www.mathonec.com
Chapitre 4

EQUATIONS AUTONOMES

On revient maintenant à des équations non-linéaires, supposées autonomes c’est à dire de la

forme

du
= f (u); (4.1)
dt

où f est une fonction de classe C 1 sur un ouvert U d’un espace de Banach X dans X. On
dit qu’une telle fonction f est un champ de vecteurs. L’étude de l’équation di¤érentielle 4.1
revient à l’étude des champs de vecteurs.

4.1 Champs de vecteurs

Dé…nition 4.1 On appelle intégrale première d’un champ de vecteur f sur U une fonction
E 2 C 1 (U; R) telle que
dE(u) f (u) = 0

qul que soit u 2 U:

On dit aussi que E est une intégrale première de l’équation di¤ érentielle 4.1.
Par dé…nition, quel que soit u 2 C 1 (J; U ) solution de 4.1, si E est une intégrale première
de 4.1 alors E u est indépendant de t.

68

www.mathonec.com
Exemple 4.1 Le champ de vecteurs associé à la loi de Newton en mécanique est de la forme

f : (x; ) 2 Rn Rn ! ( ; F (x)):

Lorsque F <<dérive d’un potentiel V >> c’est à dire

F = grad V;

alors l’énergie totale

1
E(x; ) = jj jj2 + V (x)
2

est une intégrale première.

La connaissance des intégrales premières est cruciale pour espérer résoudre explicitement
l’équation di¤ érentielle; et elle est toujours utile pour montrer des propriétés qualitatives de
cette équation.

Dé…nition 4.2 Un champ de vecteurs f est dit complet si et seulement si touts les solutions
maximales de 4.1 sont globales, c’est à dire dé…nies sur R tout entier.

Attention! La régularité d’un champ de vecteurs n’a rien à voir avec le fait qu’il soit
complet ou non, comme le montre l’exemple très simple de l’équation de Riccati, où X = R;

f (u) = u2

(donc f est analytique!) et les solutions non triviales <<explosent >> en temp …ni.
Pour démontrer qu’un champ de vecteurs est complet, on peut faire appel au

Théorème 4.1 Supposont que pour tout intervalle ouvert borné J, toute solution u 2 C 1 (J; U )
de

du
= f (u)
dt

soit à valeurs dans un compact KJ de U . Alors le champ f est complet.

69

www.mathonec.com
 Démonstration :
C’est une conséquence du théorème des bouts (théorème II.2), un peu subtile car le compact
KJ peut a priori dépendre dangereusement du diamètre de J. En fait, on va appliquer le
théorème des bouts à un système (équivalent) augmenté.
Soit u 2 C(J; U ) une solution maximale de

du
= f (u):
dt

Alors l’application
s 2 J ! (u; t) = (u(s); s)

est une solution maximale du système augmenté :

du
= f (u);
ds

dt
= 1:
ds

Supposons que = sup J et = inf J soient …nis. Choisissons T > max( ; ):

D’après le théorème II.2, la solution maximale (u; t) du système augmenté doit sortir du
compact KJ [ T; T ]: Comme u reste dans KJ , c’est que t sort de l’intervalle [ T; T ] et donc
de J; ce qui est absurde puisque t = s 2 J:
Donc l’une au moins des bornes de J est in…nie. Et en fait les deux le sont. On s’en
convaincre aisément en reprenant la démonstration du théorème des bouts : si on avait par
exemple J =] 1; [ avec < +1, alors, pour tout < et T > max( ; ); la solution
(u; t) 2 C 1 (] ; [) du système augmenté devrait sortir du compact K] ; [ [ T; T ] lorsque
s % ; ce qui est impossible.
En dimension …nie, véri…er les hypothèses de ce théorème revient simplement à trouver des
estimations a priori, c’est à dire une majoration des solutions sur les intervalles bornés.

70

www.mathonec.com
Exemple 4.2 Le champ de vecteurs de la mécanique

f : (x; ) 2 Rn Rn ! ( ; grad V (x))

est complet pourvu que V soit une fonction positive.En e¤ et, rappelons que

1
E(x; ) = jj jj2 + V (x)
2

est une intégrale première. Si V est positive, si t ! (x(t); (t)) est une solution de l’équation
di¤ érentielle associée à f , alors
(t) = x0 (t);

d’où
Zt
x(t) = x(0) + (s)ds;
0

et

1
jj (s)jj2 E(x(s); (s)) = E(x(0); (0));
2

et donc
p
jjx(t)jj jjx(0)jj + T 2E(x(0); (0))

pour t 2] T; T [:
Ainsi la restriction de (x; ) à tout intervalle borné est à valeurs dans un ensemble fermé
borné de R2n :
Donc on peut appliquer le théorème IV.1.

Dé…nition 4.3 Les trajectoires des solutions maximales de 4.1, c’est à dire les courbes dans
X de la forme fu(t); t 2 Jg où u 2 C 1 (J; U ) est une solution maximale de 4.1, sont appelées
courbes intégrales du champ de vecteurs f:

71

www.mathonec.com
Remarque 4.1 D’après l’unicité dans le théorème de Cauchy-Lipschitz, les courbes intégrales
d’un champ de vecteurs sont deux à deux disjointes!

4.2 Flot

t0 t1
Les ‡ots locaux et d’une équation di¤érentielle autonome 4.1 se déduisent l’un de l’autre
t0 t1
par translation. En e¤et, les applications t ! (t; ) et t ! (t+t1 t0 ; ) sont par dé…nition
deux solutions de 4.1 et elle prennentla même valeur en t = t0 , donc elle sont égales.
C’est pourquoi on …xe t0 = 0 une fois pour toutes. On notera désormais

0
t( )= (t; ):

4.2.1 Orbites

Si le champ de vecteurs f est complet, le ‡ot est global, c’est à dire que pour tout 2 U,
l’appliquation t ! t( ) est dé…nie sur R tout entier. De plus, l’ensemble

f t 2 C 1 (U ; U ); t 2 Rg

muni de la loi de composition est un groupe commutatif. Ene¤et, il contient IU = 0 et

en utilisant à nouveau l’unicité du théorème de Cauchy-Lipschitz on montre que

s t = t s = t+s

quels que soient s et t. Cette structure algébrique justi…e la dé…nition suivante.

Dé…nition 4.4 Pour tout 2 U onappelle orbite de (pour l’équation di¤ érentielle 4.1 la
courbe intégrale de f passant par , c’est à dire l’ensemble
f t ( ); t 2 Jgoù J est l’intervalle d’existence de au point
Bien sûr, les orbites sont deux à deux disjointes.

72

www.mathonec.com
Remarque 4.2 Si l’on dispose d’une intégrale première E, toute orbite est evidemment incluse
dans un ensemble de niveau de E:

proposition 4.1 L’orbite d’un point est réduite au singleton f g si et seulement si f ( ) = 0:

Si ce n’est pas le cas et s’il existe t0 6= s0 tel que

t0 ( )= s0 ( )

alors t ! t( ) est une solution globale périodique et l’orbite de est une courbe fermée simple

Dé…nition 4.5 On appelle point stationnaire ou point …xe ou point d’équilibre un point 2U
dont l’orbite est réduite à f g; ce qui équivaut à ce que le champ de vecteurs s’annule au point
: on dit aussi que est un point singulier du champ de vecteurs.

Dé…nition 4.6 On appelle cycle une orbite fermée.

Dé…nition 4.7 On appelle orbite hétérocline une courbe intégrale globale f t ( ); t 2 Rgreliant
deux points …xes di¤ érents en +1 et 1.
On appelle orbite homocline une orbite f t ( ); t 2 Rgreliant un même point …xe en +1
et 1:

4.2.2 Redressement du ‡ot

Dé…nition 4.8 On appelle section locale du champ f en un point non singulier y (c’est à
dire tel que f (y) 6= 0), un ouvert S d’une hypersurface contenant y tel que pour tout u 2 S;
f (u) 2
= Tu S (l’espace tangent à S au point u) : on dit que aussi que f (u) est transverse à Tu S.

73

www.mathonec.com
 En particulier, il existe toujours des sections planes, c’est à dire incluses dans un hyperplan
a¢ ne. En e¤ et, si f (y) 6= 0, il existe un hyperplan (vectoriel) fermé H tel que

Rf (y) H = X:

(En dimension in…nie, ceci repose sur le théorème de Hahn-Banach) Par continuité de f ,
l’ensemble fu 2 U ; f (u) 2
= Hg est ouvert. Donc il existe r > 0 tel que pour tout u 2 B(y; r) (la
boule ouverte de centre y et de rayon r); f (u) 2
= H: Par conséquent, l’intersection de l’hyperplan
a¢ ne y +H (dont l’espace tangent en tout point est précisément H) avec B(y; r) est une section
locale en y:

Théorème 4.2 (de redressement du ‡ot)

Si
Rf (y) H = X;

alors il existe > 0 et r > 0 tels que

Sr := y + Hr

avec
Hr := fh 2 H; jjhjjX < rg

soit une section locale en y et l’appliquation

:] ; [ Hr ! ß:= (] ; [ Sr )

(t; h) ! t (y + h);

soit un di¤ éomorphisme. En outre, pour tout 2 ß, il existe un unique t 2] ; [ tel que

t( ) 2 Sr :L’ouvert ßest appelé une boîte à ‡ot.

Démonstration :

74

www.mathonec.com
 ON a déjà vu lea première partie. Il s’agit ensuite d’appliquer le théorème d’inversion locale.
La fonction est de classe C 1 pour et r assez petits. De plus, la di¤érentielle de au point
(0; 0) est donnée par

d (0; 0) (s; k) = sf (y) + k

pour tout (s; k) 2 R H:

C’est une bijection de R H sur X puisque

Rf (y) H = X;

et elle est évidemment continue. Donc sa réciproque est aussi continue d’après le théorème de
l’application ouverte : on peut aussi le voir de façon plus <<élémentaire >> (bien qu’il y ait
derrière le théorème de Hahn-Banach), en écrivant H = ker ' avec ' 2 X 0 tel que

'(f (y)) = jjf (y)jj2 ;

de sorte que

'(x)
sf (y) + k = x , s = f (y):
'(f (y))

Par conséquent, est un di¤éomorphisme local, et donc un di¤éomorphisme de ] ; [ Hr

sur son image, quitte à réduire et r:
Ceci implique que pour tout 2ß
, il existe un unique (t; h) 2] ; [ Hr tel que

t (y + h) = ;

d’où

t( ) = y + h 2 Sr :

Inversement, si s 2] ; [ est tel que

s( ) = y + h 0 2 Sr ;

75

www.mathonec.com
alors
= s (y + h0 )

et donc s = t et h = h0 :

Remarque 4.3 Le di¤ éomorphisme inverse 1 opère un redressement du ‡ot. En e¤ et,

l’image réciproque d’un orbite f t ( ); t 2] ; [g avec

= y + h 2 Sr

n’est autre que le segment de droite f(t; h); t 2] ; [g:

4.3 Portrait de phases

Dé…nition 4.9 On appelle portrait de phases d’un champ de vecteurs f , ou de l’équation dif-
férentielle associée 4.1, la partition de U en orbites.Obtient des informations sur le portrait de
phase ne nécessite pas de phase de savoir résoudre explicitement l’équation di¤ érentielle !Le
point de départ consiste à repérer les orbite stasionnaire, ce qui revient à résoudre l’équation
algébrique
f (v) = 0:

On verra plus loin (chapitre V) comment trouver l’allure des orbites au voisinage des points
stationnaire (où précisément le théorème de redressement du ‡ot ne s’applique pas!)
Il y a d’autre orbites remarquables à rechercher. La recherche de cycles et d’orbites homocline
ou hétérocline est en général non triviale. On verra dans la suite du cours quelques outils pour
cela.

76

www.mathonec.com
4.4

Remarque 4.4 Si l’on dispose d’une intégrale première E, les orbites sont évidemment incluses
dans des ensembles a niveau de E. En dimension …nie n, si l’on dispose de n 1 intégrales
premières indépendantes, les intersections de leurs ensembles ensembles de niveaux forment des
courbes, qui fournissent le portrait de phase complet. (C’est une situation idéale et rare!)
Quoi qu’il en soit, en dimension …nie on obtient des informations intéressantes sur l’allure
des orbites par une simple étude du signe des composantes du champ de vecteurs.
Les portraits de phase en dimension 1 sont essentiellement triviaux : ils reposent seulement
sur le tableau de variations de la fonction f . En dimension supérieure à 3 ils sont di¢ ciles à
représenter et d’une grande complexité ... (On en verra une raison profonde avec le théorème
de poincaré-Bendixson.)
Il est important de remarquer que le portrait de phases dépend en fait seulement de la
direction du champ de vecteurs.

proposition 4.2 Le portrait de phases d’un champ de vecteurs f . coincide avec le portrait de
phases de tout champ de vecteurs de la forme f ,où 2 C 1 (U ; R+ ):

Démenstration :
C’est une simple question de changement de paramétrage. Soit v 2 U et u 2 C 1 (]s ; s+ [)
la solution maximale de

u0 (s) = (u(s))f (u(s)); u(0) = :

Soit alors

Zs
T (s) = (u( ))d ;
0

qui est une fonction strictement croissante puisque est par hypothèse à valeurs positives,
et soit S la fonction réciproque de T , dé…nie sur

77

www.mathonec.com
 ]t ; t+ [; où
Zs
t := (u( ))d :
0

On a

d dS 1
u(S(t)) = u0 (S(t)) = (u(S(t)))f (u(S(t))) = f (u(S(t)))
dt dt (u(S(t)))

et
u(S(0)) = u(0) = :

Donc
U =u S

est solution du problème de Cauchy :

U 0 (t) = f (U (t));

U (0) = :

De plus, U est maximale. Supposons en e¤et que l’on puisse prolonger U à ]t ; t+ +"]; " > 0:
Prolongeons alors S à ]t ; t+ + "] par

Zt
d
S(t) = > 0:
(U ( ))
0

Et on trouve que U T est solution du même problème de Cauchy que u sur l’intervalle
]s ; s+ + ]; contenant strictement ]s ; s+ [: Cela contrdit le fait que u soit maximale. C’est
donc que U est aussi une solution maximale (pour le champ f ). Autrement dit, l’orbite de
sous le champ f , c’est à dire la courbe

fu(s); s 2]s ; s+ [g = fu S(t); t 2]t ; t+ [g

est exactement l’orbite de sous le champ f:

78

www.mathonec.com
Remarque 4.5 Quel que soit le champ de vecteurs f 2 C 1 (U ); il existe 2 C 1 (U ; R+ ) tel
que
g := f

soit complet. En e¤ et, il su¢ t de poser

(u) := 1=(1 + jjf (u)jj2 )

car alors
g= f

est borné et le théorème IV.1 montre que ce champ g est complet.

4.5 Ensemble ! limite

Dans ce qui suit f est un champ de vecteurs de classe C 1 sur l’espace de Banach X tout
entier pour simpli…er.

4.5.1 Propriétés générales

Dé…nition 4.10 Soit 2 X: Si l’intervalle J de dé…nition de la courbe intégrale de f passant

par est tel que sup J = +1, on appelle ensemble ! limite de l’ensemble L! ( ) les valeurs
d’adhérence de f t ( )g lorsque t ! 1:
(Ces noms ont été choisi parce que et ! sont respectivement la première la dernière lettre
de l’alphabet grec).

Exemple 4.3 -Si appartient à une orbite hétérocline, alors par dé…nition il existe et

+ 2 U tels que
lim t( )= :
t! 1

Par conséquent,
L ( ) = L! ( ) = :

79

www.mathonec.com
proposition 4.3 Les ensembles ! limite et limite d’un point sont fermés et invaririants
par le ‡ot. Ils sont de plus inclus dans l’ensemble de niveau f!; E(!) = E( )g si E est une
intégrale première. Plus généralement, si F : X ! R est une fonction continue et monotone le
long des courbes intégrales, les ensembles !-limite et -limite sont inclus dans les ensembles de
niveau de F:

Démonstration:
Ces propriétés sont presque contenues dans la dé…nition. Si ! 2 L! (v) , il existe une suit
tn tendant vers +1 telle que tn (v) tend vers !: Alors, quel que soit t,

t (!) = lim t+tn (v)

n!+1

appartient aussi à L! (v).

Donc cet ensemble est bien invariant par le ‡ot. D’autre part, le complémentaire de
l’ensemble !-limite est par dé…nition l’ensemble des ! tels qu’il existe "0 > 0 et t0 > 0 avec

jj t (v) !jj "0

pour tout t t0 . Or pour un tel !, on a évidement

1
jj t (v) zjj "0
2

pour tout t t0 et pour z 2 B(!; 12 "0 ):

Donc le complémentaire de l’ensemble !-limite est ouvert. Si F : X ! R est continue et
telle que t ! F ( t ( )) soit simplement monotone, alors il existe a 2 [ 1; +1] tel que

lim F ( t ( )) = a:
t!+1

80

www.mathonec.com
 Si y 2 L! ( ), il existe une suite (tn )n2N tendant vers +1 telle que

lim tn ( ) = y:
n!+1

Par continuité de F , on en déduit

a = lim F ( tn ( )) = F (y)
n!+1

(en particulier, jaj < 1). Ceci prouve que L! ( ) fy; F (y) = ag:

proposition 4.4 Soit 2 X: Si l’orbite { (v)} est dé…nie pour tout t 2 R+ et reletivement
compacte, alors L! (v) est compact non vide et connexe.

Démonstration:
La compacité de L! (v) est quasi-immédiate: c’est un fermé inclus dans le compact

f t (v); t 2 R+ g

:
Il est non vide puisque l’orbite f t (v); t 2 R+ g admet au moins une valeur d’adhérence. La
connexite de L! (v) demande un peu plus de travail. Raisonnons par l’absurde et supposons
qu’il existe deux ouverts disjoints U1 et U2 tel que

L! (v) U1 [ U2 ; L! (v) \ Ui 6= ;

Par dé…nition de L! (v), pour tout n 2 N, il existe donc une suite ( n) strictement croissante
telle que

2n
(v) 2 U2 ; 2n+1
(v) 2 U1 :

comme l’enssemble f t ( ); 2n t 2n+1 g est connexe (image d’un connexe par une appli-
cation continue !), il ne peut être inclus dans U1 [ U2 : Donc il existe tn 2] 2n ; 2n+1 [ tel que

tn ( )2
= U1 [ U2 : Or la suite tn ( ) admet une valeur d’adhérece a 2 L! ( ); et cette limite
appartient au

81

www.mathonec.com
 complémentaire de U1 [ U2 (car c’est un fermé). Ceci contredit l’inclusion L! ( ) U1 [ U2 :
C’est donc que de tels ouverts U1 et U2 n’existent pas, ce qui montre que L! ( ) est connexe.

Remarque 4.6 Inversement, si L! ( ) est compact et non vide, l’orbite f t ( ); t 2 R+ g

est relativement compacte, puisqu’elle admet au moins une valeur d’adhérence, et toutes ses
valeurs d’adhérence sont inclues dans le compact L! ( ):

4.5.2 Equations dans le plan

Dans R2 , les sections locales "planes" sont des segments de droite, en particulier orientables.
De plus, on dispose du

Théorème 4.3 (de Jordan)

Une courbe fermée simple divise le plan en deux régions connexes, l’une bornée et l’autre
non bornée.

Grâce à ces deux ingrédients, on peut démontrer le

Théorème 4.4 Pour un champ de vecteurs dans R2 , tout ensemble !-limite compact non vide
et sans point …xe est un cycle.
Il existe d’autres versions de ce théorème (sans l’hypothèse <<sans point …xe >> notam-
ment). La démonstration n’est pas di¢ cile mais un peu longue. Elle utilise les résultats prélim-
inaires que voici.

Lemme 4.1 Si S est une section locale, si (tn ) est une suite croissante telle que

yn = tn ( ) 2 S;

alors la suite (yn ) est aussi monotone le long de S, c’est à dire

det(yn+1 yn ; y n yn 1) 0:

82

www.mathonec.com
 Démonstration :
Il su¢ t de faire la preuve pour trois points y0 ; y1 ; y2 : Fixons donc

y0 = t0 ( ) et y1 = t1 ( )

appartenant à S avec t0 < t1 et y0 6= y1 : On peut supposer que l’orbite de ne coupe pas S

dans l’intervalle I := [y0 ; y1 ] : si tel était le cas, on remplacerait y0 par

ye0 = t0 (y1 )

le premier point appartenant à S en partant de y1 dans le sens des t < t1 , car montrer que
y0 ; y1 et
y2 = t( ); t > t1

sont ordonnés revient à montrer que ye0 ; y1 et y2 le sont.

Soit C la courbe formée de la réunion de I et de f s( ); s 2 [t0 ; t1 ]g: C’est une courbe
fermée, simple grâce à la précaution que l’on a prise et au fait qu’une orbite ne s’intersecte pas
elle-même. Donc d’après le théorème IV.3, elle divise le plan en deux composantes connexes.
Soit D celle qui est bornée. Puisque S est une section, le champ f pointe vers le même
demi-plan le long de S,et en particulier le long de I: En particulier, il pointe soit vers D soit
vers R2 nD: Quitte à renverser le temps (et échanger y0 et y1 ), on peut supposer qu’on est dans
le second cas. Alors D est négativement invariant par le ‡ot, c’est à dire que tout y 2 D,

t (y) 2 D pour t < 0. En e¤et, le ‡ot ne peut sortir ni par I ni par { t ( ); t 2 [t0 ; t1 ]g: Ceci
implique en particulier que t (y1 ) appatient à R2 nD pour tout t > t1 : De plus, SnI est constitué
de deux intervalles I0 et I1 contenant respectivement y0 et y1 dans leur bord. Or on peut joindre
tout point de I0 assez proche de y0 à " (y0 ) avec " > 0 assez petit, point appartenant à D,
sans passer par le bord de D, on en déduit que I0 est dans D. Donc, si y2 = t2 ( ) 2 S et
t > t2 ; y2 appartient nécessairement à I1 :

Corollaire 4.1 Si S est une section locale, pour tout 2 U; L! ( ) \ S contient au plus un

83

www.mathonec.com
point.

Démonstration :
Supposons que L! ( ) \ S contienne deux points distincts, y1 et y2 : Soient alors des boîtes
à ‡ot disjointes, V1 3 y1 et V2 3 y2 :Comme ces points sont dans L! (z); l’orbite z repasse une
in…nité de fois dans chacune de ces boîtes, et donc aussi par chacun des intervalles I1 = V1 S
et I2 = V2 S. Plus précisément, il existe une suite (tn ), croissante et tendant vers +1 avec n,
telle que t2n+1 ( ) 2 I1 et t2n ( ) 2 I2 : Comme I1 et I2 sont disjoints, cesi contredit le lemme
IV.1.

84

www.mathonec.com
Chapitre 5

STABILITE DES SOLUTIONS

STATIONNAIRES

On va s’intéresser ici à deux types de stabilité : l’une concernant le comportement asymptotique

des
solutions pour les données initiales proches d’une solution stationnaire ; l’autre concernant
le portrait de phase autour d’un point …xe lorsqu’on perturbe l’équation (en faisant varier
des paramètres par exemple). Dans le premier cas, on parle de stabilité asymptotique (ou de
stabilité tout court), que l’on peut étudier grâce à la théorie de Lyapunov. Dans le second
cas, on parle de stabilité structurelle, ce qui amène aux notions de variétés stable / instable /
centrale et à des problèmes de bifurcations.

5.1 Théorie de Lyapunov

Dé…nition 5.1 Un point singulier d’un champ de vecteurs f 2 C 1 (U ) est dit stable s’il existe
un voisinage V0 de dans U tel que
i): le ‡ot t (!) est dé…ni pour tout t 0 et pour tout ! 2 V0 ;
ii):pour tout voisinage W de dans U , il existe un voisinage V V0 tel que t (!) appartient
à W pour tout t 0 et pour tout ! 2 V:

85

www.mathonec.com
 Si de plus, on a
iii):pour tout ! 2 V0 ;

lim t (!) =
t!+1

alors le point …xe est dite asymptotiquement stable.

On dit d’unpoint non stable qu’il est instable.
Concernant la stabilité des équilibres pour les équations linéaires, la proposition III.5 se
reformule ainsi.

Théorème 5.1 Soit A 2 Mn (R): Le point 0 est stable pour l’équation

du
= Au
dt

si et seulement si
i):les valeurs propres de A sont toutes de partie réelle négative ou nulle,
ii):les valeurs propres imaginaires pures de A sont semi stable
Le point d’équilibre est asymptotiquement stable si et seulement si
iii): les valeurs propres de A sont toutes de parties réelle strictement négative.
Pour les équations non-linéaire, la stabilité des équilibres est plus délicate à étudier (en
dimension …nie et a fortiori en dimension in…nie). La théorie de Lyapunov repose sur l’existece
de fonctions ( de Lyapunov !) qui permettent de <<contrôler >> , dans une certaine mesure,
le comportement asymptotique des solutions pour des données initiales proches de l’équilibre.

Théorème 5.2 (Lyapunov n 1)

Soit un point d’équilibre du champ f 2 C 1 (U ); U étant un ouvert dans un espace de
Banach X: On suppose qu’il existe un voisinage V de dans U et une fonction F 2 C 2 (V ),
telle que
le point est un minimum local de F , d’où nécessairement

DF ( ) = 0

86

www.mathonec.com
et de plus il existe > 0 tel que
D2 F ( ) Id;

et la fonction F est décroissante le long des orbites issues de V , c’est à dire

DF (u) f (u) 0 8u 2 V:

Alors est un point d’équilibre stable.

Démonstration :
Soit 0 > 0 tel que f soit bornée, disons par M , et lipschitzienne dans la boule B( ; 2 0 )
V . La démonstration du théorème de Cauchy-Lipschitz montre que le ‡ot t (!) est dé…ni pour
tout ! 2 B( ; 0) et t 2 [0; ] avec := 2M :
0

D’après l’hypothèse sur F et la formule de Taylor à l’ordre 2, il existe 0 tel que

F ( + h) F( ) jjhjj2
4

pour tout h 2 X avec jjhjj : Par suite,

2
inf F (u) F( ) + = m:
jju jj= 4

Par continuité de F;

Vn := fu 2 B( ; ); F (u) < mg

est un voisinage (ouvert) de dans U: De plus, si ! 2 Vn ;

F ( t (!)) F (!) < m

pour t 2 [0; ] et donc jj t (!) jj reste strictement inférieure à 0:

En particulier, (!) appartient à Vn B( ; 0 ): Donc on peut prolonger t (!) à t 2 [ ; 2 ];

et de proche en proche à tout intervalle [k ; (k + 1) ](k 2 N);tant en restant dans Vn . en
conclusion, pour tout ! 2 Vn ; t (!) est dé…ni pour tout t 0 et à valeurs dans Vn :

87

www.mathonec.com
Dé…nition 5.2 La fonction F du théorème V.2 est appelée une fonction de Lyapunov.

Exemple 5.1 Une intégrale première ayant un minimum locale (au sens stricte de l’énoncé du
théorème V.2) en est une fonction de Lyapunov !
Dans les modèles physiques, les fonctions de Lyapunov sont souvent reliée à une énergie.

Remarque 5.1 En dimension …nie, on peut a¤ aiblir l’hypothèse sur F , en demandant que
soit un minimum local strict sans nécessairement avoir d’estimation de D2 F ( ):
Si l’on dispose d’une fonction de Lyapunov dotée d’une propriété supplémentaire adéquate,
on peut obtenir la stabilté asymptotique de l’équilibre.

Théorème 5.3 (Lyapunov n 2)

Soit unpoint d’équilibre du champ f 2 C 1 (U ); U étant un ouvert dans un espace de Banach
X. On suppose qu’il existe un voisinage V de dans U et une fonction de Lyapunov F 2 C 2 (V )
telle que
il existe > 0 tel que

DF (u) f (u) (F (u) F ( )) 8u 2 V:

Alors est un point d’équilibre asymptotiquement stable.

Démonstration :
On sait déjà d’après la preuve du théorème V.2 qu’il existe un voisinage W de inclus dans
V tel que

t (W ) W

pour tout t 0. L’hypothèse renforcée sur F permet facilement de voir que

lim t (!) =
t!+1

88

www.mathonec.com
 pour tout ! 2 W . En e¤et,

d
(F ( t (!)) F ( )) = DF ( t (!)) f ( t (!)) (F ( t (!)) F ( ))
dt

implique par une intégration immédiate :

t
(F ( t (!)) F ( )) (F (!) F ( ))e :

Or

F ( t (!)) F( ) jj t (!) jj2

4

par construction de W:

Dé…nition 5.3 La fonction F du théorème V.3 est applée une fonction de Lyapunov forte.

Remarque 5.2 i): Bien sûr, une intégrale première n’est pas une fonction de Lyapunov forte.
Cependant, on peut parfois trouver une fonction de Lyapunov forte comme intégrale première
d’un système "approché". C’est le cas par exemple pour l’équation du pendule avec amortisse-
ment:

d2 x dx
+k + sin x = 0; k > 0:
dt2 dt

Soit E l’intégrale première connue pour k = 0 :

1
E(x; x0 ) = (x0 )2 cos x:
2

Alors pour une solution du système amorti on a :

d
E(x(t); x0 (t)) = x0 (t)(x00 (t) + sin x(t)) = k(x0 (t))2 < 0
dt

pour x0 (t) 6= 0: La fonction E est donc une fonction de Lyapunov aux points de la forme
(2k ; 0) ; ce n’est pas tout à fait une fonction de Lyapunov forte, mais elle parmet de montrer
que ces points sont asymptotiquement stables.

89

www.mathonec.com
 ii): Des systèmes admettant une fonction de Lyapunov forte évidente sont les systèmes
gradient, de la forme

du
= grad V (u):
dt

Si est un minimum local strict de V alors V est une fonction de Lyapunov forte au point
:
Un résultat remarquable de la théorie de Lyapunov, en dimension …nie, est l’équivalence
de la stabilité asymptotique pour l’équation non linéaire et pour sa version linéarisée au point
d’équilibre. Il peut s’énoncer ainsi.

Théorème 5.4 (Lyapunov n 3) [3]

Soit un point d’équilibre du champ f 2 C 1 (U ); U Rn :
On suppose que la matrice jaccobienne Df ( ) a toutes ces valeurs propres de partie réelle
strictement négative.
Alors est asymptotiqument stable.
On dit aussi en abrégé que la stabilité spectrale (c’est à dire que le spectre du linéarisé est
inclus dans f 2 C; Re < 0g) implique la stabilité non linéaire. C’est un résultat spéci…que
de la dimension …nie, qui équivaut d’après le théorème V.1 à l’assertion suivante. Si est
un point d’équilibre de f et si 0 est un point d’équilibre asymptotiquement stable du système
linéarisé autour de ce point d’équilibre :

du
= Df ( )u
dt

alors est un point d’équilibre asymptotiquement stable du système non linéaire.

Pour démontrer le théorème V.4, on aura besoin du

Lemme 5.1 Soit A 2 Mn (R) dont les valeurs propres sont toutes de partie réelle strictement
négative. Alors l’équation

du
= Au
dt

90

www.mathonec.com
 admet une fonction de lyapunov forte qui est une forme quadratique;

Démonstration :
On commence par "diagonaliser A à " près", c’est à dire à trouver, pour tout " > 0,
une matrice P" 2 GLn (C) et une matrice C" 2 Mn (C) de norme (subordonnée à la norme
hermitienne sur Cn ) inférieure à " telle que

P" 1 AP" = D + C"

avec D diagonale de coe¢ cients { 1 ; :::; ng les valeurs pripres de A. Pour cela, il su¢ t
de trigonaliser A de façon particulière : si a est l’endomorphisme de matrice A dans la base
canonique de Cn , on fabrique une base (e1 ; :::; en ) dans laquelle la matrice de a soit une matrice
triangulaire supérieure D + B, c’est à dire que pour tout j 2 f1; :::; ng :

j 1
X
a(ej ) = j ej + bij ej ;
i=1

puis on dilate les vecteurs de base en eej = je ; de sorte que

j

j 1
X
j i
a(e
ej ) ej
je = bij eei :
i=1

Autrement dit, on a
1
P AP = D + C

avec P la matrice de vecteurs colonnes ee1 ; :::; een et C la matrice triangulaire supérieure de
coe¢ cients cij := bij j i pour i j 1 (et cij = 0 pour i j). Cette matrice est de norme
inférieure à " pourvu que < min(1; " =jjBjj):
Soit alors

q(u) = u Q Qu = jjQujj2 ; Q = P" 1 :

On a

91

www.mathonec.com
 q(0) = 0; Dq(u) h = u Q Qh + h Q Qu; D2 q(u) = 2Q Q:

Comme Q Q est hermitienne dé…nie positive, la fonction q admet bien un minimum local
strict en 0 (on a la minoration évidente Q Q In avec = minjjhjj=1 jjQhjj2 ). De plus,

Dq(u)Au = 2 Re(u Q QAu) = 2 Re(u Q DQu) + 2 Re(u Q CQu):

D’après l’inégalité de Cauchy-Schwartz, on a

Re(u Q CQu) jjQujjjjCQujj "jjQujj2

par construction de C. Donc, si = maxj Re j; on a

Dq(u)Au 2( ")q(u);

ce qui implique que q est une fonction de Lyapunov forte pourvu qu’on ait choisi " < :
Démonstration [Théorème V.4]
On peut appliquer le lemme V.1 à la matrice jaccobienne A = Df ( ):
On a ainsi une forme quadratique

q(u) = jjQujj2 ;

telle que

q(0) = 0; Dq(0) = 0; Dq(0) In

et

Dq(u)Au q(u)

avec et > 0: Alors

F (!) = q(! )

92

www.mathonec.com
dé…nit une fonction de Lyapunov forte pour le champ f ; en e¤et, on a de façon évidente :

F ( ) = 0; Df ( )(! ) + B( ; !) (! ;! )

où B( ; !) est l’application bilinéaire dé…nie par :

Z1
1
B( ; !) = (1 )D2 f ( + (! ))d :
2
0

En particulier, il existe c > 0 tel que pour jj! jj 1:

2c
jjB( ; !) (! ;! )jj cjj! jj2 jjQ(! )jj2 ;

Par suite

DF (!)f (!) = 2 Re((! ) Q Qf (!))

= Dq(! )A(! ) + 2 Re((! ) Q QB( ; !) (! ;! ))

d’où par l’inégalité de Cauchy-Schwartz

4c 4c
DF (!)f (!) F (!) + jjQjjjjQ(! )jj3 ( jjQjj2 jj! jj)F (!):

Pour
jj! jj
8cjjQjj2

on a donc

DF (!)f (!) F (!):

2

Donc F est une fonction de Lyapunov forte pour f au point , et l’équilibre est asymp-
totiquement stable d’après le théorème V.3.
Pour un "gros" système, il peut parfois être plus facile de trouver une fonction de Lyapunov

93

www.mathonec.com
que de montrer que le système linéarisé a ses valeurs propres sont de partie réelle négative. En
dimension in…nie, le spectre du linéarisé est encore plus di¢ cile à localiser et il ne su¢ t pas à
prouver la stabilité. La recherche de fonctions de Lyapunov est donc très importante.

5.2 Points …xes hyperboliques

La classi…cation des points …xes pour les équations di¤érentielles linéaires dans le plan montre
que les centres (correspondant à une matrice A 2M2 (R) avec des valeurs propres imaginaires
pures) jouent un rôle particulier, au sens où ils ne sont pas structurellement stables : une
modi…cation arbitrairement petite de la matrice

0 1 0 1
0
@ A en @ A
0

avec j j << 1 modi…e dramatiquement le portrait de phase au voisinage de 0! Ces points

sont les prototypes de points …xes non hyperboliques.

Dé…nition 5.4 Un point singulier du champ f est dite hyperbolique si le spectre de l’opérateur
Df ( ) n’intersecte pas l’axe imaginaire.
Les points …xes hyperboliques sont précisément ceux qui sont structurellement stables, au
sens topologique du terme.

C’est l’objet du théorème suivant, que nous admettrons;

Théorème 5.5 (Hartman-Grobman)

Si est un point …xe hyperbolique de l’équation di¤ érentielle

du
= f (u);
dt

avec f 2 C 1 (U ); U un ouvert de Rn ; il existe un voisinage V de dans U , un voisinage O

de 0 dans Rn et un homéomorphisme h de V sur O tel que le ‡ot t (!) véri…e :

h( t (!)) = etDf ( ) h(!)

94

www.mathonec.com
 pour tout ! 2 V et t 2 R tel que t (!) soit bien dé…ni.
Ce théorème signi…e que l’on peut déformer continûment les orbites au voisinage de pour
obtenir les orbites du système linéarisé

du
= Df ( )(u );
dt

et réciproquement. Il ne dit rien sur les propriétés de nature di¤ érentielle (tangentes, con-
vexité, etc.) de ces courbes. En particulier, il ne fait pas la distinction entre les di¤ érentes
sortes noeuds en dimension 2.
Le théorème dit de la variété stable si-après ne fait pas non plus de distinction entre ces
points, mais il donne des informations sur les propriétés di¤ érentielles des orbites, souvent
utiles dans les applications.

Théorème 5.6 (de la variété stable locale)

Si est un point …xe hyperbolique de l’équation di¤ érentielle

du
= f (u)
dt

(avec f 2 C 1 (U ); U un ouvert de Rn ); dont le ‡ot est noté t (!), il existe un voisinage V

de dans U tel que les ensembles

s
Wloc = f! 2 V ; t ! t (!) 2 Cb1 (R+ )g;

u
Wloc = f! 2 V ; t ! t (!) 2 Cb1 (R )g

soient des variétés (à bord ) contenant , tangentes respectivement à E s et E u , les sous

espaces stable et instable de Df ( ), au point s et W u sont des graphes
: autrement dit, Wloc loc

de fonctions di¤ érentiables sur des voisinages de dans les sous espaces a¢ nes +E s et +E u
respectivement, dont la di¤ érentielle s’annule en .
s est positivement invariante par le ‡ot, c’est à dire
De plus, Wloc s s pour tout
t (Wloc ) Wloc
t 0, tandis que u u pour tout t
t (Wloc ) Wloc 0, et on a

95

www.mathonec.com
 s u
Wloc = f! 2 V ; lim t (!) = g; Wloc = f! 2 V ; lim t (!) = g:
t!+1 t! 1

s et W u sont applés respectivement variété stable locale et

Dé…nition 5.5 Les ensembles Wloc loc

variété instable locale.

Nous allons voir une démonstration du théorème de la variété stable (locale) à l’aide de la
méthode de Lyapunov-Schmidt.

Méthode de Lyapunov-Schmidt
Cette méthode permet de résoudre un problème non linéaire du type

F (x; ) = 0

lorsque le théorème des fonctions implicites ne s’applique pas, c’est à dire lorsque la dif-
férentielle de F par rapport à x n’est pas inversible. Plus précisément, considérons F de la
forme

F (x; ) = Bx G(x; )

pour (x; ) 2 X Rp , où X est un espace de Banach, B est un opérateur cntinu dee X dans
un autre espace de Banach Z et G 2 C 1 (X Rp ; Z) est telle que

G(0; 0) = 0; @x G(0; 0) = 0

La méthode de Lyapunov-Schmidt est fondée sur une conséquence du

Théorème 5.7 (de l’application ouverte)

Si B est une application linéaire continue surjective d’un espace de Banach X dans un
espace de Banach Y , alors B est ouverte, c’est à dire pour tout ouvert de X; B( ) est un
ouvert de Y: Par suite, si B est bijective, son inverse B 1 est continu.

96

www.mathonec.com
5.3 Variétés invariantes

Plus généralement, on peut trouver des variétés invariantes par le ‡ot d’une équation di¤éren-
tielle en séparant de façon arbitraire le spectre de l’équation linéarisé : ceci ne demande plus
d’hyperbolicité du point …xe !
De plus, on peut construire des << variétés globales >>, dé…nies comme des graphes de
fonctions dé…nies sur des sous espaces entiers (invariants par l’équation linéarisée), à condition
que la partie non linéaire du champ de vecteurs soit globalement assez petite.
C’est l’objet du résultat suivant, que l’on admettera mais dont on montrera comment il
implique en particulier le théorème V.6.

Théorème 5.8 Soit f 2 C 1 (Rn ); 2 Rn tel que

f ( ) = 0;

et E ; E + des sous espaces supplémentaires dans Rn , invariants par

A := Df ( )

tels que

(AjE ) f ; Re g; (AjE + ) f ; Re g;

où < (sans rien supposer a priori sur leur signe). En notant

g(x) = f (x) Df ( ) (x );

si jjgjjC 1 (Rn ) est assez petit, il existe un unique h 2 C 1 (E + ; E ) tel que

b

h(0) = 0; Dh(0) = 0; sup jjjDh( )jjj < +1;

2E +

et le graphe de h :

97

www.mathonec.com
 W + := f + + h( ); 2 E + g

soit invariant par le ‡ot de l’équation di¤ érentielle

u0 = f (u):

De plus, pour tout

!= + + h( ) 2 W +

et pour tout 0 > il existe C > 0 tel que

0t
jj t (!) jj Ce jj jj:

Dans cet énoncé, si les deux paramètres et sont strictement positifs, cela implique en
particulier que le point …xe est hyperbolique, et la variété obtenue est précisément la variété
instable.
Voyons comment en déduire précisément le théorème V.6. Pour …xer les idées, disons qu’il
s’applique pour jjgjjC 1 (Rn ) ". En général, on ne peut espérer cette inégalité globale. Cependant
b

on peut s’y ramener, quitte à modi…er f , et donc aussi g en dehors d’un voisinage de . En
e¤ et, il existe une boule B( ; ) dans Rn telle que

sup jjjDg(x)jjj "=3; d0ou sup jjg(x)jj "=3

x2B( ; ) x2B( ; )

par le théorème des acroissements …nis. Soit alors 2 C 1 (Rn ; R+ ); valant identiquement
1 dans B( ; =3); identiquement 0 en dehors de B( ; ) telle que

sup jj (x)jj = 1; sup jjjD (x)jjj 2= :

x2Rn x2Rn

(L’existence d’un tel , appelé fonction de troncature est laissée en exercice.)

Considérons fe dé…nie par

fe(x) = Df ( ) (x ) + ge(x); ge := g:

98

www.mathonec.com
 Evidemment on a fe = f dans B( ; =3); et donc en particulier

Dfe( ) = Df ( ):

De plus

2 " "
jjjDe
g (x)jjj jjjD (x)jjj sup jjg(x)jj + jjjDg(x)jjj + ="
x2B( ; ) 3 3

Donc le théorème V.8 s’applique à l’équation di¤ érentielle

u0 = fe(u);

dont le ‡ot coincide avec celui de

u0 = f (u)

dans B( ; =3):

Dé…nition 5.6 Lorsque = 0 dans l’énoncé du théorème V.8, la variété obtenue W + est
appelée variété centrale-instable (au point ) et généralement notée W cu : On obtient la variété
centrale-stable W cs par renversement du temps (t ! t): On dé…nit alors la variété centrale
(au point ) par
W c := W cu \ W cs :

Quite à diminuer " dans l’énoncé du théorème V.8, on peut démontrer que W c est le graphe
d’une fonction de classe C 1 au dessus de

E c := E cu \ E cs

où E cu et E cs désignent naturellement les espaces tangents respectivement à W cu et W cs au

point . D’après le théorème V.8, si ! 2 Wc , pour tout > 0 il existe C > 0 tel que

8t 2 R; jj t (!) jj Ce t :

99

www.mathonec.com
 autrement dit, les solutions du problème de Cauchy avec donnée initiale dans la variété
centrale Wc ne peuvent croître que de façon modérée (tout en restant à valeurs dans Wc du fait
de son invariance par le ‡ot).
Attension, on rappelle que le théorème V.8 ne s’applique qu’aux champs de vecteurs dont
1 n
la partie non linéaire est assez petite dans b (R ). En fait, on peut démontrer (notamment à
l’aide d’une troncature comme ci-dessus), l’existence d’une variété centrale locale.

Théorème 5.9 (de la variété centrale locale)

Si eest un point …xe de l’équation di¤ érentielle

du
= f (u)
dt

(avec f 2 C 1 (U ); U un ouvert de Rn ), dont le ‡ot est noté t (!), si E c , le sous espace

centrale de Df ( ), c’est à dire la somme des sous espaces caractéristiques associés aux valeurs
propres imaginaires pures de Df ( ), est non triviale, il existe un voisinage V de dans U et
c , tangente à E c au point , c’est à dire que
une variété locale Wloc

c
Wloc = f! = + + h( ); 2 E c g;

où est un voisinage de 0 dans E c et h : ! Es E u (somme des sous espaces carac-

téristiques associés au valeurs propres de parties réelles respectivement strictement négatives et
strictement positives) est telle que

h(0) = 0; dh(0) = 0;

véri…ant les deux propriétés suivantes :

c ; 8t 2 R;
i): Invariance par le ‡ot : 8! 2 Wloc

c
t (!) 2V ) t (!) 2 Wloc :

ii): Lieu des orbites globales proches de : 8! 2 V; si t (!) et dé…ni et à valeurs dans V;
c
8t 2 R; alors ! 2 Wloc

100

www.mathonec.com
 c : peut exister plusieurs
Attension, ce théorèmee ne prétend pas montrer l’unicité de Wloc
variétés centrales locales, tangentes à E c et véri…ant les propriétés i) ii). Un exemple très
simple est celui du système
x0 = x2 ;

y0 = y;

dont on peut calculer explicitement les solutions et constater que les orbites sont incluses
dans les courbes ayant une équation de la forme

y = C exp(1=x);

où C est une constante arbitraire. Le portrait de phase montre ansi une in…nité de courbes
tangentes à l’axe Ox et véri…ant les propriétés i) ii) : ce sont précisément les graphes de toutes
les fonctions de la forme

0 x 0;
x!
C exp(1=x) x<0 ,

5.4 Introduction aux bifurcations

On va s’intéresser ici aux bifurcations dites de co-dimension 1, c’est à dire concernant des
champs de vecteurs dépendant d’un paramètre scalaire et dont le portrait de phase (local)
change brutalement lorsque le paramètre traverse une valeur critique. D’après le théorème de
Hartman-Grobman, ceci ne peut se produire qu’au voisinage d’un point …xe non hyperbolique.
Nous allons donc considérer

f :U Rn R ! Rn

101

www.mathonec.com
 (u; ) ! f (u; );

de classe 2 , tel qu’il existe (u0 ; 0) 2U pour lequel

f (u0 ; 0) =0

et le spectre de Df ((u0 ; 0 )(où D désigne la di¤érentielle partielle par rapport à u) intersecte

l’axe imaginaire.
Plus précisément, nous allons nous placer successivemant dans les cas les plus simples :
i): La seule valeur propre imaginaire pure de Df ( (u0 ; 0) est 0;
ii):La jaccobienne Df ((u0 ; 0) admet exactement deux valeurs propres imaginaires pures
(et donc forcément conjuguées l’une à l’autre).
Sous les hypothèses ‹‹génériques››, le premier cas correspond à la bifurcation élémentaire
dite col-noeud, et le second à une bifurcation de Hopf.

5.4.1 Bifurcation col-noeud

Si Df (u0 ; 0) a seulement 0 comme valeur propre imaginaire pure, alors le système augmenté
:

du
= f (u; ); (5.1)
dt

d
= 0;
dt

:admet une variété centrale locale en (u0 ; 0) de dimension 2; Le système réduit sur cette
variété est de la forme

102

www.mathonec.com
 dx
= g(x; ); (5.2)
dt

d
= 0;
dt

où l’inconnue x est scalaire,

g(0; 0) =0

et
@x g(0; 0) = 0:

Supposons en outre que

@ g(0; 0) 6= 0; @x g((0; 0) 6= 0:

Pour …xer les idées, on supposera ces deux quantités positives. Dans ce cas, pour < 0, la
fonction g( ; ) ne s’annule pas au voisinage de 0, pour = 0 elle s’annule en deux points de
signes opposés. Plus précisément, au voisinage de (0; 0) dans le plan f(0; 0) 2 R2 g, l’ensemble
des zéros de g est une courbe, tangente à l’axe Ox en (0; 0 ), avec la concavité tournée vers
> 0. le portrait de phase du système 5.2 est particulièrement simple à tracer, car sur les
droites d’équation =constante, il se ramène au portrait de phase avec le paramètre ( ) en
abscisse, voir la …gure V.1.
Le nom de cette bifurcation provient du portrait de phase de l’équation scalaire

dx
= g(x; )
dt

augmentée d’une équation hyperbolique, par exemple

dy
= y:
dt

Le portrait de phase du système plan ainsi obtenu est qualitativement le même que celui
du système modèle :

103

www.mathonec.com
 dx
= x2 + 0
dt

dy
= y;
dt

5.4.2 Bifurcation de Hopf

On revient maintenant à un champ f tel que

f (u0 ; 0) =0

et l’on suppose que Df ((u0 ; 0) admet exactement deux valeurs propres imaginaires pures, i!
avec ! 6= 0. Ainsi, 0 n’est pas valeur propre, et donc Df (u0 ; 0) est inversible. Le théorème des
fonctions implicites montre donc que l’ensemble des zéros de f voisins de (u0 ; 0) est paramétré
par : il existe une fonction u ; dé…nie et de classe 2 au voisinage de 0; telle que

u ( 0) = u0

et

f( ; ) = 0 , = u ( );

A nouveau le système augmenté 5.1 admet une variété centrale locale en (u0 ; 0 ). Elle est
ici de dimension 3, et contient nécessairement la courbe de points …xes f(u ( ); )g: le système
réduit sur cette variété se ramène à un système plan dépendant du paramètre . C’est pourquoi
on suppose dans ce qui suit n = 2

Théorème 5.10 (bifurcation de Hopf )

On suppose que Df (u ( ); ) admet deux valeurs propres complexes conjuguées valant i!
en 0 et dont la partie réelle change de signe en 0. Alors trois cas sont possibles au voisinage
de 0 :

104

www.mathonec.com
 i); le point …xe u ( ) est un centre : le portrait de phase du champ de vecteurs (non linéaire
!) f ( ; ) est constitué de cycles concentriques;
ii):on a une bifurcation surcritique : pour < 0, le champ de vecteurs f ( ; ) n’admet
aucun cycle autour de u ( ), pour > 0,
il admet exactement un cycle autour (et assez
p
proche) de u ( ), le diamètre de ce cycle est en O( 0 ):

iii): on a une bifurcation souscritique : pour > 0, le champ de vecteurs f ( ; ) n’admet

aucun cycle autour de u ( ), et pour < 0,
il admet exactement un cycle autour (et assez
p
proche) de u ( ); le diamètre de ce cycle est en O( 0 ):

La démonstration de ce théorème repose sur la théorie des formes normales, qui sort du
cadre de ce cours. Nous nous contenterons d’observer ce qui se passe sur la forme réduite de
Poincaré-Andronov :

x0 = y+ x(x2 + y 2 );

y0 = x + y(x2 + y 2 ):

La valeur critique du paramètre est ici évidemment = 0: En passant en coordonnées

polaires, ce système se ramène à

r0 = r r3 ;

0
= 1;

c’est à dire …nilement à l’équation scalaire

dr
= r r3 :
d
p
Si < 0 cette équation a un seul point …xe, 0, et si > 0; elle a trois : 0 et : attention,
p p p
et ne sont pas des points …xes du système du départ, mais est précisément le
rayon d’un cycle.

105

www.mathonec.com
 106

www.mathonec.com
Bibliographie

[1] V.Arnold. Equations differentielles ordinaires. Edition Mir,Moscow, 1974.

Champs de vecteurs, groupes à un paramètre, di¤ éomorphismes,‡ots, systèmes linéaires,
stabilités des positions d’équilibre, théorie des oscillations, équations di¤ érentielles sur les
variétés. Traduit du russe par Djilali Embarek.

[2] V.Arnold. Chapitres supplementaires de la theorie des equations differen-

tielles ordinaires. "Mir", Moscow,1984. Translated from the russian by Djilali Embarek,
reprint of the 1980 edition.

[3] Michel Crouzeix and Alain L.Mignot. Analyse numerique des equations dif-
ferentielles.Collection mathématique appliquées pour la maitrise. Masson, Paris, 1984.

107

www.mathonec.com
 108

www.mathonec.com