Symplektische Modellordnungsreduktion von Port-Hamiltonschen Systemen

Grundinformationen

• Paper-ID: 2203.07751
• Titel: Symplektische Modellordnungsreduktion von Port-Hamiltonschen Systemen
• Autoren: Silke Glas, Mir Mamunuzzaman, Hongliang Mu, Hans Zwart
• Klassifizierung: math.OC (Optimierung und Steuerung), math.CA (Klassische Analysis)
• Veröffentlichungsdatum: März 2022 (arXiv-Preprint)
• Paper-Link: https://arxiv.org/abs/2203.07751

Zusammenfassung

Dieses Paper präsentiert eine neuartige strukturerhaltende Modellordnungsreduktionsmethode (MOR) für linear zeitinvariante Port-Hamiltonsche (pH) Systeme. Das Forschungsziel besteht darin, ein reduziertes pH-System zu konstruieren, das weiterhin in der physikalischen Domäne des ursprünglichen vollständigen Modells interpretierbar bleibt. Konkret: Wenn das ursprüngliche hochdimensionale pH-System eine Schaltung ist, soll das reduzierte Modell weiterhin als elektronische Schaltung interpretierbar sein. Für Masse-Feder-Dämpfer-Systeme (MSD) existieren bereits MOR-Methoden, die diese spezifische Strukturerhaltung garantieren. Dieses Paper zeigt, dass die neue strukturerhaltende MOR-Methode basierend auf symplektischen MOR-Verfahren im Fall von MSD-Systemen die bekannte Second-Order-Arnoldi-Methode (SO-Arnoldi) wiederherstellt. Für Schaltungs-pH-Modelle und andere Modelle mit ähnlicher Blockstruktur erzeugt die Methode jedoch neuartige Modellordnungsreduktionstechniken.

Forschungshintergrund und Motivation

Problembeschreibung

Port-Hamiltonsche Systeme sind eine Verallgemeinerung des klassischen Hamiltonschen Rahmens und äußerst nützlich in Ingenieurproblemen, da sie die klassische Hamiltonsche Methode mit dem leistungsstarken Rahmen der Modellierung offener physikalischer Systeme verbinden. Bei der Reglerentwurf hat der Regler jedoch typischerweise die gleiche Ordnung wie das ursprüngliche System, was die Notwendigkeit nach reduzierten Modellen motiviert, die das grundlegende Verhalten des ursprünglichen Systems erfassen.

Einschränkungen bestehender Methoden

1. Strukturverlustproblem: Die direkte Anwendung traditioneller MOR-Methoden führt zum Verlust der pH-Struktur
2. Fehlende physikalische Interpretierbarkeit: Bestehende Methoden können nicht garantieren, dass reduzierte Modelle in der physikalischen Domäne interpretierbar sind
3. Zerstörung der Blockstruktur: Für pH-Systeme mit spezieller Blockstruktur können bestehende Methoden diese Struktur nicht bewahren

Forschungsmotivation

Wenn physikalische Modelle in pH-Form umformuliert werden, weisen die Matrizen J, R und H typischerweise spezielle Blockstrukturen auf. Das Ziel dieses Papers ist es, diese Blockstruktur während des Reduktionsprozesses zu bewahren und sicherzustellen, dass das reduzierte Modell weiterhin als physikalisches System desselben Typs wie das vollständige Modell interpretierbar ist.

Kernbeiträge

1. Einführung der Symplektisch-pH-Methode: Präsentation einer neuartigen pH-System-MOR-Methode basierend auf symplektischen MOR-Verfahren, anwendbar auf Blockstrukturform (2), sodass das reduzierte Modell die gleiche Form eines pH-Systems beibehält
2. Theoretische Verbindung: Beweis, dass die vorgeschlagene Methode im Spezialfall von MSD-Systemen die Second-Order-Arnoldi-Methode (SO-Arnoldi) wiederherstellt
3. Numerische Verifikation: Demonstration numerischer Ergebnisse an elektronischen Schaltungen (LRCR-Schaltung) mit Vergleich zu anderen bekannten Methoden aus der Literatur
4. Stabilitätsvorteil: Feststellung, dass die SO-Arnoldi-Methode in bestimmten Fällen instabile reduzierte Modelle erzeugt, während die vorliegende Methode Stabilität bewahrt

Methodische Details

Aufgabendefinition

Gegeben ist ein linear zeitinvariantes pH-System:

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

wobei x(t) ∈ ℂ²ⁿ der Zustand, u(t) ∈ ℂᵖ die Eingabe und y(t) ∈ ℂᵖ die Ausgabe ist. Das Ziel besteht darin, ein System der gleichen Form mit kleinerer Zustandsraumdimension zu finden, sodass die beiden Übertragungsfunktionen an vordefinierten Frequenzpunkten gleich sind.

Theoretischer Kernrahmen

Symplektische Geometrie-Grundlagen

Definition 1 (Symplektische Form): Sei V ein geraddimensionaler Vektorraum über ℂ. Eine symplektische Form Ω: V × V → ℂ erfüllt:

1. Antisymmetrie: Ω(v₁, v₂) = -Ω(v₂, v₁)
2. Nicht-Degeneriertheit: Für alle v₁ ∈ V gilt Ω(v₁, v₂) = 0 genau dann, wenn v₂ = 0

Definition 2 (Symplektische Abbildung): Eine Matrix Q heißt (J₂ₙ, J₂ₖ)-symplektisch, wenn:

Q*J₂ₙQ = J₂ₖ

Symplektische MOR für den verlustfreien Fall

Theorem 4: Für verlustfreie pH-Systeme (R = 0), wenn Q ∈ ℂ²ⁿˣ²ᵏ eine (J₂ₙ⁻¹, J₂ₖ⁻¹)-symplektische Matrix ist und erfüllt:

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

dann stimmt das reduzierte pH-System mit dem ursprünglichen System an den gegebenen Frequenzen überein.

Symplektische MOR für den verlustbehafteten Fall

Theorem 5: Für verlustbehaftete pH-Systeme wird die erweiterte Matrix konstruiert:

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

wobei Q₂ (J₂ₙ, J₂ₖ)-symplektisch ist und Q₁* die linke Inverse von Q₂ ist. Unter angemessenen Bedingungen erfüllt das reduzierte System die Interpolationsbedingungen.

Blockstruktur-erhaltender Algorithmus

Algorithmus 1 (Symplektisch-pH MOR):

1. Konstruktion von V₁ basierend auf Interpolationsdaten
2. Berechnung von V₂ := HV₁
3. Aufteilung von V₁ und V₂ nach Blockstruktur
4. Konstruktion von Ψ₁, Ψ₂, Φ₁, Φ₂ durch Formel (27)
5. Aufbau von Q₁ und Q₂
6. Anwendung der Petrov-Galerkin-Projektion

Konstruktionsformeln:

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

Technische Innovationspunkte

1. Symplektische Strukturerhaltung: Bewahrung der fundamentalen geometrischen Struktur Hamiltonscher Systeme durch symplektische Abbildungen
2. Erhaltung der Blockdiagonalstruktur: Spezielle Konstruktion von Q₁ und Q₂ gewährleistet die Bewahrung der physikalisch bedeutsamen Blockstruktur
3. Interpolationsbedingungen: Garantiert, dass das reduzierte System an vordefinierten Frequenzpunkten exakt mit dem ursprünglichen System übereinstimmt
4. Einheitlicher Rahmen: Enthält die SO-Arnoldi-Methode für MSD-Systeme als Spezialfall

Experimentelle Einrichtung

Numerisches Beispiel: LRCR-Schaltung

Betrachtet wird die in Abbildung 1 gezeigte LRCR-Schaltung mit n=50 Zellen und folgenden Parametern:

• Kapazitäten: Cᵢ = 1 pF
• Induktivitäten: Lᵢ = 1 nH
• Widerstände: Rᵢ,₁ = 1 mΩ, Rᵢ,₂ = 1 kΩ

Bewertungsmetriken

• H₂-Fehler:

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

• H∞-Fehler:

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

Vergleichsmethoden

1. pH-Arnoldi: Strukturerhaltende Methode basierend auf Theorem 2
2. SO-Arnoldi: Second-Order-Arnoldi-Methode
3. Symplektisch-pH: Die in diesem Paper vorgeschlagene Methode

Wahl der Interpolationspunkte

Gleichmäßig verteilte Interpolationspunkte auf logarithmischer Skala im Intervall 10⁻², 10², beispielsweise 6 Interpolationspunkte: 10⁻²i, -10⁻²i, 10⁰i, -10⁰i, 10²i, -10²i

Experimentelle Ergebnisse

Hauptergebnisse

Aus Abbildung 3 können folgende Beobachtungen gemacht werden:

1. Symplektisch-pH vs pH-Arnoldi: Symplektisch-pH liefert in den meisten Fällen niedrigere H₂- und H∞-Fehler
2. Symplektisch-pH vs SO-Arnoldi:
   • Ähnliche Leistung bei Reduktionsdimension < 24
   • Drastischer Leistungsabfall der SO-Arnoldi bei Reduktionsdimension = 32
   • SO-Arnoldi erzeugt instabile Systeme bei Reduktionsdimensionen 28 und 36

Stabilitätsanalyse

Abbildung 4 zeigt die von der SO-Arnoldi-Methode bei r=28 und r=36 erzeugten instabilen Pole, während die Symplektisch-pH-Methode bei allen getesteten Dimensionen Stabilität bewahrt.

Theoretische Verifikation: MSD-Systemäquivalenz

Theorem 7 beweist, dass die Symplektisch-pH-Methode für MSD-Systeme durch eine unitäre Transformation äquivalent zur SO-Arnoldi-Methode ist: Es existiert eine unitäre Matrix Θ, sodass zred(t) = Θẑred(t)

Verwandte Arbeiten

Entwicklung der pH-System-MOR

1. Frühe Arbeiten: Methoden für passive Systeme basierend auf positiv reellen Übertragungsfunktionen
2. Momentenabgleich: Arbeiten von Polyuga und van der Schaft
3. Dirac-Struktur-Approximation: Alternative Wege zur Strukturerhaltung
4. Spektralzerlegungsmethoden: Methoden basierend auf Popov-Funktionen
5. Optimierungsmethoden: Methoden basierend auf Parameteroptimierung und Riemannscher Optimierung

Entwicklung der Symplektischen MOR

1. Autonome Hamiltonsche Systeme: Bahnbrechende Arbeiten von Peng und Mohseni
2. Nichtlineare Methoden: Auf Mannigfaltigkeiten basierende symplektische Reduktion
3. Datengesteuerte Methoden: Operatorinferenz und Systemidentifikation

Schlussfolgerungen und Diskussion

Hauptschlussfolgerungen

1. Erfolgreiche Entwicklung einer MOR-Methode, die pH-Struktur und Blockstruktur bewahrt
2. Theoretische Vereinigung verschiedener MOR-Methoden für MSD-Systeme
3. Numerische Experimente verifizieren die Effektivität und Stabilitätsvorteile der Methode

Einschränkungen

1. Interpolationsmethode: Derzeit nur Abgleich von Übertragungsfunktionswerten, keine Ableitungsabgleichung
2. Blockstruktur-Beschränkung: Methode speziell für bestimmte Blockstrukturen konzipiert
3. Rechenkomplexität: Rechenaufwand für die Konstruktion symplektischer Matrizen bedarf weiterer Analyse

Zukünftige Richtungen

1. Ableitungsabgleichung: Erweiterung auf Abgleichung von Ableitungen der Übertragungsfunktion zur Verbesserung der Anpassungsqualität
2. Adaptive Auswahl: Entwicklung adaptiver Auswahlstrategien für Interpolationspunkte
3. Nichtlineare Erweiterung: Erweiterung auf nichtlineare pH-Systeme
4. Großskalige Anwendungen: Optimierung des Algorithmus für größere Systeme

Tiefgreifende Bewertung

Stärken

1. Theoretische Strenge: Solide mathematische Grundlagen basierend auf symplektischer Geometrie
2. Strukturerhaltung: Gleichzeitige Bewahrung der pH-Struktur und physikalisch bedeutsamen Blockstruktur
3. Einheitlichkeit: Vereinigung bestehender Methoden als Spezialfälle in einem Rahmen
4. Praktikabilität: Methode mit klarer algorithmischer Implementierung und numerischer Verifikation

Mängel

1. Anwendungsbereich: Methode speziell für bestimmte Blockstrukturen konzipiert, begrenzte Allgemeingültigkeit
2. Recheneffizienz: Rechenkomplexitätsanalyse der Konstruktion symplektischer Matrizen nicht ausreichend
3. Großskalige Verifikation: Mangel an numerischer Verifikation für größere Systeme

Auswirkungen

1. Theoretischer Beitrag: Neue theoretische Perspektive für pH-System-MOR
2. Praktischer Wert: Direkter praktischer Wert in Ingenieuranwendungen wie Schaltungsanalyse
3. Methodologische Bedeutung: Demonstriert die Wichtigkeit geometrischer Methoden in der MOR

Anwendungsszenarien

1. Schaltungsanalyse: Ordnungsreduktion großer elektronischer Schaltungen
2. Mechanische Systeme: Reglerentwurf für Masse-Feder-Dämpfer-Systeme
3. Multiphysik-Systeme: Komplexe physikalische Systeme mit Energieerhaltungseigenschaften

Literaturverzeichnis

Das Paper enthält 31 Literaturquellen, die wichtige Arbeiten aus mehreren Bereichen abdecken, darunter MOR-Theorie, pH-Systemmodellierung, symplektische Geometrie und Regelungstheorie, und bieten eine solide theoretische Grundlage für die Forschung.