Sistemi dinamici a tempo continuo: variabili di stato Flashcards
Sistemi dinamici: caratteristiche
Legame dinamico tra u e y: y dipende non solo da u ma anche dallo stato x
Per determinare l’uscita, occorre conoscere, oltre all’ingresso, la condizione iniziale dello stato
Classificazione sistemi
1) Sistema statico o dinamico, in base alla dipendenza di y da x
2) Lineare: f e g sono funzioni lineari di x e u
Non lineare: f e g sono funzioni non lineari di x e u
3) Tempo invariante: né f né g sono funzioni esplicite del tempo
Tempo variante: f o g sono funzioni esplicite del tempo
4) Ordine del sistema, in base al numero di variabili di stato
5) SISO: u e y sono scalari
MIMO: u e/o y sono vettoriali
6) Strettamente proprio: g non dipende direttamente da u
Proprio: g è funzione di u
Rappresentazione di stato
Equazione di stato: x’(t)=f(x(t),u(t),t)
Trasformazione di uscita: y(t)=g(x(t),u(t),t)
Rappresentazione sistemi LTI
x'(t) = Ax(t) + Bu(t) y(t) = Cx(t) + Du(t)
x app R^n
y app R^p
u app R^m
A, B, C, D sono matrici a coefficienti costanti
A nxn, B nxm, C pxn, d pxm
Scelta delle variabili di stato: criterio generale
x1=y
x2=y’
…
xn=y(n-1)
- > Equazioni di stato: x1’=x2, x2’=x3,…, xn’=y(n)=phi
Dove phi è l’equazione differenziale di ordine n che descrive il sistema
Variabili di stato associate in questo modo a fenomeni di accumulo
Movimento di un sistema dinamico
Noti u(t), Vt > = 0 e x0=x(t0) condizione inziale, Si chiama movimento dello stato la corrispondente soluzione dell'equazione di stato Si chiama movimento dell'uscita la soluzione della trasformazione di uscita ottenuta sostituendo x e u
In sistemi TI si può sempre porre t0=0, x0=x(0)
Movimento di equilibrio (sistemi TI)
Fissati x0=x_ costante e u(t)=u_ costante Vt > = 0
Se esiste un corrispondente movimento dello stato x(t)=x_ cost. e un movimento dell’uscita y(t)=y_ cost., x_ e y_ sono lo stato e l’uscita di equilibrio associati a u(t)=u_ e x0=x_
Calcolo dell’equilibrio (sistemi TI)
Se esiste x_ di equilibrio, deve essere soluzione di:
0=f(x_,u_), perché x’=0
L’uscita è y_=g(x_,u_)
Equilibrio nei sistemi LTI
0=Ax_+Bu_ - > Ax_=-Bu_
- > x_ = -A^-1 B u_
y_=Cx_+Du_=(-CA^-1B+D)u_
Esistenza di x_ dipende dalla matrice A:
- Se A è invertibile - > x_ esiste ed è unico
- Se A non è invertibile - > Ax_=-Bu_ può avere nessuna soluzione (zero stati di equilibrio) o infinite soluzioni (infiniti stati di equilibrio)
Studio del movimento per sistemi LTI
Dato il sistema LTI:
x’(t)=Ax(t)+Bu(t)
y(t)=Cx(t)+Du(t)
Dati u(t) Vt > = 0 e x(0) = x0 Calcolare x(t) e y(t) Vt > = 0
Si usano formule di Lagrange:
Movimento dello stato:
x(t) = e^Atx0 + int(0,t) (e^A(t-tau)Bu(tau)dtau)
Movimento dell’uscita:
y(t) = Ce^Atx0 + Cint(0,t) (e^A(t-tau)Bu(tau)dtau) + Du(t)
Entrambe le equazioni sono formate da due parti:
- Movimento libero: dipende solo dalla condizione iniziale
- Movimento forzato: dipende solo dall’ingresso
Principio di sovrapposizione degli effetti per sistemi dinamici LTI
Se dati x(0)=x0’ e u(t)=u’(t) si ottiene x’(t) e y’(t)
e dati x(0)=x0’’ e u(t)=u’‘(t) si ottiene x’‘(t) e y’‘(t)
Allora dati x(0)=ax0’+bx0’’ e u(t)=au’(t)+bu’‘(t) si ottiene
x(t)=ax’(t)+bx’‘(t) e y(t)=ay’(t)+by’‘(t)
In particolare, il movimento libero si ottiene da una condizione iniziale generica con ingresso nullo, il movimento forzato da un ingresso generico con condizione inziale nulla
Linearizzazione
Dato un sistema non lineare tempo invariante, il sistema linearizzato approssima linearmente il sistema di partenza in un intorno di una condizione di equilibrio x_,u_,y_
Il sistema non lineare x'(t)=f(x(t),u(t)) y(t)=g(x(t),u(t)) diventa: δx'(t)=Alin δx(t) + Blin δu(t) δy(t) = Clin δx(t) + Dlin δu(t)
Dove
δx(t) = x(t)-x_
δu(t) = u(t)-u_
δy(t) = y(t)-y_
E se il sistema è SISO: Alin = dfi/dxj |x_,u_ Blin = dfi/du |x_,u_ Clin = dg/dxi |x_,u_ Dlin = dg/du |x_,u_
Esponenziale di matrice
Definizione non operativa:
e^At = Σ(k=0, +inf) (At)^k/k!
Se A è diagonale:
e^At = diag(e^λit) , dove λi sono gli autovalori
Se A è diagonalizzabile:
e^At = T^-1 e^Adt T
dove T^-1 è formata dagli autovalori di A, e Ad è la diagonalizzata di A
- > Modi del tipo e^λt
Se A non è diagonalizzabile:
- > Modi del tipo e^λt e t^k*e^λt , k > = 1 legato alla differenza tra molteplicità algebrica e geometrica
Se A ha autovalori complessi e coniugati:
- > Modi del tipo e^(σ +- jω)t - > andamento oscillatorio
Oscillazioni crescenti all’infinito per σ > 0;
Oscillazioni smorzate per σ < 0
Oscillazioni di ampiezza costante per σ = 0
Se A ha autovalori complessi coniugati non regolari:
- > compaiono inoltre modi del tipo t^k*e^(σ +- jω)t
Rappresentazioni equivalenti sistemi LTI
Se A è diagonalizzabile, esiste T matrice di trasformazione tale che T A T^-1 = Ad diagonale
Mediante questa matrice si può effettuare un cambio di coordinate equivalente dal punto di vista input/output, cioè trasformare il sistema originale in un sistema che con lo stesso ingresso produce la stessa uscita.
Si ottiene ponendo:
A* = TAT^-1 = Ad B* = TB C* = CT^-1 D* = D
Il movimento dello stato invece:
x(0) - > x(t)
x*(0) - > Tx(t)
Stabilità: definizione generale
La stabilità è la proprietà per cui un movimento del sistema, a fronte di una perturbazione dello stato (delle condizioni iniziali), tende a ritornare nella condizione di partenza