Sistema dinamico lineare stazionario

In teoria dei sistemi, un sistema dinamico lineare stazionario, anche detto sistema lineare tempo-invariante o sistema LTI, è un sistema dinamico lineare tempo-invariante, soggetto cioè al principio di sovrapposizione degli effetti e tale che il suo comportamento sia costante nel tempo. Si tratta di un modello matematico che riveste particolare importanza in numerose applicazioni, in particolare in elettronica e nella teoria del controllo.

Descrizione

modifica

Un sistema stazionario (o tempo invariante) è un sistema i cui parametri non dipendono dal tempo. Il processo fisico di cui il sistema è il modello matematico, pertanto, è un sistema di equazioni differenziali, derivate rispetto al tempo, a coefficienti costanti:

 

dove  ,  ,   e   sono vettori colonna. Il vettore   rappresenta le variabili di stato in funzione del tempo  , che in generale non possono essere fissate né osservate direttamente, il vettore   rappresenta le variabili di stato all'istante iniziale  ,   sono gli ingressi, cioè le variabili su cui si agisce per modificare l'andamento o traiettoria dello stato, e   sono le uscite, cioè le variabili misurate da cui si deduce, a seconda delle caratteristiche di osservabilità del sistema, il valore o la stima dello stato. Ci possono essere particolari variabili di ingresso, dette disturbi o rumori, su cui non si può agire in alcun modo. Il termine   è inoltre la derivata in   di  , e le funzioni   e   non dipendono direttamente da  .

Un sistema è inoltre lineare quando dipende linearmente dalle variabili di stato e dalle variabili di ingresso:

 

dove  ,  ,   e   sono matrici di dimensioni opportune che premoltiplicano   e  . In generale esse possono variare nel tempo, ma non nel caso di un sistema stazionario:

 

Tra le caratteristiche dei sistemi LTI più studiate vi sono la stabilità, le proprietà di raggiungibilità e di osservabilità: se sono verificate allora per il sistema di controllo (cioè il sistema ottenuto retroazionando il sistema dinamico LTI con un controllore LTI) esiste sempre un controllore che rende il sistema di controllo asintoticamente stabile.

La risposta in frequenza dei sistemi LTI può essere studiata a partire dalle caratteristiche della funzione di trasferimento, una funzione complessa della quale il comportamento dei poli è sintomatico della stabilità del sistema che descrive.

Un sistema lineare stazionario è particolarmente importante perché, oltre ad offrire innumerevoli risultati pratici e teorici, si usa spesso per linearizzare anche sistemi non lineari o non stazionari in modo da facilitare il calcolo e le applicazioni. Nel caso di variabili continue i sistemi lineari e stazionari sono descritti da equazioni algebriche nel dominio del tempo se statici, altrimenti si hanno equazioni differenziali ordinarie se dinamici. Inoltre, i sistemi lineari e stazionari possono essere studiati anche nel dominio della frequenza.

Nel caso generale e con la sola dipendenza da una variabile temporale, sia   una qualsiasi sollecitazione di ingresso. Sia   un operatore che riassume tutte le operazioni che il sistema può compiere sulla sollecitazione di ingresso  . Allora la relazione che lega ingresso e uscita di un sistema è in generale:

 

I sistemi lineari sono soggetti al principio di sovrapposizione, ovvero un sistema è lineare se valgono le seguenti proprietà:

 

dove   è un numero arbitrario. I sistemi tempo invarianti, anche detti stazionari o statici, sono inoltre quei sistemi per i quali la risposta dipende solo dai valori istantanei dell'ingresso:

 

anche nel caso in cui i parametri del sistema sono indipendenti dal tempo.

Esistono anche sistemi statici in elettronica digitale e sono chiamati combinatori. In contrapposizione esistono sistemi dinamici lineari nei quali l'uscita è dipendente sia dai valori istantanei dell'ingresso che dalla storia passata del segnale in ingresso. Allo stesso modo in elettronica digitale esistono sistemi dinamici che sono chiamati sequenziali. In elettronica, tra i sistemi lineari sono notevolmente importanti elementi circuitali quali resistori, condensatori, induttori, mentre tra i sistemi non lineari vi sono il diodo e i transistor.

Sistemi a tempo continuo

modifica

L'uscita   di un sistema dinamico lineare tempo-invariante a tempo continuo soggetto a un segnale in ingresso   è descritta dalla convoluzione:

 

dove   è la risposta del sistema all'impulso, ovvero quando l'ingresso   è una funzione a delta di Dirac. L'uscita   è quindi proporzionale alla media dell'ingresso   pesata dalla funzione  , traslata di un tempo  .

Se la funzione   è nulla quando   allora   dipende soltanto dai valori assunti da   precedentemente al tempo  , ed il sistema è detto causale.

Per mostrare come la risposta all'impulso determini completamente il comportamento del sistema LTI, sia   l'azione del sistema al tempo  . Per l'invarianza temporale si ha:

 

da cui:

 

in modo che si ottiene:

 

Funzione di trasferimento

modifica

Un'autofunzione   di un operatore lineare   è una funzione che viene trasformata dall'operatore nella stessa funzione moltiplicata per un numero  , detto autovalore:

 

Per un sistema LTI a tempo continuo le autofunzioni sono le funzioni esponenziali  , con   e   in  . Infatti, sia   l'ingresso e   la risposta del sistema alla delta di Dirac. L'uscita è data da:

 

La trasformata di Laplace:

 

è la funzione di trasferimento del sistema, che permette così di ottenere gli autovalori a partire dalla risposta all'impulso di Dirac. Per ogni   e   in   l'uscita è dunque il prodotto dell'ingresso   per una costante dipendente solo dal parametro  , autovalore del sistema LTI relativo all'autovettore   (elemento di uno spazio vettoriale funzionale). Di particolare interesse è il caso in cui l'ingresso è un esponenziale complesso  , con   e  . La funzione di trasferimento è data in tal caso dalla trasformata di Fourier:

 

Mentre la trasformata di Laplace è utilizzata per segnali che sono nulli prima di un certo tempo  , solitamente lo zero, la trasformata di Fourier consente di trattare funzioni di durata infinita, con la richiesta (a differenza della trasformata di Laplace in sistemi stabili) di essere a quadrato sommabili.

Grazie alle proprietà della convoluzione, nel dominio della trasformata l'integrale si riduce ad una moltiplicazione:

 

Tale fatto consente di trasformare le equazioni differenziali ed integrali che solitamente governano i sistemi dinamici LTI in equazioni algebriche.

Risposta nel dominio della frequenza

modifica
 
Descrizione di un sistema LTI nel dominio del tempo (in blu) e nel dominio delle frequenze (la trasformata di Laplace è mostrata in rosso).

I sistemi lineari e stazionari possono essere studiati nel dominio della frequenza analizzandone la risposta a ingressi sinusoidali puri, la cui frequenza non viene cambiata in seguito alla trasformazione lineare compiuta dal sistema (ad esempio, la derivazione o l'integrazione del segnale). Questo permette di rappresentare un segnale periodico come combinazione lineare di segnali sinusoidali tramite la serie di Fourier. Nel caso di funzioni non periodiche si utilizzano la trasformata di Fourier o la trasformata di Laplace.

Lo studio dei sistemi lineari e stazionari nel dominio della frequenza passa attraverso i metodi simbolico e/o il metodo operatoriale, utili anche per studiare sistemi in cascata. Lo scopo è determinare una funzione di trasferimento che determina completamente la risposta del sistema.

Il legame tra la risposta nel dominio del tempo e la risposta nel dominio della frequenza è di importanza notevole. Tali relazioni si ricavano esattamente solo in casi semplici, e sono in particolare i legami tra gli ingressi impulsivi o unitari e le varie frequenza di taglio o di risonanza, e i valori dell'ampiezza e fase in termini di frequenza. Ancora più in particolare, ci sono connessioni semplici tra i tempi di salita dei segnali, la banda passante e la fase di un sistema. Nel caso di sistemi del primo ordine sono connessioni esatte, nel caso di sistemi del secondo ordine o superiori al secondo sono utili per un'approssimazione.

Sistemi a tempo discreto

modifica
  Lo stesso argomento in dettaglio: Sistema dinamico lineare stazionario discreto.

Un sistema a tempo discreto trasforma la successione in ingresso   in un'altra successione  , data dalla convoluzione discreta con la risposta   alla delta di Kronecker:

 

Gli elementi di   possono dipendere da ogni elemento di  . Solitamente   dipende maggiormente dagli elementi in prossimità del tempo  .

La maggior parte dei segnali a tempo discreto sono ottenuti da un segnale a tempo continuo considerandone il valore assunto in precisi istanti di tempo, solitamente separati da un intervallo temporale fisso  . La procedura che permette di ottenere un segnale discreto a partire da uno continuo è detta campionamento, ed è alla base della conversione analogico-digitale (ADC). Essa trasforma una funzione continua   nel segnale discreto:

 

con   la frequenza di campionamento. Il teorema del campionamento pone un limite alla massima frequenza del segnale continuo, che non può essere superiore ad   se si vuole evitare perdita di informazione (fenomeno di aliasing).

Come nel caso di sistemi a tempo continuo, se   è l'operatore di trasformazione al tempo  :

 

la successione:

 

caratterizza completamente il sistema. Per mostrare questo, considerando l'invarianza temporale:

 

e dato che vale l'identità:

 

si ha:

 

L'operatore   restituisce un'uscita proporzionale alla media pesata di   con funzione peso data da  . Se   per valori di   negativi il sistema è causale.

Funzione di trasferimento

modifica

Gli esponenziali del tipo  , con  , sono autofunzioni di un operatore lineare tempo-invariante. Infatti, detto   il periodo di campionamento e  , con   e   in  , si supponga   l'ingresso del sistema. Se   è la risposta impulsiva, si ha:

 

La funzione:

 

dipende solo dal parametro   ed è l'autovalore associato all'autovettore (autofunzione)   del sistema LTI.

La trasformata zeta:

 

è la funzione di trasferimento del sistema. Di particolare interesse è il caso in cui le autofunzioni siano sinusoidi pure  , con  , che possono essere scritte come  , dove  . Per funzioni di questo tipo la funzione di trasferimento è data dalla trasformata di Fourier a tempo discreto:

 

Grazie alle proprietà della convoluzione, nel dominio della trasformata si ottiene una moltiplicazione:

 

che analogamente al caso continuo risulta di notevole utilità nell'analisi dei sistemi LTI.

Descrizione matriciale

modifica

Un sistema LTI è descritto da un'equazione del tipo:

 

in cui  ,  ,   e   non sono funzione del tempo,  ,   e  ; inoltre la matrice   ha dimensione  ,   ha dimensione  ,   ha dimensione   e   ha dimensione  .

Soluzione dell'equazione matriciale

modifica

Volendo risolvere la precedente equazione bisogna distinguere i seguenti casi:

  •   ammette soltanto autovalori reali con molteplicità algebrica uguale alla molteplicità geometrica per ogni autovalore.
  •   ammette soltanto autovalori complessi coniugati.
  •   ammette sia autovalori reali che complessi coniugati.

Autovalori reali con molteplicità algebrica pari alla molteplicità geometrica

modifica

Si consideri la trasformazione di coordinate:

 

con   matrice di dimensione  , le cui colonne sono gli autovettori di   linearmente indipendenti che generano ciascun autospazio relativo ad ogni autovalore, e   vettore di dimensione  . Si ha che  , dove   è la matrice inversa di  , mentre:

 

Dalla teoria della diagonalizzazione delle matrici si ha che  , dove   è la matrice diagonale in cui sulla diagonale principale vi sono gli autovalori di   ripetuti, eventualmente ciascuno con la propria molteplicità. Pertanto si ottiene la seguente equazione differenziale matriciale:

 

In particolare, se gli autovalori di   sono reali e distinti sulla matrice diagonale   vi saranno gli   autovalori distinti di  .

Moltiplicando ambo i membri dell'equazione per la matrice esponenziale  , che sulla diagonale principale possiede gli esponenziali   (con   gli autovalori di  ), si ha la seguente equazione differenziale:

 

Integrando si ricava, scegliendo come primitive quelle che si annullano in   e moltiplicando l'equazione per  :

 

Quindi si ha:

 

da cui si ricava il vettore  . Pertanto la soluzione dell'equazione differenziale matriciale è:

 

Essendo   costante rispetto a  , si ottiene:

 

Si nota che la risposta libera nello stato, ottenuta ponendo  , è:

 

cioè basta moltiplicare la matrice   degli autovettori di  , la matrice esponenziale  , l'inversa di   ed il vettore di stato  .

La risposta forzata nello stato è invece ottenuta ponendo  , cioè:

 

Inoltre, la risposta libera nell'uscita per   è:

 

mentre la risposta forzata nell'uscita per   è:

 

Autovalori complessi coniugati

modifica

Volendo analizzare il caso in cui   ammette autovalori complessi coniugati, si supponga che essa sia una matrice di dimensione 2 e siano   e   i due autovalori complessi coniugati di  . Siano inoltre   e   i due autovettori complessi coniugati corrispondenti. Applicando la definizione di autovalore e di autovettore si ha la seguente equazione algebrica:

 

dove   è la matrice identica di dimensione 2. Si possono separare parte reale e parte immaginaria nella forma:

 

Affinché l'equazione sia vera è necessario che parte reale e parte immaginaria si annullino entrambe, pertanto si ha il sistema:

 

che può essere posto nella forma:

 

Se si pone   uguale alla matrice le cui colonne sono la parte reale e immaginaria dei due autovettori complessi coniugati si ha che:

 

Ragionando come nel caso degli autovalori reali e distinti si ottiene:

 

e quindi in tal caso la soluzione dell'equazione differenziale matriciale è:

 

Sviluppando in serie di Taylor la matrice esponenziale:

 

per l'identità di Eulero si ha che:

 

Per cui, sostituendo si ha:

 

Autovalori reali e complessi coniugati

modifica

Si supponga che la matrice   di ordine n ammetta   autovalori reali distinti   a cui corrispondono   autovettori distinti   Allora si hanno le seguenti equazioni:

 

Supponendo inoltre che   ammetta   coppie di autovalori complessi coniugati la cui p-esima coppia è   e  , a cui corrisponde la coppia di autovettori complessi coniugati   e  , allora per quanto visto nel caso precedente si ha per la  -esima coppia:

 

Posto   uguale alla matrice le cui colonne sono i   autovettori corrispondenti agli autovalori reali e le parti reali e immaginarie delle   coppie di autovettori complessi coniugati, cioè:

 

allora dalle precedenti equazioni si ha la matrice diagonale a blocchi:

 

pertanto:

 

Proprietà dei sistemi LTI

modifica

Lo stato   di un sistema LTI può essere esplicitato in funzione dell'ingresso   applicando la trasformata di Laplace all'equazione differenziale che lo definisce:

 

Da cui, trasformando e ipotizzando che  , si ha:

 

da cui:

 

e quindi:

 

essendo   e   le trasformate di   e  ,   la matrice unità di dimensione  , e   la matrice inversa di  . Lo stato può essere ricavato antitrasformando:

 

Poiché l'uscita del sistema è data da  , trasformando si ha:

 

cioè:

 

La matrice   è la matrice di trasferimento o funzione di trasferimento del sistema.

Esempio

modifica
 

Nel caso del circuito elettrico lineare mostrato in figura il vettore di stato   è costituito dalla corrente   che passa attraverso l'induttore di induttanza   e dalla tensione   ai capi del condensatore di capacità  , dove l'ingresso   è la tensione del generatore mentre il vettore delle uscite   è dato, ad esempio, dalle correnti che passano attraverso il resistore di resistenza   e resistore di resistenza  . Applicando le equazioni costitutive dei bipoli nonché le equazioni topologiche o leggi di Kirchhoff si ha:

 

Pertanto, sostituendo l'ultima relazione nelle precedenti e ponendo:

 

in tal caso si ha che:

 
 
 
 

Si supponga per esempio di voler determinare l'andamento della seconda variabile di stato a partire da un dato istante  , ipotizzando che il valore iniziale della stessa fosse nullo e l'andamento dell'ingresso coincida con un impulso di Dirac centrato in  . Nel dominio di Laplace l'ingresso ha dunque valore identicamente unitario, quindi avremo:

 
 
 

Pertanto:

 

Antitrasformando per passare al dominio del tempo:

 

Dove:

 
 

Dalla terza equazione topologico-costitutiva del modello matematico del sistema si ricava anche la prima variabile di stato:

 

Stabilità

modifica

Di un sistema LTI possono essere studiati diversi tipi di stabilità, come la stabilità interna o quella esterna. Facendo riferimento ai sistemi causali, ovvero nei sistemi in cui le uscite non dipendono dai valori futuri degli ingressi, la funzione di trasferimento ha un polinomio a denominatore di grado non inferiore al grado del polinomio a numeratore. Se gli zeri dei denominatori, che sono i poli della funzione di trasferimento, appartengono al semipiano a parte reale positiva del piano complesso, il sistema è instabile e la risposta all'impulso tende ad un valore infinito al crescere del tempo.

Se invece i poli della funzione di trasferimento appartengono al semipiano a parte reale negativa del piano complesso, il sistema è asintoticamente stabile e la risposta impulsiva tende asintoticamente a zero al crescere del tempo. Se, infine, i poli della funzione di trasferimento appartengono alla retta verticale a parte reale nulla del piano complesso ed hanno molteplicità singola, il sistema è semplicemente stabile e la risposta all'impulso è maggiorata in valore assoluto da un certo valore al crescere del tempo.

Per determinare come variano le posizioni dei poli e degli zeri al variare di particolari parametri, che generalmente rappresentano i guadagni ed altre caratteristiche associate al compensatore dinamico che si vuole progettare per stabilizzare il sistema, si usano particolari grafici, quali ad esempio il diagramma di Bode, il diagramma di Nyquist e il luogo delle radici.

Raggiungibilità

modifica

Un sistema lineare tempo invariante è raggiungibile se per ogni stato iniziale   lo stato generico   è raggiungibile, cioè se per ogni stato iniziale   esiste un ingresso   che permette al sistema di raggiungere lo stato generico  .

Il criterio di Kalman stabilisce che un sistema LTI di dimensione   è completamente raggiungibile se e solo se:

 

dove   indica il rango di  , che se è pari a   rende il determinante diverso da zero. Qualora risulti che   allora vi sono autovalori raggiungibili, quindi modificabili (in numero pari a  ) e autovalori non raggiungibili, detti autovalori fissi (in numero pari a  ). Il sistema si dice in tal caso non completamente raggiungibile.

Ad esempio, in un sistema con un ingresso ( ) ed una variabile di stato ( ) le matrici   e   si riducono a scalari e l'equazione di stato relativa è:

 

avendo indicato con   la derivata prima di   rispetto al tempo. Anche la matrice di raggiungibilità (o matrice di Kalman) è uno scalare  : se   il sistema è completamente raggiungibile poiché  , mentre se è nullo il sistema non è completamente raggiungibile e l'equazione di stato diventa:

 

Un'equazione differenziale di questo tipo è un sistema autonomo.

Un sistema lineare tempo invariante è quindi raggiungibile quando tutti i suoi stati sono raggiungibili, ovvero quando la matrice di raggiungibilità ha rango massimo. Un modo per stabilirlo è il test di Popov-Belevitch-Hautus, anche detto PBH test di raggiungibilità, il quale stabilisce che il rango della matrice ottenuta affiancando la matrice degli ingressi sullo stato   alla matrice  , a cui è sottratta la matrice della dinamica del sistema lineare  , deve essere pari al numero totale degli stati   al variare di  :

 

Esso deriva dalla trasformata di Laplace dell'equazione:

 

cioè:

 

da cui, essendo sempre possibile moltiplicare uno scalare per la matrice identità:

 

essendo   e   le trasformate di   e  . Lo stato del sistema, nel dominio di Laplace, è quindi definito come:

 

Stabilizzabilità

modifica

Un sistema lineare tempo invariante è stabilizzabile se esiste una matrice di retroazione dallo stato che rende asintoticamente stabile il sistema complessivo. Questo è possibile se e solo se:

  • Il sistema è completamente raggiungibile.
  • Il sistema non è completamente raggiungibile e gli autovalori non raggiungibili sono asintoticamente stabili.

Gli autovalori si dicono asintoticamente stabili se hanno parte reale negativa (nei sistemi a tempo continuo) o se hanno modulo minore di 1 (nei sistemi a tempo discreto). In particolare, la completa raggiungibilità di un sistema garantisce la stabilizzabilità in quanto il sottosistema non raggiungibile non esiste e quindi con la retroazione dallo stato è possibile allocare arbitrariamente gli autovalori della matrice  .

Osservabilità e rilevabilità

modifica

Un sistema si dice completamente osservabile se e solo se ogni stato non nullo genera un'uscita libera non identicamente nulla, ovvero quando la matrice di osservabilità ha rango massimo, ed in tal caso è verificato il PBH test di osservabilità. L'uscita di un sistema è somma del contributo dell'evoluzione libera, dipendente dallo stato iniziale, ed evoluzione forzata, dipendente unicamente dall'ingresso.

Un sistema LTI di dimensione   è completamente osservabile se e solo se:

 

oppure se è soddisfatto il PBH test di osservabilità:

 

Se i test di osservabilità di cui sopra falliscono, non necessariamente non si può osservare il sistema. Con riferimento alla decomposizione di Kalman rispetto alla osservabilità, se gli autovalori della parte inosservabile   si trovano già in una parte del piano complesso delimitata da un'ascissa passante per  , detto  , allora il sistema è detto   rilevabile.

Risposta al gradino di Heaviside

modifica

Se   è il gradino di Heaviside unitario di ingresso e   è l'insieme delle operazioni che il sistema effettua su tale ingresso, si definisce la risposta unitaria:

 

Ricordando che la funzione delta di Dirac è la derivata del gradino di heaviside:

 

e che quindi la risposta impulsiva è legata alla risposta unitaria da:

 

la risposta unitaria è data dall'integrale:

 

Nella rappresentazione dinamica dei segnali ogni segnale deterministico può essere rappresentato mediante gradini di Heaviside:

 

allora la risposta ad un tale segnale rappresentato secondo la funzione di Heaviside è:

 

Dopo aver ottenuto la funzione di trasferimento del sistema tramite l'applicazione della trasformata di Laplace, applicare un ingresso a gradino corrisponde a moltiplicare la F. di T. per  

Bibliografia

modifica
  • (EN) Phillips, C.l., Parr, J.M., & Riskin, E.A, Signals, systems and Transforms, Prentice Hall, 2007, ISBN 0-13-041207-4.
  • (EN) Hespanha,J.P., Linear System Theory, Princeton university press, 2009, ISBN 0-691-14021-9.
  • E. Fornasini, G. Marchesini, Appunti di Teoria dei Sistemi, Edizioni Libreria Progetto, Padova, 2003.
  • A. Ruberti, S. Monaco, Teoria dei Sistemi - Appunti dalle lezioni, Pitagora Editrice, Bologna, 1998.
  • O. M. Grasselli, Proprietà strutturali dei sistemi lineari e stazionari, Pitagora Editrice, Bologna, 1978.

Voci correlate

modifica

Collegamenti esterni

modifica
Controllo di autoritàLCCN (ENsh85077184 · GND (DE4213494-8 · J9U (ENHE987007529340805171
  Portale Controlli automatici: accedi alle voci di Wikipedia che trattano di Controlli automatici