Riduzione dell'ordine del modello simplettico di sistemi port-hamiltoniani

Informazioni Fondamentali

• ID Articolo: 2203.07751
• Titolo: Symplectic model order reduction of port-Hamiltonian systems
• Autori: Silke Glas, Mir Mamunuzzaman, Hongliang Mu, Hans Zwart
• Classificazione: math.OC (Ottimizzazione e Controllo), math.CA (Analisi Classica)
• Data di Pubblicazione: Marzo 2022 (preprint arXiv)
• Link Articolo: https://arxiv.org/abs/2203.07751

Riassunto

Questo articolo propone un innovativo metodo di riduzione dell'ordine del modello (MOR) che preserva la struttura per sistemi port-hamiltoniani (pH) lineari e invarianti nel tempo. L'obiettivo della ricerca è costruire un sistema pH ridotto che rimanga interpretabile nel dominio fisico del modello originale a ordine completo. In particolare, se il sistema pH ad alta dimensionalità iniziale è un circuito, desideriamo che il modello ridotto rimanga interpretabile come circuito elettronico. Per i sistemi massa-molla-smorzamento (MSD), esistono già metodi MOR che garantiscono la preservazione di questa struttura specifica. L'articolo dimostra che il nuovo metodo MOR che preserva la struttura, basato sul metodo MOR simplettico, recupera il noto metodo Arnoldi del secondo ordine nel caso dei sistemi MSD. Tuttavia, per modelli pH con circuiti e strutture a blocchi simili, il metodo produce innovative tecniche di riduzione dell'ordine del modello.

Contesto e Motivazione della Ricerca

Descrizione del Problema

I sistemi port-hamiltoniani rappresentano una generalizzazione del framework hamiltoniano classico, estremamente utile nei problemi ingegneristici poiché combinano i metodi hamiltoniani classici con il potente framework per la modellazione di sistemi fisici aperti. Tuttavia, nella progettazione dei controllori, il controllore tipicamente possiede lo stesso ordine del sistema originale, il che motiva la necessità di modelli ridotti che catturino il comportamento fondamentale del sistema originale.

Limitazioni dei Metodi Esistenti

1. Perdita di Struttura: L'applicazione diretta di metodi MOR tradizionali comporta la perdita della struttura pH
2. Mancanza di Interpretabilità Fisica: I metodi esistenti non garantiscono l'interpretabilità del modello ridotto nel dominio fisico
3. Distruzione della Struttura a Blocchi: Per sistemi pH con strutture a blocchi speciali, i metodi esistenti non preservano questa struttura

Motivazione della Ricerca

Quando si riformulano modelli fisici nella forma pH, le matrici J, R e H tipicamente possiedono strutture a blocchi speciali. L'obiettivo di questo articolo è preservare questa struttura a blocchi durante il processo di riduzione, garantendo che il modello ridotto rimanga interpretabile come un sistema fisico dello stesso tipo del modello a ordine completo.

Contributi Principali

1. Introduzione del Metodo Simplettico-pH: Propone un innovativo metodo MOR per sistemi pH basato su MOR simplettico, applicabile alla forma con struttura a blocchi (2), in modo che il modello ridotto preservi la forma di sistema pH identica
2. Connessione Teorica: Dimostra che nel caso speciale dei sistemi MSD, il metodo proposto recupera il metodo Arnoldi del secondo ordine (SO-Arnoldi)
3. Verifica Numerica: Presenta risultati numerici su circuiti elettronici (circuiti LRCR), con confronti con altri metodi noti in letteratura
4. Vantaggi di Stabilità: Scopre che il metodo SO-Arnoldi produce in alcuni casi modelli ridotti instabili, mentre il metodo proposto preserva la stabilità

Dettagli del Metodo

Definizione del Compito

Dato un sistema pH lineare e invariante nel tempo:

ẋ(t) = (J - R)Hx(t) + Bu(t)
y(t) = B*Hx(t)

dove x(t) ∈ C²ⁿ è lo stato, u(t) ∈ Cᵖ è l'ingresso, y(t) ∈ Cᵖ è l'uscita. L'obiettivo è trovare un sistema della stessa forma ma con dimensione dello spazio degli stati più piccola, tale che le due funzioni di trasferimento siano uguali in punti di frequenza predeterminati.

Framework Teorico Fondamentale

Fondamenti di Geometria Simplettica

Definizione 1 (Forma Simplettica): Sia V uno spazio vettoriale di dimensione pari su C, una forma simplettica Ω: V × V → C soddisfa:

1. Antisimmetria: Ω(v₁, v₂) = -Ω(v₂, v₁)
2. Non-degeneratezza: Per tutti v₁ ∈ V, Ω(v₁, v₂) = 0 se e solo se v₂ = 0

Definizione 2 (Applicazione Simplettica): Una matrice Q è detta (J₂ₙ, J₂ₖ)-simplettica se:

Q*J₂ₙQ = J₂ₖ

MOR Simplettico per il Caso Senza Dissipazione

Teorema 4: Per sistemi pH senza dissipazione (R = 0), se Q ∈ C²ⁿˣ²ᵏ è una matrice (J₂ₙ⁻¹, J₂ₖ⁻¹)-simplettica e soddisfa:

• ran(B) ⊆ ran(Q)
• (sₘI₂ₙ - J₂ₙH)⁻¹Bμₘ ∈ ran(Q)

allora il sistema pH ridotto è uguale al sistema originale alle frequenze date.

MOR Simplettico per il Caso Con Dissipazione

Teorema 5: Per sistemi pH con dissipazione, si costruisce la matrice estesa:

QE := [Q₁  0  ]
      [0   Q₂]

dove Q₂ è (J₂ₙ, J₂ₖ)-simplettica, Q₁* è l'inversa sinistra di Q₂. Sotto condizioni appropriate, il sistema ridotto soddisfa le condizioni di interpolazione.

Algoritmo di Preservazione della Struttura a Blocchi

Algoritmo 1 (MOR Simplettico-pH):

1. Costruire V₁ in base ai dati di interpolazione
2. Calcolare V₂ := HV₁
3. Dividere V₁ e V₂ secondo la struttura a blocchi
4. Costruire Ψ₁, Ψ₂, Φ₁, Φ₂ mediante la formula (27)
5. Costruire Q₁ e Q₂
6. Applicare la proiezione Petrov-Galerkin

Formule di costruzione:

Ψ₂ = orth(V₂₂)
Ψ₁ = V₁₂(Ψ₂*V₁₂)⁻¹
Φ₂ = V₂₁(Ψ₂*Ĵₙ*V₂₁)⁻¹Ĵₖ*
Φ₁ = V₁₁(Φ₂*V₁₁)⁻¹

Punti di Innovazione Tecnica

1. Preservazione della Struttura Simplettica: Mediante applicazioni simplettiche si preserva la struttura geometrica fondamentale del sistema hamiltoniano
2. Mantenimento della Struttura Diagonale a Blocchi: La costruzione speciale di Q₁ e Q₂ garantisce che la struttura a blocchi di significato fisico sia preservata
3. Condizioni di Interpolazione: Garantisce che il sistema ridotto corrisponda esattamente al sistema originale in punti di frequenza predeterminati
4. Framework Unificato: Include il metodo SO-Arnoldi per sistemi MSD come caso particolare

Configurazione Sperimentale

Esempio Numerico: Circuito LRCR

Si consideri il circuito LRCR mostrato in Figura 1, contenente n=50 celle, con parametri impostati come:

• Capacità: Cᵢ = 1 pF
• Induttanza: Lᵢ = 1 nH
• Resistenza: Rᵢ,₁ = 1 mΩ, Rᵢ,₂ = 1 kΩ

Metriche di Valutazione

• Errore H₂:

  H₂error = √(1/2π ∫ Trace((G(iω) - Gred(iω))*(G(iω) - Gred(iω)))dω)
  

• Errore H∞:

  H∞error = max‖G(iω) - Gred(iω)‖₂
  

Metodi di Confronto

1. pH-Arnoldi: Metodo che preserva la struttura basato sul Teorema 2
2. SO-Arnoldi: Metodo Arnoldi del secondo ordine
3. Simplettico-pH: Metodo proposto in questo articolo

Selezione dei Punti di Interpolazione

Punti di interpolazione uniformemente distribuiti su scala logaritmica nell'intervallo 10⁻², 10², ad esempio 6 punti di interpolazione: 10⁻²i, -10⁻²i, 10⁰i, -10⁰i, 10²i, -10²i

Risultati Sperimentali

Risultati Principali

Dalla Figura 3 si osserva:

1. Simplettico-pH vs pH-Arnoldi: Simplettico-pH fornisce errori H₂ e H∞ inferiori nella maggior parte dei casi
2. Simplettico-pH vs SO-Arnoldi:
   • Prestazioni simili quando la dimensione ridotta <24
   • Prestazioni di SO-Arnoldi degradate drasticamente quando la dimensione ridotta =32
   • SO-Arnoldi produce sistemi instabili quando la dimensione ridotta è 28 e 36

Analisi di Stabilità

La Figura 4 mostra i poli instabili prodotti dal metodo SO-Arnoldi per r=28 e r=36, mentre il metodo Simplettico-pH mantiene la stabilità per tutte le dimensioni testate.

Verifica Teorica: Equivalenza per Sistemi MSD

Il Teorema 7 dimostra che per sistemi MSD, il metodo Simplettico-pH è equivalente al metodo SO-Arnoldi attraverso una trasformazione unitaria: esiste una matrice unitaria Θ tale che zred(t) = Θẑred(t)

Lavori Correlati

Sviluppo di MOR per Sistemi pH

1. Lavori Iniziali: Metodi per sistemi passivi basati su funzioni di trasferimento reali positive
2. Matching dei Momenti: Lavori di Polyuga e van der Schaft
3. Approssimazione della Struttura di Dirac: Approcci alternativi per la preservazione della struttura
4. Metodi di Decomposizione Spettrale: Metodi basati sulla funzione di Popov
5. Metodi di Ottimizzazione: Metodi basati su ottimizzazione parametrica e ottimizzazione riemanniana

Sviluppo di MOR Simplettico

1. Sistemi Hamiltoniani Autonomi: Lavori pioneristici di Peng e Mohseni
2. Metodi Non-lineari: Riduzione simplettica basata su varietà
3. Metodi Guidati dai Dati: Inferenza di operatori e identificazione di sistemi

Conclusioni e Discussione

Conclusioni Principali

1. Sviluppo riuscito di un metodo MOR che preserva la struttura pH e la struttura a blocchi
2. Unificazione teorica di diversi metodi MOR per sistemi MSD
3. Verifica sperimentale numerica dell'efficacia del metodo e dei vantaggi di stabilità

Limitazioni

1. Metodo di Interpolazione: Attualmente corrisponde solo ai valori della funzione di trasferimento, senza considerare il matching delle derivate
2. Limitazioni della Struttura a Blocchi: Il metodo è specificamente progettato per strutture a blocchi particolari
3. Complessità Computazionale: La complessità computazionale della costruzione di matrici simplettiche richiede ulteriore analisi

Direzioni Future

1. Matching delle Derivate: Estensione al matching delle derivate della funzione di trasferimento per migliorare la qualità dell'approssimazione
2. Selezione Adattiva: Sviluppo di strategie di selezione adattiva per i punti di interpolazione
3. Estensione Non-lineare: Estensione a sistemi pH non-lineari
4. Applicazioni su Larga Scala: Ottimizzazione dell'algoritmo per gestire sistemi di dimensioni maggiori

Valutazione Approfondita

Punti di Forza

1. Rigore Teorico: Fondamenti matematici solidi basati sulla geometria simplettica
2. Preservazione della Struttura: Preservazione simultanea della struttura pH e della struttura a blocchi di significato fisico
3. Unificazione: Unificazione dei metodi esistenti come casi particolari in un unico framework
4. Praticità: Il metodo possiede un'implementazione algoritmica esplicita e verifica numerica

Insufficienze

1. Ambito di Applicabilità: Il metodo è specificamente progettato per strutture a blocchi particolari, con generalità limitata
2. Efficienza Computazionale: L'analisi della complessità computazionale della costruzione di matrici simplettiche non è sufficientemente approfondita
3. Verifica su Larga Scala: Mancanza di verifica numerica su sistemi di dimensioni maggiori

Impatto

1. Contributo Teorico: Fornisce una nuova prospettiva teorica per MOR di sistemi pH
2. Valore Pratico: Possiede valore diretto in applicazioni ingegneristiche come l'analisi di circuiti
3. Significato Metodologico: Dimostra l'importanza dei metodi geometrici in MOR

Scenari Applicabili

1. Analisi di Circuiti: Modellazione ridotta di circuiti elettronici su larga scala
2. Sistemi Meccanici: Progettazione di controllo per sistemi massa-molla-smorzamento
3. Sistemi Multi-fisici: Sistemi fisici complessi con caratteristiche di conservazione dell'energia

Bibliografia

L'articolo contiene 31 riferimenti bibliografici che coprono molteplici campi inclusa la teoria MOR, la modellazione di sistemi pH, la geometria simplettica e la teoria del controllo, fornendo una base teorica solida per la ricerca.