Per <strong>il</strong> sistema <strong>di</strong>namico così ottenuto valgono le tre ipotesi del Teorema 2.30 e quin<strong>di</strong><strong>il</strong> punto fisso x e è asintoticamente stab<strong>il</strong>e. Se x e è asintoticamente stab<strong>il</strong>e <strong>per</strong> dx/dτ =−f(x), allora sarà instab<strong>il</strong>e <strong>per</strong> <strong>il</strong> sistema <strong>di</strong>namico <strong>di</strong> partenza.Teorema 2.31 □Una funzione che sod<strong>di</strong>sfi ai primi due requisiti nell’enunciato del teorema precedenteè detta <strong>di</strong> Liapunov <strong>per</strong> <strong>il</strong> punto fisso x e , se sod<strong>di</strong>sfa anche all’ulteriore terzo requisito èdetta <strong>di</strong> Liapunov in senso stretto. È imme<strong>di</strong>ato che <strong>il</strong> teorema precedente <strong>per</strong>mette <strong>di</strong><strong>di</strong>mostrare la stab<strong>il</strong>ità <strong>di</strong> un punto fisso che sia un minimo relativo proprio <strong>di</strong> un integraleprimo u, infatti l’essere minimo relativo proprio implica che u sod<strong>di</strong>sfa alla prima ipotesidel teorema, <strong>il</strong> fatto che u è un integrale primo implica che è sod<strong>di</strong>sfatta anche la secondaipotesi.Esempio 2.32. Si <strong>di</strong>scutono alcuni sistemi <strong>di</strong>namici stu<strong>di</strong>ati negli esempi precedenti alla luce del teorema<strong>di</strong> stab<strong>il</strong>ità <strong>di</strong> Liapunov. È fac<strong>il</strong>e <strong>di</strong>mostrare che (1, 1) è punto fisso stab<strong>il</strong>e del sistema <strong>di</strong> Lotka–Volterra(2.10): l’integrale primo ha un minimo realtivo proprio in (0, 0), ma u(0, 0) ≠ 0; ciò non è grave, èsufficiente considerare w(x, y) = u(x, y) − u(0, 0) = α(x − log x) + y − log y − α − 1. Ovviamente w ha unminimo relativo proprio in (1, 1), in tale punto assume valore 0 e inoltre L f w(x, y) = 0 <strong>per</strong>ché anche w èun integrale primo del sistema <strong>di</strong>namico <strong>di</strong> Lotka–volterra; in conclusione w è una funzione <strong>di</strong> Liapunove quin<strong>di</strong> (1, 1) è un punto critico stab<strong>il</strong>e.È imme<strong>di</strong>ato <strong>di</strong>mostrare che (0, 0) è punto fisso stab<strong>il</strong>e del sistema <strong>di</strong>namico (2.16) usando l’integraleprimo come funzione <strong>di</strong> Liapunov. Per quanto riguarda i punti fissi del tipo (¯x, 0), con ¯x > 0, si può usarecome funzione <strong>di</strong> Liapunov <strong>il</strong> paraboloide centrato nel punto fisso w(x, y) := (x − ¯x) 2 + y 2 ; si verificafac<strong>il</strong>mente che le ipotesi del Teorema 2.30 sono sod<strong>di</strong>sfatte, in particolareL f w(x, y) = ∇w(x, y) · f(x, y) = 2(x − ¯x)e x y 2 − 2ye x xy = −2¯xy 2 e x ≤ 0∀x, y ∈ Rdove si è usato che ¯x è positivo; è un ut<strong>il</strong>e esercizio tentare <strong>di</strong> capire <strong>per</strong> quale motivo la derivata <strong>di</strong> Liedel paraboloide rispetto al campo f del sistema <strong>di</strong>namico (2.16) risulta minore o uguale a zero, si ricor<strong>di</strong>che le orbite <strong>di</strong> fase del sistema <strong>di</strong>namico sono archi <strong>di</strong> circonferenza che tendono all’asse x da sotto e dasopra l’asse.Si considera, ora, <strong>il</strong> sistema <strong>di</strong>namico (2.18) stu<strong>di</strong>ato in dettaglio nell’Esempio 2.19; l’origine (0, 0) èun punto fisso del sistema <strong>di</strong>namico, ma l’integrale primo u(x, y) = (x 2 − 1) 2 − 1 ha un punto <strong>di</strong> massimorelativo proprio in (0, 0) e u(0, 0) = 1, quin<strong>di</strong> non è una funzione <strong>di</strong> Liapunov. Si osserva, <strong>per</strong>ò, che seu è integrale primo allora sarà integrale primo anche w(x, y) := −[u(x, y) − 1], cioè L f w(x, y) = 0, dovef è <strong>il</strong> campo delle <strong>di</strong>rezioni del sistema <strong>di</strong>namico in questione. D’altro canto w ha un minimo relativoproprio in (0, 0) e w(0, 0) = 0; in coclusione w è una funzione <strong>di</strong> Liapunov e quin<strong>di</strong> (0, 0) è stab<strong>il</strong>e.Esempio 2.33. Si considera <strong>il</strong> sistema <strong>di</strong>namico planare{ ẋ = −x 3 + xyẏ = −y 3 − x 2 (2.24)Il sistema <strong>di</strong>namico è nella forma (2.6) con campo delle <strong>di</strong>rezioni f(x, y) = (−x 3 + xy, −y 3 − x 2 ). Sideterminano i punti fissi e se ne stu<strong>di</strong>a la stab<strong>il</strong>ità. Punti fissi:{ x(−x 2 + y) = 0f(x, y) = 0 ⇒y 3 + x 2 ⇒ x = 0 e y = 0= 0Quin<strong>di</strong> l’origine è l’unico punto fisso. Per stu<strong>di</strong>arne la stab<strong>il</strong>ità si considera dapprima sistema linearizzatonell’intorno dell’origine; si calcolano le derivate del campo delle <strong>di</strong>rezioni∂f 1∂x = −3x2 + y,∂f 1∂y = x, ∂f 2∂x = −2x2 ,∂f 2∂y = 3y2fismat05.tex – 20 Apr<strong>il</strong>e 2006 – 13:12 pagina 44
e si stu<strong>di</strong>a la matrice associata al sistema linearizzato( ) 0 0A =⇒ det(A − λ1I) = λ 2 = 0 ⇒ λ 1 = λ 2 = 00 0Il Teorema 2.23 non <strong>per</strong>mette <strong>di</strong> trarre conclusioni a proposito della stab<strong>il</strong>ità dell’origine.Si cerca, allora, una funzione <strong>di</strong> Liapunov; tipicamente conviene provare con funzioni <strong>di</strong> secondogrado contenenti solo le potenze <strong>di</strong> or<strong>di</strong>ne 2 (paraboloide), cioè w(x, y) = ax 2 + by 2 , con a, b > 0; conquesta scelta, infatti, si è sicuri che l’origine sia un punto <strong>di</strong> minimo relativo proprio <strong>per</strong> w. AlloraL f w(x, y) = (2ax, 2by) · (−x 3 + xy, −y 3 − x 2 ) = −2ax 4 + 2ax 2 y − 2by 4 − 2bx 2 yIn questo caso se si sceglie a = b > 0 <strong>il</strong> gioco è fatto, infatti:L f w(x, y) = −2ax 4 + 2ax 2 y − 2ay 4 − 2ax 2 y = −2a(x 4 + y 4 )È imme<strong>di</strong>ato osservare che w è una funzione <strong>di</strong> Liapunov in senso stretto <strong>per</strong> <strong>il</strong> punto fisso (0, 0) e quin<strong>di</strong>,in virtù del Teorema 2.30, (0, 0) è un punto fisso asintoticamente stab<strong>il</strong>e.Esempio 2.34. Moto armonico semplice. Si considera l’equazione <strong>di</strong>fferenziale lineare del secondo or<strong>di</strong>neche descrive <strong>il</strong> moto armonico semplice <strong>di</strong> una palla <strong>di</strong> massa m sottoposta all’azione <strong>di</strong> una molla <strong>di</strong>costante elastica κ:m¨z + κz = 0 (2.25)con m, κ ∈ R + .In primo lujogo si osserva che l’equazione (2.25) può essere scritta nella forma (2.6) introducendo ledue nuove variab<strong>il</strong>i q = z e p = mż, rispettivamente la posizione e l’impulso della palla. Infatti{ q = zp = mż⇒ { ˙q = żṗ = m¨z = −κz⇒{ ˙q = p/mṗ = −κq(2.26)Il sistema (2.26) è un sistema <strong>di</strong>namico nella forma (2.6) con x = (q, p) ∈ R 2 e campo delle <strong>di</strong>rezionif(q, p) = (p/m, −κq).Si scrive ora l’energia meccanica totale che, come è stato già <strong>di</strong>mostrato nell’Esempio 2.16, è unintegrale primo. Si osserva che detta u(q) = κq 2 /2, si ha du/dq = κq, allora l’energia meccanica totale èh(q, p) = 12m p2 + 1 2 κq2 (2.27)I punti <strong>di</strong> equ<strong>il</strong>ibrio sod<strong>di</strong>sfano <strong>il</strong> sistema <strong>di</strong> equazioni f(q, p) = 0, ovvero{ {f1 (q, p) = 0 p/m = 0f 2 (q, p) = 0 ⇒ ⇒ P = (0, 0)−κq = 0è l’unico punto <strong>di</strong> equ<strong>il</strong>ibrio.Per stu<strong>di</strong>are la stab<strong>il</strong>ità del punto fisso si linearizza <strong>il</strong> sistema in un intorno dell’origine. In realtà <strong>il</strong>sistema (2.26) è lineare e può essere scritto nella forma ẋ = Ax con A la matrice 2 × 2 reale( )0 1/mA =−κ 0Si scrive, quin<strong>di</strong>, l’equazione secolare e si determinano gli autovalori della matrice A:( )−λ 1/mdet(A − λ1I) = det= λ 2 + κ √ κ−κ −λm = 0 ⇒ λ 1,2 = ±imGli autovalori <strong>di</strong> A sono complessi coniugati con parte reale nulla. Il Teorema 2.23 non <strong>per</strong>mette <strong>di</strong>concludere nulla sulla stab<strong>il</strong>ità dell’origine, ma, ricordando che <strong>il</strong> sistema sotto esame è un sistema lineare,fismat05.tex – 20 Apr<strong>il</strong>e 2006 – 13:12 pagina 45
- Page 1 and 2: Esercizi e appunti per il corso di
- Page 3 and 4: ⇒ x(t) = Ce tcon C ∈ RImpondend
- Page 5 and 6: ✻x✟✻✟ ✟ ✟ ✟x✻x 0x
- Page 7 and 8: viene scelto vicino a x 2 il sistem
- Page 9 and 10: Esercizio 1.4. Supponendo uniforme
- Page 11 and 12: Questo è parzialmente vero nel cas
- Page 13 and 14: L’integrale, infatti, ha senso pe
- Page 15 and 16: sensato, perché si ricorda che le
- Page 17 and 18: dimostrare che il tempo t 1 − t 0
- Page 19 and 20: con f : R n → R una funzione asse
- Page 21 and 22: Come nel caso unidimensionale verif
- Page 23 and 24: - asintoticamente stabile se e solo
- Page 25 and 26: Esempio 2.9. Sulla base di argoment
- Page 27 and 28: x 2✻✲✻ ✲✛ ✻❄✛ ❄
- Page 29 and 30: Per esempio nel punto ¯x = (0, 1)
- Page 31 and 32: livello chiusa, questa osservazione
- Page 33 and 34: (x 2 + y 2 )/2. Si può verificare
- Page 35 and 36: - Per a = 1, la curva di livello pa
- Page 37 and 38: Teorema 2.21 Si consideri il sistem
- Page 39 and 40: ciò fa intuire che in qualche sens
- Page 41 and 42: Teorema 2.28 (Stabilità dei sistem
- Page 43: La palla B δ (x e ) è proprio que
- Page 47 and 48: −ω 2 sin θ, dove ω = √ g/l
- Page 49 and 50: Da questa proprietà segue che w è
- Page 51 and 52: ✲✲✛✛P 3p✻✲✲✛✛P 2
- Page 53 and 54: −πp✲✻✲✲ ✲✛0 π q✛
- Page 55 and 56: è soddisfatta, perché L f ′w(q,
- Page 57 and 58: y✻✛− √ 2✲✲ ✲✲✛
- Page 59 and 60: Esercizio 2.16. Per i seguenti sist
- Page 61 and 62: degli assi cartesiani dello spazio
- Page 63 and 64: Dal momento che I 1 < I 2 < I 3 si
- Page 65 and 66: Il problema (3.1) può essere ricon
- Page 67 and 68: dell’analisi avendo come riferime
- Page 69 and 70: Γ e = {(q, p) ∈ R 2 : q e 1 ≤
- Page 71 and 72: Osservato che dall’equazione dell
- Page 73 and 74: assunto nel minimo, cioè a zero. I
- Page 75 and 76: La tesi, allora, segue in virtù de
- Page 77 and 78: u(q)✻q 3 q 4 q 1 q 2✲ qu 0p
- Page 79 and 80: p✻q0(q − q 0 ) 1/2 p (q − q 0
- Page 81 and 82: Si può osservare che ˙θ ha segno
- Page 83 and 84: Troncando lo sviluppo delle potenze
- Page 85 and 86: deve specificare il valore della ca
- Page 87 and 88: valore della soluzione sull’asse
- Page 89 and 90: x ∈ [a, b], u x (a, t) = u 1,a (t
- Page 91 and 92: Inoltre si è usata l’identità n
- Page 93 and 94: con q = (x, y, z) ∈ R 3 , t ∈ R
- Page 95 and 96:
Le equazioni del moto possono esser
- Page 97 and 98:
Sostituendo queste espressioni nell
- Page 99 and 100:
1. u x + u y = u;2. 2u x − 3u y =
- Page 101 and 102:
In realtà lo studio è limitato al
- Page 103 and 104:
3. u(x, y) = log √ x 2 + y 2 , (x
- Page 105 and 106:
nella regione√2E := {(ξ, η) ∈
- Page 107 and 108:
ammette l’unica soluzioneu(x, t)
- Page 109 and 110:
Esercizio 6.43. Una corda semi-illi
- Page 111 and 112:
Esercizio 6.50. Come l’Esercizio
- Page 113 and 114:
Esercizio 6.59. Per effetto di una
- Page 115 and 116:
Esercizio 6.66. Si risolva l’equa
- Page 117 and 118:
4. u 0 (x) = exp{−α|x|} con α
- Page 119:
6.3. Equazione di Laplace: funzioni