ATLAS Ouvrir dans l'atlas →

Théorème 6Équations différentielles

Différentiabilité du flot et équation variationnelle

Soit X:U→ℝnX : U \to {\mathbb R^n} de classe C1C^1 sur un ouvert UU, φ\varphi son flot et DD son domaine (Domaine du flot). Pour x∈Ux \in U, la fonction Ax(s)=dX(φs(x))A_x(s) = dX(\varphi_s(x)) est continue sur ]α(x),β(x)[]\alpha(x), \beta(x)[, comme composée de deux applications continues.

L'équation variationnelle le long de la courbe issue de xx est le système linéaire matriciel Z′=dX(φs(x)) Z,Z(0)=I.Z' = dX(\varphi_s(x))\,Z, \qquad Z(0) = I. Résolu colonne par colonne, il a une unique solution ZxZ_x, définie sur ]α(x),β(x)[]\alpha(x), \beta(x)[ tout entier (Systèmes linéaires, (i)). On munit ℝn\mathbb R^n de ∥⋅∥∞\|\cdot\|_\infty, et les matrices de la norme subordonnée ∣∣∣⋅∣∣∣|||\cdot|||.

Énoncé

(i) Différentiabilité. Pour tout (t,x)∈D(t, x) \in D, l'application φt\varphi_t est différentiable en xx, et dφt(x)=Zx(t).d\varphi_t(x) = Z_x(t).

(ii) Classe C1C^1. L'application (t,x)↦Zx(t)(t, x) \mapsto Z_x(t) est continue sur DD, et (t,x)↦φt(x)(t, x) \mapsto \varphi_t(x) est de classe C^1 sur DD, de dérivée partielle X(φt(x))X(\varphi_t(x)) par rapport à tt et de différentielle Zx(t)Z_x(t) par rapport à xx.

(iii) Difféomorphismes. Pour tout réel tt, φt\varphi_t est un C^1-difféomorphisme de DtD_t sur D−tD_{-t}, de réciproque φ−t\varphi_{-t}. En particulier, Zx(t)Z_x(t) est inversible, d'inverse Zφt(x)(−t)Z_{\varphi_t(x)}(-t).

Démonstration

Deux inégalités sur les intégrales. Pour v:[0,t]→ℝnv : [0, t] \to \mathbb R^n continue, on a ∥∫0tv∥⩽∫0t∥v∥\|\int_0^t v\| \leqslant \int_0^t \|v\|, puisque ∣∫0tvi∣⩽∫0t∣vi∣⩽∫0t∥v∥|\int_0^t v_i| \leqslant \int_0^t |v_i| \leqslant \int_0^t \|v\| pour chaque coordonnée (Propriétés de l'intégrale). Pour une fonction matricielle continue MM et ∥w∥⩽1\|w\| \leqslant 1, (∫0tM)w=∫0tMw\big(\int_0^t M\big)w = \int_0^t Mw donne alors ∣∣∣∫0tM∣∣∣⩽∫0t∣∣∣M∣∣∣|||\int_0^t M||| \leqslant \int_0^t |||M|||.

Les constantes. Fixons x0∈Ux_0 \in U et un segment [a,b][a, b] comme au nœud Domaine du flot, (ii), avec les mêmes TT, Γ\Gamma, δ\delta, KK et LL, et posons ρ=δe−LT\rho = \delta e^{-LT}. La fonction ω(r)=sup⁡{∣∣∣dX(y)−dX(y′)∣∣∣:y,y′∈K, ∥y−y′∥⩽r}\begin{aligned} \omega(r) = \sup\big\{&|||dX(y) - dX(y')||| : \\ &y, y' \in K,\ \|y - y'\| \leqslant r\big\} \end{aligned} est bien définie, majorée par 2L2L, et croissante.

Elle tend vers 00 quand r→0r \to 0 : par le théorème de Heine (Théorème de Heine), appliqué à dXdX, continue sur le compact KK, il existe pour tout ε>0\varepsilon > 0 un η>0\eta > 0 tel que ω(r)⩽ε\omega(r) \leqslant \varepsilon pour r<ηr < \eta.

Bornes sur ZZ. Pour ∥x−x0∥<ρ\|x - x_0\| < \rho et s∈[a,b]s \in [a, b], φs(x)∈K\varphi_s(x) \in K, donc ∣∣∣Ax(s)∣∣∣⩽L|||A_x(s)||| \leqslant L. De Zx(t)=I+∫0tAxZxZ_x(t) = I + \int_0^t A_xZ_x (Cauchy-Lipschitz, (ii), forme intégrale) on tire, pour t∈[0,b]t \in [0, b], ∣∣∣Zx(t)∣∣∣⩽1+L∫0t∣∣∣Zx∣∣∣|||Z_x(t)||| \leqslant 1 + L\int_0^t |||Z_x|||, et Grönwall (Grönwall, (i)) donne ∣∣∣Zx(t)∣∣∣⩽eLT|||Z_x(t)||| \leqslant e^{LT}.

Étape 1 : le reste de Taylor. Soit hh tel que 0<∥h∥<ρ0 < \|h\| < \rho, x=x0+hx = x_0 + h et r=∥h∥eLT<δr = \|h\|e^{LT} < \delta. Pour s∈[0,b]s \in [0, b], posons p=φs(x0)p = \varphi_s(x_0) et q=φs(x)q = \varphi_s(x), de sorte que ∥q−p∥⩽r\|q - p\| \leqslant r (Domaine du flot, (ii)).

L'application g(y)=X(y)−dX(p)yg(y) = X(y) - dX(p)y est différentiable sur la boule ouverte V=B(p,min⁡(2r,δ))V = B(p, \min(2r, \delta)), convexe et contenue dans KK, de différentielle dX(y)−dX(p)dX(y) - dX(p) (Différentielle, (iii)), de norme au plus ω(2r)\omega(2r). Comme q∈Vq \in V, l'inégalité des accroissements finis (Accroissements finis, (ii)), appliquée à gg entre pp et qq, donne ∥X(q)−X(p)−dX(p)(q−p)∥⩽ω(2r) r.\begin{aligned} &\|X(q) - X(p) - dX(p)(q - p)\| \\ &\quad \leqslant \omega(2r)\,r. \end{aligned}

Étape 2 : Grönwall. Posons u(t)=φt(x)−φt(x0)−Zx0(t)hu(t) = \varphi_t(x) - \varphi_t(x_0) - Z_{x_0}(t)h et notons R(s)R(s) le vecteur que l'on vient de majorer. Par les formes intégrales de φt(x)\varphi_t(x), de φt(x0)\varphi_t(x_0) et de Zx0Z_{x_0}, et puisque Ax0(s)u(s)+R(s)=X(q)−X(p)−Ax0(s)Zx0(s)hA_{x_0}(s)u(s) + R(s) = X(q) - X(p) - A_{x_0}(s)Z_{x_0}(s)h, u(t)=∫0t(Ax0(s)u(s)+R(s)) ds.u(t) = \int_0^t \big(A_{x_0}(s)u(s) + R(s)\big)\,ds.

Il vient ∥u(t)∥⩽Tω(2r)r+L∫0t∥u∥\|u(t)\| \leqslant T\omega(2r)r + L\int_0^t \|u\| pour t∈[0,b]t \in [0, b], et le lemme de Grönwall donne ∥u(t)∥⩽Te2LTω(2∥h∥eLT) ∥h∥.\|u(t)\| \leqslant Te^{2LT}\omega(2\|h\|e^{LT})\,\|h\|.

Conclusion pour t⩾0t \geqslant 0. Posons ε(h)=u(t)/∥h∥\varepsilon(h) = u(t)/\|h\| pour h≠0h \neq 0 et ε(0)=0\varepsilon(0) = 0. Alors φt(x0+h)=φt(x0)+Zx0(t)h+∥h∥ ε(h)\varphi_t(x_0 + h) = \varphi_t(x_0) + Z_{x_0}(t)h + \|h\|\,\varepsilon(h), et ∥ε(h)∥⩽Te2LTω(2∥h∥eLT)\|\varepsilon(h)\| \leqslant Te^{2LT}\omega(2\|h\|e^{LT}) tend vers 00 avec hh : φt\varphi_t est différentiable en x0x_0, de différentielle Zx0(t)Z_{x_0}(t), pour tout t∈[0,b]t \in [0, b], et la majoration du reste ne dépend pas de tt.

Les temps négatifs. Le champ −X-X a pour flot ψs=φ−s\psi_s = \varphi_{-s} (Domaine du flot, (i)), et W(s)=Zx(−s)W(s) = Z_x(-s) vérifie W′(s)=−Ax(−s)W(s)=d(−X)(ψs(x)) W(s)W'(s) = -A_x(-s)W(s) = d(-X)(\psi_s(x))\,W(s) et W(0)=IW(0) = I : par unicité, c'est la solution de l'équation variationnelle de −X-X.

Le segment [−b,−a][-b, -a] donne pour −X-X les mêmes constantes ; le cas précédent, appliqué à −X-X, donne dφt(x0)=dψ−t(x0)=Zx0(t)d\varphi_t(x_0) = d\psi_{-t}(x_0) = Z_{x_0}(t) pour t∈[a,0]t \in [a, 0].

(i) est ainsi démontré, puisque tout (t,x0)∈D(t, x_0) \in D a son temps tt dans un tel segment (Domaine du flot, preuve de (iii)).

(ii) Continuité de ZZ. Pour ∥h∥<ρ\|h\| < \rho, x=x0+hx = x_0 + h et t∈[0,b]t \in [0, b], la différence V(t)=Zx(t)−Zx0(t)V(t) = Z_x(t) - Z_{x_0}(t) vérifie V(t)=∫0t(AxV+(Ax−Ax0)Zx0) ds.V(t) = \int_0^t \big(A_xV + (A_x - A_{x_0})Z_{x_0}\big)\,ds.

Les points φs(x)\varphi_s(x) et φs(x0)\varphi_s(x_0) sont dans KK à distance au plus r=∥h∥eLTr = \|h\|e^{LT}, d'où ∣∣∣Ax(s)−Ax0(s)∣∣∣⩽ω(r)|||A_x(s) - A_{x_0}(s)||| \leqslant \omega(r). Avec ∣∣∣Zx0∣∣∣⩽eLT|||Z_{x_0}||| \leqslant e^{LT}, on obtient ∣∣∣V(t)∣∣∣⩽Tω(r)eLT+L∫0t∣∣∣V∣∣∣|||V(t)||| \leqslant T\omega(r)e^{LT} + L\int_0^t |||V|||, puis, par Grönwall, ∣∣∣V(t)∣∣∣⩽Te2LTω(∥h∥eLT)|||V(t)||| \leqslant Te^{2LT}\omega(\|h\|e^{LT}). Pour t∈[a,0]t \in [a, 0], la même majoration s'obtient en passant à −X-X, comme plus haut.

Tout (t1,x0)∈D(t_1, x_0) \in D admet un tel segment [a,b][a, b] avec t1∈ ]a,b[t_1 \in \,]a, b[, comme on l'a vu dans la preuve du nœud Domaine du flot, (iii). Pour t∈[a,b]t \in [a, b] et ∥x−x0∥<ρ\|x - x_0\| < \rho, on a donc ∣∣∣Zx(t)−Zx0(t1)∣∣∣⩽Te2LTω(∥x−x0∥eLT)+∣∣∣Zx0(t)−Zx0(t1)∣∣∣,\begin{aligned} &|||Z_x(t) - Z_{x_0}(t_1)||| \\ &\quad \leqslant Te^{2LT}\omega(\|x - x_0\|e^{LT}) \\ &\qquad + |||Z_{x_0}(t) - Z_{x_0}(t_1)|||, \end{aligned} qui tend vers 00 quand (t,x)→(t1,x0)(t, x) \to (t_1, x_0) : (t,x)↦Zx(t)(t, x) \mapsto Z_x(t) est continue sur DD.

Classe C1C^1. Sur l'ouvert DD (Domaine du flot, (iii)), la dérivée partielle de φ\varphi par rapport à tt est X(φt(x))X(\varphi_t(x)), puisque t↦φt(x)t \mapsto \varphi_t(x) est solution ; elle est continue, comme composée de XX et de φ\varphi, continue sur DD.

La dérivée partielle par rapport à xjx_j est Zx(t)ejZ_x(t)e_j, par (i) et Différentielle, (ii) ; elle est continue par ce qui précède. Le nœud Classe Cᵏ, (i), conclut.

(iii) Pour tt fixé, l'application affine x↦(t,x)x \mapsto (t, x) envoie l'ouvert DtD_t dans DD ; par composition (Classe Cᵏ, (iii)), φt\varphi_t est de classe C1C^1 sur DtD_t, et de même φ−t\varphi_{-t} sur D−tD_{-t}. Avec Domaine du flot, (iv), φt\varphi_t est un C1C^1-difféomorphisme de DtD_t sur D−tD_{-t}.

Enfin, la règle de la chaîne (Différentielle, (iv)) appliquée à φ−t∘φt=id\varphi_{-t} \circ \varphi_t = \mathrm{id} et à φt∘φ−t=id\varphi_t \circ \varphi_{-t} = \mathrm{id} montre que Zφt(x)(−t)Z_{\varphi_t(x)}(-t) et Zx(t)Z_x(t) sont inverses l'une de l'autre.

□
Remarque

Le champ X(y1,y2)=(1,y12)X(y_1, y_2) = (1, y_1^2) sur ℝ2\mathbb R^2 a pour flot φt(y)=(y1+t, y2+(y1+t)3−y133).\varphi_t(y) = \Big(y_1 + t,\ y_2 + \frac{(y_1 + t)^3 - y_1^3}{3}\Big). Le long de la courbe issue de 00, φs(0)=(s,s33)\varphi_s(0) = (s, \frac{s^3}3), et dX(φs(0))dX(\varphi_s(0)) a pour seul coefficient non nul 2s2s, en bas à gauche.

L'équation variationnelle a pour solution la matrice Z0(s)Z_0(s) de lignes (1,0)(1, 0) et (s2,1)(s^2, 1), qui est bien la matrice jacobienne de φs\varphi_s en 00. Au temps 11, on calcule exactement le reste : φ1(h)−φ1(0)−Z0(1)h=(0,h12)\varphi_1(h) - \varphi_1(0) - Z_0(1)h = (0, h_1^2), qui est d'ordre 22 en hh, comme (i) le prévoit.

y1 y2 h12 h 0 φ1(0)
Le cercle gris, de centre 00 et de rayon 12\frac12, est emporté par le flot au temps 11 : son image exacte φ1(cercle)\varphi_1(\text{cercle}), en rouge, et l'image φ1(0)+Z0(1)h\varphi_1(0) + Z_0(1)h par l'approximation linéaire, en bleu, sont deux courbes autour de φ1(0)=(1,13)\varphi_1(0) = (1, \frac13). Pour un même hh, l'écart entre le point rouge et le point bleu est vertical et vaut h12⩽14h_1^2 \leqslant \frac14, d'ordre 22 en hh, comme on le voit pour h=(12,0)h = (\frac12, 0) ; la courbe verte est la trajectoire du centre.

La différentielle Z0(1)Z_0(1) est un cisaillement : elle déforme le cercle en ellipse, et le terme h12h_1^2 courbe encore cette ellipse.

Remarque

Pour y′=y2y' = y^2, la dérivée de φt(x)=x1−tx\varphi_t(x) = \frac{x}{1 - tx} par rapport à xx est 1(1−tx)2\frac{1}{(1 - tx)^2} ; on vérifie qu'elle résout Z′=2φt(x) ZZ' = 2\varphi_t(x)\,Z, Z(0)=1Z(0) = 1. Elle tend vers +∞+\infty quand tt s'approche de l'instant d'explosion 1x\frac1x : deux solutions voisines s'écartent d'autant plus vite que l'on s'approche de l'explosion, comme le laisse attendre l'explosion de la solution au nœud Solutions maximales.

Remarque

L'hypothèse C1C^1 ne peut pas être affaiblie en une hypothèse lipschitzienne. Le champ X(y)=∣y∣X(y) = |y| sur ℝ\mathbb R est 11-lipschitzien, et son flot est φt(x)=xet\varphi_t(x) = xe^t pour x⩾0x \geqslant 0 et xe−txe^{-t} pour x<0x < 0.

Pour t≠0t \neq 0, φt\varphi_t a en 00 une dérivée à droite ete^t et une dérivée à gauche e−te^{-t} : il n'est pas dérivable. Le flot reste lipschitzien en xx, de rapport e∣t∣e^{|t|}, par le lemme de Grönwall (Grönwall, (ii)), qui ne demande qu'une hypothèse lipschitzienne ; c'est le reste de Taylor de l'étape 1 qui demande la différentiabilité de XX, et la majoration uniforme de ce reste, la continuité de dXdX.

Remarque

La différentielle du flot transporte le champ : Zx(t)X(x)=X(φt(x))Z_x(t)X(x) = X(\varphi_t(x)). En effet, s↦φt(φs(x))=φt+s(x)s \mapsto \varphi_t(\varphi_s(x)) = \varphi_{t+s}(x) a pour dérivée en s=0s = 0, d'un côté dφt(x)X(x)d\varphi_t(x)X(x) par la règle de la chaîne, de l'autre X(φt(x))X(\varphi_t(x)). Pour y′=y2y' = y^2, on retrouve x2(1−tx)2=φt(x)2\frac{x^2}{(1 - tx)^2} = \varphi_t(x)^2.

Ouvrir dans l'atlas interactif →