Zonotope-Based Elastic Tube Model Predictive Control

Riepilogo Tecnico: Controllo Predittivo a Modello con Tubo Elastico Basato su Zono-topi

Enunciato del Problema
Il Controllo Predittivo a Modello (MPC) basato su tubi è un quadro robusto per sistemi lineari vincolati soggetti a disturbi additivi, che garantisce che tutte le traiettorie future rimangano all'interno di un "tubo" di insiemi. Sebbene efficace, la trattabilità numerica dell'MPC basato su tubi è spesso limitata dalla complessità delle parametrizzazioni degli insiemi. Gli approcci poliedrici tradizionali soffrono di una crescita quadratica delle variabili decisionali rispetto alla lunghezza dell'orizzonte o richiedono descrizioni complesse che aumentano le dimensioni del problema. I tubi omotetici (scalati) ed elastici (controllati tramite offset) offrono domini di attrazione ampliati, ma spesso introducono bilinearità o richiedono la pre-selezione di variabili ausiliarie per gestire le inclusioni di insiemi, portando a una selezione dei parametri subottimale o a oneri computazionali elevati, specialmente in sistemi ad alta dimensionalità.

Metodologia
Il documento propone un nuovo quadro di MPC a Tubo Elastico che utilizza zonotopi scalati per parametrizzare le sezioni trasversali del tubo. La metodologia fondamentale comprende:

1. Rappresentazione a Zonotopo Scalato: Il tubo è definito come una sequenza di zonotopi $\langle c_k, G\Delta_k \rangle$, dove $c_k$ è il centro, $G$ è una matrice generatrice fissa e $\Delta_k = \text{diag}(\delta_k)$ rappresenta fattori di scala variabili nel tempo. Ciò consente al tubo di allungarsi e restringersi elasticamente lungo direzioni specifiche.
2. Condizioni di Inclusione Lineari: Viene derivata una nuova classe di condizioni di inclusione per zonotopi scalati. A differenza dei metodi precedenti che si basano su termini bilineari o richiedono la pre-selezione di variabili ausiliarie (come $\Gamma$) per gestire la scalatura, l'approccio proposto utilizza una matrice $\Phi$ per codificare linearmente il contenimento delle somme di Minkowski di zonotopi.
   • Il Teorema 1 formula l'ottimizzazione MPC come un Programma Quadratico (QP) in cui i vincoli di raggiungibilità in un passo e di ammissibilità di stato/in ingresso sono espressi come uguaglianze e disuguaglianze lineari che coinvolgono i fattori di scala $\delta_k$ e la matrice ausiliaria $\Phi_k$.
   • Crucialmente, i centri degli zonotopi $c_k$ sono trattati come variabili decisionali indipendenti, aumentando la flessibilità senza introdurre non convessità.
3. Alternative di Pre-calcolo: Per ridurre ulteriormente la complessità di runtime, il documento introduce due varianti:
   • Basata su $\Phi_0$: Una matrice costante $\Phi_0$ viene pre-calcolata offline, sostituendo l'ottimizzazione di $\Phi_k$ ad ogni passo con un vincolo lineare sui fattori di scala.
   • Basata su $\Gamma$: Viene utilizzata una coppia pre-calcolata $(\Gamma, \gamma)$, simile agli approcci poliedrici precedenti ma adattata per gli zonotopi.
4. Insiemi Positivamente Invarianti Robusti (RPI): La matrice generatrice $G$ è derivata da un'approssimazione di un insieme RPI per la dinamica dell'errore. Il documento migliora questo aspetto trattando il centro dell'insieme RPI come una variabile decisionale nel problema di minimizzazione, riducendo il numero di variabili decisionali rispetto ai lavori precedenti.
5. Stabilità e Fattibilità: La fattibilità ricorsiva e la stabilità asintotica sono garantite costruendo un insieme terminale tramite ricorrenza di insiemi e imponendo condizioni di contrazione simili a quelle di Lyapunov sul costo terminale e sui fattori di scala.

Contributi Chiave

• Codifica Lineare Innovativa: Il documento introduce una codifica lineare per il contenimento degli zonotopi (Corollario 1) che generalizza le condizioni standard. Ciò elimina la necessità di pre-selezionare variabili ausiliarie per aggirare l'ottimizzazione non lineare, un problema comune nella scalatura poliedrica.
• Efficienza Computazionale: La formulazione proposta rimane lineare sia nelle variabili ausiliarie che nei fattori di scala, facilitando l'integrazione in framework di ottimizzazione su larga scala. Evita le bilinearità tipicamente associate ai fattori di scala in altri metodi di MPC a tubo.
• Riduzione della Complessità: Sfruttando gli zonotopi, il metodo riduce significativamente il numero di vincoli e fattori di scala richiesti rispetto ai tubi elastici poliedrici, specialmente all'aumentare della dimensionalità del sistema.
• Analisi Esplicita del Trade-off: Il documento fornisce un'analisi completa della complessità che confronta tubi rigidi, omotetici ed elastici nei contesti poliedrici e zonotopici, illustrando il compromesso tra il dominio di attrazione ampliato e lo sforzo computazionale.

Risultati
L'approccio proposto è stato validato attraverso esperimenti numerici su due benchmark:

1. Doppio Integratore: I confronti tra tubi elastici poliedrici, tubi elastici zonotopici con $\Phi_0$ pre-calcolato e tubi elastici zonotopici con $\Phi$ ottimizzato hanno mostrato che le varianti zonotopiche producono domini di attrazione più ampi. La variante con $\Phi$ completamente ottimizzato ha offerto un leggero vantaggio rispetto alla variante con $\Phi_0$ pre-calcolato, ma a un costo computazionale superiore.
2. Esperimento con Molle Accoppiate (CSE): Questo sistema ad alta dimensionalità (fino a 20 stati) ha dimostrato la scalabilità del metodo. Le implementazioni poliedriche non sono riuscite a trovare soluzioni fattibili per sistemi con 3 o più molle a causa dei limiti computazionali. Al contrario, tutte le varianti zonotopiche hanno risolto con successo il problema nell'intero intervallo di parametri.
   • Runtime: I metodi zonotopici hanno mostrato una complessità computazionale lineare (o quadratica per la variante $\Phi$), mentre i metodi poliedrici sono diventati intrattabili.
   • Centri: Trattare i centri degli zonotopi come variabili decisionali ha comportato un sovraccarico computazionale solo marginale, fornendo al contempo garanzie di robustezza cruciali.

Significato
Il documento afferma che l'integrazione di zonotopi scalati elasticamente nell'MPC a tubo riduce significativamente il runtime e migliora la scalabilità per sistemi ad alta dimensionalità rispetto ai metodi poliedrici tradizionali. Sgravando la necessità di una specifica a priori dei vincoli di contenimento degli insiemi ed evitando le bilinearità, l'approccio raggiunge un equilibrio più efficiente tra dominio di attrazione e sforzo computazionale. Il lavoro estende la ricerca preliminare sull'MPC a tubo elastico fornendo un'alternativa basata sul pre-calcolo computazionalmente efficiente e stabilendo condizioni rigorose per la fattibilità ricorsiva e la stabilità. Gli autori notano che, sebbene il metodo sia robusto, il lavoro futuro potrebbe affrontare la sensibilità all'insieme RPI iniziale e potenziali estensioni all'MPC basato sui dati.