Voir la source de Équation différentielle ordinaire

En [[mathématiques]], une '''équation différentielle ordinaire''' (parfois simplement appelée ''équation différentielle'' et abrégée en '''EDO''') est une [[équation différentielle]] dont la ou les fonctions inconnues ne dépendent que d'une seule variable; elle se présente sous la forme d'une relation entre ces [[fonction (mathématiques)|fonctions]] inconnues et leurs [[dérivée]]s successives.

Le terme ''ordinaire'' est utilisé par opposition au terme ''équation différentielle partielle'' (plus communément [[équation aux dérivées partielles]], ou ''EDP'') où la ou les fonctions inconnues peuvent dépendre de plusieurs variables. ''Dans la suite de l'article, le terme équation différentielle est utilisé pour signifier équation différentielle ordinaire''.

L'ordre d'une équation différentielle correspond au degré maximal de dérivation auquel l'une des fonctions inconnues a été soumise. Il existe une forme de référence à laquelle on essaie de ramener les équations différentielles ordinaires par divers procédés mathématiques :
:<math>X'(t)=F(t,X(t))</math>,
équation d'ordre 1 où <math>X</math> est la fonction inconnue, et <math>t</math> sa variable.

Les équations différentielles représentent un objet d'étude de toute première importance, aussi bien en [[mathématiques pures]] qu'en [[mathématiques appliquées]].
Elles sont utilisées pour construire des [[Modèle mathématique|modèles mathématiques]] de processus d'évolution [[physique]]s et [[biologie|biologiques]], par exemple pour l'étude de la [[radioactivité]], la [[mécanique céleste]] ou la [[dynamique des populations]]... La variable <math>t</math> représente alors souvent le temps, même si d'autres choix de modélisation sont possibles.

Les objectifs principaux de la théorie des équations ordinaires sont la résolution explicite complète quand elle est possible, la résolution approchée par des procédés d'[[analyse numérique]], ou encore l'étude qualitative des solutions. Ce dernier domaine s'est progressivement étoffé, et constitue l'un des composants principaux d'une vaste branche des mathématiques contemporaines : l'étude des [[systèmes dynamiques]].

== Définition ==
=== Premiers exemples ===

Même si ce n'est pas la discipline qui a fait naître les équations différentielles, la [[dynamique des populations]] en illustre de façon simple des exemples parmi les plus accessibles. Ainsi, l'étude d'une population isolée dans un milieu produisant de la nourriture en abondance conduit au modèle suivant pour l'effectif <math>x</math> en fonction du temps <math>t</math> :
:<math>x'(t)=K\,x(t)</math>,
c'est-à-dire que l'accroissement de population <math>x'(t)</math> est, à chaque instant, proportionnel à la taille de la population <math>x(t)</math>. Les solutions de cette équation font apparaître un phénomène de [[croissance exponentielle]].

[[Fichier:Lotka_Volterra_Phase_plane_100.svg|gauche|vignette|Les courbes d'évolution des populations pour les équations de Lotka-Volterra.]]
Un système plus complexe, formé de deux espèces, proie et prédateur, conduit aux [[équations de Lotka-Volterra]]
:<math>\begin{cases} x'(t)&=A\,x(t)-B\,x(t)\,y(t) \\ y'(t)&=-C\,y(t)+D\,x(t)\,y(t).\end{cases}</math>
L'effectif des proies est <math>x(t)</math>, celui des prédateurs <math>y(t)</math>. On retombe sur le cas précédent si <math>y</math> est nul. La quantité <math>x(t)y(t)</math> est une probabilité de rencontre, qui influe négativement sur une population (les proies), positivement sur l'autre (les prédateurs). À chaque instant, connaissant les populations en présence, on peut décrire la tendance. Ces deux équations sont couplées c'est-à-dire qu'il faut les résoudre ensemble. Mathématiquement, il faut les concevoir comme une seule équation dont l'inconnue est le couple <math>(x(t),y(t))</math> . Si l'effectif initial des populations est connu, l'évolution ultérieure est parfaitement déterminée. Elle se fait le long d'une des courbes d'évolution figurées ci-contre, qui laissent apparaître un comportement cyclique.

Une des plus célèbres équations différentielles est la [[Lois du mouvement de Newton#Deuxième loi de Newton ou principe fondamental de la dynamique de translation|relation fondamentale de la dynamique]], de [[Isaac Newton|Newton]] : <math>f=ma</math>, où <math>m</math> est la masse d'une particule, <math>f</math> la force exercée sur celle-ci et <math>a</math> l'accélération qui en résulte. Dans le cas d'un mouvement rectiligne, si la force subie est fonction de la position (par exemple dans le cas d'un [[ressort]]) on obtient une équation de la forme
:<math>\displaystyle m\frac{\mathrm{d}^2 x}{\mathrm{d}t^2}=f(x).</math>
Cette fois, pour déterminer parfaitement le mouvement, il faut se donner position et vitesse initiales.

=== Équation différentielle, processus d'évolution et déterminisme ===

Les caractéristiques d'un système dynamique régi par une équation différentielle sont les suivantes :
* les états ''a priori'' possibles pour le système forment un espace de dimension finie, c’est-à-dire peuvent être décrits par un nombre fini de variables. Cet espace est l'[[espace des phases]]. Par exemple, pour décrire le mouvement d'une particule dans l'espace usuel, il faut trois variables. Pour le mouvement d'un solide, six sont nécessaires ;
* les lois qui gouvernent l'évolution temporelle sont des fonctions au moins [[dérivée|dérivables]] ;
* l'évolution du système est [[déterminisme|déterministe]] : connaissant les conditions initiales, c'est-à-dire l'état du système au temps présent, on peut en déduire l'état du système à n'importe quel temps du futur ou du passé.

L'aspect déterministe des équations différentielles a des implications particulièrement fortes, et se concrétise mathématiquement par le [[#Conditions initiales, théorème de Cauchy-Lipschitz|théorème de Cauchy-Lipschitz]].

=== Définition générale ===

Soit ''E'' un [[espace vectoriel normé]]. Une '''équation différentielle ordinaire''' (appelée simplement ''équation différentielle'' dans la suite de l'article) est une [[équation]] de la forme<center>
<math>F(x,y,y', \dots, y^{(n)})=0</math></center>
où ''F'' est une [[Continuité (mathématiques)|fonction continue]] sur un ouvert ''U'' de ℝ × ''E''{{exp|''n ''+ 1}}, appelé domaine.

L'ordre de cette équation différentielle est l'ordre ''n'' de la plus haute dérivée y apparaissant. Soient ''y'' une fonction de ''x'' définie d'un intervalle ''I'' dans ''E'' et ''y<nowiki>'</nowiki>'', ''y"'', …, ''y''{{exp|(''n'')}} les [[dérivée]]s successives de la fonction ''y''. Cette fonction ''y'' est dite solution si elle est [[Classe de régularité|de classe C{{exp|''n''}}]] et si
<center><math>\forall x \in I, \qquad F(x,y(x),y'(x), \dots, y^{(n)}(x))=0.</math></center>

Résoudre une équation différentielle revient à trouver les fonctions solution ''y''. Par exemple, l'équation différentielle ''y" + y ''= 0 a une solution générale de la forme : ''y''(''x'') = ''A'' cos ''x + B'' sin ''x'',
où ''A, B'' sont des constantes complexes (qu'on peut déterminer si on ajoute des conditions initiales).

Dans une équation différentielle, la fonction ''y'' peut être par exemple à valeurs dans un [[espace vectoriel de dimension finie]], ainsi si ''y'' a pour composantes ''y<sub>1</sub>'' et ''y<sub>2</sub>'' :
:<math>\begin{cases} y'_1&=y_1-2xy_2 +x^2\\
y'_2&=xy_1-2y_2.\end{cases}</math>

L'usage en physique est de parler alors de système d'équations différentielles couplées. Mais le point de vue fécond en mathématiques est de n'y voir qu'une seule équation, pour une [[fonction à valeurs vectorielles]].

On peut encore élargir la définition, en considérant des équations différentielles sur des [[variété différentielle|variétés différentielles]].

=== Équation différentielle sous forme résolue ===

Une équation différentielle d'ordre ''n'' est mise sous forme résolue quand on peut exprimer la dérivée la plus forte en fonction de ''x'' et des dérivées précédentes
<center><math>y^{(n)}=G(x,y,y', \dots, y^{(n-1)}).</math></center>

=== Réduction à 1 de l'ordre d'une équation ===

Une équation différentielle d'ordre ''n''<center>
<math>F(x,y,y',\dots,y^{(n)})=0</math></center>peut se lire aussi comme une équation du premier ordre à valeurs dans ''E{{exp|n}}'', de fonction inconnue {{nobr|''v''(''x'') {{=}} (''y''{{ind|0}}(''x''), …, ''y''{{ind|''n ''– 1}}(''x'')).}} L'équation se réécrit en effet, en notant ''y ''= ''y''{{ind|0}} :
<center><math>y_0'=y_1,\ldots,y_{n-2}'=y_{n-1},F(x,y_0,y_1,\dots,y_{n-1},y_{n-1}')=0</math></center>
ou encore, en définissant ''f ''par ''f''(''x'', ''v''{{ind|0}}, …, ''v''{{ind|''n ''– 1}}, ''w''{{ind|0}}, …, ''w''{{ind|''n ''– 1}}) = (''w''{{ind|0}} – ''v''{{ind|1}}, …, ''w''{{ind|''n ''– 2}} – ''v''{{ind|''n ''– 1}}, ''F''(''x'', ''v''{{ind|0}}, …, ''v''{{ind|''n ''– 1}}, ''w''{{ind|''n ''– 1}})) :

<center><math>f(x,v,v')=0.</math></center>

Si l'équation d'ordre ''n ''était sous forme résolue
<center><math>y^{(n)}=G(x,y,y', \dots, y^{(n-1)}),</math></center>
l'équation équivalente d'ordre 1 le sera aussi :
<center><math>v'=g(x,v)</math></center>
avec ''g''(''x'', ''v''{{ind|0}}, …, ''v''{{ind|''n ''– 1}}) = (''v''{{ind|1}}, …, ''v''{{ind|''n ''– 2}}, ''G''(''x'', ''v''{{ind|0}}, …, ''v''{{ind|''n ''– 1}})).

De plus, dans les deux cas (forme implicite ou forme résolue), si l'équation d'ordre ''n ''était [[#Équation différentielle autonome|autonome]], celle d'ordre 1 le sera aussi (c'est-à-dire que si ''F ''ou ''G ''ne dépend pas de la variable ''x ''alors ''f ''ou ''g ''non plus) et si l'équation était [[#Équation différentielle linéaire|linéaire]], elle le reste. Par exemple, l'[[Équation différentielle linéaire#Réduction à 1 de l'ordre|équation différentielle linéaire]] d'ordre 2, résolue et autonome
<center><math>\frac {\mathrm d^2 y}{\mathrm dx^2} = y</math></center>
se transforme en équation du premier ordre à valeurs dans ℝ{{2}} : la fonction inconnue de la nouvelle équation différentielle est une fonction ''x ''↦ ''v''(''x'') = (''y''(''x''), ''z''(''x'')) de ℝ dans ℝ{{2}} et l'équation s'écrit :
<center><math>\frac {\mathrm d v}{\mathrm dx} = g(v)\text{ avec }g(y,z)=(z,y).</math></center>

== Solutions ==
=== Durée de vie ===

Si ''y'' est solution d'une équation différentielle sur l'intervalle I, on peut considérer sa [[Fonction (mathématiques)|restriction]] à un intervalle J inclus dans I. Celle-ci restera solution de l'équation différentielle. Une solution est encore appelée courbe intégrale.

Il est souvent judicieux de ne considérer que les solutions maximales, encore appelées courbes intégrales maximales, c'est-à-dire celles qui ne sont les restrictions d'aucune autre. L'intervalle de définition est nommé intervalle maximal.

Il ne faut pas croire pour autant que les solutions maximales sont définies sur ℝ entier. Il est tout à fait possible qu'elles aient une durée de vie finie dans le futur ou dans le passé. Il en est ainsi des solutions de l'équation ''y' = y''{{2}}, par exemple.

Cependant, si une solution reste confinée dans un domaine [[Compacité (mathématiques)|compact]], alors elle a une durée de vie infinie.

==== Exemple ====

L'équation différentielle scalaire d'ordre 1 sous forme résolue : ''y' = G''(''x'', ''y''), admet une interprétation géométrique simple dans le plan ramené à un repère d'axes (''Ox''), (''Oy''). On peut représenter, attaché en chaque point de coordonnées ''x'', ''y'', le vecteur de composantes 1 et ''G''(''x'', ''y''), ce qui constitue un [[champ de vecteurs]] du plan. Les courbes solutions sont des représentations graphiques de fonctions {{nobr|1=''y = f''(''x''),}} continûment dérivables, dont la tangente en chaque point est donnée par le champ de vecteurs.

==== Forme résolue et forme implicite ====

Les équations différentielles qui peuvent se mettre sous forme résolue jouissent de bonnes propriétés théoriques, avec – sous certaines hypothèses – un [[théorème d'existence]] et d'unicité de solutions : le [[#Conditions initiales, théorème de Cauchy-Lipschitz|théorème de Cauchy-Lipschitz]].

Dans le cas contraire, on dit que l'équation différentielle est sous forme implicite. On essaie, sur les domaines les plus grands possibles, de mettre l'équation différentielle sous forme résolue. Puis on doit procéder au raccordement des solutions obtenues. Le traitement des équations différentielles de ce type sera évoqué en fin d'article.

=== Conditions initiales, théorème de Cauchy-Lipschitz ===

{{Article détaillé| Théorème de Cauchy-Lipschitz}}

Une [[condition initiale]] (ou condition de Cauchy), pour une équation d'ordre ''n'' d'inconnue ''y'' est la donnée d'une valeur ''x''<sub>0</sub> et de ''n'' vecteurs ''Y''<sub>0</sub>,..., ''Y''<sub>''n''-1</sub>. La fonction solution ''y'' satisfait ces conditions initiales si
:<math>y(x_0)=Y_0, \qquad y'(x_0)=Y_1, \qquad \dots \qquad y^{(n-1)}(x_0)=Y_{n-1}.</math>

Un [[problème de Cauchy]] est la donnée d'une équation différentielle avec un jeu de conditions initiales.

Pour une équation différentielle ''sous forme résolue'', moyennant une certaine hypothèse de régularité (caractère [[Propriété locale|localement]] [[application lipschitzienne|lipschitzien]] à ''x'' fixé, par rapport au bloc des autres variables), le [[théorème de Cauchy-Lipschitz]] énonce que, pour chaque condition initiale :
* il existe une solution qui la satisfait et qui est définie sur un intervalle de la forme ]''x''{{ind|0}} – α, ''x''{{ind|0}} + α[ ;
* il existe une unique solution maximale qui la satisfait.

=== Conditions aux limites ===

Un autre problème classique est celui des [[conditions aux limites]], pour lequel on prescrit les valeurs d'une fonction solution en plusieurs points, voire les valeurs des limites d'une fonction solution aux bornes du domaine. Ainsi le problème :
:<math>\begin{cases} y''+y=0\\ y(0)=y(2\pi)=0.\end{cases}</math>

Un tel problème (parfois appelé [[problème de Dirichlet]]) peut très bien n'avoir aucune solution ou au contraire une infinité de fonctions solutions.

Dans l'exemple donné, les solutions sont les fonctions de la forme ''x'' ↦ ''k'' sin(''x''), pour toute constante ''k''.

=== Résolution explicite ===

{{Article détaillé|Résolution d'équations différentielles par quadrature}}

La résolution par quadrature qui consiste à obtenir une forme explicite des solutions de l'équation différentielle à l'aide des fonctions usuelles et de l'opérateur de primitivation, est rarement possible. Un petit nombre d'équations possédant des formes particulières peuvent être ramenées par [[Changement de variable (simplification algébrique)|changements de variables]] successifs à l'équation la plus simple de toutes : l'équation <math>y'=f</math>, qui est une simple primitivation.

Parmi les équations différentielles pouvant être entièrement résolues sont les [[Équation différentielle linéaire d'ordre un|équations linéaires scalaires d'ordre 1]], les [[Équation différentielle d'ordre un à variables séparées|équations à variables séparées]], les [[Équation différentielle homogène#Équation différentielle du premier ordre, homogène de degré n|équations homogènes du premier ordre]], l'[[Équation différentielle de Bernoulli|équation de Bernoulli]], les équations différentielles vectorielles à coefficients constants.

D'autres peuvent être entièrement résolues dès lors qu'une solution particulière est connue, ainsi l'[[équation différentielle linéaire d'ordre 2]], l'[[Équation de Riccati|équation différentielle de Riccati]].

À défaut de résolution par quadrature, il est parfois possible de rechercher une expression au moins locale des solutions, sous forme de [[série entière]]. Cette démarche, employée de façon systématique sur certaines classes d'équations différentielles linéaires, porte le nom de [[méthode de Frobenius]].

== Propriétés de continuité des solutions ==
=== Continuité par rapport aux conditions initiales et aux paramètres ===
{{Article détaillé|Flot (mathématiques)}}
La donnée des conditions initiales ''x''<sub>0</sub>, ''Y''<sub>0</sub>, … , ''Y''<sub>''n''–1</sub> définit une unique fonction solution qu'on peut noter ''S''(''x''<sub>0</sub>, ''Y''<sub>0</sub>, … , ''Y''<sub>''n''–1</sub>, ''x''). On définit ainsi une fonction globale ''S'' qui prend le nom de [[Flot (mathématiques)|flot]]  (parfois appelé coulée ou encore courant) et qui rend compte de la façon dont les solutions varient avec les conditions initiales. Son domaine d'existence est un [[topologie|ouvert]].

En se plaçant dans les hypothèses du [[théorème de Cauchy-Lipschitz]], les solutions dépendent continûment des conditions initiales, c'est-à-dire que la fonction ''S'' est une fonction continue de l'ensemble de ses variables.

Si on fait dépendre continûment le système d'un paramètre λ, il y a également continuité de ''S'' par rapport à ce paramètre. En effet ajouter un paramètre peut se ramener à modifier le système. Il suffit de rajouter une composante Λ à la fonction cherchée, et lui demander de vérifier l'équation Λ' = 0 et la condition initiale Λ(''x''{{ind|0}}) = λ.

=== Propriétés globales ===

Soit ''y'' une solution particulière de l'équation différentielle, avec pour conditions initiales ''x''<sub>0</sub>, ''Y''<sub>0</sub>, … , ''Y''<sub>''n''–1</sub>. La propriété de continuité permet de donner le comportement des solutions correspondant à des conditions initiales voisines.

* Si l'on restreint la solution à un [[Intervalle (mathématiques)#Connexité et compacité|segment]] [''x<sub>i</sub>'', ''x<sub>f</sub>''] contenant ''x<sub>0</sub>'', les solutions de conditions initiales voisines forment un tube de solutions autour de la solution ''y''.<br />Plus précisément, pour tout ε > 0, il existe η > 0 tel que si ''z'' est solution avec des conditions initiales ''x<sub>0</sub>'', ''Z''<sub>0</sub>, … , ''Z''<sub>''n''–1</sub> et les ''Z<sub>i</sub>'' η-proches des ''Y<sub>i</sub>'', alors la solution ''z'' est tracée dans un [[voisinage tubulaire]] de ''y'', de rayon ε.<br />Par conséquent, si on prend une suite ''z<sub>n</sub>'' de telles solutions, dont les conditions initiales tendent vers celles de ''y'', la suite ''z<sub>n</sub>'' [[convergence uniforme|converge uniformément]] vers ''y''.
* Si l'on étudie la solution sur tout son domaine d'existence, une telle propriété n'est plus vérifiée.

=== Stabilité des solutions ===
{{Article détaillé|Stabilité de Lyapunov}}
{{Quoi|date=janvier 2015|La solution (ℝ, 0)}} de l'équation différentielle ''x' = f''(''t'', ''x'') est stable s'il existe une [[fonction de Liapounov]].

=== Effet papillon, chaos ===
{{Article détaillé|Théorie du chaos}}
Les propriétés de continuité précédentes sont à manier avec précaution, puisqu'elles n'apportent pas d'information quantifiée. Dans la pratique, on observe dans de nombreux systèmes une [[sensibilité aux conditions initiales|sensibilité]] extrême sur le long terme à de petites variations initiales, phénomène popularisé par [[Edward Lorenz]] sous le nom d'[[effet papillon]]. Pour rendre compte de façon satisfaisante de l'évolution d'un [[système physique]] sur un temps très long, il faudrait pousser les mesures de conditions initiales jusqu'à une précision inenvisageable. Ainsi il faudrait englober dans le calcul des [[prévision numérique du temps|prévisions météorologiques]] de très long terme jusqu'aux battements des ailes de papillon.

Les systèmes régis par des équations différentielles, bien qu'étant en principe déterministes, peuvent arborer des comportements extrêmement complexes et paraissant désordonnés, chaotiques. [[Henri Poincaré]] fut le premier à éclaircir cette notion de [[théorie du chaos|chaos]] déterministe. Ses idées tarderont à être reprises, mais servent maintenant de fondement à la théorie des [[systèmes dynamiques]].

== Cas particuliers ==
=== Équation différentielle autonome ===

[[Fichier:Champ tangentes.jpg|vignette|Tracé d'une courbe intégrale suivant un champ de vecteur ; cas ''y' = y''.]]
{{Article détaillé|Équation différentielle autonome|Champ vectoriel}}
Un cas particulier important est celui où la variable n'apparaît pas dans l'équation fonctionnelle, alors qualifiée d'autonome : ainsi l'équation ''y' = f''(''y'').

Les lois de la physique s'appliquent en général à des fonctions du temps, et se présentent sous forme d'équations différentielles autonomes, ce qui manifeste l'invariance de ces lois dans le temps. Ainsi si un système autonome revient à sa position initiale au bout d'un intervalle de temps ''T'', il connaît dès lors une évolution périodique de période ''T''.

L'étude des équations autonomes du premier ordre est équivalente à la recherche des [[courbe intégrale|courbes intégrales]] pour les [[champs de vecteurs]] qui sont tangentes au champ de vecteurs en chaque point. Le [[théorème de Cauchy-Lipschitz]] montre que ces courbes solutions ne se coupent pas et remplissent l'espace. Le [[théorème de Poincaré-Bendixson]] en décrit le comportement asymptotique, dans le plan : une courbe qui ne part pas à l'infini à tendance à s'enrouler autour d'un ensemble [[attracteur]] simple (réunion de courbes). Mais à partir de la dimension 3, il peut y avoir des phénomènes plus complexes, avec des attracteurs étranges<ref>{{ouvrage|langue=en|auteur1=[[Vladimir Arnold]]|auteur2={{Lien|Yulij Ilyashenko}}|titre=Ordinary differential equations and smooth dynamical systems|sous-titre=Ordinary differential equations|lieu=Heidelberg|éditeur=Springer|année=1997}}, {{p.|37}}</ref>.

=== Équation différentielle linéaire ===
{{Article détaillé|Équation différentielle linéaire}}
Une équation différentielle est dite linéaire quand l'expression de l'équation est [[équation linéaire|linéaire]] (ou plus généralement affine) relativement au bloc de variables <math>(y,y',...y^{(n)})</math>. Une équation différentielle linéaire scalaire d'ordre {{math|''n''}} et d'inconnue {{math|''y''}} est donc de la forme
: <math>a_0 y + a_1 y' + a_2 y'' + ... + a_n y^{(n)}= a_{n+1} \,</math>
où <math>a_0,a_1,\dots,a_n,a_{n+1}</math> sont des fonctions numériques.

Une équation différentielle linéaire vectorielle d'ordre ''n'' a le même aspect, en remplaçant les <math>a_i</math> par des applications linéaires (ou souvent des matrices) fonctions de ''x''. Une telle équation est parfois aussi appelée [[système différentiel]] linéaire.

Particularités des équations différentielles linéaires sous forme résolue :
* les solutions ont une durée de vie infinie ;
* on peut superposer (faire des [[Combinaison linéaire|combinaisons linéaires]]) de solutions d'équations différentielles linéaires ;
* quand l'équation est homogène ({{math|''a''{{ind|''n'' + 1}} {{=}} 0}}), son ensemble de solutions est un espace vectoriel de dimension ''n'' fois la dimension de ''E'' ;
* il suffit donc d'exhiber un nombre suffisant de solutions indépendantes de l'équation homogène pour la résoudre. On peut tester l'indépendance de solutions à l'aide du [[wronskien]] ;
* l'ensemble des solutions de l'équation générale est un [[espace affine]] : la solution générale est formée de la somme de cette solution particulière avec la solution générale de l'équation linéaire homogène associée ;
* la [[méthode de variation des constantes]] permet, une fois résolue l'équation homogène, de résoudre l'équation complète ;
* dans le cas d'équations à coefficients constants, on dispose de formules de résolution explicites à l'aide d'[[exponentielle d'une matrice|exponentielles]] de matrices ou d'endomorphismes, ou encore en utilisant la [[Application de la transformée de Laplace aux équations différentielles|transformation de Laplace]].

== Équation différentielle holomorphe ==
{{Article détaillé|Équation différentielle holomorphe}}
Une [[équation différentielle holomorphe]] est l'homologue, pour la variable complexe, d'une équation différentielle ordinaire. La théorie générale en est beaucoup plus complexe.

=== Résultats locaux ===

Une équation différentielle holomorphe sous forme résolue vérifie l'analogue du théorème de Cauchy-Lipschitz : existence et unicité locales d'une fonction solution, elle-même holomorphe.

En outre si l'équation dépend de paramètres de façon holomorphe, la solution aussi. Il y a aussi dépendance holomorphe en les conditions initiales.

Cependant il n'y a plus en général raccordement en une unique solution maximale.

=== Résultats globaux ===

On connaît des difficultés même pour l'équation différentielle la plus simple : le calcul de [[primitive]]s. Par exemple la construction d'une fonction telle que le [[logarithme complexe]] n'est pas univoque. On peut chercher à construire des déterminations de la fonction logarithme sur les ouverts les plus grands possibles : par exemple des plans fendus. On peut aussi construire une primitive « le long d'un chemin ». Apparaît alors le phénomène de [[monodromie]] : si le chemin fait un tour dans le sens direct autour de l'origine, la primitive est modifiée d'une constante (2iπ). Pour rendre compte de la situation, il faut faire intervenir les concepts de [[revêtement (mathématiques)|revêtement]], [[point de branchement]].

Les [[Fonction puissance|fonctions puissances]] sont également solutions d'équations différentielles simples et susceptibles de présenter de la monodromie. Ainsi l'équation ''z' = –z''{{3}} n'admet aucune solution non nulle holomorphe, ni même [[méromorphe]] sur le plan entier.

=== Cas linéaire ===

La théorie des équations différentielles holomorphes linéaires sous forme résolue est très semblables à celle des équations pour la variable réelle, tant qu'on reste sur des domaines [[simplement connexe]]s. Sinon elle donne également lieu à des problèmes de type [[point de branchement]].

== Méthodes numériques ==

La résolution des équations différentielles par quadrature (c'est-à-dire à l'aide des [[Opération élémentaire|opérations élémentaires]] et de la primitivation) n'est possible que dans un nombre de cas très restreints. Par exemple, même les équations différentielles linéaires scalaires d'ordre 2 n'admettent pas de telle formule de résolution générale. Il est donc indispensable de disposer de techniques de résolution approchée.

=== Méthode d'Euler ===
{{Article détaillé|Méthode d'Euler}}

Cette méthode, la plus ancienne et la plus simple, possède également un intérêt théorique puisqu'elle permet de prouver un résultat d'existence de solutions sous des hypothèses plus faibles que le théorème de Cauchy-Lipschitz : c'est le [[théorème de Cauchy-Peano-Arzela]].

On considère une équation différentielle d'ordre 1 sous forme résolue ''y' = f''(''x'', ''y''), avec la condition initiale ''y''(''x''<sub>0</sub>) = ''y''<sub>0</sub>.

Le principe est d'approcher la solution ''y'' sur [''a'', ''b''] par une [[fonction affine par morceaux]], en opérant une discrétisation du paramètre : on pose
:<math>x_i = a+ih</math> où <math>h=(b-a)/n</math> est le pas.

La fonction affine par morceaux joindra donc les points de coordonnées (''x<sub>i</sub>'', ''y<sub>i</sub>''), et il s'agit de proposer un algorithme pour construire les ''y<sub>i</sub>'' à partir de ''y''<sub>0</sub>. Sur chaque intervalle [''x<sub>i</sub>'', ''x''<sub>''i''+1</sub>] on prend pour pente du segment affine celle que suggère l'équation : ''f''(''x<sub>i</sub>'', ''y<sub>i</sub>'').

=== Autres méthodes ===

{{Article détaillé|Résolution numérique des équations différentielles}}

Les plus classiques sont les [[méthodes de Runge-Kutta]], la [[méthode de Newmark]], la [[méthode des différences finies]] ou la [[méthode des éléments finis]] qui est plus adaptée pour les E.D.P.

== Équation différentielle sous forme implicite ==
=== Traitement d'un exemple ===

[[Fichier:Diffeq.svg|vignette|Tracé de quelques solutions de <math>(y')^2 + xy' - y = 0 \,</math>. En bleu les solutions régulières, en vert la solution singulière, en rouge la solution hybride mentionnée dans le texte.]]

Soit l'équation différentielle implicite
:<math>(y')^2 + xy' - y = 0.</math>

Pour l'étudier on effectue un régionnement du [[Plan (mathématiques)|plan]] : on distingue les valeurs (''x'', ''y'') pour lesquelles l'équation <math>T^2 + xT - y = 0</math> admet 0, 1 ou 2 solutions. On obtient trois régions ''U'', ''V'', ''W''. La région ''V'' est la [[parabole]] d'équation <math>y=-\frac14 x^2</math>, les régions ''U'' et ''W'' sont les deux [[Ensemble ouvert|ouverts]] qu'elle délimite.

On commence par s'intéresser aux solutions qui ne sont tracées que sur un des trois domaines
# Dans la région ''U'', l'équation n'admet aucune solution.
# Il y a une solution tout entière tracée dans ''V'', c'est la solution singulière <math>y=-\frac14 x^2</math> tracée en vert ci-contre.
# Dans l'ouvert ''W'', l'équation peut être mise sous l'une des deux formes résolues
:<math>y'=\frac{-x+\sqrt{x^2+4y}}2\qquad y'=\frac{-x-\sqrt{x^2+4y}}2.</math>
Chacune de ces deux équations vérifie le [[théorème de Cauchy-Lipschitz]]. Si on se restreint à l'ouvert ''W'', il y a donc exactement deux solutions pour chaque couple de solutions initiales. Elles sont tracées en bleu sur la figure ci-contre. Dans le cas présent il s'agit d'ailleurs de droites, d'équation
:<math>y = Ax + A^2.</math>
Elles sont tangentes à la parabole d'équation <math>y=-\frac14 x^2</math>. Plus précisément, les solutions tracées sur ''W'' sont ces droites, arrêtées au point de tangence puisqu'on sort de ''W''.

On peut maintenant faire l'étude de l'équation différentielle sur le plan entier. Il existe alors des solutions « hybrides » formées en raccordant de façon C{{exp|1}} un arc de parabole (verte) avec les solutions rectilignes (bleues). Ainsi la solution représentée en rouge :
: <math>y : x \mapsto \begin{cases} x + 1, & \mbox{si } x < -2 \\ - \frac{1}{4} x^2, & \mbox{si } -2 \le x < 2 \\ -x + 1, & \mbox{si } 2 \le x.\end{cases} </math>

Un tel raccordement ne peut se faire qu'en un point de ''V''. La description de l'ensemble de toutes les solutions se ferait en discutant en fonction de la condition initiale ''x<sub>0</sub>, y<sub>0</sub>''
# Condition initiale dans ''U'' : pas de solution
# Condition initiale dans ''V'' : pour les valeurs de ''x'' supérieures à ''x<sub>0</sub>'', la solution peut être la parabole entière, ou alors on suit un arc de parabole puis on bifurque sur la tangente. De même pour les valeurs inférieures à ''x<sub>0</sub>''.
# Condition initiale dans ''W'' : il y a d'abord deux droites solutions, tangentes à la parabole. Soit on les prolonge indéfiniment, soit on les quitte pour la parabole au niveau du point de tangence. On continue alors sur la parabole, ou on repart sur une tangente un peu plus loin.

=== Généralisation ===

Pour généraliser cette étude il faut se placer dans un espace à trois dimensions, de coordonnées notées (''x'', ''y'', ''p''). À l'équation différentielle est associée la [[Surface (géométrie)|surface]] d'équation ''F''(''x'', ''y'', ''p'') = 0 (la coordonnée ''p'' permet de représenter ''y' ''). Les solutions sont des courbes tracées sur la surface. Les difficultés rencontrées viennent de ce que ces courbes sont projetées sur le plan (''x'', ''y''). L'application de [[Projection orthogonale|projection]] connaît des [[point critique (mathématiques)|points critiques]] aux points où le [[gradient]] de ''F'' est « vertical ». Ce sont ces points qui se projettent en la parabole verte.

Finalement, le cadre d'étude des ''équations différentielles implicites'' est le même que celui de la théorie des [[enveloppe (géométrie)|enveloppes]]. La parabole, solution singulière est ici l'enveloppe de la famille des droites, solutions régulières.

== Notes et références ==
{{Références}}

== Voir aussi ==
{{Autres projets
  | wikiversity = Équation différentielle
  | wikt = équation différentielle
}}

=== Bibliographie ===
* {{en}} Thomas Archibald, [[Craig Fraser]] et [[Ivor Grattan-Guinness]] (éd.), ''The History of Differential Equations 1670 - 1950'', [[Institut de recherches mathématiques d'Oberwolfach|MFO]] Report 51/2004
* [[Vladimir Arnold]], ''Équations Différentielles Ordinaires'', [[Éditions Mir|Mir]], Moscou, 1974
* [[Sylvie Benzoni|Sylvie Benzoni-Gavage]], ''Calcul différentiel et équations différentielles'', [[Dunod]], 2014
* Vladimir Damgov, Nonlinear and parametric phenomena. Applications to radiometric and mechanical systems, [[World Scientific]], Series on Nonlinear Sciences, 2004
* [[Jean-Pierre Demailly]], ''Analyse numérique et équations différentielles'', {{3e|éd.}}, [[EDP Sciences]], 2006
* [[Jean Dieudonné]], ''[[Éléments d'analyse]]'', 9 tomes, [[Gauthier-Villars]], 1960-1982
* Jean Dieudonné, ''Calcul infinitésimal'', {{2e|éd.}}, [[Hermann (éditions)|Hermann]], 1980
* Abdelhafed Elkhadiri, ''Calcul différentiel et équations différentielles - Licence de mathématiques'', Ellipses, 2015
* {{en}} [[Andrew Russell Forsyth]], ''Theory of Differential Equations'', [[Cambridge University Press]] (6 vol.), 2012
* Jean-Baptiste Hiriart-Urruty, ''Les équations différentielles pour les débutants'', H & K, 2013
* [[John H. Hubbard]], Beverly H. West, V. Gautheron (trad.), ''Équations différentielles et systèmes dynamiques'', Cassini, 1999
* {{en}} {{Lien|langue=en|fr=Edward Lindsay Ince}}, ''Ordinary Differential Equations'', [[Dover Publications|Dover]], New-York, 1926{{Commentaire biblio|Très utile pour la théorie des équations différentielles linéaires  avec des conditions aux limites.}}
* [[Camille Jordan (mathématicien)|C. Jordan]], ''Cours d'analyse de l'école polytechnique, tome 3, équations différentielles'', {{3e|éd.}}, Gauthier-Villars, 1959
* [[François Laudenbach]], ''Calcul différentiel et intégral'', Éd. École Polytechnique, 2011, {{Google Livres|Ws9A7ZoRJNEC|aperçu}}
* Ahmed Lesfari, ''Équations différentielles ordinaires et équations aux dérivées partielles - Cours et exercices corrigés'', [[Éditions Ellipses|Ellipses]], 2015
* Pierre Meunier, ''Les équations différentielles'', [[Éditions Cépaduès|Cépaduès]], 2014
* [[Lev Pontriaguine]], ''Équations différentielles ordinaires'', Mir, Moscou, 1969
* Hervé Reinhard, ''Équations différentielles : fondements et applications'', {{2e|éd.}}, Dunod, 1989
* N. Rouche et [[Jean Mawhin|J. Mawhin]], ''Équations différentielles ordinaires'' (vol. 1 : Théorie générale, vol. 2 : Stabilité et solutions périodiques), [[Éditions Masson|Masson]], Paris, 1973
* Robert Roussarie et Jean Roux, ''Des équations différentielles aux systèmes dynamiques'' (2 vol.), EDP Sciences, 2012
* Claude Wagschal, ''Fonctions holomorphes, équations différentielles'', Hermann, 2003

=== Articles connexes ===
* [[Équation différentielle]]
* [[Parcours des équations différentielles|Histoire des équations différentielles]]
* [[Équation différentielle linéaire]]
* [[Exemples d'équations différentielles]]
* [[Équation aux dérivées partielles]]
* [[Application de la transformée de Laplace aux équations différentielles]]
* [[Système oscillant à un degré de liberté]]
* ''[[Differential Equations]]''

=== Lien externe ===
* [https://www.canal-u.tv/video/les_amphis_de_france_5/les_equations_differentielles_lineaires_dans_le_champ_complexe.293 Les équations différentielles linéaires dans le champ complexe] (vidéo)

{{Palette|Équations différentielles}}
{{Portail|analyse}}

[[Catégorie:Équation différentielle|Ordinaire]]