Equilibrio e stabilità di sistemi dinamici Stabilità ...

Stabilità interna di sistemi dinamici LTI

Equilibrio e stabilità di sistemi dinamici

2


Stabilità interna di sistemi dinamici LTI TCCriteri di stabilità per sistemi dinamici LTI TCStabilità interna di sistemi dinamici LTI TDCriteri di stabilità per sistemi dinamici LTI TDEsempi di analisi della stabilità interna

Stabilità interna di sistemi dinamici LTI TC


4

Stabilità interna di sistemi dinamici LTI TC (1/5)

Dato un sistema dinamico, a dimensione finita, MIMO, a tempo continuo, lineare e stazionario (LTI), descritto dall’equazione di stato , se ne considerino due diverse evoluzioni temporali:

Un movimento “nominale” ottenuto applicando un ingresso “nominale” al sistema posto in uno stato iniziale “nominale” ⇒

soddisfa il seguente sistema di equazioni

Un movimento “perturbato” ottenuto applicando lo stesso ingresso “nominale” al sistema posto inuno stato iniziale differente (“perturbato”) ⇒


( ) ( ) ( )x t Ax t Bu t= +

( )x t( )u t

0 0( 0)x t x= =

( )x t

0 0x x≠( )u t

( )x t

( )x t0 0( ) ( ) ( ), ( 0)x t Ax t Bu t x t x= + = =

0 0( ) ( ) ( ), ( 0)x t Ax t Bu t x t x= + = =

5


La differenza fra i due diversi movimenti costituisce la perturbazione sullo stato del sistema:

L’evoluzione temporale della perturbazione sullo stato δx (t ) è soluzione dell’equazione differenziale

con condizione iniziale

( ) ( ) ( ) ( ) ( ) ( )nx t x t x t x t x t x tδ δ= − ∈ ⇒ = +

0 0 0 0 0 0( 0) ( 0) ( 0) 0x t x t x t x x xδ δ= = = − = = − = ≠

( ) ( )

[ ]( )

( ) ( ) ( )( ) ( ) ( )

( ) ( ) ( ) ( )( ) ( ) ( ) ( ) ( )

d x t d x t x tdt dt

x t x t x t

Ax t Bu t Ax t Bu tAx t Ax t A x t x t A x t

δδ

δ

−= = = − =

= + − + == − = − =

6


La soluzione δx (t ) dell’equazione differenziale

è data da

Nel caso di un sistema dinamico LTI TC, l’evoluzione temporale della perturbazione sullo stato non dipende quindi dallo stato iniziale “nominale” o dall’ingresso “nominale” , cioè non dipende dal particolare movimento “nominale” considerato ⇒ nel caso dei sistemi dinamici LTI TC, la proprietàdi stabilità riguarda l’intero sistema e non i singoli movimenti, come avviene invece nel caso dei sistemi dinamici non lineari

0 0 0 0( ) ( ), ( 0)x t A x t x t x x xδ δ δ δ= = = − =

0( ) , 0Atx t e x tδ δ= ∀ ≥

( )x tδ0x

( )u t( )x t

7


Come conseguenza dei risultati dell’analisi modale, l’evoluzione temporale della perturbazione sullo stato

è combinazione lineare dei modi propri del sistema, che in generale sono del tipo

L’evoluzione temporale dei modi propri dipende dagli autovalori λi della matrice di stato A . In particolare, i modi propri di un sistema dinamico LTI TC sono:

Esponenzialmente convergenti se Re(λi ) < 0Limitati se Re(λi ) = 0 e μi’ = 1Polinomialmente divergenti se Re(λi ) = 0 e μi’ > 1Esponenzialmente divergenti se Re(λi ) > 0

0( ) , 0Atx t e x tδ δ= ∀ ≥

( ) ( )( )e1, ( ) cos m , 1ii

ii i

ti i im t t e t

λμμ λ ϕ μ μ

′−′ ′= + ≤ ≤I

8


La perturbazione δx (t ) rimane limitata nel tempo e tende a zero asintoticamente ( ) se e soltanto setutti i modi propri sono convergenti⇒ il sistema dinamico LTI è asintoticamente stabileLa perturbazione δx (t ) rimane limitata nel tempo ma non tende a zero asintoticamente se e soltanto se nessun modo proprio è divergente ed almeno un modo proprio è limitato⇒ il sistema dinamico LTI è semplicemente stabileLa perturbazione δx (t ) non è limitata nel tempo e anzi diverge se e soltanto se almeno un modo proprio è divergente⇒ il sistema dinamico LTI è instabileSi possono così formulare alcuni criteri di stabilità

t→∞

Criteri di stabilità persistemi dinamici LTI TC


10

Criterio di asintotica stabilità per sistema LTI TC

Dato un sistema dinamico LTI a tempo continuo, condizione necessaria e sufficiente affinché risulti asintoticamente stabile è che

In tal caso la perturbazione sullo stato δx (t ), oltre a rimanere limitata, tende a zero asintoticamente per qualsiasi perturbazione iniziale δx (t0 ) ⇒ il sistema è globalmente asintoticamente stabileSe l’ingresso nominale è costante e pari ad , esiste un unico stato di equilibrio pari a (infatti la matrice A è invertibile, poiché )ed è globalmente asintoticamente stabile

( ): e ( ) 0ii Aλ∀ <

1det 0ni iA λ== ≠∏

1x A Bu−= −u

11

Criterio di instabilità per sistema LTI TC

Dato un sistema dinamico LTI a tempo continuo, condizione (soltanto) sufficiente affinché risulti instabile è che

( ): e ( ) 0ii Aλ∃ >

12

Criterio di semplice stabilità per sistema LTI TC

Dato un sistema dinamico LTI a tempo continuo, condizione (soltanto) sufficiente affinché risulti semplicemente stabile è che

( )( )( )

: e ( ) 0

: e ( ) 0

: e ( ) 0 1

i

k

k k

i A

k A

k A

λ

λ

λ μ

∀ ≤

∃ =

∀ = ⇒ =

13

Caso critico per lo studio della stabilità (1/3)

Dato un sistema dinamico LTI a tempo continuo tale che

esso può risultare semplicemente stabile oppureinstabile, a seconda che modi polinomialmente divergenti siano effettivamente presenti nella evoluzione temporale della perturbazione δx (t )

( )( )

: e ( ) 0

: e ( ) 0, 1i

k k

i A

k A

λ

λ μ

∀ ≤

∃ = >

14


Esempio #1: si consideri il sistema dinamico LTI TC

avente matrice di stato , caratterizzato

dall’autovalore con molteplicità . La perturbazione sullo stato δx (t ) è soluzione della equazione differenziale

δx (t ) è costante e quindi è limitata ma non tende a zero asintoticamente⇒ il sistema risulta semplicemente stabile

0 00 0

A ⎡ ⎤= ⎢ ⎥⎣ ⎦0λ = 2μ =

1 1 1

2 2 2

( ) 0 ( ) (0)0( ) ( )0 ( ) 0 ( ) (0)

x t x t xx t A x t

x t x t x

δ δ δδ δ

δ δ δ

= =⎧ ⎧⎪ ⎪⎡ ⎤= = ⇒ ⇒⎨ ⎨⎢ ⎥ = =⎣ ⎦ ⎪ ⎪⎩ ⎩

15


Esempio #2: si consideri il sistema dinamico LTI TC


ancora dall’autovalore con molteplicità . La perturbazione sullo stato δx (t ) è soluzione della equazione differenziale

δx (t ) diverge polinomialmente, se δx1(0) ≠ 0⇒ il sistema risulta instabile

0 01 0

A ⎡ ⎤= ⎢ ⎥⎣ ⎦0λ = 2μ =

1

12

1 1 1

1 12 2 2

( ) 00 0( ) ( ) ( )( )1 0( ) 0 ( ) (0)( ) ( ) ( ) (0) (0)

x tx t A x t x tx tx t x t xx t x t x t t x x

δδ δ δδδ δ δδ δ δ δ δ

⎡ ⎤ ⎡ ⎤⎡ ⎤= = = ⇒⎢ ⎥ ⎢ ⎥⎢ ⎥⎣ ⎦ ⎣ ⎦⎣ ⎦= =⎧ ⎧

⇒⎨ ⎨= = +⎩ ⎩

16

Prospetto riassuntivo per sistemi dinamici LTI TC

Autovalori λi (A ) del sistema

Semplice stabilità

Oltre a quelli che convergono esponenzialmente, è presente almeno uno limitato

Instabilitàoppure

Semplicestabilità

Oltre a quelli che convergono, almeno uno diverge polinomialmente oppure è limitato

InstabilitàAlmeno uno divergeesponenzialmente

Asintotica stabilità

Convergono tutti esponenzialmente

Proprietà di stabilità del sistema

Modi propri del sistema

: e ( ) 0( )ii Aλ∀ <

: e ( ) 0( )ii Aλ∃ >

: e ( ) 0: e ( ) 0: e ( ) 0

1

( )( )( )

i

k

kk

i Ak Ak A

λλλμ

∀ ≤∃ =∀ =

⇒ =

: e ( ) 0: e ( ) 0,

1

( )( )

i

kk

i Ak A

λλμ

∀ ≤∃ =

>

Stabilità interna di sistemi dinamici LTI TD


18

Stabilità interna di sistemi dinamici LTI TD (1/5)

Dato un sistema dinamico, a dimensione finita, MIMO, a tempo discreto, lineare e stazionario (LTI), descritto dall’equazione di stato se ne considerino due diverse evoluzioni temporali:

Un movimento “nominale” ottenuto applicando un ingresso “nominale” al sistema posto in uno stato iniziale “nominale” ⇒


Un movimento “perturbato” ottenuto applicando lo stesso ingresso “nominale” al sistema posto inuno stato iniziale differente (“perturbato”) ⇒


( 1) ( ) ( ),x k Ax k Bu k+ = +

( )x k( )u k

0 0( 0)x k x= =

( )x k

0 0x x≠( )u k

( )x k

( )x k0 0( 1) ( ) ( ), ( 0)x k Ax k Bu k x k x+ = + = =

0 0( 1) ( ) ( ), ( 0)x k Ax k Bu k x k x+ = + = =

19


La differenza fra i due diversi movimenti costituisce la perturbazione sullo stato del sistema:

L’evoluzione temporale della perturbazione sullo stato δx (k ) è soluzione dell’equazione alle differenze

con condizione iniziale

( ) ( ) ( ) ( ) ( ) ( )nx k x k x k x k x k x kδ δ= − ∈ ⇒ = +

0 0 0 0 0 0( 0) ( 0) ( 0) 0x k x k x k x x xδ δ= = = − = = − = ≠

[ ]( )

( 1) ( 1) ( 1)( ) ( ) ( ) ( )( ) ( ) ( ) ( ) ( )

x k x k x kAx k Bu k Ax k Bu kAx k Ax k A x k x k A x k

δ

δ

+ = + − + == + − + == − = − =

20


La soluzione δx (k ) dell’equazione alle differenze

è data da

Nel caso di un sistema dinamico LTI TD, l’evoluzione temporale della perturbazione sullo stato non dipende quindi dallo stato iniziale “nominale” o dall’ingresso “nominale” , cioè non dipende dal particolare movimento “nominale” considerato ⇒ nel caso dei sistemi dinamici LTI TD, la proprietàdi stabilità riguarda l’intero sistema e non i singoli movimenti, come avviene invece nel caso dei sistemi dinamici non lineari

0 0 0 0( 1) ( ), ( 0)x k A x k x k x x xδ δ δ δ+ = = = − =

0( ) , 0kx k A x kδ δ= ∀ ≥

( )x kδ0x

( )u k( )x k

21


Come conseguenza dei risultati dell’analisi modale, l’evoluzione temporale della perturbazione sullo stato

è combinazione lineare dei modi propri del sistema, che in generale sono del tipo

L’evoluzione temporale dei modi propri dipende dagli autovalori λi della matrice di stato A . In particolare, i modi propri di un sistema dinamico LTI TD sono:

Geometricamente convergenti se |λi | < 1Limitati se |λi | = 1 e μi’ = 1Polinomialmente divergenti se |λi | = 1 e μi’ > 1Geometricamente divergenti se |λi | > 1

0( ) , 0kx k A x kδ δ= ∀ ≥

( )( )1, ( ) cos arg , 1i

ii i

ki i i im k k kμμ λ λ ϕ μ μ

′−′ ′= + ≤ ≤

22


La perturbazione δx (k) rimane limitata nel tempo e tende a zero asintoticamente ( ) se e soltanto setutti i modi propri sono convergenti⇒ il sistema dinamico LTI è asintoticamente stabileLa perturbazione δx (k) rimane limitata nel tempo ma non tende a zero asintoticamente se e soltanto se nessun modo proprio è divergente ed almeno un modo proprio è limitato⇒ il sistema dinamico LTI è semplicemente stabileLa perturbazione δx (k) non è limitata nel tempo e anzi diverge se e soltanto se almeno un modo proprio è divergente⇒ il sistema dinamico LTI è instabileSi possono così formulare alcuni criteri di stabilità

k→∞

Criteri di stabilità persistemi dinamici LTI TD


24

Criterio di asintotica stabilità per sistema LTI TD

Dato un sistema dinamico LTI a tempo discreto, condizione necessaria e sufficiente affinché risulti asintoticamente stabile è che

In tal caso la perturbazione sullo stato δx (k ), oltre a rimanere limitata, tende a zero asintoticamente per qualsiasi perturbazione iniziale δx (k0) ⇒ il sistema è globalmente asintoticamente stabileSe l’ingresso nominale è costante e pari ad , esiste un unico stato di equilibrio (infatti la matrice I−A è invertibile: )ed è globalmente asintoticamente stabile

: ( ) 1ii Aλ∀ <

1det( ) (1 ) 0ni iI A λ=− = − ≠∏

1( )x I A Bu−= −u

25

Criterio di instabilità per sistema LTI TD

Dato un sistema dinamico LTI a tempo discreto, condizione (soltanto) sufficiente affinché risulti instabile è che

: ( ) 1ii Aλ∃ >

26

Criterio di semplice stabilità per sistema LTI TD

Dato un sistema dinamico LTI a tempo discreto, condizione (soltanto) sufficiente affinché risulti semplicemente stabile è che

: ( ) 1

: ( ) 1

: ( ) 1 1

i

k

k k

i Ak A

k A

λ

λ

λ μ

∀ ≤

∃ =

∀ = ⇒ =

27


Dato un sistema dinamico LTI a tempo discreto tale che

esso può risultare semplicemente stabile oppureinstabile, a seconda che modi polinomialmente divergenti siano effettivamente presenti nella evoluzione temporale della perturbazione δx (k )

: ( ) 1

: ( ) 1, 1i

k k

i Ak A

λ

λ μ

∀ ≤

∃ = >

28


Esempio #1: si consideri il sistema dinamico LTI TD


dall’autovalore con molteplicità . La perturbazione sullo stato δx (k ) è soluzione della equazione alle differenze

δx (k ) è costante e quindi è limitata ma non tende a zero asintoticamente⇒ il sistema risulta semplicemente stabile

1 00 1

A ⎡ ⎤= ⎢ ⎥⎣ ⎦1λ = 2μ =

1

2

1 1 1 1

2 2 2 2

( )( 1) ( ) ( ) ( )

( 1) ( ) ( ) (0)( 1) ( ) ( ) (0)

x kx k A x k x k x kx k x k x k xx k x k x k x

δδ δ δ δδ δ δ δδ δ δ δ

⎡ ⎤+ = = = ⇒⎢ ⎥

⎣ ⎦+ = =⎧ ⎧

⇒⎨ ⎨+ = =⎩ ⎩

29


Esempio #2: si consideri il sistema dinamico LTI TD


ancora dall’autovalore con molteplicità .La perturbazione sullo stato δx (k ) è soluzione della equazione alle differenze

δx (k ) diverge polinomialmente, se δx1(0) ≠ 0

1 01 1

A ⎡ ⎤= ⎢ ⎥⎣ ⎦

1 1

12 2

1 1 1 1

1 12 2 2 2

( ) ( )1 0( 1) ( ) ( ) ( ) ( )1 1( 1) ( ) ( ) (0)( 1) ( ) ( ) ( ) (0) (0)

x k x kx k A x k x k x k x kx k x k x k xx k x k x k x k k x x

δ δδ δ δ δ δδ δ δ δδ δ δ δ δ δ

⎡ ⎤ ⎡ ⎤⎡ ⎤+ = = = ⇒⎢ ⎥ ⎢ ⎥⎢ ⎥ +⎣ ⎦ ⎣ ⎦ ⎣ ⎦+ = =⎧ ⎧

⇒⎨ ⎨+ = + = +⎩ ⎩

1λ = 2μ =

⇒ il sistema risulta instabile

30

Prospetto riassuntivo per sistemi dinamici LTI TD

Autovalori λi (A ) del sistema

Semplice stabilità

Oltre a quelli che convergono geometricamente, è presente almeno uno limitato

Instabilitàoppure

Semplicestabilità

Oltre a quelli che convergono, almeno uno diverge polinomialmente oppure è limitato

InstabilitàAlmeno uno divergegeometricamente

Asintotica stabilità

Convergono tutti geometricamente

Proprietà di stabilità del sistema

Modi propri del sistema

: | ( )| 1ii Aλ∀ <

: | ( )| 1ii Aλ∃ >

: | ( )| 1: | ( )| 1: | ( )| 1

1

i

k

kk

i Ak Ak A

λλλ

μ

∀ ≤∃ =∀ =

⇒ =

: | ( )| 1: | ( )| 1,

1

i

kk

i Ak A

λλμ

∀ ≤∃ =

>

Esempi di analisi della stabilità interna


32

Esempio #1 di analisi della stabilità

Dato il sistema dinamico LTI TC descritto dal modello

analizzarne la proprietà di stabilità interna, sapendo che gli autovalori λi (A ) della matrice di stato A sono

Per analizzare la stabilità interna dei sistemi dinamici LTI a tempo continuo, occorre considerare la parte reale degli autovalori

Tutti gli autovalori hanno Re(λi (A )) < 0⇒ il sistema è globalmente asintoticamente stabile

( ) ( ) ( )( ) ( ) ( )

x t Ax t Bu ty t C x t Du t

= += +

{ } { }( ) 2, 0.4, 0.2, 0.1i Aλ = − − − −

( ){ } { }e ( ) 2, 0.4, 0.2, 0.1i Aλ = − − − −

33





L’autovalore λ4 = 0.5 ha Re(λ4) = 0.5 > 0⇒ il sistema è instabile

( ) ( ) ( )( ) ( ) ( )


= += +

{ } { }( ) 1, 0.5, 0, 0.5i Aλ = − −

( ){ } { }e ( ) 1, 0.5, 0, 0.5i Aλ = − −

34





La coppia λ3,4 = 0.4± 0.2j ha Re(λ3,4) = 0.4 > 0⇒ il sistema è instabile

( ) ( ) ( )( ) ( ) ( )


= += +

{ } { }( ) 0.5 0.1 , 0.4 0.2i A j jλ = − ± ±

( ){ } { }e ( ) 0.5, 0.5, 0.4, 0.4i Aλ = − −

35





Tutti gli autovalori hanno Re(λi (A )) ≤ 0 ed un solo autovalore λ4 = 0 ha Re(λ4)=0⇒ il sistema è semplicemente stabile

( ) ( ) ( )( ) ( ) ( )


= += +

{ } { }( ) 10, 5, 0.1, 0i Aλ = − − −

( ){ } { }e ( ) 10, 5, 0.1, 0i Aλ = − − −

36





Tutti gli autovalori hanno Re(λi (A ))≤ 0 e tutti quellicon Re(λk (A)) = 0 sono distinti, cioè hanno μk = 1⇒ il sistema è semplicemente stabile

( ) ( ) ( )( ) ( ) ( )


= += +

{ } { }( ) 1, , 0i A jλ = − ±

( ){ } { }e ( ) 1, 0, 0, 0i Aλ = −

37


Dato il sistema dinamico LTI TD descritto dal modello


Per analizzare la stabilità interna dei sistemi dinamici LTI a tempo discreto, occorre considerare il modulo degli autovalori

L’autovalore λ1 = −2 ha |λ1| = 2 > 1⇒ il sistema è instabile

( 1) ( ) ( )( ) ( ) ( )

x k Ax k Bu ky k C x k Du k+ = +

= +

{ } { }( ) 2, 0.4, 0.2, 0.1i Aλ = − − − −

{ } { }( ) 2, 0.4, 0.2, 0.1i Aλ =

38





Tutti gli autovalori hanno |λi (A )| ≤ 1 ed un solo autovalore λ1 =−1 ha |λ1| = 1⇒ il sistema è semplicemente stabile

{ } { }( ) 1, 0.5, 0, 0.5i Aλ = − −

( 1) ( ) ( )( ) ( ) ( )


= +

{ } { }( ) 1, 0.5, 0, 0.5i Aλ =

39





Tutti gli autovalori hanno |λi (A )| < 1⇒ il sistema è globalmente asintoticamente stabile

{ } { }( ) 0.5 0.1 , 0.4 0.2i A j jλ = − ± ±

( 1) ( ) ( )( ) ( ) ( )


= +

{ } { }2 2 2 20.5 0.1 0.4 0.20.26 0.26 0.2 0.2( ) , , ,i Aλ + += ==

40





Gli autovalori λ1 =−10 e λ2 =−5 hanno |λ1| = 10 > 1 e |λ2| = 5 > 1⇒ il sistema è instabile

{ } { }( ) 10, 5, 0.1, 0i Aλ = − − −

( 1) ( ) ( )( ) ( ) ( )


= +

{ } { }( ) 10, 5, 0.1, 0i Aλ =

41





Tutti gli autovalori hanno |λi (A )| ≤ 1 e tutti quelli con |λk (A )|= 1 sono distinti, cioè hanno μk = 1⇒ il sistema è semplicemente stabile

{ } { }( ) 1, , 0i A jλ = − ±

( 1) ( ) ( )( ) ( ) ( )


= +

{ } { }( ) 1,1,1, 0i Aλ =

42

Esempio #11 di analisi della stabilità (1/2)

Dato il sistema dinamico LTI TC con matrice di stato

analizzarne la proprietà di stabilità internaPer analizzare la stabilità interna dei sistemi dinamici LTI a tempo continuo, occorre considerare la parte reale degli autovaloriIn questo caso, la matrice A è diagonale a blocchi ⇒ gli autovalori λi (A ) sono gli autovalori dei blocchi quadrati posti sulla diagonale principale di A

0 1 0 04 4 0 00 0 0.2 00 0 0 0.3

A

⎡ ⎤⎢ ⎥− −= ⎢ ⎥−⎢ ⎥

−⎣ ⎦

43


Tutti gli autovalori hanno Re(λi (A )) < 0⇒ il sistema è globalmente asintoticamente stabile

0 1 0 04 4 0 00 0 0.2 00 0 0 0.3

A

⎡ ⎤⎢ ⎥− −= ⎢ ⎥−⎢ ⎥

−⎣ ⎦

A1

{ } { } { }1( ) ( ) 0.2, 0.3i iA Aλ λ= − −∪

1 12 21. .( ) det( ) 4 4 ( 2)

4 4p c A I A λλ λ λ λλ

−= − = = + + = ++

{ } { }( ) 2, 2, 0.2, 0.3i Aλ⇒ = − − − −

( ){ } { }e ( ) 2, 2, 0.2, 0.3i Aλ = − − − −

44


Dato il sistema dinamico LTI TD con matrice di stato

analizzarne la proprietà di stabilità internaPer analizzare la stabilità interna dei sistemi dinamici LTI a tempo discreto, occorre considerare il modulo degli autovaloriIn questo caso, la matrice A è diagonale a blocchi ⇒ gli autovalori λi (A ) sono gli autovalori dei blocchi quadrati posti sulla diagonale principale di A

0 1 0 04 4 0 00 0 0.2 00 0 0 0.3

A

⎡ ⎤⎢ ⎥− −= ⎢ ⎥−⎢ ⎥

−⎣ ⎦

45


Gli autovalori λ1 =λ2 =−2 hanno |λ1| = |λ2| = 2 > 1⇒ il sistema è instabile

0 1 0 04 4 0 00 0 0.2 00 0 0 0.3

A

⎡ ⎤⎢ ⎥− −= ⎢ ⎥−⎢ ⎥

−⎣ ⎦

A1

1 12 21. .( ) det( ) 4 4 ( 2)

4 4p c A I A λλ λ λ λλ

−= − = = + + = ++

{ } { }( ) 2, 2, 0.2, 0.3i Aλ =

{ } { } { }1( ) ( ) 0.2, 0.3i iA Aλ λ= − −∪

{ } { }( ) 2, 2, 0.2, 0.3i Aλ⇒ = − − − −

Equilibrio e stabilità di sistemi dinamici Stabilità ...

Documents

Transcript of Equilibrio e stabilità di sistemi dinamici Stabilità ...