|
*[[Fondamenti di automatica/Sistemi lineari tempoinvarianti|Sistemi lineari tempoinvarianti]]
==*[[Fondamenti di automatica/Proprietà e prestazioni|Proprietà deie sistemi==prestazioni]]
===Stabilità===
Esistono criteri di stabilità esterna che si riferiscono alle uscite del sistema e criteri di stabilità interna che si riferiscono allo stato del sistema;
\vedilibro{rif:c}{604, sezione 14.5: Stability of nonlinear systems}
\vedilibro{rif:b}{52, sezione 2.7: Stabilità}
la stabilità interna implica la stabilità esterna, ma in generale la stabilità esterna non implica la stabilità interna
(anche se sotto alcune ipotesi questo è vero
\vedilibro{rif:b}{97} )
Si dice \emph{perturbazione} una variazione di movimento corrispondente a variazioni di ingresso o di condizioni iniziali
un movimento $x(t_{0} \cdots t)$ si dice \emph{stabile}
\vedilibro{rif:b}{52}
se per ogni $\epsilon$ positivo esiste $\delta$ positivo tale che
per tutti gli stati iniziali $x(t_{1})$ che soddisfano la relazione
$||x(t_{0}) - x(t_{1})|| \leq \delta$
risulta
<math>
||x(t_{0} \cdots t) - x(t_{1} \cdots t)|| \leq \epsilon
</math>
per ogni $t \geq 0$;
Se $\delta$ non dipende da $x(t_{0})$ il sistema è \emph{uniformemente stabile},
se $\delta$ può essere grande a piacere il movimento è \emph{globalmente stabile}
\vedilibro{rif:b}{56, sezione 2.7.2: Regione di attrazione}
Per un sistema lineare, se questo ha un movimento asintoticamente stabile rispetto a perturbazioni, allora anche un altro qualsiasi movimento è globalmente ed uniformemente stabile
\vedilibro{rif:b}{73, sezione 3.3.2: Stabilità e movimento libero}
un movimento si dice \emph{instabile}
\vedilibro{rif:b}{53}
se non è stabile
un movimento si dice \emph{asintoticamente stabile}
\vedilibro{rif:b}{54}
se è stabile ed inoltre
<math>
\lim_{t \rightarrow \infty} ||x(t_{0} \cdots t) - x(t_{1} \cdots t)|| = 0
</math>
per ogni $t \geq 0$
Un sistema lineare stazionario
\begin{itemize}
\item
è \emph{stabile} se e solo se tutti i movimenti liberi dello stato sono limitati,
ovvero se tutti i suoi autovalori sono negativi e quelli nulli hanno moltelicità 1
\vedilibro{rif:b}{74, sezione 3.3.3: Stabilità e autovalori}
\item
è \emph{asintoticamente stabile} se e solo se tutti i movimenti liberi dello stato tendono a $0$ per $t \rightarrow \infty$,
ovvero se e solo se tutti i suoi autovalori hanno parte reale strettamente negativa
\item
è \emph{instabile} se e solo se almeno un movimento libero dello stato non è limitato,
ovvero se almeno un autovalore ha parte reale positiva
(oppure se ci sono autovalori doppi nulli)
\end{itemize}
Un sistema asintoticamente stabile ha proprietà utili
\vedilibro{rif:b}{80, sezione 3.3.5: proprietà dei sistemi asintoticamente stabili} :
il movimento dello stato per tempi sufficentemente grandi coincide con il movimento forzato,
in quanto il movimento libero tende ad essere nullo;
la risposta ad un qualunque ingresso di durata finita tende ad annullarsi;
l'uscita è limitata ed il sistema è stabile BIBO
Definiamo la \emph{matrice di transizione dello stato} $\Phi$ tale che
<math>
x(t) = \Phi(t,t_{0}) x(t_{0})
</math>
(che per un sistema stazionario è $e^{At} = \invLaplaceTrasf{(sI-A)^{-1}}$)
La matrice di transizione dello stato soddisfa le proprietà seguenti:
\begin{eqnarray}
\Phi'(t, t_{0}) = A \Phi(t, t_{0}) \\
\Phi(t_{0}, t_{0}) = I \\
\Phi(t_{2}, t_{0}) = \Phi(t_{2},t_{1}) \Phi(t_{1},t_{0}) \\
\Phi^{-1}(t, t_{0}) = \Phi (t_{0}, t)
\end{eqnarray}
In generale un sistema (strettamente proprio) è \emph{stabile BIBS} (bounded input - bounded state) se
<math>
||x(t)|| = \left|| \int_{t_{0}}^{t} \Phi{t, \tau}B(\tau)u(\tau)d\tau \right|| < \epsilon
</math>
mentre è \emph{stabile BIBO} (bounded input - bounded output) se
<math>
||y(t)|| = \left|| C(t)\int_{t_{0}}^{t} \Phi{t, \tau}B(\tau)u(\tau)d\tau \right|| < \epsilon
</math>
La definizione di stabilità secondo Lyapunov è applicabile a sistemi in generale anche non lineari
\vedilibro{rif:c}{646, sezione 14.2: Lyapunov's stability criterion}
====Teorema di Lyapunov====
dato $x' = f(x)$ (con $x$ vettore di $n_{x}$ componenti e $f(\cdot)$ campo vettoriale)
in cui $x = 0 \cdots 0$ è punto di equilibrio (a meno di traslazioni dell'origine),
si associa $V(x) = x^{T}Px$
(funzione di Lyapunov, quadratica e tale che P è una matrice simmetrica definita positiva, con tutti gli autovalori positivi),
allora $V'(x) = 2x^{T}Px' = x^{T}(A^{T}P + PA)x$,
che è negativa se e solo se $(A^{T}P + PA)_{i} < 0$,
nel qual caso il sistema è asintoticamente stabile
Gli autovalori di $A$ hanno parte reale negativa se e solo se esiste $P > 0$ simmetrica tale che $A^{T}P + PA + Q = 0$ con $Q>0$ (in genere $Q=I$),
inoltre se $P$ esiste è unica e il sistema è asintoticamente stabile
====Stabilità marginale e instabilità====
Un sistema è \emph{marginalmente stabile}
se la sua risposta ad un ingresso finito è sempre limitata ma finita;
è \emph{instabile}
se la sua risposta ad un ingresso limitato è illimitata;
è \emph{stabilizzabile}
se tutte le sue vabiabili instabili sono controllabili
===Osservabilità===
Uno stato di un sistema si dice non osservabile se in un tempo finito il movimento libero da esso generato è nullo,
\vedilibro{rif:b}{90, sezione 3.5.3: Osservabilità}
\vedilibro{rif:c}{500, State variable Analysis and Design: Observability}
Un sistema è completamente osservabile se è privo di stati non osservabili,
se e solo se è massimo il rango della \emph{matrice di osservabilità} definita come
<math>
M_{o} =
[C^{T} \phantom{1} A^{T}C^{T} \phantom{1} A^{T2}C^{T} \cdots A^{T n-1}C^{T}]
\in \Re^{n \times pn}
</math>
\`E condizione necessaria per l'osservabilità che la matrice $C$ non abbia elementi nulli: $C_{1..n_{x}} \not= 0$
Se il sistema ha una sola uscita allora la matrice di osservabilità è quadrata e quindi è sufficente che $det(M_{c}) \not= 0$
La proprietà di osservabilità di un sistema è associata alla coppia $(A,C)$, inoltre questa coincide con la controllabilità di $(A^{T}, C^{T})$
Se il sistema è in coordinate modali, è necessario e sufficente che ogni variabile di stato non sia mai nulla (???) o che la matrice $C$ trasformata $CT$ non abbia elementi nulli perché il sistema sia osservabile,
oppure se il sistema è in forma di Jordan, la matrice $C$ trasformata non deve essere nulla nella prima riga corrispondente al blocco di Jordan (????)
\vedilibro{rif:b}{92}
La proprietà di osservabilità coincide per i sistemi lineari con la prorietà di \emph{non ricostruibilità} consistente nell'impossibilità di distinguere lo stato finale (anziché iniziale) di un sistema da quello nullo mediante l'analisi di un transitorio libero dell'uscita di qualunque durata
Dato un sistema non completamente osservabile è sempre possibile isolare la sua parte non osservabile (????)
===Controllabilità===
Consiste nel poter cambiare lo stato del sistema da un valore arbitrario al valore nullo in un tempo finito mediante un opportunuo ingresso
\vedilibro{rif:b}{87, sezione 3.5.2: Raggiungibilità}
\vedilibro{rif:c}{494, State variable Analysis and Design: Controllability}
Un sistema è completamente controllabile se è privo di stati non controllabili,
se e solo se è massimo il rango della \emph{matrice di controllabilità} definita come
<math>
M_{c} =
[B \phantom{1} AB \phantom{1} A^{2}B \cdots A^{n-1}B]
\in \Re^{n \times pn}
</math>
Se il sistema ha un solo ingresso allora la matrice di controllabilità è quadrata e quindi è sufficente che $det(M_{c}) \not= 0$
Dato un sistema non completamente controllabile è sempre possibile isolare la sua parte non controllabile (????)
===Causalità===
Un sistema è causale se:
\begin{itemize}
\item
la sua funzione di trasferimento ha tanti zeri quanti poli oppure un numero minore di zeri rispetto ai poli.
\item
la sua risposta impulsiva è nulla per tempi negativi
\end{itemize}
Un sistema causale è fisicamente realizzabile;
sistemi non causali sono utili solo al livello accademico
===Risposta al gradino===
Si studia la risposta del sistema ad un ingresso gradino unitario
\vedilibro{rif:k}{385, sezione 7-4: Unit step response and time-domain specifications}
per valutare le prorietà della risposta transitoria del sistema, in quanto
la risposta ad un gradino unitario può essere utilizzata per approssimare la risposta ad un ingresso qualunque che sia esprimibile come somma di gradini
che siano sufficentemente piccoli rispetto alla rapidità di variazione della risposta
(i gradini devono essere separati l'uno dall'altro da un intervallo di tempo pari almeno al tempo di assestamento del sistema)
La risposta transitoria è per definizione la risposta del sistema che diventa trascurabile per tempi molto grandi
e corrisponde alla soluzione dell'equazione differenziale omogenea
(movimento libero del sistema)
Le caratteristiche della risposta al gradino che vengono valutate generalmente sono:
\begin{itemize}
\item
sovraelongazione massima $M_{p}$ (oppure $M_{p}\%$ se in percentuale)
\item
tempo di ritardo $t_{d}$
\item
tempo di salita $t_{r}$
\item
tempo di assestamento $t_{s}$
\end{itemize}
====Sovraelongazione massima====
Si definisce \emph{sovraelongazione massima} (o maximum overshoot) $M_{p}$ la differenza tra il valore di regime della risposta e il valore massimo che essa raggiunge in fase transitoria;
la \emph{sovraelongazione massima percentuale} $M_{p}\%$ è misurata in percentuale rispetto al valore di regime.
Valori di sovraelongazione bassi sono desiderabili per un sistema ($ < 5\% $)
====Tempo di assestamento====
Si chiama \emph{tempo di assestamento} (o settling time) $t_{s}$ il tempo necessario alla risposta a portarsi definitivamente a valori vicini al valore di regime;
In genere si misura il tempo dopo cui la risposta non differisce di più del 10\% dal valore di regime
Un valore del tempo di assestamento più basso possibile è desiderabile per un sistema
====Tempo di salita====
Si chiama \emph{tempo di salita} (o rise time) $t_{r}$ il tempo necessario alla risposta per portarsi dal 10\% al 90\% del valore di regime
====Tempo di ritardo====
Si chiama \emph{tempo di ritardo} (o delay time) $t_{d}$ il tempo necessario alla risposta per raggiungere il 10\% del valore di regime
===Comportamento a regime===
Se il sistema è un sistema in retroazione unitaria che ha come scopo che la sua uscita segua il segnale di comando in ingresso il più accuratamente possibile, allora
si studia la risposta del sistema a ciascuno delgli ingressi canonici
(gradino, rampa e parabola)
\vedilibro{rif:k}{371}
a transitorio esaurito
e si valuta la stabilità esterna del sistema,
ovvero se il sistema segue l'ingresso oppure diverge da esso
Se il sistema ha altri scopi
(ad esempio mantenere la derivata del segnale in uscita piccola \dots)
o diversa struttura
allora l'errore a regime deve essere definito in maniera differente
L'errore a regime può essere valutato con vari metodi, ad esempio con applicando il teorema del valore finale
====Errore a regime====
Si chiama \emph{errore a regime}
\vedilibro{rif:k}{365, sezione 7-3: Steady-state error}
(o steady-state error)
$\epsilon_{ss}$ o semplicemente $\epsilon$
è la differenza tra il segnale in uscita e quello in ingresso in condizioni di regime (per il tempo che tende all'infinito)
L'errore ad un ingresso $u(t)$ su può valutare considerando la funzione di trasferimento ingresso-errore in ciclo chiuso (o funzione di sensitività) del sistema
$S(s) = 1 / (1 + G(s))$
moltiplicata per la trasformata di Laplace dell'ingresso $U(s)$
e quindi applicando il teorema del valore finale, per cui
<math>
\epsilon_{ss} = \lim_{s \rightarrow 0} \frac{1}{1 + G(s)} U(s)
</math>
Si definiscono le \emph{costanti di errore agli ingressi canonici}
\vedilibro{rif:k}{371}
come:
\begin{itemize}
\item
$K_{p} = \lim_{s \rightarrow 0} G(s)$: costante di errore al gradino
\item
$K_{v} = \lim_{s \rightarrow 0} sG(s)$: costante di errore alla rampa
\item
$K_{a} = \lim_{s \rightarrow 0} s^{2}G(s)$: costante di errore al alla parabola
\end{itemize}
====Grado dei sistemi====
Un sistema che è stabile in retroazione ha un errore a regime agli ingressi canonici nullo, finito o infinito dipendentemente dal suo numero di poli nulli $g$ (detto \emph{grado del sistema}) in anello aperto.
\vedilibro{rif:k}{375, tabella 7-1}
\begin{itemize}
\item
Un \emph{sistema di ordine zero} ha errore a regime finito non nullo
$1/(1+\lim_{s \rightarrow 0} G(s))$
al gradino unitario ed errore infinito per ingressi a rampa o parabolici
\item
Un \emph{sistema di ordine uno} ha errore a regime nullo al gradino, finito alla rampa
$1/\lim_{s \rightarrow 0} sG(s)$
e infinito alla parabola
\item
Un \emph{sistema di ordine due} ha errore a regime nullo per rampa e gradino ed errore finito alla parabola
$1/\lim_{s \rightarrow 0} s^{2}G(s)$
\end{itemize}
===Comportamento in frequenza===
====Guadagno}
====Banda}
\vedilibro{rif:b}{188, sezione 6.8: Azione filtrante dei sistemi dinamici}
====Margine di guadagno e margine di fase}
\vedilibro{rif:b}{319, sezione 11.7: Margine di guadagno e margine di fase}
\vedilibro{rif:k}{605, sezione 9-14: Gain margin and phase margin}
Si definisce \emph{pulsazione di taglio}
\vedilibro{rif:k}{608, sezione 9-14-1: Gain Margin}
\vedilibro{rif:b}{319, sezione 11.7.1: Margine di guadagno}
(o phase-crossover frequency) $\omega_{\pi}$ di un sistema
la pulsazione per cui la fase della risposta è in ritardo di $\pi$ radianti
(supponendo che questa sia unica, altrimenti si sceglie quella per cui il margine di guadagno è minore)
<math>
\angle G(j\omega_{\pi}) = -\pi
</math>
Il guadagno in corrispondenza della pulsazione di taglio è detto
\emph{margine di guadagno}
$GM$ (o gain-margin)
<math>
GM = -20 \log_{10} |G(j\omega_{\pi})| \phantom{5} \textrm{dB}
</math>
Per un sistema stabile il margine di guadagno è positivo.
Il margine di guadagno indica di quanto è possibile aumentare il guadagno d'anello di un sistema retroazionato prima che esso diventi instabile
Si definisce \emph{pulsazione critica}
\vedilibro{rif:b}{321, sezione 11.7.2: Margine di fase}
\vedilibro{rif:k}{610, sezione 9-14-2: Phase margin}
(o gain-crossover frequancy) $\omega_{c}$ di un sistema
(ed in sua corrispondenza la frequenza critica $f_{c}$)
la pulsazione per cui il modulo della risposta del sistema è unitario (0 dB)
(supponendo che sia unica, altrimenti si sceglie quella per cui il margine di fase è minore)
<math>
|G(j\omega_{c})| = 1
</math>
La differenza tra la fase in corrispondenza della frequenza critica e $-\pi$ è detto
\emph{margine di fase}
(o phase-margin) $PM$
<math>
PM = \angle G(j\omega_{c}) + \pi \phantom{5} \textrm{rad}
</math>
ed è definito in questo modo solo per sistemi a sfasamento minimo
\vedilibro{rif:k}{612, phase margin of nonminimum-phase systems}
Il margine di fase indica quale è lo sfasamento massimo che dei ritardi di tempo possono inserire nell'anello di reazione prima che il sistema diventi instabile
==Metodi di analisi==
| 