Dissipativity is an input-output (IO) characterization of nonlinear systems that enables compositional robust control through Vidyasagar's Network Dissipativity Theorem. However, determining the dissipativity of a system is an involved and, often, model-specific process. We present a general method to determine the local dissipativity properties of nonlinear, control affine systems. We simultaneously search for the optimal IO characterization of a system and synthesize a continuous piecewise affine (CPA) storage function via a convex optimization problem. To do so, we reformulate the relationship between the Hamilton-Jacobi inequality and the dissipation inequality as an linear matrix inequality (LMI) and develop novel LMI bounds for a triangulation. Further, we develop a method to synthesize a combined quadratic and CPA storage function to expand the systems the optimization problem is applicable to. Finally, we demonstrate that our method will always find a feasible IO characterization and a CPA or quadratic storage function given that the system is strictly locally dissipative.
academic
Analisi della Dissipatività Locale di Sistemi Nonlineari
La dissipatività è una caratteristica ingresso-uscita (IO) di sistemi nonlineari che consente il controllo robusto composizionale attraverso il teorema di dissipatività di rete di Vidyasagar. Tuttavia, la determinazione della dissipatività di un sistema è un processo complesso e generalmente dipendente dal modello specifico. Questo articolo propone un metodo generale per determinare le proprietà di dissipatività locale QSR (Quadratic Supply Rate) di sistemi nonlineari affini nel controllo. Il metodo ricerca simultaneamente le proprietà IO ottimali del sistema e sintetizza funzioni di immagazzinamento affini a tratti continue (CPA) attraverso un problema di ottimizzazione convessa. A tal fine, gli autori riformulano la relazione tra le disuguaglianze di Hamilton-Jacobi e le disuguaglianze di dissipatività come disuguaglianze matriciali lineari (LMI), sviluppando nuovi limiti LMI per triangolazioni. Inoltre, sviluppano metodi per sintetizzare funzioni di immagazzinamento combinate quadratiche e CPA, estendendo l'applicabilità del problema di ottimizzazione. Infine, dimostrano che il metodo trova sempre proprietà IO fattibili e funzioni di immagazzinamento CPA o quadratiche quando il sistema è strettamente dissipativo localmente.
La teoria della stabilità ingresso-uscita (IO) è cruciale nel controllo robusto, adottando un approccio composizionale per l'analisi della stabilità—la stabilità IO in anello chiuso della rete dipende dalle proprietà IO in anello aperto grossolane dei singoli sottosistemi. Tuttavia, per i sistemi nonlineari, la determinazione delle proprietà IO è un processo complesso e generalmente dipendente dal modello specifico.
Per i sistemi lineari, le proprietà IO possono essere stabilite attraverso analisi nel dominio della frequenza, risolvendo LMI derivati da HJI, o analisi guidata dai dati
Attraverso il lemma KYP, la dissipatività è equivalente all'esistenza di una funzione di immagazzinamento quadratica semidefinita positiva, risolvibile attraverso ottimizzazione LMI
Tuttavia, le funzioni di immagazzinamento per sistemi nonlineari mancano di forme funzionali fisse e non esiste un metodo sistematizzato
Problemi dei Metodi Nonlineari Esistenti:
Metodi di Programmazione Dinamica6,14: Metodi iterativi che richiedono buona inizializzazione e non sempre garantiscono convergenza
Metodi di Trasformazione26,33: Validi solo in regioni di spazio degli stati finite dove la trasformazione è efficace, richiedono tecniche specifiche
Metodi SOS25,27,1: Applicabili solo a sistemi polinomiali, dipendono da ottimizzazione non convessa che può essere eccessivamente conservativa
Metodi CPA Precedenti17: HJI e i suoi termini di errore sono polinomi nelle variabili di progettazione, portando a ottimizzazione non convessa conservativa, escludendo termini lineari e affini nel controllo
È necessario un metodo generale e sistematizzato per determinare le proprietà di dissipatività locale di sistemi nonlineari affini nel controllo, che dovrebbe:
Essere applicabile a una vasta classe di sistemi nonlineari
Realizzarsi attraverso ottimizzazione convessa
Fornire garanzie di fattibilità
Sintetizzare simultaneamente funzioni di immagazzinamento e proprietà IO
Framework di Ottimizzazione Convessa: Propone un metodo che ricerca simultaneamente i parametri QSR e la funzione di immagazzinamento CPA attraverso un problema di ottimizzazione convessa, garantendo che HJI valga nell'intera regione dello spazio degli stati
Nuovi Limiti LMI di Errore (Sezione 3):
Sviluppa nuovi limiti LMI di errore per imporre vincoli LMI su funzioni CPA su triangolazioni (Teorema 10)
Sviluppa metodi per forzare disuguaglianze su palle chiuse intorno all'origine, consentendo di considerare funzioni di immagazzinamento quadratiche (Teorema 11)
Riformulazione LMI di HJI (Sezione 4):
Riformula HJI come LMI, convessificando il problema di ottimizzazione precedente
Propone condizioni affinché la funzione di immagazzinamento CPA soddisfi HJI (Teorema 12)
Sviluppa metodi di sintesi per funzioni di immagazzinamento combinate quadratiche e CPA (Teorema 13)
Garanzie Teoriche (Sezione 5):
Dimostra che per sistemi nonlineari affini nel controllo strettamente dissipativi localmente, il problema di ottimizzazione verifica sempre proprietà di dissipatività con funzioni di immagazzinamento CPA o quadratiche (Teoremi 15, 16)
Verifica Numerica (Sezione 6):
Dimostra l'efficacia del metodo su tre sistemi, inclusi confronti con limiti analitici e metodi esistenti
Modello del Sistema: Considerare la mappatura G(x0):L2em→L2ep, dove y=G(x0)u è definita da:
{x˙=f(x)+Gˉ(x)uy=h(x)+Jˉ(x)u
dove Gˉ(x)=B+G(x), Jˉ(x)=D+J(x), e f(0)=0,G(0)=0,h(0)=0,J(0)=0.
Obiettivo: Determinare la dissipatività QSR locale del sistema su un insieme invariante positivamente robusto Ω⊂Rn, cioè l'esistenza di una funzione di immagazzinamento V:Ω→R+ tale che:
Triangolazione: Partizionare la regione dello spazio degli stati Ω in un insieme finito di n-simplessi T={σi}i=1mT
Definizione della Funzione CPA: Una funzione CPA su Ω è univocamente determinata dai suoi valori sui vertici della triangolazione EΩ. Per il simplesso σi=co{xi,j}j=0n, il valore della funzione è:
W(x)=x⊤Xi−1Wˉi
dove la j-esima riga di Xi è xi,j−xi,0, e il j-esimo elemento di Wˉi è Wxi,j−Wxi,0.
Teorema 12 (Funzione di Immagazzinamento CPA Pura, B=0, D=0):
Se la funzione CPA V={Vx}x∈ET soddisfa:
Vx≥0,∀x∈Ex (non-negatività)
∥∇Vi∥1≤li (limite del gradiente)
M(xi,j)+E(xi,j)⪯0 (condizione QSR con limite di errore)
allora V soddisfa la disuguaglianza di dissipatività HJI su X∘.
Teorema 13 (Funzione di Immagazzinamento Combinata Quadratica-CPA):
Per sistemi con termini affini nel controllo non nulli, definire la funzione di immagazzinamento combinata:
Contributo Metodologico: Propone il primo metodo sistematizzato che sintetizza simultaneamente funzioni di immagazzinamento CPA e determina la dissipatività QSR locale attraverso ottimizzazione convessa
Garanzie Teoriche: Dimostra che per sistemi strettamente dissipativi localmente, il metodo garantisce di trovare una soluzione fattibile (Teoremi 15-16)
Innovazioni Tecniche:
Nuovi limiti di errore LMI che consentono di implementare vincoli CPA attraverso condizioni finite sui vertici
Funzione di immagazzinamento combinata quadratica-CPA che estende l'applicabilità del metodo
Valore Pratico: Gli esperimenti numerici verificano l'efficacia del metodo su diverse dimensioni e tipi di sistemi
Contributo Rivoluzionario: Prima convessificazione del problema di sintesi della funzione CPA, progresso importante nel campo
Limite di Errore LMI: La progettazione strutturata del Teorema 10 utilizza abilmente la struttura matriciale, evitando la non convessità del lavoro precedente
Funzione di Immagazzinamento Combinata: Soluzione elegante al problema dei termini affini nel controllo all'origine
12 Hill & Moylan (1976): "The stability of nonlinear dissipative systems" - Lavoro fondamentale della teoria della dissipatività
8,9 Giesl & Hafstein (2012, 2014): Lavori pioneristici sulla sintesi di funzioni di Lyapunov CPA, base importante di questo articolo
25 Summers et al. (2013): "Quantitative local L2-gain and reachability analysis" - Rappresentante del metodo SOS
17 Lavaei & Bridgeman (2022): Lavoro precedente sulla funzione di immagazzinamento CPA, miglioramento diretto di questo articolo
28 Van Der Schaft (1992): "L2-gain analysis of nonlinear systems" - Letteratura classica sull'analisi del guadagno L2 di sistemi nonlineari
Riepilogo: Questo è un lavoro teorico di alta qualità che fornisce importanti contributi all'analisi della dissipatività di sistemi nonlineari. La convessificazione della sintesi della funzione di immagazzinamento CPA è un'innovazione chiave, e le garanzie teoriche aumentano l'affidabilità del metodo. Le principali limitazioni riguardano la complessità computazionale dei sistemi ad alta dimensione e una certa conservatività rispetto ai metodi manuali. Questo lavoro fornisce uno strumento sistematizzato e prezioso per il controllo robusto nonlineare, con impatto previsto significativo nel campo della teoria del controllo.