Voir la source de Flot (mathématiques)

{{autre4|le terme mathématique||Flot (homonymie)}}

Le '''flot''', '''coulée''' ou encore '''courant''' est, en [[mathématiques]], un concept fondamental utilisé en [[géométrie différentielle]]. 

La notion de flot permet notamment de modéliser le déplacement dans le temps des éléments d'un fluide. Pour ce faire, on crée une application α qui, à  chaque point ''x'' de l'espace concerné par l'écoulement, associe un autre point α(''x'',''t''), correspondant à la position qu'aurait une particule du fluide à l'instant ''t'', si elle avait été située en ''x'' à l'instant 0.

[[Fichier:PenduleEspaceDesPhases.png|thumb|Flot associé à l'équation différentielle d'un pendule. Les abscisses représentent la position et les ordonnées la vitesse.]]
[[File:Line-Integral.gif|thumb|Une courbe intégrale d'un champ de vecteurs.]]

Le flot est associé à la notion de [[champ de vecteurs]], c'est-à-dire à une application ''f'', qui, à un point ''x'' d'un [[Ouvert (topologie)|ouvert]] Ω d'un [[espace de Banach]] ''E'', associe un vecteur de ''E''. Un tel champ définit une [[équation différentielle]] du type  α'(''t'') = ''f''(α(''t'')), (c'est-à-dire, dans l'exemple du fluide, que la fonction ''f'' associe au point x à l'instant t, ses direction et vitesse de déplacement). Si la fonction ''f'' est [[Propriété locale|localement]] [[fonction lipschitzienne|lipschitzienne]], pour chaque point ''x'' de Ω, il existe une [[Théorème de Cauchy-Lipschitz#Vocabulaire spécifique|solution maximale]] α{{ind|''x''}} du [[problème de Cauchy]] constitué de cette équation différentielle et de la condition dite de Cauchy α{{ind|''x''}}(0) = ''x''. Vue comme une [[Fonction de plusieurs variables|fonction de deux variables]], ''t'' et ''x'', l'application α est appelée le flot du champ ''f'' de vecteurs. Cette définition se généralise dans le cas d'un champ de vecteurs temporel (c'est-à-dire dépendant d'une variable ''t'' qui prend ses valeurs dans '''R''') et dépendant d'un paramètre λ. Le flot et le champ de vecteurs deviennent des fonctions de trois variables : ''t'', ''x'' et λ.

Si le champ de vecteurs ''f'' est régulier, le flot est le support de plusieurs théorèmes, piliers de la théorie des équations différentielles. Si la fonction ''f'' est de classe C{{exp|''p''}}, le flot l'est aussi. Ce résultat est parfois considéré comme une forme avancée du [[théorème de Cauchy-Lipschitz]]. Si la fonction ''f'' ne dépend pas du temps, le [[théorème du redressement|théorème du redressement du flot]] indique que, localement, le champ de vecteurs est ''équivalent'' à un champ constant et cette équivalence transforme le flot en une fonction qui à (''x'', ''t'') associe ''x + tv'', où ''v'' est l'unique image du champ constant.

Le flot est utilisé dans diverses branches des mathématiques. En analyse qualitative des équations différentielles, il est le cadre d'expression de théorèmes, comme [[théorème de Poincaré-Bendixson|celui de Poincaré-Bendixson]]. On trouve la notion de flot de manière générale dans l'étude d'un [[système dynamique]] continu. En [[topologie algébrique]], il est utilisé pour démontrer le [[théorème de la boule chevelue]] ou encore [[théorème du point fixe de Brouwer|celui du point fixe de Brouwer]] ; des applications plus avancées définissent des flots caractéristiques de la géométrie des objets étudiés, tels que le [[flot de Ricci]], outil de base utilisé par [[Grigori Perelman]] pour démontrer la [[conjecture de Poincaré]]. L'usage du flot dépasse le cadre strict des mathématiques ; ainsi, le flot de Ricci est à l'origine d'un des modes d'expression des équations de la [[relativité générale]] en physique.

== Définitions ==
=== Introduction ===
Dans la suite de l'article, ''E'' désigne un espace de Banach. Le symbole (1) désigne l'équation différentielle ''x' = f''(''x''), où ''f'' est fonction définie sur un ouvert Ω de ''E'' et à valeurs dans ''E'' et (2) l'équation ''x' = f''(''t'', ''x''). Dans le cas de l'équation (2), Ω désigne un ouvert de '''R'''×''E'' et ''f'' est toujours une fonction définie sur Ω et à valeurs dans ''E''. La condition de Cauchy appelée ''C'' dans l'article, est celle que vérifie une solution ''s'' de l'équation. Dans le cas de (1), elle signifie que ''s''(0) = ''x''<sub>0</sub>, dans le cas (2) les notations choisies sont : ''s''(''t''<sub>0</sub>) = ''x''<sub>0</sub><ref group="Note">Les différentes définitions du vocabulaire utilisé ici sont données dans l'article « [[Théorème de Cauchy-Lipschitz]] ».</ref>. Il existe un troisième cas, plus général et qui permet en particulier l'étude des singularités de certaines équations différentielles. L'équation (3) est : ''x' = f''(''t'', λ, ''x''). Dans ce cas, Ω est un ouvert de '''R'''×''F''×''E'', où ''F'' est encore un espace de Banach. Dans les trois cas, la fonction ''f'' est supposée au moins continue, et localement lipschitzienne par rapport à ''x''. 

Le flot permet de formaliser un vocabulaire adapté pour l'étude de deux questions concernant les équations différentielles : la sensibilité à la [[condition initiale]] et le comportement asymptotique. Ces questions sont au cœur d'une branche mathématique appelée [[système dynamique]].

Si la condition initiale, dite [[théorème de Cauchy-Lipschitz|condition de Cauchy]], est un peu modifiée, il se peut que la solution de l'équation soit de plus en plus éloignée de la courbe intégrale originelle, si la variable ''t'' augmente. Tel est le cas, par exemple, pour les [[théorie du chaos|systèmes chaotiques]]. L'étude de la régularité du flot apporte des premiers éléments de réponse.

Dans certaines conditions et pour une condition initiale donnée, la solution se prolonge à l'infini et « se stabilise » plus ou moins. Elle peut converger vers une valeur, s'approcher de plus en plus d'un comportement périodique, diverger ou encore adopter un comportement différent et appelé chaotique. L'étude de ces différents comportements est l'objet de la seconde question.

=== Sensibilité à la condition initiale ===
[[Fichier:Poincaré-Bendixon-(1).jpg|thumb|Pour l'étude de la sensibilité initiale, on étudie les ''flots locaux'' ou ''globaux'', qui permettent d'énoncer une version plus forte du [[théorème de Cauchy-Lipschitz]].]]
L'étude de la sensibilité à la condition initiale impose un vocabulaire et une représentation géométrique un peu différents de ceux utilisés dans les approches plus élémentaires. Pour comprendre l'origine de cette différence, le plus simple est de considérer le cas d'une [[équation différentielle autonome]], c'est-à-dire de type (1), d'imaginer que Ω est un plan d'eau et que '''R''' représente le temps. le plan d'eau est agité par un courant, représenté par la fonction ''f'', appelée champ de vecteurs. En dimension 2, on représente ce champ de vecteurs en associant à certains points ''x'' de Ω une représentation graphique du vecteur ''f''(''x''), à l'image de la figure de droite. Une courbe intégrale satisfaisant à la condition de Cauchy ''C'' peut s'imaginer comme la trajectoire d'un bouchon placé dans l'eau à l'instant initial 0 et à la position ''x''<sub>0</sub>. Pour connaître d'un seul coup toutes les solutions de l'équation différentielle, il suffit de connaître le mouvement de la surface de l'eau, appelé ''flot'', ''coulée'' ou encore ''courant''<ref>{{Harvsp|Leborgne|1982|p=228}}.</ref>. Dans le cas général de l'équation (2), c'est-à-dire celui où l'équation n'est pas nécessairement autonome, on dispose de la définition suivante :

: Soit ''U'' un voisinage de ''x''<sub>0</sub> tel que {''t''<sub>0</sub>}×''U'' soit inclus dans Ω et ''J'' un [[Intervalle (mathématiques)|intervalle ouvert]] contenant ''t''<sub>0</sub>. Un '''flot local''' est la donnée d'une application β<sub>t<sub>0</sub></sub> de ''J''×''U'' dans ''E'' tel que pour tout ''x'' de ''U'', l'application de ''J'' dans ''E'', qui à ''t'' associe ''β''<sub>''t''<sub>0</sub></sub>(''t'', ''x'') soit une [[Théorème de Cauchy-Lipschitz|courbe intégrale]] de condition de Cauchy (''t''<sub>0</sub>, ''x'') (c'est-à-dire que ''β''<sub>''t''<sub>0</sub></sub>(''t''<sub>0</sub>, ''x'') = ''x'')<ref>{{Harvsp|Berger|Gostiaux|p=40}}.</ref>.

Le [[théorème de Cauchy-Lipschitz]] garantit l'existence de [[théorème de Cauchy-Lipschitz|courbes intégrales maximales]]. C'est-à-dire que, pour une condition de Cauchy donnée, il existe un intervalle maximal, une unique solution de l'équation définie sur cet intervalle et vérifiant la condition de Cauchy. Ce résultat permet d'enrichir les définitions associées au flot.

: Soit ''x'' un point de ''E'' et ''J{{ind|x}}'' l'intervalle égal au [[Ensemble de définition|domaine de définition]] de la solution maximale satisfaisant la condition de Cauchy (''t''<sub>0</sub>, ''x''), qui peut être vide si, par exemple, (''t''<sub>0</sub>, ''x'') n'est pas élément de Ω. On note ''D''(''f'') et {{refnec|date=mai 2017|on appelle '''domaine de définition du flot global''' l'ensemble des couples (''t'', ''x'') de '''R'''×''E'' tel que ''t'' soit un élément de ''J{{ind|x}}''}}.

:{{refconf|date=mai 2017|Le '''flot global''' de ''f'' est l'application α<sub>''t''<sub>0</sub></sub>}} de ''D''(''f'') dans ''E'', telle que l'application qui à ''t'' associe α<sub>''t''<sub>0</sub></sub>(''t'', ''x'') soit la solution maximale de (2) associée à la condition de Cauchy (''t''<sub>0</sub>, ''x'')<ref>{{Harvsp|Berger|Gostiaux|p=46}}.</ref>.

Dans toute la suite de l'article α désigne le flot global du champ vectoriel ''f''. Si α ne comporte pas d'indice, cela signifie que la fonction qui à ''t'' associe α(''t'', ''x'') est la courbe intégrale maximale ''s'' vérifiant la condition de Cauchy : ''s''(0) = ''x''. Sinon, α<sub>''t''<sub>0</sub></sub> désigne le flot global vérifiant α<sub>''t''<sub>0</sub></sub>(''t''<sub>0</sub>, ''x'') = ''x''.

=== Comportement asymptotique ===
La deuxième question concerne le comportement asymptotique<ref name=Sari5>{{Harvsp|Sari|2006|p=5}}.</ref> du flot, autrement dit ce qui se passe une fois le système stabilisé, si jamais il se stabilise. Pour la traiter, un vocabulaire spécifique existe ; il suppose en général que l'équation différentielle considérée est autonome<ref group="Note">C'est le choix de {{Harvsp|Sari|2006}}.</ref>, c'est-à-dire de type (1) {{Supra|Introduction}}.

: L''''orbite'''<ref name=Sari4/> ou la '''trajectoire''' d'un point ''x'' de Ω est l'[[Image d'une application|image]] d'une courbe intégrale maximale passant par ''x''.

Le théorème de Cauchy-Lipschitz montre que les orbites forment une [[Partition (mathématiques)|partition]] de Ω.

: Un '''ensemble invariant par le flot'''<ref name=Sari4>{{Harvsp|Sari|2006|p=4}}.</ref> est une réunion d'orbites.

Pour reprendre la [[métaphore]] de l'introduction, un ensemble ''A'' est [[invariant]] par le flot si un bouchon qui est dans ''A'' y a toujours été et y restera toujours. {{Citation|Pour les systèmes dynamiques, les parties invariantes jouent le même rôle que les [[Connexité (mathématiques)|connexes]] en [[topologie|topologie élémentaire]]<ref name=Cantat/>}} : on peut limiter l'étude du comportement asymptotique à une telle zone.

: Un point de ''x'' de Ω est dit '''point d'équilibre''', '''point fixe''' ou '''point stationnaire'''<ref name=Sari4/> si son orbite est réduite au [[Singleton (mathématiques)|singleton]] {''x''} ou encore si ''x'' est un [[zéro d'une fonction|zéro de la fonction]] ''f''. S'il existe un [[Voisinage (mathématiques)|voisinage]] de ''x'' tel que toute courbe intégrale maximale qui rencontre ce voisinage converge vers ''x'', le point est dit '''attractif'''.

Un point d'équilibre ''x''{{ind|e}} est un point sans courant, ou encore un point tel que, si l'on y place le bouchon, il reste indéfiniment immobile. Deux comportements différents peuvent se produire, le point peut être attractif ou non. Il peut exister une zone pas trop petite (d'intérieur non vide) tel que tout point dans cette zone finit par s'immobiliser en ''x''{{ind|e}}. Ce point semble « attirer » les trajectoires vers lui, pour cette raison, on parle de point attractif. Sinon, tout point se rapprochant de ''x''{{ind|e}} finit par s'éloigner, ce point d'équilibre est alors instable et il faut positionner le bouchon exactement sur le point pour qu'il y reste immobile. Dans le cas général, d'autres comportements que la convergence sont possibles :

:Soit ''s'' une courbe intégrale définie sur '''R'''.
:*Lorsque ''s'' admet une limite en {{math|+∞}} et en {{math|–∞}}, l'orbite est dite '''hétérocline''' si les deux limites sont distinctes et '''homocline''' sinon.
:*Si ''s'' est [[Fonction périodique|périodique]], l'orbite est dite '''périodique'''.
:*L''''ensemble ω-limite'''<ref name=Cantat>{{Harvsp|Cantat|1995}}.</ref> de l'orbite est l'ensemble des [[Valeur d'adhérence#Définitions|valeurs d'adhérence de la fonction]] ''s'' en {{math|+∞}}. On le note ω(''x''), où ''x'' est un point quelconque de l'orbite. On définit de même l''''ensemble α-limite''' α(''x'') comme l'ensemble des valeurs d'adhérence de ''s'' en {{math|–∞}}.
Ces deux ensembles sont invariants par le flot. Si l'orbite est périodique, ses deux ensembles limites lui sont égaux.

Une notion connexe est celle d'attracteur. L'<nowiki/>'''attracteur futur''' est le plus petit ensemble contenant tous les ensembles ω(''x'') si ''x'' décrit Ω, à l'exception, peut-être d'un [[Ensemble négligeable|ensemble de mesure nulle]]. L''''attracteur passé''' correspond à la même définition, mais cette fois-ci avec les ensembles α-limites<ref name=Sari5/>.

== Sensibilité à la condition initiale ==
=== Résumé ===
{{Article détaillé|Théorème de Cauchy-Lipschitz}}
Dans sa version élémentaire, le théorème de Cauchy-Lipschitz montre l'existence et l'unicité d'une courbe intégrale maximale. On peut alors définir le flot global α comme la fonction qui, à un triplet (''t''{{ind|1}}, ''t''{{ind|0}}, ''x''{{ind|0}}) avec (''t''{{ind|0}}, ''x''{{ind|0}}) élément de Ω, associe (lorsqu'elle est définie) la valeur en ''t''{{ind|1}} de la solution maximale du problème de Cauchy ''x{{'}}'' = ''f''(''t'', ''x''), ''x''(''t''{{ind|0}}) = ''x''{{ind|0}}. On sait déjà que par rapport à ''t''{{ind|1}}, α est de classe C{{exp|''p''+1}} si ''f'' est de classe C{{exp|''p''}}, mais on va préciser son comportement par rapport au triplet (''t''{{ind|1}}, ''t''{{ind|0}}, ''x''{{ind|0}}) :

* ''Le flot global est défini sur un ouvert et localement lipschitzien<ref>{{Harvsp|Berger|Gostiaux|p=50}}, le démontrent dans le cas particulier où (localement) ''f'' est lipschitzienne par rapport à (''t'', ''x''), et pas seulement par rapport à ''x'' pour ''t'' fixé.</ref>.''
* ''Si la fonction f est de classe'' C''{{exp|p}}, il en est de même du flot global<ref name=Lang373>{{Harvsp|Lang|1993|p=373}}, pour une équation autonome.</ref>.

Pour certaines études, comme celles analysant des singularités de courbes intégrales, il est utile d'ajouter un paramètre à l'équation, qui prend la forme : ''x' = f''(''t'', ''x'', λ). Le flot global dépend alors du paramètre λ, choisi dans un Banach ''F''. La question est celle de la régularité du flot global en fonction du paramètre. Le théorème décrivant cette situation est parfois dénommé {{Citation|théorème de Cauchy-Lipschitz : cas non autonome avec paramètre<ref>{{Harvsp|Leborgne|1982|p=219}}.</ref>}}. Si la fonction ''f'' est de classe C{{exp|''p''}}, le flot, considéré comme une application sur un ouvert de '''R'''×'''R'''×''E''×''F'', est aussi de classe C{{exp|''p''}}.

=== Continuité du flot===
Dans un premier temps, on étudie le cas de l'équation (2), c'est-à-dire celui d'une équation non nécessairement [[équation différentielle autonome|autonome]]<ref>Le cas autonome est traité dans {{Harvsp|Berger|Gostiaux|p=42}}. [[Théorème de Cauchy-Lipschitz#Équation différentielle non autonome|De même que dans le théorème de Cauchy-Lipschitz classique, un jeu d'écriture permet]] {{Harv|Berger|Gostiaux|p=49}} d'en déduire le cas non autonome, mais seulement pour ''t''{{ind|0}} fixé et sous réserve que (localement) ''f'' soit lipschitzienne par rapport à (''t'', ''x''), et pas seulement par rapport à ''x'' pour ''t'' fixé.</ref>, mais sans paramètre. L'objectif est d'établir la continuité locale d'un flot, c'est-à-dire d'analyser ce qui se passe si la condition de Cauchy est peu modifiée.

On commence par préciser le comportement de α(''t'', ''t''{{ind|0}}, ''x''{{ind|0}}) par rapport à (''t''{{ind|0}}, ''x''{{ind|0}}) dans le cas où ''t'' est proche de ''t''{{ind|0}}.

En affinant la [[Théorème de Cauchy-Lipschitz#Solution locale|technique de point fixe utilisée pour démontrer le théorème de Cauchy-Lipschitz classique]], on parvient au résultat suivant :

{{Théorème|nom=Lemme|Soit (''t''{{ind|0}}, ''x''{{ind|0}}) un élément de Ω. Au voisinage de (''t''{{ind|0}}, ''t''{{ind|0}}, ''x''{{ind|0}}), le flot global associé à l'équation (2) est défini et lipschitzien<ref>{{Note autre projet|wikiversité|Calcul différentiel/Équations différentielles#Dépendance des conditions initiales|« Dépendance des conditions initiales » du chapitre sur les équations différentielles|début=Ce « lemme de continuité locale » est démontré dans le paragraphe}}</ref>.
|style=display:table}}

[[Fichier:Flot-discontinu.jpg|thumb|Un flot discontinu ne peut pas se produire si le champ de vecteurs le définissant est continu et défini sur un ouvert.]]

Pour rendre global le résultat précédent de continuité « au voisinage de la diagonale », il reste encore un peu de travail. Il s'agit surtout de montrer que la configuration de droite ne peut pas se produire. En reprenant la métaphore de l'introduction, on suppose que l'étang contient un rocher. Sur la zone rouge, le flot contourne le rocher, sur la zone bleue il passe tout droit. Le flot possède alors des points de discontinuité. Pour s'en rendre compte, on considère l'évolution d'un voisinage d'un point symbolisé par une zone grise circulaire à un instant 0, sur la figure de droite. Ce point est situé à la frontière des deux zones. À l'instant ''t''<sub>0</sub>, ce voisinage est découpé en deux parties éloignées l'une de l'autre. Cette configuration peut se produire, par exemple si Ω n'est pas un ouvert. En revanche, dans les hypothèses étudiées, cette configuration est impossible :

{{Théorème|Le flot global associé à l'équation (2) est défini sur un ouvert et localement lipschitzien<ref>{{Note autre projet|wikiversité|Calcul différentiel/Équations différentielles#Dépendance des conditions initiales|« Dépendance des conditions initiales » du chapitre sur les équations différentielles|début=Ce « théorème de continuité du flot » est démontré dans le paragraphe}}</ref>.
|style=display:table}}

Une remarque est utile pour le démontrer. Elle est formulée ici uniquement pour les équations autonomes<ref>Pour la formulation générale, voir {{Harvsp|Paulin|2008|page=257}}.</ref>. Le point β(''s'', ''x'') désigne la position, à l'instant ''s'', du point qui était en ''x'' à l'instant 0. Le point β(''t'', β(''s'', ''x'')) désigne la position à l'instant ''t'', du point qui était en β(''s'', ''x''), à l'instant 0. Il désigne aussi la position à l'instant ''s'' + ''t'' du point qui était en position ''x'' à l'instant ''s''. Ce qui signifie, si les différentes valeurs sont dans les bons domaines de définitions :
<center><math>\beta(t,\beta(s,x)) = \beta(t+s,x)\;</math></center>

Ainsi, si les domaines de définitions s'y prêtent, il est possible de combiner les flots locaux, pour en construire un plus vaste. Cette technique permet de munir les flots locaux d'une structure de [[Théorème de Hille-Yosida#Semi-groupes|semi-groupe]]. Si le flot global est toujours défini sur '''R''', les solutions intégrales maximales partageant une portion de domaine de définition forment un groupe, image de '''R''' par un morphisme de groupes.

=== Équation différentielle avec paramètre ===
On souhaite maintenant étudier le cas d'un flot local associé à l'équation :
<center><math>(3)\quad\frac{\mathrm dx}{\mathrm dt}=f\left(t,x,\lambda\right)</math>.</center>

On suppose que λ est élément d'un espace de Banach ''F'', Ω est un ouvert de '''R'''×''E''×''F'' et ''f'' est une fonction continue de Ω dans ''E'', localement lipschitzienne ''par rapport à'' (λ, ''x'').

Cette situation, en apparence plus générale que celle de l'équation (2), en est en fait un cas particulier, en remplaçant simplement, dans cette dernière :
*l'espace de Banach ''E'' (dans lequel chaque courbe intégrale prend ses valeurs) par ''E''×''F'' ;
*le champ de vecteurs (d'un ouvert de '''R'''×''E'' dans ''E'') par la fonction :<center><math>g:\Omega\to E\times F,\quad y=\left(t,x,\lambda\right)\mapsto(f(y),0)</math> ;</center>
*l'équation ''x{{'}}'' = ''f''(''t'', ''x'') par (''x{{'}}'', λ{{'}}) = ''g''(''t'', ''x'', λ) ;
*la « condition initiale » (''t''{{ind|0}}, ''x''{{ind|0}}) par (''t''{{ind|0}}, ''x''{{ind|0}}, λ{{ind|0}}).

Le flot α dépend maintenant de trois variables : ''t'', ''t''{{ind|0}}, ''x''{{ind|0}} et la « variable » supplémentaire λ, mais compte tenu de la définition de ''g'', la valeur de λ le long d'une courbe intégrale reste constamment égale à sa valeur initiale λ{{ind|0}}. Cet artifice permet d'étendre directement aux équations à paramètres le résultat du paragraphe précédent : soit ''f'' une fonction continue, d'un ouvert de '''R'''×''F''×''E'' dans ''E'', localement lipschitzienne par rapport à ses variables dans ''F''×''E''. Le flot global associé à l'équation (3) est défini sur un ouvert de '''R'''×'''R'''×''F''×''E'' et localement lipschitzien<ref>{{Harvsp|Berger|Gostiaux|p=50}}, pour ''t''{{ind|0}} fixé.</ref>.

Mais ce corollaire n'est d'aucune utilité pour la suite, contrairement au théorème suivant, qui généralise le lemme du paragraphe précédent et se démontre (directement) de la même manière :
{{Théorème|Soient <math>\Lambda</math> un [[espace topologique]], <math>\Omega</math> un ouvert de <math>\R\times\Lambda\times E</math> et <math>f:\Omega\to E</math> une application continue et localement lipschitzienne ''par rapport à sa dernière variable'' <math>x</math>.

Alors, pour tout <math>\left(t_0,\lambda_0,x_0\right)\in\Omega</math>, le flot correspondant est défini et continu au voisinage de <math>\left(t_0,t_0,\lambda_0,x_0\right)</math>.
|style=display:table}}

Le jeu d'écriture lui-même sera utile en revanche, pour étendre aux équations de type (3) un lemme (ci-dessous) sur la régularité C{{1}}.

=== Régularité du flot au voisinage de la diagonale===
Les démonstrations associées à la régularité du flot ont pendant longtemps été un peu complexes. De manière indépendante, {{refsou|date=mai 2017|Pugh}} et Robbin<ref>{{Article|langue=en|auteur=Joel W. Robbin|titre=On the existence theorem for differential equations|revue=[[Proceedings of the American Mathematical Society|Proc. Amer. Math. Soc.]]|vol=19|année=1968|p.=1005-1006|url=http://www.ams.org/journals/proc/1968-019-04/S0002-9939-1968-0227583-0}}.</ref>{{,}}<ref>Pour une variante proche de celle de Robbin mais utilisant le [[théorème d'inversion locale]], voir {{Lafontaine1}}, 2010, {{p.|44-45}}.</ref>{{,}}<ref group=Note>{{Harvsp|Lang|1993|p=371}}, mentionne Pugh et Robbin, mais sans référence. Sa démonstration s'écarte de celle de {{Harvsp|Robbin|1968}} par le choix des espaces de Banach sur lesquels il applique le théorème des fonctions implicites, ce qui la rend plus acrobatique (voir en particulier {{Harvsp|Lang|1993|p=374-376}}).</ref> ont trouvé une démonstration élémentaire utilisant le [[théorème des fonctions implicites]]. Cette fois, l'équation considérée est celle notée (1), c'est-à-dire qu'on ne traite que d'une équation différentielle autonome. On en déduit une partie du résultat analogue pour les équations non autonomes.

Soient ''x''<sub>0</sub> un point de Ω et ''b'' un réel strictement positif tel que la courbe intégrale β{{ind|0}} qui vaut ''x''{{ind|0}} en 0 soit définie au moins sur [–''b'', ''b'']. Choisissons un ouvert Ω' contenant ''x''{{ind|0}} et un ouvert Ω" contenant Im(β{{ind|0}}) – ''x''{{ind|0}}, tels que Ω' + Ω" soit inclus dans Ω. Notons enfin ''G'' l'espace de Banach des fonctions de classe C{{1}} de [–''b'', ''b''] dans ''E'' et nulles en 0, muni de la norme de la [[convergence uniforme]] de la dérivée, ''V'' l'ouvert de celles à valeurs dans Ω"<ref group="Note">''V'' est même un ouvert pour la topologie de la [[convergence simple]], comme intersection des φ{{ind|''t''}}{{-1}}(Ω"), où φ{{ind|''t''}} : ''G'' → ''E'' est l'application d'évaluation en ''t''.</ref>, et ''H'' l'espace de Banach des fonctions continues de [–''b'', ''b''] dans ''E'', muni de la norme de la convergence uniforme.

On considère l'application ''T'' de Ω'×''V'' dans ''H'' définie par :
<center><math>\forall x \in\Omega'\quad\forall \sigma \in V\quad\forall t \in [-b,b]\quad T(x,\sigma)(t)=\sigma'(t)-f(x+\sigma(t))</math>.</center>

Dire qu'un couple (''x'', σ) est un [[zéro d'une fonction|zéro]] de la fonction ''T'' revient à dire que ''x'' + σ est une solution sur [–''b'', ''b''] de l'équation (1) avec la condition de Cauchy ''x'' + σ(0) = ''x''. En particulier, (''x''{{ind|0}}, β{{ind|0}} – ''x''{{ind|0}}) est un zéro de ''T''.

On démontre que si ''f'' est de classe C{{1}} alors ''T'' aussi et pour ''b'' suffisamment petit, ''T'' satisfait les hypothèses du théorème des fonctions implicites au point (''x''{{ind|0}}, β{{ind|0}} – ''x''{{ind|0}}). Cela signifie que la [[fonction implicite]] donnée par l'équation ''T''(''x'', σ) = 0 est de classe C{{1}} au voisinage de ''x''{{ind|0}}. Or cette fonction est celle qui associe à ''x'' la fonction β('''∙''', ''x'') – ''x'' définie sur [–''b'', ''b''], ce qui permet de démontrer un premier résultat<ref>Pour une démonstration détaillée (sous des hypothèses plus faibles : équation non autonome, avec ''f'', fonction de (''t'', ''x''), seulement continue et de différentielle ''[[Dérivée partielle|partielle par rapport à x]]'' continue), voir le {{Note autre projet|wikiversité|Calcul différentiel/Équations différentielles#Dépendance des conditions initiales|preuve du « théorème de régularité locale », dans le chapitre sur les équations différentielles|début=point 1 de la}}</ref> :

{{Théorème|nom=Lemme|
Si ''f'' est de classe C{{1}}, le flot global l'est aussi au voisinage de (0, ''x''{{ind|0}}), pour tout point ''x''{{ind|0}} de Ω.
|style=display:table}}

Ce lemme s'étend aux équations différentielles de type (3), par le jeu d'écriture ci-dessus pour le paramètre et [[Théorème de Cauchy-Lipschitz#Équation différentielle non autonome|un autre, analogue]], pour la variable temporelle.

Or dès que ''f'' est de classe C{{1}}, la deuxième différentielle partielle de β vérifie une équation différentielle de type (3) :
<center><math>D_1D_2\beta\left(t,x\right)=D_2f\left(t,\beta\left(t,x\right)\right)\circ D_2\beta\left(t,x\right)</math><ref>{{Note autre projet|wikiversité|Calcul différentiel/Équations différentielles#Dépendance des conditions initiales|preuve du « théorème de régularité locale », dans le chapitre sur les équations différentielles|début=Point 2 de la}}</ref>.</center>
Ceci permet de montrer par récurrence que pour toute équation de type (3), si ''f'' est de classe C{{exp|''p''}} alors le flot α(''t'', ''t''{{ind|0}}, λ, ''x''{{ind|0}}) l'est aussi, par rapport à (''t'', λ, ''x''{{ind|0}}). Nous admettrons la régularité par rapport à ''t''{{ind|0}}<ref>{{Note autre projet|wikiversité|Calcul différentiel/Équations différentielles#Dépendance des conditions initiales|preuve du « théorème de régularité locale », dans le chapitre sur les équations différentielles|début=Points 2 et 3 de la}}</ref>, qui utilise entre autres, pour l'initialisation de la récurrence, le théorème de continuité [[#Équation différentielle avec paramètre|ci-dessus]]. Moyennant quoi, on peut énoncer :

{{Théorème|Soit ''f'' une fonction de classe C{{exp|''p''}} définie sur un ouvert Ω de '''R'''×''F''×''E'' et (''t''<sub>0</sub>, λ<sub>0</sub>, ''x''<sub>0</sub>) un élément de Ω. Au voisinage de (''t''<sub>0</sub>, ''t''<sub>0</sub>, λ<sub>0</sub>, ''x''<sub>0</sub>), le flot global associé à l'équation (3) est de classe C{{exp|''p''}}<ref name=Lang373/>.
|style=display:table}}

=== Régularité globale du flot===
Le théorème de Cauchy-Lipschitz, sous sa forme élémentaire, garantit à la fois l'existence et l'unicité du flot global α, et sa régularité par rapport à la première variable. Le paragraphe précédent apporte aussi des informations. Il établit que si ''t'' est proche de 0 et si ''t'' et ''x'' varient peu, le flot est de classe C{{exp|''p''}}, si ''f'' l'est. En revanche, la régularité globale du flot α suppose aussi que l'expression α(''t'', ''x'') soit de classe C{{exp|''p''}} en ''x'', même si ''t'' n'est pas proche de 0. On a supposé ici que (''t'', ''x'') est un point du domaine de définition de α.

La même technique que dans le paragraphe « Continuité du flot » {{supra|Continuité du flot}} permet de rendre global le résultat précédent de régularité du flot « au voisinage de la diagonale » pour une équation de type (2) — donc aussi, par le « jeu d'écriture » {{supra|Équation différentielle avec paramètre}}, pour une équation de type (3) :
* '''Théorème''' : ''Pour une équation différentielle non autonome à paramètre, si le champ de vecteurs f est de classe ''C''{{exp|p}}, le flot global l'est aussi<ref>{{Harvsp|Berger|Gostiaux|p=47}}, pour ''t''{{ind|0}} fixé.</ref>.''

== Comportement asymptotique ==
=== Comportement à la limite du domaine ===
[[Fichier:Flot-borné.jpg|thumb|L'équation différentielle ''x' = –x''{{2}} n'admet pas de solution sur '''R''' tout entier si l'image de 0 n'est pas nulle.]]
Pour une valeur de ''x'' prise dans ''E'', la courbe intégrale maximale ''t'' ↦ ''s''(''t'') = α<sub>''t''<sub>0</sub></sub>(''t'', ''x'') est définie sur un intervalle ouvert. Cet intervalle n'est pas toujours égal à '''R'''. Deux « obstacles » peuvent se présenter. Un premier cas est celui où la courbe a pour extrémité un point frontière de l'ouvert Ω. Si tel n'est pas le cas, il n'existe qu'un seul obstacle possible : la dérivée de ''s'' n'est pas bornée. Cela résulte de la proposition suivante<ref>Reformulation de la proposition 7.25 de {{Harvsp|Paulin|2008|page=262}}, qui correspond au 10.5.5 de {{Ouvrage |langue=en |auteur1=[[Jean Dieudonné|J. Dieudonné]] |titre=Foundations of Modern Analysis |éditeur= |année=1960 |passage=284-285 |lire en ligne=https://archive.org/stream/FoundationsOfModernAnalysis_578/Dieudonne-FoundationsOfModernAnalysis#page/n300/mode/1up}}. {{Note autre projet|wikiversité|Calcul différentiel/Exercices/Équations différentielles non linéaires#Exercice 2|exercice corrigé de la leçon « Calcul différentiel »|début=Voir aussi cet}}</ref> :

* ''Soit b la borne supérieure du domaine de définition d'une courbe intégrale maximale s. Si b est fini et si s'''{{'}}''' est bornée au voisinage de b, alors s admet au point b une limite c et le couple ''(''b'', ''c'')'' appartient à la [[Frontière (topologie)|frontière]] de ''Ω.

Un exemple<ref>{{Harvsp|Berger|Gostiaux|p=52}}.</ref>{{,}}<ref group="Note">On pourrait exhiber des exemples plus représentatifs car dans celui-ci, non seulement la dérivée ''s''', mais la fonction ''s'' elle-même n'est pas bornée au voisinage du point limite. Cela est dû au fait que — dans les conditions d'application du théorème de Cauchy-Lipschitz — toute solution d'une équation autonome à valeurs scalaires est monotone.</ref> est donné par l'équation ''x' = –x''{{2}}. La solution s'écrit :
<center><math>\alpha(t,x)=\frac x{xt+1}</math>.</center>
Les courbes intégrales dont l'[[ordonnée à l'origine]] est différente de 0 ne se prolongent pas au-delà de –1/''x''.

=== Ensemble ω-limite ===
[[Fichier:Attracteur étrange de Lorenz.png|thumb|Le [[système dynamique de Lorenz]] montre un exemple où les ensembles ω-limites sont en général de structure complexe.]] 
On suppose ici que l'équation différentielle étudiée est autonome. Une méthode pour étudier le comportement asymptotique est d'analyser l'ensemble ω-limite d'une orbite. Cet ensemble possède toujours les propriétés suivantes :
* ''Soit x un élément de ''Ω'', l'ensemble ''ω''-limite ''ω(''x'')'' est un fermé invariant par le flot.''

En particulier, si ''y'' est un point de la trajectoire passant par ''x'', alors ω(''x'') = ω(''y''). Si l'orbite de ''x'' n'est pas à valeurs dans un ensemble borné, la proposition peut s'avérer moins forte qu'il n'y parait : l'ensemble ω(''x'') peut être vide. Tel est le cas si la courbe intégrale ''s'' ayant pour condition de Cauchy ''s''(0) = ''x'' est définie par ''s''(''t'') = ''tv'' où ''v'' est un vecteur non nul de ''E''. Si l'orbite est [[Compacité (mathématiques)|compacte]], on a par contre :
* ''Avec les notations précédentes, si la trajectoire de ''x'' est à valeurs dans un compact ''K'', ω(''x'') est un [[Connexité (mathématiques)|connexe]] compact non vide<ref>Une démonstration, détaillant celle de {{Harvsp|Cantat|1995}}, est proposée dans la [[Théorème de Poincaré-Bendixson#Flot|section « Flot » de l'article « Théorème de Poincaré-Bendixson »]].</ref>.''

La trajectoire peut néanmoins être complexe, le [[système dynamique de Lorenz]] est un exemple exhibant des ensembles ω-limites compacts dont la [[dimension de Hausdorff]] n'est pas égale à 1.

=== Théorème de Poincaré-Bendixson ===
{{Article détaillé|Théorème de Poincaré-Bendixson}}
Il existe un cas où le flot dispose d'un comportement asymptotique relativement simple, il se produit si l'équation différentielle associée est autonome, si l'espace ''E'' est le plan et si la courbe intégrale considérée est compacte. On dispose alors du théorème suivant, si α est le flot et ''x'' un point de Ω.

* ''Dans le cadre des hypothèses du paragraphe, la fonction qui à t associe ''α(''t'', ''x'')'' est définie sur '''''R'''''. Si cette fonction n'est pas convergente, l'ensemble ''ω''-limite est l'image d'une orbite cyclique<ref>{{Lien web|langue=en|auteur=Oliver Knill|url=http://www.math.harvard.edu/archive/118r_spring_05/handouts/poincare.pdf|titre=The Poincaré-Bendixon theorem|site=[[Université Harvard]]|année=2005}}.</ref>.''

Ce théorème se généralise mal. Si ''E'' est de dimension strictement supérieure à 2, le système dynamique de Lorenz montre que le résultat du théorème n'est plus généralement vrai.

== Analyse locale ==
{{Article détaillé|Théorème du redressement}}

==Notes et références==
===Ouvrages cités===
* {{BergerGostiaux1}} {{Commentaire biblio|Ce livre se limite au cas où ''E'' est un espace vectoriel de dimension finie. Il traite presque de l'intégralité des informations de l'article sur la sensibilité à la condition initiale, à l'exception du caractère C{{exp|''p''}} du flot local.}}
*{{Article|auteur=Serge Cantat|url=http://math.univ-lyon1.fr/~benzoni/CantatJME3.pdf|titre=Théorème de Poincaré-Bendixson|éditeur=[[ENS Lyon]]|revue=Le journal de maths des élèves|vol=1|année=1995|numéro=3}}
*{{Ouvrage |langue=en |auteur1=[[Serge Lang]] |titre=Real and Functional Analysis |éditeur= |collection=[[Graduate Texts in Mathematics|GTM]] |numéro dans collection=142 |année=1993 |numéro d'édition=3 |pages totales=580 |isbn= |lire en ligne=https://archive.org/stream/springer_10.1007-978-1-4612-0897-6/10.1007-978-1-4612-0897-6}}
* [[Paul Malliavin]], ''Géométrie différentielle intrinsèque'', [[Hermann (éditions)|Hermann]], 1972 {{ISBN|2-7056-5696-0}} {{Commentaire biblio|Cette référence, encore sur la sensibilité à la condition initiale, est plus complète, mais plus ardue que la précédente. L'orientation est définitivement géométrique.}}
*{{Ouvrage |langue=fr |auteur1=Daniel Leborgne |titre=Calcul différentiel et géométrie |éditeur=[[PUF]] |lieu=Paris |année=1982 |pages totales=262 |isbn=2-13-037495-6}} {{Commentaire biblio|Ce livre suppose déjà connu le théorème dans le cas des Banach et le généralise aux variétés différentielles. Il ne traite pas l'analyse du comportement asymptotique d'un flot.}}
*{{Lien web|auteur=Frédéric Paulin|titre=Topologie, analyse et calcul différentiel|url=https://www.math.u-psud.fr/~paulin/notescours/cours_analyseI.pdf|éditeur=École Normale Supérieure|année=2008}}<!--http://archive.wikiwix.com/cache/?url=http%3A%2F%2Fwww.mmfai.ens.fr%2Fcours%2Fcours_analyseIpaulin0809.pdf-->
*{{Lien web|auteur=Tewfik Sari|url=http://www.math.uha.fr/sari/papers/cimpa05.pdf|titre=Introduction aux systèmes dynamiques et applications à un modèle cosmologique|année=2006|site=[[Université de Haute-Alsace]]}}, publié dans ''Géométries et dynamiques'', Hermann, 2008 {{ISBN|978-2-70566772-6}}, [http://archive.wikiwix.com/cache/?url=http%3A%2F%2Fwww.umpa.ens-lyon.fr%2F~zeghib%2Fsari.pdf {{p.|259-274}}]

=== Notes ===
{{Références|groupe="Note"}}

=== Références ===
{{Références|colonnes=2}}

==Bibliographie complémentaire==
{{Ouvrage |auteur1=[[François Laudenbach]] |titre=Calcul différentiel et intégral |éditeur=Éditions de l'École polytechnique |année=2001 |isbn= |lire en ligne={{Google Livres|Ws9A7ZoRJNEC}}}}

{{Palette|Variété différentielle}}
{{Portail|analyse|géométrie}}

[[Catégorie:Équation différentielle]]
[[Catégorie:Champ de vecteurs]]