Fondamenti di automatica/Sistemi lineari tempoinvarianti
Descrizione con equazioni differenziali
modificaSi può creare un modello matematico di molti sistemi per mezzo delle equazioni differenziali [1] Trattiamo solo sistemi modellizzabili con equazioni differenziali lineari a coefficienti costanti [2]
In generale un sistema lineare tempoinvariante SISO rappresentato con equazioni differenziali avrà uno stato interno descritto da equazioni del tipo dove è la funzione di stato, i coefficienti sono numeri reali e è la funzione di ingresso. L'uscita del sistema è data in generale come funzione lineare dello stato e dell'ingresso
Il polinomio caratteristico dell'equazione è
La soluzione generale dell'equazione sarà composta da una componente data dalla soluzione dell'equazione omogenea associata (dove ) più una componente detta soluzione particolare che sarà dello stesso tipo di
Per ottenere una soluzione unica, è necessario conoscere tante condizioni () su o sulle sue derivate pari al grado dell'equazione
Sistemi di equazioni differenziali
modificaÈ possibile scomporre sempre una o più equazioni differenziali in un sistema di equazioni differenziali del primo ordine [3] del tipo dove e sono scalari reali che ha soluzione [4]
Diagramma analogico
modificaÈ possibile rappresentare un sistema in una forma a blocchi detta diagramma analogico [5] utilizzando blocchi che rappresentano sistemi semplici:
- Amplificatore con guadagno :
- Integratore
- Sommatore
Esistono due (???) strutture standard
Se si considera l'equazione differenziale che descrive il sistema espressa come
è possibile utilizzare un sommatore per ottenere , integratori in serie per ottenere e , dei guadagni e per ottenere i termini a secondo membro dell'equazione; si porta nel sommatore le uscite degli integratori moltiplicate per i relativi coefficienti e si manda l'ingresso nel sommatore (anch'esso moltiplicato per i suoi coefficienti)
L'uscita del sistema si ottiene sommando le uscite degli integratori, appositamente moltiplicate
Un sistema rappresentato in questo modo, portato in variabili di stato, viene ad essere espresso in forma canonica di controllabilità. Per ricavare rapidamente il diagramma dall'equazione differenziale, si esprimono le derivate con , quindi si divide per più volte per portarle tutte a denominatore, si raccoglie e si ottiene un'equazione che rappresenta la struttura del circuito
Variabili di stato
modificaUn sistema MIMO lineare tempoinvariante si descrive facilmente utilizzando le matrici
Un sistema può essere descritto da un'equazione differenziale o da un sistema di equazioni differenziali in cui le funzioni incognite sono dette variabili di stato e sono il minimo numero di variabili tali che la conoscenza di queste ad un istante è sufficiente a determinare lo stato del sistema per ogni istante successivo
i valori rappresentano le condizioni iniziali del sistema
Un sistema lineare può essere espresso utilizzando le matrici
dove le matrici sono tali che , , , ; la matrice è detta matrice di stato.
Lo stato in funzione del tempo (o movimento dello stato) [6] [7] è in generale
dove è la matrice di transizione dello stato di un sistema lineare.
L'uscita del sistema è
Gli stati di equilibrio di un sistema lineare in variabili di stato [8] sono le soluzioni dell'equazione quando l'ingresso è costante.
Se la matrice è invertibile allora lo stato di equilibrio è unico e corrisponde allo stato e all'uscita
la matrice rappresenta il guadagno statico del sistema [9] che nel caso di sistemi SISO si indica con e costituisce il rapporto tra l'uscita e l'ingresso quanto tutte le variabili del sistema, ingresso e stato, sono costanti
Soluzione del sistema in variabili di stato
modificaL'unico passo non banale nella soluzione di un sistema in variabili di stato consiste nella valutazione della matrice di transizione dello stato [10] che è possibile risolvere utilizzando la trasformata di Laplace o portando la matrice in forma diagonale.
Se consideriamo il sistema a ingresso nullo e trasformiamo entrambi i membri, si ottiene , ovvero ; e quindi antitrasformando si trova che
Se consideriamo il sistema con ingresso non nullo, con lo stesso metodo, otteniamo
Utilizzando le trasformazioni lineari e portando la matrice di stato in forma diagonale o a blocchi di Jordan, la valutazione di è più semplice.
Matrici
modificaAlcuni richiami sulle matrici ...[11]
Gli autovalori di una matrice sono le radici dell'equazione
radici multiple corrispondono ad autovalori con molteplicità algebrica maggiore di uno.
Ad ogni autovalore corrisponde un autovettore tale che
Per ogni autovalore di molteplicità algebrica esistono autovettori generalizzati che si determinano ripetendo volte l'equazione precedente sostituendo a 0 l'opposto dell'ultimo autovettore trovato.
Si dice molteplicità geometrica di un autovalore il numero di autovettori linearmente indipendenti ad esso associati corrispondenti alla dimensione del nullo di A () (?)
Il rango di una matrice è pari al numero di righe linearmente indipendenti, valgono le seguenti proprietà:
- è uguale al numero di colonne o righe linearmente indipendenti
- è uguale al rango della matrice trasposta di :
- è minore del numero di colonne e del numero di righe della matrice
- se la matrice è quadrata e il suo rango è massimo (), allora è invertibile
- il rango è inveriante rispetto a pre o post moltiplicazioni per matrici non singolari (come le matrici in una similitudine)
- il rango del prodotto di due matrici è minore dei ranghi di ciascuna delle due
Cambiamento di coordinate
modificaLe matrici , , , (matrici rappresentative del sistema) rappresentano una trasformazione lineare (similitudine) dallo spazio di e a quello di [12] [13]
È possibile cambiare sistema di riferimento applicando una trasformazione lineare alla matrici rappresentative del sistema [14]
Applicando una trasformazione al sistema [15] si ottiene
dove è una matrice quadrata e è un nuovo vettore di stato tale che .
Una matrice quadrata rappresenta la stessa trasformazione lineare di una matrice se non è singolare [16] , due matrici simili hanno stesso polinomio caratteristico e stessi autovalori.
È possibile scegliere in modo da rendere in una forma che consenta di effettuare meglio calcoli, [17] diagonale [18] (si indica con ) o in forma di Jordan [19] ()
La matrice è diagonale e gli elementi sulla diagonale sono gli autovalori di
La matrice a blocchi di Jordan è diagonale, tranne degli 1 al di sopra degli elementi della diagonale (che corrispondono ad autovettori linearmente dipendenti), e gli elementi sulla diagonale sono gli autovalori di e gli sulla diagonale superiore sono in corrispondenza degli autovettori generalizzati che sono linearmente dipendenti (???).
Per portare una matrice in forma diagonale o di Jordan si usa la trasformazione
dove sono gli autovettori (intesi come vettori-colonna) della matrice (si può costruire anche con gli autovettori sinistri messi per righe)
Esponenziale di una matrice
modificaPer ogni matrice quadrata e ogni scalare è definita la matrice esponenziale [20]
estendendo la definizione di esponenziale come serie valida per gli scalari.
Per il teorema di Cayley-Hamilton [21], qualunque matrice quadrata soddisfa la sua equazione caratteristica (il suo polinomio caratteristico eguagliato a zero)
e quindi la serie necessaria per il calcolo dell'esponenziale si riduce ad una somma finita, in quanto tutte le potenze da in poi possono essere espresse come combinazione lineare di
dove i coefficienti si possono ricavare utilizzando il determinante di Vandermonte
È possibile calcolare l'esponenziale di una matrice avvalendosi della trasformata di Laplace
Oppure se la matrice è in forma diagonale e sono i suoi autovalori (e quindi gli elementi della diagonale principale) allora
Se la matrice è in forma di Jordan, per ogni blocco di Jordan che si ottiene con ognuno degli autovalori il suo esponenziale è
Forma canonica di controllo
modificaDato un sistema in variabili di stato di definisce forma canonica di controllo [22] il sistema simile al sistema dato in cui
Dove si suppone che (eventualmente si aggiungono coefficienti nulli)
Questa forma è tale che per qualunque valore di il sistema risulta controllabile; può essere ricavata direttamente dal diagramma analogico applicato all'equazione differenziale che descrive il sistema (????)
Forma canonica di osservabilità
modificaDato un sistema in variabili di stato di definisce forma canonica di osservabilità [23] il sistema simile al sistema dato in cui
Dove si suppone che (eventualmente si aggiungono coefficienti nulli).
Questa forma è tale che per qualunque valore di il sistema risulta osservabile.
Funzione di trasferimento
modificaLa funzione di trasferimento ingresso-uscita di un sistema lineare SISO [24] è definita come la funzione tale che
dove e rappresentano la trasformata di Laplace dell'ingresso e dell'uscita ; è in generale esprimibile da un rapporto di polinomi in , ma sono possibili varie forme [25]
- Funzione di trasferimento in forma polinomiale
- Funzione di trasferimento in forma poli-zeri
- Funzione di trasferimento in forma di Bode
I simboli utilizzati rappresentano:
- sono detti poli della funzione di trasferimento;
- sono gli zeri;
- è il numero dei poli e il numero degli zeri;
- e sono gli inversi dei poli e degli zeri non nulli (detti anche costanti di tempo) e ;
- è il grado del sistema, il numero di poli nulli;
- rappresenta qui il guadagno generalizzato
[26] [27] del sistema, ovvero il guadagno per frequenze minori delle frequenze dei poli e degli zeri del sistema (se nel sistema non si considerano i poli a frequenza nulla).
I poli sono gli opposti delle del polinomio caratteristico ed influenzano la stabilità del sistema, poli negativi conducono a radici positive.
Il denominatore della funzione di trasferimento rappresenta il polinomio caratteristico del sistema.
Gli zeri di un sistema rappresentano gli ingressi per cui il sistema ha uscita nulla, questa è detta proprietà bloccante degli zeri [28]
Poiché i coefficienti e sono reali, i poli e gli zeri di un sistema (che sono considerati talvolta le radici dei polinomi a numeratore e a denominatore della funzione di trasferimento, talvolta i loro opposti) sono reali oppure complessi coniugati.
Un sistema è detto a sfasamento minimo se non ha poli o zeri negativi (instabili) o immaginari puri (eccetto siano nulli); in questo caso la variazione di fase totale tra la pulsazione nulla e infinita è , il valore della funzione di trasferimento non è mai zero o infinito per pulsazioni limitate non nulle.
La funzione di trasferimento è invariante rispetto a trasformazioni lineari [29] e dipende dalla sola parte raggiungibile ed osservabile del sistema [30]
La trasformata di Laplace
modificaÈ definita per funzioni reali la trasformata di Laplace [31]
come funzione della variabile complessa secondo la formula
deve essere definita almeno per
Il minimo valore di , indicato con per cui la trasformata esiste è detto ascissa di convergenza.
L'operazione trasforma una funzione in un altro spazio le cui basi sono le funzioni ovvero
- sono quindi trasformabili tutte le funzioni che possono essere espresse come somma di funzioni sinusoidali esponenzialmente smorzate, che includono le funzioni esponenziali semplici, lineari e sinusoidali
L'operazione inversa è dove è un qualsiasi valore maggiore di
Teorema del valore finale
modificaSe ha trasformata razionale (o anche solo se esiste il limite a primo membro) con grado del denominatore maggiore del grado del numeratore e poli negativi o nulli, allora [32]
Teorema del valore iniziale
modificaSe ha trasformata razionale (o anche solo se esiste) con grado del denominatore maggiore del grado del numeratore, allora [33] deve essere minore di 0 altrimenti non è definita in 0 e il teorema non è applicabile
Risposta impulsiva
modificaLa risposta impulsiva di un sistema è definita come l'uscita del sistema quando si ha in ingresso una delta di Dirac su un singolo ingresso [34]; si definisce risposta all'impulso dello stato e risposta all'impulso dell'uscita
Si definisce anche la risposta al gradino , la cui derivata è la risposta impulsiva, come la risposta del sistema ad un ingresso a gradino unitario dove per analogia con gli scalari
Risposta in frequenza
modificaSe si restringe la funzione di trasferimento di un sistema sui valori di immaginari puri positivi, si ottiene la risposta in frequenza del sistema [35] che rappresenta l'uscita del sistema quando ha in ingresso una sinusoide.
Essa è definita per tutti i valori della pulsazione positivi che non siano poli immaginari puri di
Quindi se , allora si ha che
La risposta in frequenza è una funzione a immagine complessa ; il modulo della risposta in frequenza rappresenta la variazione del modulo dell'uscita in funzione della pulsazione dell'ingresso; la fase della risposta in frequenza rappresenta la variazione della fase dell'uscita in funzione della pulsazione dell'ingresso
si può rappresentare in varie forme corrispondenti a quelle della funzione di trasferimento, le più comuni sono:
- Risposta in frequenza in forma poli-zeri
- Risposta in frequenza in forma di Bode
- Risposta in frequenza in modulo e fase
se si vuole valutare la risposta in frequenza per pulsazioni lontane dai poli e dagli zeri, è possibile trascurare poli e zeri a pulsazioni elevate (tali che ) e semplificare l'espressione di poli e zeri a pulsazioni inferiori (tali che ) ottenendo la risposta in frequenza semplificata che è valida per pulsazioni ;
Se si vuole valutare il modulo a frequenze maggiori di un numero pari di poli e zeri (); ovvero in corrispondenza di un tratto a guadagno costante, è possibile semplificare anche dall'espressione della risposta in ampiezza, ottenendo la forma ulteriormente semplificata
La trasformata di Fourier
modificaLa trasformata di Fourier [36] è definita come la trasformata di laplace per funzioni reali ma prendendo come base le funzioni sinusoidali
l'operazione inversa è
Passaggio tra le varie rappresentazioni
modificaDalle equazioni differenziali alle variabili di stato
modificaSe il sistema è descritto da un'equazione differenziale e se si sceglie come variabili di stato
il sistema viene descritto in forma canonica di osservabilità [38]
Dalle variabili di stato alla funzione di trasferimento
modificaSe si considera il sistema espresso in variabili di stato in forma canonica e si applica la trasformata di Laplace a entrambi i membri delle due equazioni, con alcuni passaggi ( poiché non è identicamente nullo per ogni ) si ottiene: [39]
La funzione di trasferimento rappresenta l'uscita divisa per l'ingresso con stato iniziale nullo, quindi in questo caso, poiché il sistema in variabili di stato è in generale MIMO, la funzione di trasferimento è una matrice con termini che sono funzioni razionali fratte con al denominatore il polinomio caratteristico della matrice
Conoscendo le matrici è possibile ricavare dalla funzione di trasferimento ingresso-uscita anche la funzione di trasferimento ingresso-stato
Dalle equazioni differenziali alla funzione di trasferimento
modificaDalle variabili di stato alla risposta impulsiva
modificaConoscendo le matrici rappresentative del sistema , la risposta impulsiva [41] si ricava in modo simile alla funzione di trasferimento è possibile ricavare la risposta impulsiva anche antitrasformando la funzione di trasferimento
Note
modifica- ↑ Automatic Control Systems di Benjamin C. Kuo; pag 25, sezione 2-3: Differential equations
- ↑ Calcolo, volume terzo - Analisi 2 di Tom M. Apostol; pag 34, sezione 2.4: Equazioni differenziali lineari di ordine
- ↑ Automatic Control Systems di Benjamin C. Kuo; pag 26, sezione 2-3-3: First-order differential equations
- ↑ Calcolo, volume terzo - Analisi 2 di Tom M. Apostol; pag 108, sezione 3.16: Sistemi lineari non omogenei a coefficienti costanti
- ↑ Controls systems engineering di I. J. Nagrath, M. Gopal; Wiley International edition, 1982, seconda edizione; pag 462, figura 12.7
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag. 60, sezione 3.2.1: Formula di Lagrange
- ↑ Controls systems engineering di I. J. Nagrath, M. Gopal; Wiley International edition, 1982, seconda edizione; pag. 481, sezione 12.5: Solutions of state equations
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag. 70, sezione 3.2.5: Equilibrio
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag. 70
- ↑ Controls systems engineering di I. J. Nagrath, M. Gopal; Wiley International edition, 1982, seconda edizione; pag. 485, Computation of state transition matrix
- ↑ Automatic Control Systems di Benjamin C. Kuo; Prentice Hall, settima edizione 1995; pag. 245, sezione 5-7: Characteristic equation, eigenvalues, eigenvectors
- ↑ Calcolo, volume secondo - Geometria di Tom M. Apostol; Bollati Boringhieri, ristampa del febbraio 1996; pag. 125, Capitolo 4: Trasformazioni lineari e matrici
- ↑ Controls systems engineering di I. J. Nagrath, M. Gopal; Wiley International edition, 1982, seconda edizione; pag. 472, sezione 12.4: Diagonalization
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag. 61, sezione 3.2.2: Rappresentazioni equivalenti
- ↑ Controls systems engineering di I. J. Nagrath, M. Gopal; Wiley International edition, 1982, seconda edizione
- ↑ Calcolo, volume secondo - Geometria di Tom M. Apostol; Bollati Boringhieri, ristampa del febbraio 1996; pag. 210, Sezione 6-9: Matrici simili
- ↑ Calcolo, volume secondo - Geometria di Tom M. Apostol; Bollati Boringhieri, ristampa del febbraio 1996; pag. 224, sezione 7.9: Diagonalizzazione di matrici Hermitiane o anti-Hermitiane
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag.63, sezione 3.2.3: Matrice della dinamica diagonalizzabile
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag.65, sezione 3.2.3: Matrice della dinamica diagonalizzabile
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag. 591, appendice A.5: Esponenziale
- ↑ Calcolo, volume terzo - Analisi 2 di Tom M. Apostol; Bollati Boringhieri, settima impressione del settembre 1991; pag.98, sezione 3.11: Il teorema di Cayley-Hamilton
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag.125, sezione 4.5.1: Forma canonica di raggiungibilità
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag.125, sezione 4.5.2: Forma canonica di osservabilità
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag.99, capitolo 4: Funzione di trasferimento
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag. 106, sezione 4.3: Rappresentazioni e parametri della funzione di trasferimento
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag. 70, sezione 3.2.5: Equilibrio
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag. 106, sezione 4.3.1: Guadagno
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag. 155, sezione 6.2.2: Proprietà bloccante degli zeri
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag. 102, sezione 4.2.3: Invarianza della funzione di trasferimento
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag. 102, sezione 4.2.4: Cancellazioni e stabilità
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag. 598, appendice B.3: Trasformata di Laplace
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag.603
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag. 603
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag. 68, sezione 3.2.4: Risposta all'impulso e movimento forzato
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag. 154, capitolo 6: Risposta in frequenza
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag. 617, appendice B.5: Trasformata di Fourier
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag. 129, figura 4.16: Rappresentazioni dei sistemi dinamici e relazioni corrispondenti
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag. 127, sezione 4.5.4: Interpretazione delle forme canoniche
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag. 99, sezione 4.2.1: Definizione e proprietà della FdT
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag. 129, sezione 4.5.4
- ↑ Fondamenti di controlli automatici di Paolo Bolzern, Riccardo Scattolini, Nicola Schiavoni; McGraw-Hill, prima edizione del marzo 1998; pag. 68, sezione 3.2.4: Risposta all'impulso e movimento forzato