2025-11-10T02:47:53.458764

KLAP: KYP lemma based low-rank approximation for $\mathcal{H}_2$-optimal passivation

Nicodemus, Voigt, Gugercin et al.
We present a novel passivity enforcement (passivation) method, called KLAP, for linear time-invariant systems based on the Kalman-Yakubovich-Popov (KYP) lemma and the closely related Lur'e equations. The passivation problem in our framework corresponds to finding a perturbation to a given non-passive system that renders the system passive while minimizing the $\mathcal{H}_2$ or frequency-weighted $\mathcal{H}_2$ distance between the original non-passive and the resulting passive system. We show that this problem can be formulated as an unconstrained optimization problem whose objective function can be differentiated efficiently even in large-scale settings. We show that any minimizer of the unconstrained problem yields the same passive system. Furthermore, we prove that, in the absence of a feedthrough term, every local minimizer is also a global minimizer. For cases involving a non-trivial feedthrough term, we analyze global minimizers in relation to the extremal solutions of the Lur'e equations, which can serve as tools for identifying local minima. To solve the resulting numerical optimization problem efficiently, we propose an initialization strategy based on modifying the feedthrough term and a restart strategy when it is likely that the optimization has converged to a non-global local minimum. Numerical examples illustrate the effectiveness of the proposed method.
academic

KLAP: Approssimazione a basso rango basata sul lemma KYP per la passivazione ottimale in H2\mathcal{H}_2

Informazioni Fondamentali

  • ID Articolo: 2501.05178
  • Titolo: KLAP: KYP lemma based low-rank approximation for H2\mathcal{H}_2-optimal passivation
  • Autori: Jonas Nicodemus, Matthias Voigt, Serkan Gugercin, Benjamin Unger
  • Classificazione: math.OC (Ottimizzazione e Controllo Matematico)
  • Data di Pubblicazione: 14 ottobre 2025 (preprint arXiv)
  • Link Articolo: https://arxiv.org/abs/2501.05178

Riassunto

Il presente articolo propone un nuovo metodo di imposizione della passività denominato KLAP, destinato alla passivazione di sistemi lineari tempo-invarianti basato sul lemma di Kalman-Yakubovich-Popov (KYP) e sulle relative equazioni di Lur'e. Nel contesto di questo framework, il problema della passivazione corrisponde alla ricerca di una perturbazione applicata a un sistema non passivo assegnato, al fine di renderlo passivo, minimizzando contemporaneamente la distanza in norma H2\mathcal{H}_2 o H2\mathcal{H}_2 ponderata in frequenza tra il sistema originale non passivo e il sistema passivo risultante. La ricerca dimostra che il problema può essere formulato come un problema di ottimizzazione non vincolata, la cui funzione obiettivo è efficacemente differenziabile anche in contesti di grandi dimensioni. È provato che qualsiasi minimizzatore del problema non vincolato produce il medesimo sistema passivo, e in assenza di termini di feedforward, ogni minimizzatore locale è anche un minimizzatore globale.

Contesto di Ricerca e Motivazione

Importanza del Problema

  1. Necessità di Modellazione di Sistemi Fisici: Numerosi campi fisici quali sistemi circuitali, sistemi di potenza, sistemi meccanici ed elasticità porosa richiedono garanzie di passività per ottenere risultati di simulazione fisicamente significativi
  2. Requisiti di Interconnessione in Rete: I sistemi passivi fungono da blocchi costitutivi di modelli di rete di grandi dimensioni, e l'interconnessione conservativa di sistemi passivi produce un sistema globalmente passivo
  3. Sfide di Modellazione Pratica: Sebbene i processi fisici sottostanti siano intrinsecamente passivi, i modelli ottenuti mediante metodi di riduzione dell'ordine non strutturati o tecniche di identificazione di sistemi guidate dai dati risultano frequentemente non passivi

Limitazioni dei Metodi Esistenti

I metodi di passivazione esistenti si dividono principalmente in tre categorie:

  1. Metodi LMI basati sul Lemma KYP: Il costo computazionale cresce rapidamente con le dimensioni del sistema, a causa della necessità dell'esistenza della matrice di Lyapunov
  2. Metodi basati sulle Caratteristiche Spettrali della Matrice Hamiltoniana: Mancano garanzie di convergenza e potrebbero richiedere iterazioni multiple
  3. Metodi basati su Frequenze Discrete: Garantiscono la passività solo in intervalli di frequenza specifici

Motivazione della Ricerca

Il presente lavoro mira allo sviluppo di un metodo di passivazione efficiente, capace di:

  • Gestire sistemi di grandi dimensioni
  • Fornire garanzie di convergenza
  • Trovare soluzioni ottimali nel senso della norma H2\mathcal{H}_2

Contributi Fondamentali

  1. Parametrizzazione Esplicita: Sfruttando l'esistenza della soluzione a rango minimo della disuguaglianza KYP, si ottiene una parametrizzazione esplicita di qualsiasi sistema passivo con nmnm variabili decisionali
  2. Riformulazione dell'Ottimizzazione Non Vincolata: Il problema di ottimizzazione convessa vincolata viene riformulato come problema di ottimizzazione non convessa non vincolata, stabilendo metodi di risolubilità, unicità e calcolo del gradiente
  3. Teoria dell'Ottimalità Globale: È provato che nel caso di termini di feedforward antisimmetrici (D+DT=0D + D^T = 0), qualsiasi minimizzatore locale è anche un minimizzatore globale
  4. Rilevamento dell'Ottimalità Locale: Viene fornito un nuovo criterio per verificare, utilizzando la soluzione estremale della disuguaglianza KYP, se un minimizzatore locale è un minimizzatore globale
  5. Strategie Algoritmiche Pratiche: Sono proposte strategie di inizializzazione basate sulla modifica del termine di feedforward e strategie di riavvio

Dettagli Metodologici

Definizione del Compito

Dato un sistema dinamico lineare tempo-invariante: Σ:{x˙(t)=Ax(t)+Bu(t)y(t)=Cx(t)+Du(t)\Sigma : \begin{cases} \dot{x}(t) = Ax(t) + Bu(t) \\ y(t) = Cx(t) + Du(t) \end{cases}

L'obiettivo è trovare il sistema modificato: Σ^(C^):{x˙(t)=Ax(t)+Bu(t)y(t)=C^x(t)+Du(t)\hat{\Sigma}(\hat{C}) : \begin{cases} \dot{x}(t) = Ax(t) + Bu(t) \\ y(t) = \hat{C}x(t) + Du(t) \end{cases}

tale che Σ^(C^)\hat{\Sigma}(\hat{C}) sia passivo e minimizzi la distanza in norma H2\mathcal{H}_2 dal sistema originale.

Fondamenti Teorici Principali

Lemma KYP e Parametrizzazione della Passività

Basandosi sul lemma KYP, un sistema è passivo se e solo se esistono matrici L∈Rn×mL \in \mathbb{R}^{n \times m} e M∈Rm×mM \in \mathbb{R}^{m \times m} tali che: C=BTL−1(−LLT)+MLTC = B^T\mathcal{L}^{-1}(-LL^T) + ML^TD+DT=MMTD + D^T = MM^T

dove L\mathcal{L} è l'operatore di Lyapunov: L(X)=ATX+XA\mathcal{L}(X) = A^TX + XA.

Funzione Obiettivo e Gradiente

La funzione obiettivo può essere espressa come: J(L)=tr((C−C^(L))P(CT−C^(L)T))J(L) = \text{tr}((C - \hat{C}(L))P(C^T - \hat{C}(L)^T))

dove PP è il Gramiano di controllabilità. Il gradiente è: ∇J(L)=2XL−2P(CT−C^(L)T)M\nabla J(L) = 2XL - 2P(C^T - \hat{C}(L)^T)M

Architettura dell'Algoritmo

Flusso dell'Algoritmo KLAP

  1. Inizializzazione: Ottenimento dell'iniziale L0L_0 utilizzando l'Algoritmo 1
  2. Ottimizzazione: Risoluzione del problema non vincolato mediante L-BFGS
  3. Rilevamento della Globalità: Verifica degli autovalori di Y∗=A−B(D+DT)−1M(L∗)TY^* = A - B(D+D^T)^{-1}M(L^*)^T
  4. Strategia di Riavvio: Se viene rilevato un ottimo locale, esecuzione di un passo di gradiente e riavvio

Strategia di Inizializzazione

Perturbazione del termine di feedforward DD per rendere il sistema passivo:

  • Calcolo di λmin⁡=min⁡ωλmin⁡(Φ(iω))\lambda_{\min} = \min_\omega \lambda_{\min}(\Phi(i\omega))
  • Impostazione di Dpert=D−(λmin⁡/2−ϵ)ImD_{\text{pert}} = D - (\lambda_{\min}/2 - \epsilon)I_m
  • Risoluzione dell'equazione di Riccati algebrica corrispondente per l'inizializzazione

Configurazione Sperimentale

Sistemi di Prova

  1. Problema di Riferimento ACC: Sistema di piccole dimensioni (n=4,m=1n=4, m=1)
  2. Braccio di Lettore CD: Sistema di medie dimensioni (n=120,m=2n=120, m=2)
  3. Collegamento Interconnessione Smartphone ad Alta Velocità: Sistema di grandi dimensioni (n=800,m=4n=800, m=4)

Metodi di Confronto

  • LMI: Metodo LMI standard basato sul lemma KYP
  • LMI-TP: Metodo LMI con parametrizzazione della traccia
  • Metodo Hamiltoniano: Metodo basato sulla perturbazione degli autovalori Hamiltoniani

Indicatori di Valutazione

  • Errore in norma H2\mathcal{H}_2: \|G - \hat{G}(\cdot; \hat{C})\|_{\mathcal{H}_2
  • Tempo di calcolo e numero di iterazioni
  • Tasso di successo della convergenza all'ottimo globale

Risultati Sperimentali

Risultati Principali

Confronto dell'Efficienza Computazionale

ModelloMetodoIterazioniTempo Totale(s)Tempo per Iterazione(s)Errore H2\mathcal{H}_2
ACCKLAP122.29×10⁻⁴1.91×10⁻⁵8.71×10⁻¹
ACCLMI134.61×10⁻³3.54×10⁻⁴8.71×10⁻¹
ACCLMI-TP113.59×10⁻²3.26×10⁻³8.71×10⁻¹
Lettore CDKLAP305.44×10⁻¹1.81×10⁻²1.06×10⁶
Lettore CDLMI-TP1166.04×10²5.21×10⁰1.00×10⁶
SmartphoneKLAP22081.46×10²6.63×10⁻²8.32×10⁵

Risultati Chiave

  1. Efficienza Computazionale: KLAP è da 1 a 2 ordini di grandezza più veloce dei metodi LMI tradizionali
  2. Convergenza Globale: In assenza di termini di feedforward, tutti gli ottimi locali sono ottimi globali
  3. Efficacia della Strategia di Riavvio: La strategia di riavvio recupera con successo da ottimi locali non globali
  4. Applicabilità a Grandi Dimensioni: Rimane efficace anche su sistemi di dimensione 800

Analisi di Casi Studio

Problema di Riferimento ACC

  • Senza termine di feedforward: tutte le inizializzazioni convergono all'ottimo globale
  • Con termine di feedforward: il 40% delle inizializzazioni casuali converge a un ottimo locale non globale
  • Dopo l'utilizzo della strategia di riavvio: tutte le inizializzazioni convergono all'ottimo globale

Collegamento Interconnessione Smartphone ad Alta Velocità

  • Miglioramento dell'errore in norma H2\mathcal{H}_2 di circa il 31% rispetto al metodo di riferimento
  • Mediante trasformazione diagonale, il tempo di risoluzione di una singola equazione di Lyapunov è ridotto da 550ms a 4ms

Lavori Correlati

Classificazione dei Metodi di Passivazione

  1. Metodi basati sul Lemma KYP: Producono problemi di ottimizzazione convessa ma con elevato costo computazionale
  2. Metodi basati sullo Spettro Hamiltoniano: Mancano garanzie di convergenza
  3. Metodi basati su Campionamento in Frequenza: Efficaci solo in intervalli di frequenza specifici

Vantaggi del Presente Lavoro

  • Evita la risoluzione di LMI su larga scala
  • Fornisce garanzie teoriche di convergenza
  • Applicabile a sistemi di grandi dimensioni
  • Possiede criteri espliciti per l'ottimalità globale

Conclusioni e Discussione

Conclusioni Principali

  1. Il metodo KLAP trasforma con successo il problema di ottimizzazione vincolata in un problema non vincolato
  2. Garantisce l'ottimalità globale nel caso di termini di feedforward antisimmetrici
  3. Fornisce meccanismi efficaci di rilevamento dell'ottimalità locale e di riavvio
  4. Dimostra efficienza computazionale superiore su molteplici test di riferimento

Limitazioni

  1. Per termini di feedforward non banali, possono esistere molteplici ottimi locali
  2. Richiede che il sistema soddisfi l'ipotesi di stabilità asintotica
  3. Attualmente focalizzato principalmente sull'ottimizzazione in norma H2\mathcal{H}_2

Direzioni Future

  1. Estensione al lemma della realtà limitata per trovare il sistema contrattile più vicino
  2. Applicazione a sistemi parametrizzati e equazioni differenziali algebriche
  3. Ricerca del problema di passivazione ottimale in norma H∞\mathcal{H}_\infty

Valutazione Approfondita

Punti di Forza

  1. Contributi Teorici Solidi: Fornisce un'analisi teorica completa, inclusa l'esistenza, l'unicità e l'ottimalità globale
  2. Forte Innovazione Metodologica: Sfrutta abilmente la decomposizione a basso rango del lemma KYP, evitando i colli di bottiglia computazionali dei metodi tradizionali
  3. Praticità Eccezionale: L'algoritmo è facile da implementare e applicabile a sistemi di grandi dimensioni
  4. Sperimentazione Completa: La validazione su sistemi di riferimento di diverse dimensioni dimostra l'efficacia del metodo

Carenze

  1. Problema dell'Ottimalità Locale: Per termini di feedforward generali, il metodo può ancora rimanere intrappolato in ottimi locali
  2. Dipendenza dall'Inizializzazione: Le prestazioni del metodo dipendono in una certa misura dalla qualità dell'inizializzazione
  3. Limitazioni dell'Analisi Teorica: L'analisi per il caso D+DT⊁0D + D^T \not\succ 0 non è sufficientemente completa

Impatto

  1. Valore Accademico: Fornisce una nuova prospettiva teorica e un metodo di risoluzione per il problema della passivazione
  2. Valore Pratico: Particolarmente adatto alle esigenze di passivazione di sistemi ingegneristici di grandi dimensioni
  3. Riproducibilità: Il codice e i dati sono stati resi pubblici, facilitando la verifica e l'applicazione

Scenari di Applicazione

  1. Passivazione di sistemi lineari di grandi dimensioni
  2. Recupero della passività dopo la riduzione dell'ordine del modello
  3. Post-elaborazione dopo l'identificazione di sistemi guidata dai dati
  4. Progettazione di sistemi interconnessi in rete

Bibliografia

L'articolo cita 58 riferimenti correlati, principalmente riguardanti:

  • Fondamenti della teoria dei sistemi dissipativi Willems, 1972
  • Lemma KYP e teoria della realtà positiva Anderson & Vongpanitlerd, 1973
  • Rassegna dei metodi di passivazione Grivet-Talocia & Gustavsen, 2016
  • Metodi di ottimizzazione numerica Boyd et al., 1989