ATLAS Ouvrir dans l'atlas →

Proposition 1Équations différentielles

Domaine du flot, continuité et réciproque de φₜ

Soit UU un ouvert de ℝn{\mathbb R^n}, X:U→ℝnX : U \to \mathbb R^n de classe C1C^1 et φ\varphi son flot (Flot), la courbe intégrale issue de xx étant définie sur ]α(x),β(x)[]\alpha(x), \beta(x)[. Le domaine du flot est D={(t,x)∈ℝ×U:α(x)<t<β(x)},D = \{(t, x) \in \mathbb R \times U : \alpha(x) < t < \beta(x)\}, et, pour tt réel, Dt={x∈U:(t,x)∈D}D_t = \{x \in U : (t, x) \in D\} est l'ensemble des points où φt\varphi_t est défini.

On munit ℝn\mathbb R^n de ∥⋅∥∞\|\cdot\|_\infty, et les matrices de la norme subordonnée ∣∣∣⋅∣∣∣|||\cdot|||. Pour une partie AA de ℝn\mathbb R^n, AδA^\delta est la réunion des boules fermées de rayon δ\delta centrées en un point de AA (Sortie des compacts).

Énoncé

(i) Renversement du temps. Le champ −X-X est de classe C1C^1, et son flot est ψt(x)=φ−t(x)\psi_t(x) = \varphi_{-t}(x), défini pour −β(x)<t<−α(x)-\beta(x) < t < -\alpha(x).

(ii) Estimation uniforme. Soit x0∈Ux_0 \in U et [a,b][a, b] un segment contenant 00 et contenu dans ]α(x0),β(x0)[]\alpha(x_0), \beta(x_0)[. Posons T=max⁡(−a,b)T = \max(-a, b) et Γ={φt(x0):a⩽t⩽b}\Gamma = \{\varphi_t(x_0) : a \leqslant t \leqslant b\}. Il existe δ>0\delta > 0 tel que K=ΓδK = \Gamma^\delta soit un compact contenu dans UU ; soit LL le maximum de ∣∣∣dX∣∣∣|||dX||| sur KK.

Alors, pour x∈Ux \in U tel que ∥x−x0∥<δe−LT\|x - x_0\| < \delta e^{-LT} et pour t∈[a,b]t \in [a, b], φt(x)\varphi_t(x) est défini, appartient à KK, et ∥φt(x)−φt(x0)∥⩽∥x−x0∥ eL∣t∣<δ.\|\varphi_t(x) - \varphi_t(x_0)\| \leqslant \|x - x_0\|\,e^{L|t|} < \delta.

(iii) Ouverture et continuité. DD est un ouvert de ℝ×ℝn\mathbb R \times \mathbb R^n, et (t,x)↦φt(x)(t, x) \mapsto \varphi_t(x) est continue sur DD.

(iv) Réciproque. Pour tout réel tt, DtD_t est ouvert, et φt\varphi_t est une bijection de DtD_t sur D−tD_{-t}, de réciproque φ−t\varphi_{-t}.

Démonstration

(i) Si y:J→Uy : J \to U est solution de y′=X(y)y' = X(y), la fonction z(t)=y(−t)z(t) = y(-t), définie sur −J-J, vérifie z′(t)=−y′(−t)=−X(z(t))z'(t) = -y'(-t) = -X(z(t)) par la règle de la chaîne (Différentielle, (ii) et (iv)), et réciproquement.

Le changement de tt en −t-t échange donc les solutions des deux équations et respecte l'inclusion des intervalles, donc la maximalité : la solution maximale de y′=−X(y)y' = -X(y) issue de xx est t↦φ−t(x)t \mapsto \varphi_{-t}(x), sur ]−β(x),−α(x)[]{-\beta(x)}, -\alpha(x)[. Enfin −X-X est de classe C1C^1, de différentielle −dX-dX.

(ii) Les constantes. Γ\Gamma est compact, image du segment [a,b][a, b] (Compacts de ℝ, (i)) par l'application continue t↦φt(x0)t \mapsto \varphi_t(x_0) (Parties compactes, (iii)). Le nœud Sortie des compacts, (i), fournit δ\delta, et LL existe puisque y↦∣∣∣dX(y)∣∣∣y \mapsto |||dX(y)||| est continue sur le compact KK (Compacts de ℝ, (iii)).

Les temps positifs. Le compact K=ΓδK = \Gamma^\delta contient les boules Bˉ(φt(x0),δ)\bar B(\varphi_t(x_0), \delta) pour t∈[0,b]t \in [0, b], et LL est le maximum de ∣∣∣dX∣∣∣|||dX||| sur KK. Comme b⩽Tb \leqslant T, la condition ∥x−x0∥<δe−LT\|x - x_0\| < \delta e^{-LT} entraîne ∥x−x0∥<δe−Lb\|x - x_0\| < \delta e^{-Lb}, et le nœud Flot, (ii), appliqué sur [0,b][0, b] avec ce compact, donne exactement l'estimation cherchée pour t∈[0,b]t \in [0, b], avec φt(x)∈K\varphi_t(x) \in K.

Les temps négatifs. Le champ −X-X a pour flot ψt=φ−t\psi_t = \varphi_{-t}, par (i). Le segment [−b,−a][-b, -a] contient 00 et est contenu dans ]−β(x0),−α(x0)[]{-\beta(x_0)}, -\alpha(x_0)[ ; il donne les mêmes TT, Γ\Gamma, δ\delta et KK, et le même LL puisque ∣∣∣d(−X)∣∣∣=∣∣∣dX∣∣∣|||d(-X)||| = |||dX|||.

Le nœud Flot, (ii), appliqué à −X-X sur [0,−a][0, -a] avec le même compact KK, qui contient les boules Bˉ(ψt(x0),δ)\bar B(\psi_t(x_0), \delta) pour t∈[0,−a]t \in [0, -a], est l'estimation cherchée pour t∈[a,0]t \in [a, 0].

(iii) Soit (t1,x0)∈D(t_1, x_0) \in D. Comme α(x0)<min⁡(t1,0)\alpha(x_0) < \min(t_1, 0) et max⁡(t1,0)<β(x0)\max(t_1, 0) < \beta(x_0), il existe η>0\eta > 0 tel que le segment [a,b]=[min⁡(t1,0)−η,max⁡(t1,0)+η][a, b] = [\min(t_1, 0) - \eta, \max(t_1, 0) + \eta] soit contenu dans ]α(x0),β(x0)[]\alpha(x_0), \beta(x_0)[ ; il contient 00, et t1∈ ]a,b[t_1 \in \,]a, b[. Soit ρ=δe−LT\rho = \delta e^{-LT}, avec les constantes de (ii).

Par (ii), ]a,b[×B(x0,ρ)]a, b[ \times B(x_0, \rho) est contenu dans DD ; c'est un ouvert qui contient (t1,x0)(t_1, x_0), donc DD est ouvert. Pour (t,x)(t, x) dans cet ouvert, ∥φt(x)−φt1(x0)∥⩽∥x−x0∥ eLT+∥φt(x0)−φt1(x0)∥,\begin{aligned} &\|\varphi_t(x) - \varphi_{t_1}(x_0)\| \leqslant \|x - x_0\|\,e^{LT} \\ &\quad + \|\varphi_t(x_0) - \varphi_{t_1}(x_0)\|, \end{aligned} qui tend vers 00 quand (t,x)→(t1,x0)(t, x) \to (t_1, x_0), par continuité de la solution issue de x0x_0.

(iv) DtD_t est l'image réciproque de l'ouvert DD par l'application continue x↦(t,x)x \mapsto (t, x). Soit x∈Dtx \in D_t et y=φt(x)y = \varphi_t(x). La fonction s↦φs+t(x)s \mapsto \varphi_{s+t}(x), définie sur ]α(x)−t,β(x)−t[]\alpha(x) - t, \beta(x) - t[, est solution de l'équation autonome y′=X(y)y' = X(y) et vaut yy en 00.

Par Solutions maximales, (i), elle coïncide avec s↦φs(y)s \mapsto \varphi_s(y) sur l'intersection de leurs intervalles ; la fonction égale à l'une ou à l'autre sur la réunion, qui est un intervalle contenant 00, est une solution qui prolonge la solution maximale issue de yy. Donc ]α(x)−t,β(x)−t[  ⊂  ]α(y),β(y)[]\alpha(x) - t, \beta(x) - t[ \;\subset\; ]\alpha(y), \beta(y)[.

Comme 0∈ ]α(x),β(x)[0 \in \,]\alpha(x), \beta(x)[, le réel −t-t appartient au premier intervalle : y∈D−ty \in D_{-t}, et φ−t(y)=φ−t+t(x)=x\varphi_{-t}(y) = \varphi_{-t+t}(x) = x. Ainsi φ−t∘φt\varphi_{-t} \circ \varphi_t est l'identité de DtD_t ; en changeant tt en −t-t, φt∘φ−t\varphi_t \circ \varphi_{-t} est l'identité de D−tD_{-t}, et φt\varphi_t est une bijection de DtD_t sur D−tD_{-t}.

□
Remarque

Pour l'équation y′=y2y' = y^2 sur U=ℝU = \mathbb R, la solution issue de xx est φt(x)=x1−tx\varphi_t(x) = \frac{x}{1 - tx}, définie tant que tx<1tx < 1 : le domaine du flot est D={(t,x):tx<1}D = \{(t, x) : tx < 1\}.

t x D1/2 D−1/2 1/2 −1/2 2 −2 t x = 1 D D
En gris, la région tx⩾1tx \geqslant 1, hors du domaine DD ; les deux branches de l'hyperbole tx=1tx = 1, en rouge, n'appartiennent pas à DD, qui est ouvert. DD n'est pas un produit : la verticale t=12t = \frac12 le coupe selon D1/2= ]−∞,2[D_{1/2} = \,]{-\infty}, 2[, en bleu, que φ1/2\varphi_{1/2} envoie bijectivement sur D−1/2= ]−2,+∞[D_{-1/2} = \,]{-2}, +\infty[, en vert.

On vérifie en effet que φ1/2(x)=2x2−x\varphi_{1/2}(x) = \frac{2x}{2 - x} croît de −2-2 à +∞+\infty quand xx parcourt ]−∞,2[]{-\infty}, 2[, et que φ−1/2(φ1/2(x))=x\varphi_{-1/2}\big(\varphi_{1/2}(x)\big) = x. L'erreur la plus répandue est de voir dans φt\varphi_t une bijection de UU : elle n'est définie que sur DtD_t, et son image est D−tD_{-t}, qui en diffère en général.

Remarque

Le point (iii) dit que le temps d'existence ne peut pas chuter brusquement : si β(x0)>t\beta(x_0) > t, alors β(x)>t\beta(x) > t pour xx voisin de x0x_0. Il peut en revanche augmenter brusquement, et β\beta n'est pas continue en général.

Pour le champ constant X=(1,0)X = (1, 0) sur U=ℝ2U = \mathbb R^2 privé de la demi-droite {(1,s):s⩾0}\{(1, s) : s \geqslant 0\}, on a β(0,s)=1\beta(0, s) = 1 pour s⩾0s \geqslant 0, mais β(0,s)=+∞\beta(0, s) = +\infty pour s<0s < 0 : la courbe issue d'un point situé sous la demi-droite passe sans la rencontrer.

Ouvrir dans l'atlas interactif →