ATLAS Ouvrir dans l'atlas →

Théorème 3Équations différentielles

Solutions maximales et critère d'explosion

On suppose désormais FF continue sur l'ouvert Ω\Omega et localement lipschitzienne en yy : pour tout cylindre C=[t1−a,t1+a]×Bˉ(y1,r)C = [t_1 - a, t_1 + a] \times \bar B(y_1, r) contenu dans Ω\Omega, il existe LL tel que ∥F(t,y)−F(t,z)∥⩽L∥y−z∥\|F(t, y) - F(t, z)\| \leqslant L\|y - z\| pour (t,y),(t,z)∈C(t, y), (t, z) \in C.

C'est le cas si FF est de classe C1C^1. La fonction continue ∣∣∣dF∣∣∣|||dF||| est alors bornée par un réel LL sur le compact CC (L'espace ℝⁿ, (iii), Compacts de ℝ, (iii)). Pour tt fixé, z↦F(t,z)z \mapsto F(t, z) est différentiable sur la boule ouverte B(y1,r)B(y_1, r), convexe, de différentielle k↦dF(t,z)(0,k)k \mapsto dF(t, z)(0, k) (Différentielle, (iii) et (iv)), de norme au plus LL. Elle est donc LL-lipschitzienne sur cette boule (Accroissements finis, (ii)), puis sur la boule fermée, par continuité.

Une solution est maximale si elle n'est la restriction d'aucune solution définie sur un intervalle strictement plus grand.

Énoncé

(i) Deux solutions qui coïncident en un point coïncident sur l'intersection de leurs intervalles de définition.

(ii) Par tout point (t0,y0)∈Ω(t_0, y_0) \in \Omega passe une unique solution maximale, définie sur un intervalle ouvert ]α,β[]\alpha, \beta[, et toute solution passant par ce point en est une restriction.

(iii) Explosion. Si Ω=I×ℝn\Omega = I \times {\mathbb R^n} avec II un intervalle ouvert, et si β<sup⁡I\beta < \sup I, alors \|y(t)\| \to +\infty quand t→βt \to \beta. De même, si α>inf⁡I\alpha > \inf I, alors ∥y(t)∥→+∞\|y(t)\| \to +\infty quand t→αt \to \alpha.

(iv) Si ∥F(t,y)∥⩽c(t)(1+∥y∥)\|F(t, y)\| \leqslant c(t)(1 + \|y\|) avec cc continue sur II, les solutions maximales sont définies sur II tout entier.

Démonstration

(i) L'ensemble des tt où elles coïncident est fermé dans l'intervalle commun, par continuité. Il est ouvert dans cet intervalle, par l'unicité locale (Cauchy-Lipschitz, (ii)), appliquée sur un cylindre contenu dans l'ouvert Ω\Omega et centré au point commun, sur lequel FF est bornée, par continuité et compacité, et lipschitzienne en yy par hypothèse. Par connexité, il est égal à l'intervalle entier.

(ii) Par (i), la réunion des graphes de toutes les solutions passant par (t0,y0)(t_0, y_0) est le graphe d'une fonction yy, définie sur la réunion J∗J^* de leurs intervalles. Aucun choix n'intervient : on prend la réunion de toutes.

La fonction yy est une solution. Soit en effet t∈J∗t \in J^*. Si t>t0t > t_0, ou bien J∗J^* contient un s>ts > t, et l'intervalle d'une solution définie en ss contient [t0,s][t_0, s], voisinage de tt dans J∗J^* ; ou bien tt est la borne droite de J∗J^*, et l'intervalle d'une solution définie en tt contient [t0,t][t_0, t]. Le cas t<t0t < t_0 est symétrique, et, pour t=t0t = t_0, la solution fournie par Cauchy-Lipschitz est définie sur un voisinage [t0−h,t0+h][t_0 - h, t_0 + h] de t0t_0. Dans tous les cas, yy coïncide au voisinage de tt dans J∗J^* avec une solution : elle est dérivable en tt, de dérivée F(t,y(t))F(t, y(t)).

Toute solution passant par (t0,y0)(t_0, y_0) est une restriction de yy, par construction ; yy est donc maximale, et c'est la seule solution maximale passant par ce point, puisqu'une autre, restriction de yy, ne pourrait être maximale qu'en lui étant égale. Enfin, si J∗J^* contenait l'une de ses extrémités, Cauchy-Lipschitz en ce point prolongerait la solution : J∗J^* est un intervalle ouvert.

(iii) Supposons β<sup⁡I\beta < \sup I et que ∥y∥\|y\| ne tende pas vers l'infini quand t→βt \to \beta. Il existe alors RR tel que, pour tout η>0\eta > 0, on ait ∥y(t)∥⩽R\|y(t)\| \leqslant R pour un t∈ ]β−η,β[t \in \,]\beta - \eta, \beta[.

Comme β∈I\beta \in I et que II est ouvert, il existe δ>0\delta > 0 tel que [β−2δ,β+δ]⊂I[\beta - 2\delta, \beta + \delta] \subset I. Le compact K=[β−2δ,β+δ]×Bˉ(0,R+1)K = [\beta - 2\delta, \beta + \delta] \times \bar B(0, R + 1) est contenu dans Ω\Omega, et FF y est bornée par un réel M>0M > 0, quitte à l'agrandir. Posons h=min⁡(δ,1/M)h = \min(\delta, 1/M), qui ne dépend que de RR et de δ\delta.

Soit t1∈ ]β−h,β[t_1 \in \,]\beta - h, \beta[ tel que ∥y(t1)∥⩽R\|y(t_1)\| \leqslant R. Comme β−δ<t1<β\beta - \delta < t_1 < \beta, le cylindre [t1−δ,t1+δ]×Bˉ(y(t1),1)[t_1 - \delta, t_1 + \delta] \times \bar B(y(t_1), 1) est contenu dans KK : FF y est bornée par MM et lipschitzienne en yy. Cauchy-Lipschitz, avec a=δa = \delta et r=1r = 1, fournit une solution zz sur [t1−h,t1+h][t_1 - h, t_1 + h] telle que z(t1)=y(t1)z(t_1) = y(t_1).

Par (i), yy et zz coïncident sur l'intersection de leurs intervalles ; la fonction égale à yy sur ]α,β[]\alpha, \beta[ et à zz sur [t1−h,t1+h][t_1 - h, t_1 + h] est une solution, puisqu'au voisinage de chaque point elle coïncide avec l'une ou l'autre. Elle est définie jusqu'en t1+h>βt_1 + h > \beta, ce qui contredit la maximalité de yy.

Du côté de α\alpha, on applique ce qui précède à z(t)=y(−t)z(t) = y(-t), qui vérifie z′(t)=−F(−t,z(t))z'(t) = -F(-t, z(t)) : le second membre F~(t,y)=−F(−t,y)\tilde F(t, y) = -F(-t, y) est continu et localement lipschitzien en yy sur (−I)×ℝn(-I) \times \mathbb R^n, et le changement de tt en −t-t échange les solutions des deux équations en respectant l'inclusion des intervalles, donc la maximalité. La solution maximale zz est définie sur ]−β,−α[]{-\beta}, -\alpha[, avec −α<sup⁡(−I)-\alpha < \sup(-I).

(iv) Sur [t0,β[[t_0, \beta[, on a ∥y(t)∥⩽∥y0∥+∫t0tc(s)(1+∥y(s)∥) ds,\begin{aligned} \|y(t)\| &\leqslant \|y_0\| \\ &\quad + \int_{t_0}^t c(s)(1 + \|y(s)\|)\,ds, \end{aligned} par l'inégalité sur les intégrales vectorielles du nœud Cauchy-Lipschitz. Si β<sup⁡I\beta < \sup I, cc est bornée par CC sur [t0,β][t_0, \beta], et Grönwall (Grönwall) appliqué à 1+∥y∥1 + \|y\| montre que ∥y∥\|y\| reste bornée, ce que (iii) interdit. Du côté de α\alpha, on passe à t↦y(−t)t \mapsto y(-t) comme en (iii), le second membre −F(−t,y)-F(-t, y) vérifiant la même majoration avec la fonction continue t↦c(−t)t \mapsto c(-t).

□
Remarque

L'équation y′=y2y' = y^2, y(0)=1y(0) = 1, a pour solution maximale y(t)=11−ty(t) = \frac{1}{1 - t}, définie sur ]−∞,1[]{-\infty}, 1[ : elle explose en temps fini, bien que F(t,y)=y2F(t, y) = y^2 soit définie et régulière partout.

t y 1 1/2 R = 2 4/5 R = 5 (0, 1) y = 1/(1 − t)
La solution maximale y=11−ty = \frac{1}{1 - t} sort de la bande ∣y∣⩽R|y| \leqslant R à l'instant 1−1R1 - \frac1R, soit 12\frac12 pour R=2R = 2 et 45\frac45 pour R=5R = 5, et n'y revient plus, puisqu'elle est croissante : elle tend vers +∞+\infty quand tt tend vers 11, comme l'annonce (iii).

L'erreur la plus répandue est de croire qu'une équation à second membre partout défini a des solutions partout définies ; c'est la croissance de FF qui en décide, comme le montre (iv).

Ouvrir dans l'atlas interactif →