ATLAS Ouvrir dans l'atlas →

Théorème 6Stabilité des équilibres

Instabilité par linéarisation

On garde les notations du nœud Linéarisation : aa est un équilibre du champ XX, de classe C1C^1 sur l'ouvert UU, A=dX(a)A = dX(a), et μ(A)\mu(A) est la plus grande partie réelle d'une valeur propre complexe de AA. Pour z,w∈ℂnz, w \in \mathbb C^n, ⟨z,w⟩\langle z, w \rangle et ∣z∣|z| sont définis comme au nœud Trigonalisation sur ℂ.

Énoncé

Si μ(A)>0\mu(A) > 0, c'est-à-dire si AA a une valeur propre complexe de partie réelle strictement positive, l'équilibre aa est instable.

Démonstration

Les constantes. Les parties réelles des valeurs propres sont en nombre fini (Trigonalisation sur ℂ, (i)). Parmi les n+1n + 1 réels kμ(A)n+2\frac{k\mu(A)}{n + 2}, 1⩽k⩽n+11 \leqslant k \leqslant n + 1, l'un au moins n'est pas une partie réelle de valeur propre ; soit ss le plus petit de ceux qui ne le sont pas. Alors 0<s<μ(A)0 < s < \mu(A), et κ=min⁡λ∣Re⁡λ−s∣>0\kappa = \min_\lambda |\operatorname{Re}\lambda - s| > 0.

Soit RR fournie par Trigonalisation sur ℂ, (ii) et (iii), avec Λ={λ:Re⁡λ>s}\Lambda = \{\lambda : \operatorname{Re}\lambda > s\} et η=κ4\eta = \frac\kappa4 : R−1AR=D+NR^{-1}AR = D + N, et les coefficients diagonaux d1,…,dnd_1, \dots, d_n de DD vérifient Re⁡dk>s\operatorname{Re} d_k > s pour k⩽qk \leqslant q et Re⁡dk<s\operatorname{Re} d_k < s pour k>qk > q, avec q⩾1q \geqslant 1, puisque toute valeur propre figure sur la diagonale. Posons σk=1\sigma_k = 1 pour k⩽qk \leqslant q, σk=−1\sigma_k = -1 sinon, et, pour z∈ℂnz \in \mathbb C^n et y∈ℝny \in \mathbb R^n, H(z)=∑k=1nσk ∣(R−1z)k∣2,V(y)=H(y−a).\begin{aligned} H(z) &= \sum_{k=1}^n \sigma_k\,\big|(R^{-1}z)_k\big|^2, \\ V(y) &= H(y - a). \end{aligned}

VV prend des valeurs positives près de aa. Soit u=Re1u = Re_1 ; comme (D+N)e1=d1e1(D + N)e_1 = d_1e_1, on a Au=d1uAu = d_1u. La matrice AA étant réelle, Auˉ=dˉ1uˉA\bar u = \bar d_1\bar u, et Re⁡dˉ1=Re⁡d1>s\operatorname{Re}\bar d_1 = \operatorname{Re} d_1 > s.

Soit w=R−1uˉ≠0w = R^{-1}\bar u \neq 0 et jj le plus grand indice tel que wj≠0w_j \neq 0. La ligne jj de (D+N)w=dˉ1w(D + N)w = \bar d_1 w se réduit, D+ND + N étant triangulaire supérieure, à djwj=dˉ1wjd_jw_j = \bar d_1w_j : donc dj=dˉ1d_j = \bar d_1, j⩽qj \leqslant q, et ww n'a de coordonnées non nulles que d'indice au plus qq. Ainsi H(uˉ)=∣w∣2>0H(\bar u) = |w|^2 > 0, et H(u)=∣e1∣2=1H(u) = |e_1|^2 = 1.

Écrivons u=u1+iu2u = u_1 + iu_2 avec u1,u2∈ℝnu_1, u_2 \in \mathbb R^n. Pour deux complexes α,β\alpha, \beta, ∣α+iβ∣2+∣α−iβ∣2=2∣α∣2+2∣β∣2|\alpha + i\beta|^2 + |\alpha - i\beta|^2 = 2|\alpha|^2 + 2|\beta|^2 ; appliqué aux coordonnées de R−1u1R^{-1}u_1 et R−1u2R^{-1}u_2, cela donne H(u)+H(uˉ)=2H(u1)+2H(u2).H(u) + H(\bar u) = 2H(u_1) + 2H(u_2). L'un des vecteurs réels u1,u2u_1, u_2, notons-le h0h_0, vérifie donc H(h0)>0H(h_0) > 0, et V(a+τh0)=τ2H(h0)>0V(a + \tau h_0) = \tau^2H(h_0) > 0 pour tout réel τ≠0\tau \neq 0.

Croissance de VV. Soit ρ\rho et ρ′\rho' comme dans la preuve de Systèmes linéaires stables, (i), de sorte que ∥h∥⩽ρ∣R−1h∣\|h\| \leqslant \rho|R^{-1}h| et ∣R−1h∣2⩽nρ′2∥h∥2|R^{-1}h|^2 \leqslant n\rho'^2\|h\|^2. Écrivons X(a+h)=Ah+∥h∥ε(h)X(a + h) = Ah + \|h\|\varepsilon(h), avec ε(h)→0\varepsilon(h) \to 0, et soit r0>0r_0 > 0 tel que B(a,r0)⊂UB(a, r_0) \subset U et 2nn ρρ′∥ε(h)∥⩽κ22n\sqrt n\,\rho\rho'\|\varepsilon(h)\| \leqslant \frac\kappa2 pour ∥h∥<r0\|h\| < r_0.

Le long d'une courbe intégrale qui reste dans B(a,r0)B(a, r_0), posons h=φt(x)−ah = \varphi_t(x) - a, w=R−1hw = R^{-1}h, w~=(σkwk)k\tilde w = (\sigma_kw_k)_k et g(t)=V(φt(x))g(t) = V(\varphi_t(x)). On a w′=(D+N)w+∥h∥ R−1ε(h),w' = (D + N)w + \|h\|\,R^{-1}\varepsilon(h), et, comme dans la preuve de Systèmes linéaires stables, (i), g′=2Re⁡⟨w~,w′⟩g' = 2\operatorname{Re}\langle \tilde w, w' \rangle.

Le terme diagonal vaut 2∑kσkRe⁡(dk)∣wk∣2=2sg+2∑kσk(Re⁡dk−s)∣wk∣2⩾2sg+2κ∣w∣2,\begin{aligned} &2\sum_k \sigma_k\operatorname{Re}(d_k)|w_k|^2 \\ &\quad = 2sg + 2\sum_k \sigma_k(\operatorname{Re} d_k - s)|w_k|^2 \\ &\quad \geqslant 2sg + 2\kappa|w|^2, \end{aligned} puisque σk(Re⁡dk−s)=∣Re⁡dk−s∣\sigma_k(\operatorname{Re} d_k - s) = |\operatorname{Re} d_k - s|. Le terme en NN est au plus 2η∣w∣2=κ2∣w∣22\eta|w|^2 = \frac\kappa2|w|^2 en valeur absolue, puisque ∣w~∣=∣w∣|\tilde w| = |w|.

Pour le reste, ∣⟨w~,z⟩∣⩽∑k∣wk∣∣zk∣⩽n∣w∣∣z∣|\langle \tilde w, z \rangle| \leqslant \sum_k |w_k||z_k| \leqslant n|w||z|, avec ∣R−1ε(h)∣⩽nρ′∥ε(h)∥|R^{-1}\varepsilon(h)| \leqslant \sqrt n\rho'\|\varepsilon(h)\| et ∥h∥⩽ρ∣w∣\|h\| \leqslant \rho|w| ; il est au plus 2nnρρ′∥ε(h)∥ ∣w∣2⩽κ2∣w∣22n\sqrt n\rho\rho'\|\varepsilon(h)\|\,|w|^2 \leqslant \frac\kappa2|w|^2. Au total, g′⩾2sg+κ∣w∣2⩾2sg.g' \geqslant 2sg + \kappa|w|^2 \geqslant 2sg.

Conclusion. Supposons aa stable, et soit δ\delta associé à ε=r0\varepsilon = r_0. Soit τ>0\tau > 0 tel que ∥τh0∥<δ\|\tau h_0\| < \delta, et x=a+τh0x = a + \tau h_0. La courbe issue de xx est définie sur [0,+∞[[0, +\infty[ et reste dans B(a,r0)B(a, r_0) ; la fonction g e−2stg\,e^{-2st} y croît (Signe de la dérivée), et g(t)⩾V(x)e2stg(t) \geqslant V(x)e^{2st} tend vers +∞+\infty, puisque V(x)>0V(x) > 0. Or g(t)⩽∣w(t)∣2⩽nρ′2r02g(t) \leqslant |w(t)|^2 \leqslant n\rho'^2r_0^2 : c'est absurde, et aa est instable.

□
Remarque

Pour un champ linéaire, le point (iii) du nœud Systèmes linéaires stables donnait le résultat par un vecteur propre. Pour un champ non linéaire, la solution issue d'un vecteur propre n'est plus connue, et c'est la fonction VV, différence des carrés des composantes « instables » et « stables », qui la remplace : le cône {V>0}\{V > 0\} repousse les trajectoires hors de toute petite boule.

Remarque

La réciproque est fausse : un équilibre peut être instable avec μ(A)=0\mu(A) = 0. Pour y′=y2y' = y^2 sur ℝ\mathbb R, A=0A = 0, et la solution issue de x>0x > 0 explose en temps fini (Solutions maximales) ; pour le champ (−y+x3, x+y3)(-y + x^3,\ x + y^3) du nœud Linéarisation, le linéarisé est une rotation, et V=x2+y2V = x^2 + y^2 vérifie V˙=2(x4+y4)⩾V2\dot V = 2(x^4 + y^4) \geqslant V^2, donc (1/V)′⩽−1(1/V)' \leqslant -1 : VV explose avant l'instant 1/V(x)1/V(x).

Remarque

Lorsque aucune valeur propre n'a de partie réelle nulle, l'équilibre est dit hyperbolique, et les deux nœuds Linéarisation et Instabilité décident entièrement de sa stabilité. On montre davantage, ce qui est admis ici : les points dont la trajectoire reste dans un petit voisinage de aa forment, près de aa, une sous-variété tangente à la somme des sous-espaces caractéristiques des valeurs propres de partie réelle négative ; c'est le théorème de la variété stable.

Ouvrir dans l'atlas interactif →