Teoria: Stabilità degli equilibri dei sistemi autonomi 2x2In un sistema autonomo y' = f(y) le soluzioni sono invarianti per traslazioni nel tempo (D) e le orbite nel piano delle fasi non si intersecano (D). Un equilibrio (f(P) = 0) è stabile se le soluzioni vicine restano vicine, asintoticamente stabile se inoltre convergono a P, instabile altrimenti. Per y' = Ay (2×2, det A ≠ 0) l'origine è asintoticamente stabile se e solo se entrambi gli autovalori hanno parte reale negativa, cioè tr A < 0 e det A > 0; i tipi sono nodo, sella, fuoco, centro.Stabilità degli equilibri dei sistemi autonomi 2x2 → , Sistemi lineari ed equazioni del secondo ordine a coefficienti costantiUn sistema lineare y' = A(t)y + g(t) con coefficienti continui ha, per ogni dato iniziale, una e una sola soluzione, definita su tutto l'intervallo (D). Le soluzioni del sistema omogeneo formano uno spazio vettoriale di dimensione n (D); quelle del sistema completo sono una soluzione particolare più l'integrale omogeneo (D). Per y'' + by' + cy = 0 si cercano e^(λt): tre casi (D). Un sistema 2×2 a coefficienti costanti si riduce a u'' − (tr A)u' + (det A)u = g, risolta con il metodo di somiglianza.Sistemi lineari ed equazioni del secondo ordine a coefficienti costanti → .
Testo
Per α ∈ ( − 5 , 3 ) \alpha\in(-5,3) α ∈ ( − 5 , 3 ) si consideri
{ x ′ = x + y y ′ = − 4 x − α y . \begin{cases}x'=x+y\\ y'=-4x-\alpha y.\end{cases} { x ′ = x + y y ′ = − 4 x − α y .
(a) Studiare la stabilità dell'origine al variare di α \alpha α .
(b) Per α = 1 \alpha=1 α = 1 determinare l'integrale generale.
(c) Per α = − 1 \alpha=-1 α = − 1 risolvere il problema di Cauchy con x ( 0 ) = 1 x(0)=1 x ( 0 ) = 1 , y ( 0 ) = 0 y(0)=0 y ( 0 ) = 0 .
Soluzione
La matrice è A = ( 1 1 − 4 − α ) A=\begin{pmatrix}1&1\\-4&-\alpha\end{pmatrix} A = ( 1 − 4 1 − α ) con τ = 1 − α \tau=1-\alpha τ = 1 − α e d = − α + 4 d=-\alpha+4 d = − α + 4 . Il discriminante è
Δ = τ 2 − 4 d = ( 1 − α ) 2 − 16 + 4 α = α 2 + 2 α − 15 = ( α + 5 ) ( α − 3 ) , \Delta=\tau^2-4d=(1-\alpha)^2-16+4\alpha=\alpha^2+2\alpha-15=(\alpha+5)(\alpha-3), Δ = τ 2 − 4 d = ( 1 − α ) 2 − 16 + 4 α = α 2 + 2 α − 15 = ( α + 5 ) ( α − 3 ) ,
negativo proprio per α ∈ ( − 5 , 3 ) \alpha\in(-5,3) α ∈ ( − 5 , 3 ) : in tutto l'intervallo gli autovalori sono complessi coniugati, con parte reale τ 2 = 1 − α 2 \frac\tau2=\frac{1-\alpha}2 2 τ = 2 1 − α . Inoltre d = 4 − α > 1 > 0 d=4-\alpha>1>0 d = 4 − α > 1 > 0 .
(a) Stabilità
α > 1 \alpha>1 α > 1 (cioè α ∈ ( 1 , 3 ) \alpha\in(1,3) α ∈ ( 1 , 3 ) ): τ < 0 \tau<0 τ < 0 , d > 0 d>0 d > 0 : asintoticamente stabile (fuoco stabile).
α = 1 \alpha=1 α = 1 : τ = 0 \tau=0 τ = 0 : autovalori ± i 3 \pm i\sqrt3 ± i 3 , centro : stabile, non asintoticamente.
α < 1 \alpha<1 α < 1 (cioè α ∈ ( − 5 , 1 ) \alpha\in(-5,1) α ∈ ( − 5 , 1 ) ): τ > 0 \tau>0 τ > 0 : instabile (fuoco instabile).
(b) α = 1 \alpha=1 α = 1
A = ( 1 1 − 4 − 1 ) A=\begin{pmatrix}1&1\\-4&-1\end{pmatrix} A = ( 1 − 4 1 − 1 ) , τ = 0 \tau=0 τ = 0 , d = 3 d=3 d = 3 . Poiché a 12 = 1 ≠ 0 a_{12}=1\ne0 a 12 = 1 = 0 e g = 0 g=0 g = 0 : x ′ ′ − 0 ⋅ x ′ + 3 x = 0 x''-0\cdot x'+3x=0 x ′′ − 0 ⋅ x ′ + 3 x = 0 , λ = ± i 3 \lambda=\pm i\sqrt3 λ = ± i 3 ,
x ( t ) = β cos 3 t + γ sin 3 t . x(t)=\beta\cos\sqrt3t+\gamma\sin\sqrt3t . x ( t ) = β cos 3 t + γ sin 3 t .
Dalla prima equazione y = x ′ − x y=x'-x y = x ′ − x e x ′ = − 3 β sin 3 t + 3 γ cos 3 t x'=-\sqrt3\beta\sin\sqrt3t+\sqrt3\gamma\cos\sqrt3t x ′ = − 3 β sin 3 t + 3 γ cos 3 t :
y ( t ) = ( 3 γ − β ) cos 3 t − ( 3 β + γ ) sin 3 t . y(t)=(\sqrt3\gamma-\beta)\cos\sqrt3t-(\sqrt3\beta+\gamma)\sin\sqrt3t . y ( t ) = ( 3 γ − β ) cos 3 t − ( 3 β + γ ) sin 3 t .
Tutte le soluzioni sono periodiche, di periodo 2 π 3 \frac{2\pi}{\sqrt3} 3 2 π (orbite chiuse).
(c) α = − 1 \alpha=-1 α = − 1
A = ( 1 1 − 4 1 ) A=\begin{pmatrix}1&1\\-4&1\end{pmatrix} A = ( 1 − 4 1 1 ) , τ = 2 \tau=2 τ = 2 , d = 5 d=5 d = 5 : λ 2 − 2 λ + 5 = 0 \lambda^2-2\lambda+5=0 λ 2 − 2 λ + 5 = 0 , λ = 1 ± 2 i \lambda=1\pm2i λ = 1 ± 2 i . Quindi x = e t ( c 1 cos 2 t + c 2 sin 2 t ) x=e^t(c_1\cos2t+c_2\sin2t) x = e t ( c 1 cos 2 t + c 2 sin 2 t ) e y = x ′ − x y=x'-x y = x ′ − x .
Condizioni iniziali. x ( 0 ) = c 1 = 1 x(0)=c_1=1 x ( 0 ) = c 1 = 1 . Poi x ′ ( 0 ) = x ( 0 ) + y ( 0 ) = 1 x'(0)=x(0)+y(0)=1 x ′ ( 0 ) = x ( 0 ) + y ( 0 ) = 1 , mentre x ′ = e t [ ( c 1 + 2 c 2 ) cos 2 t + ( c 2 − 2 c 1 ) sin 2 t ] x'=e^t[(c_1+2c_2)\cos2t+(c_2-2c_1)\sin2t] x ′ = e t [( c 1 + 2 c 2 ) cos 2 t + ( c 2 − 2 c 1 ) sin 2 t ] dà x ′ ( 0 ) = c 1 + 2 c 2 = 1 + 2 c 2 x'(0)=c_1+2c_2=1+2c_2 x ′ ( 0 ) = c 1 + 2 c 2 = 1 + 2 c 2 . Quindi c 2 = 0 c_2=0 c 2 = 0 . Allora
x ( t ) = e t cos 2 t , y ( t ) = x ′ − x = e t ( cos 2 t − 2 sin 2 t ) − e t cos 2 t = − 2 e t sin 2 t . x(t)=e^t\cos2t,\qquad y(t)=x'-x=e^t(\cos2t-2\sin2t)-e^t\cos2t=-2e^t\sin2t . x ( t ) = e t cos 2 t , y ( t ) = x ′ − x = e t ( cos 2 t − 2 sin 2 t ) − e t cos 2 t = − 2 e t sin 2 t .
Controllo. y ′ = − 2 e t sin 2 t − 4 e t cos 2 t = − 4 x − α y = − 4 e t cos 2 t − 2 e t sin 2 t y'=-2e^t\sin2t-4e^t\cos2t=-4x-\alpha y=-4e^t\cos2t-2e^t\sin2t y ′ = − 2 e t sin 2 t − 4 e t cos 2 t = − 4 x − α y = − 4 e t cos 2 t − 2 e t sin 2 t ✓ (con α = − 1 \alpha=-1 α = − 1 ). La soluzione è una spirale che si allontana (fuoco instabile, α = − 1 < 1 \alpha=-1<1 α = − 1 < 1 ).
Versione ripasso
Testo. x ′ = x + y x'=x+y x ′ = x + y , y ′ = − 4 x − α y y'=-4x-\alpha y y ′ = − 4 x − α y , α ∈ ( − 5 , 3 ) \alpha\in(-5,3) α ∈ ( − 5 , 3 ) : (a) stabilità; (b) α = 1 \alpha=1 α = 1 integrale generale; (c) α = − 1 \alpha=-1 α = − 1 Cauchy x ( 0 ) = 1 x(0)=1 x ( 0 ) = 1 , y ( 0 ) = 0 y(0)=0 y ( 0 ) = 0 .
Passaggi (Stabilità degli equilibri dei sistemi autonomi 2x2In un sistema autonomo y' = f(y) le soluzioni sono invarianti per traslazioni nel tempo (D) e le orbite nel piano delle fasi non si intersecano (D). Un equilibrio (f(P) = 0) è stabile se le soluzioni vicine restano vicine, asintoticamente stabile se inoltre convergono a P, instabile altrimenti. Per y' = Ay (2×2, det A ≠ 0) l'origine è asintoticamente stabile se e solo se entrambi gli autovalori hanno parte reale negativa, cioè tr A < 0 e det A > 0; i tipi sono nodo, sella, fuoco, centro.Stabilità degli equilibri dei sistemi autonomi 2x2 → ).
τ = 1 − α \tau=1-\alpha τ = 1 − α , d = 4 − α > 0 d=4-\alpha>0 d = 4 − α > 0 , Δ = ( α + 5 ) ( α − 3 ) < 0 \Delta=(\alpha+5)(\alpha-3)<0 Δ = ( α + 5 ) ( α − 3 ) < 0 : sempre complessi. (a) α ∈ ( 1 , 3 ) \alpha\in(1,3) α ∈ ( 1 , 3 ) as. stabile; α = 1 \alpha=1 α = 1 centro; α ∈ ( − 5 , 1 ) \alpha\in(-5,1) α ∈ ( − 5 , 1 ) instabile.
(b) x ′ ′ + 3 x = 0 x''+3x=0 x ′′ + 3 x = 0 : x = β cos 3 t + γ sin 3 t x=\beta\cos\sqrt3t+\gamma\sin\sqrt3t x = β cos 3 t + γ sin 3 t , y = x ′ − x y=x'-x y = x ′ − x .
(c) λ = 1 ± 2 i \lambda=1\pm2i λ = 1 ± 2 i ; x = e t ( c 1 cos 2 t + c 2 sin 2 t ) x=e^t(c_1\cos2t+c_2\sin2t) x = e t ( c 1 cos 2 t + c 2 sin 2 t ) , c 1 = 1 c_1=1 c 1 = 1 , x ′ ( 0 ) = 1 ⇒ c 2 = 0 x'(0)=1\Rightarrow c_2=0 x ′ ( 0 ) = 1 ⇒ c 2 = 0 : x = e t cos 2 t x=e^t\cos2t x = e t cos 2 t , y = − 2 e t sin 2 t y=-2e^t\sin2t y = − 2 e t sin 2 t .
Errori tipici. Calcolare x ′ ( 0 ) x'(0) x ′ ( 0 ) senza passare dal sistema; segno di τ \tau τ per α \alpha α (cresce con − α -\alpha − α ).