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: KYP-Lemma-basierte Niedrigrang-Approximation für H2\mathcal{H}_2-optimale Passivierung

Grundinformationen

  • Papier-ID: 2501.05178
  • Titel: KLAP: KYP lemma based low-rank approximation for H2\mathcal{H}_2-optimal passivation
  • Autoren: Jonas Nicodemus, Matthias Voigt, Serkan Gugercin, Benjamin Unger
  • Klassifikation: math.OC (Mathematische Optimierung und Regelungstheorie)
  • Veröffentlichungsdatum: 14. Oktober 2025 (arXiv-Preprint)
  • Papierlink: https://arxiv.org/abs/2501.05178

Zusammenfassung

In diesem Artikel wird eine neue Passivierungsmethode namens KLAP vorgestellt, die auf dem Kalman-Yakubovich-Popov (KYP)-Lemma und den zugehörigen Lur'e-Gleichungen für die Passivierung linearer zeitinvarianter Systeme basiert. Das Passivierungsproblem im vorgeschlagenen Rahmen entspricht der Suche nach einer Störung eines gegebenen nicht-passiven Systems, um es passiv zu machen, während gleichzeitig der H2\mathcal{H}_2- oder frequenzgewichtete H2\mathcal{H}_2-Abstand zwischen dem ursprünglichen nicht-passiven System und dem resultierenden passiven System minimiert wird. Es wird gezeigt, dass dieses Problem als uneingeschränktes Optimierungsproblem formuliert werden kann, dessen Zielfunktion auch bei großen Systemen effizient differenzierbar ist. Es wird bewiesen, dass jeder Minimierer des uneingeschränkten Problems dasselbe passive System erzeugt, und dass ohne Vorwärtskopplungsterm jeder lokale Minimierer auch ein globaler Minimierer ist.

Forschungshintergrund und Motivation

Bedeutung des Problems

  1. Anforderungen der physikalischen Systemmodellierung: Systeme in vielen physikalischen Bereichen wie Schaltkreisen, Stromnetzen, mechanischen Systemen und poröser Elastizität erfordern Passivitätsgarantien für physikalisch aussagekräftige Simulationsergebnisse
  2. Anforderungen der Netzwerkinterverbindung: Passive Systeme als Bausteine großer Netzwerkmodelle; leistungserhaltende Verbindungen passiver Systeme führen zu einem insgesamt passiven System
  3. Praktische Modellierungsherausforderungen: Obwohl physikalische Prozesse selbst passiv sind, sind Modelle, die durch unstrukturierte Modellreduktionsmethoden oder datengesteuerte Systemidentifikation gewonnen werden, häufig nicht-passiv

Einschränkungen bestehender Methoden

Bestehende Passivierungsmethoden lassen sich hauptsächlich in drei Kategorien einteilen:

  1. LMI-Methoden basierend auf dem KYP-Lemma: Die Rechenkosten wachsen schnell mit der Systemgröße aufgrund der erforderlichen Existenz der Lyapunov-Matrix
  2. Methoden basierend auf spektralen Eigenschaften der Hamilton-Matrix: Mangelnde Konvergenzgarantien, möglicherweise mehrere Iterationen erforderlich
  3. Methoden basierend auf diskreten Frequenzen: Passivität kann nur in bestimmten Frequenzbereichen garantiert werden

Forschungsmotivation

Dieser Artikel zielt darauf ab, eine effiziente Passivierungsmethode zu entwickeln, die:

  • Großskalige Systeme verarbeiten kann
  • Konvergenzgarantien bietet
  • Optimale Lösungen im Sinne der H2\mathcal{H}_2-Norm findet

Kernbeiträge

  1. Explizite Parametrisierung: Unter Ausnutzung der Existenz von Lösungen mit minimaler Rangordnung der KYP-Ungleichung wird eine explizite Parametrisierung eines beliebigen passiven Systems mit nmnm Entscheidungsvariablen erhalten
  2. Umformulierung zu uneingeschränkter Optimierung: Das konvexe Optimierungsproblem mit Nebenbedingungen wird in ein nicht-konvexes uneingeschränktes Optimierungsproblem umformuliert, mit Etablierung von Lösbarkeit, Eindeutigkeit und Gradientenberechnungsmethoden
  3. Theorie der globalen Optimalität: Es wird bewiesen, dass bei schiefsymmetrischen Vorwärtskopplungstermen (D+DT=0D + D^T = 0) jeder lokale Minimierer auch ein globaler Minimierer ist
  4. Erkennung lokaler Optimalität: Es wird ein neues Kriterium bereitgestellt, um mittels extremaler Lösungen der KYP-Ungleichung zu überprüfen, ob ein lokaler Minimierer ein globaler Minimierer ist
  5. Praktische Algorithmusstrategien: Es werden Initialisierungsstrategien basierend auf Vorwärtskopplungsmodifikation und Neustartstrategien vorgeschlagen

Methodische Details

Aufgabendefinition

Gegeben sei ein lineares zeitinvariantes dynamisches System: Σ:{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}

Das Ziel besteht darin, das modifizierte System zu finden: Σ^(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}

so dass Σ^(C^)\hat{\Sigma}(\hat{C}) passiv ist und der H2\mathcal{H}_2-Abstand zum ursprünglichen System minimiert wird.

Theoretische Grundlagen

KYP-Lemma und Passivitätsparametrisierung

Basierend auf dem KYP-Lemma ist ein System passiv genau dann, wenn Matrizen L∈Rn×mL \in \mathbb{R}^{n \times m} und M∈Rm×mM \in \mathbb{R}^{m \times m} existieren, so dass: C=BTL−1(−LLT)+MLTC = B^T\mathcal{L}^{-1}(-LL^T) + ML^TD+DT=MMTD + D^T = MM^T

wobei L\mathcal{L} der Lyapunov-Operator ist: L(X)=ATX+XA\mathcal{L}(X) = A^TX + XA.

Zielfunktion und Gradient

Die Zielfunktion kann ausgedrückt werden als: 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))

wobei PP das Steuerbarkeits-Gramian ist. Der Gradient ist: ∇J(L)=2XL−2P(CT−C^(L)T)M\nabla J(L) = 2XL - 2P(C^T - \hat{C}(L)^T)M

Algorithmusarchitektur

KLAP-Algorithmusablauf

  1. Initialisierung: Erhalten Sie das initiale L0L_0 mit Algorithmus 1
  2. Optimierung: Lösen Sie das uneingeschränkte Problem mit L-BFGS
  3. Globalitätserkennung: Überprüfen Sie die Eigenwerte von Y∗=A−B(D+DT)−1M(L∗)TY^* = A - B(D+D^T)^{-1}M(L^*)^T
  4. Neustartstrategien: Falls lokale Optimalität erkannt wird, führen Sie einen Gradientenschritt durch und starten Sie neu

Initialisierungsstrategie

Störung des Vorwärtskopplungsterms DD zur Passivierung des Systems:

  • Berechnen Sie λmin⁡=min⁡ωλmin⁡(Φ(iω))\lambda_{\min} = \min_\omega \lambda_{\min}(\Phi(i\omega))
  • Setzen Sie Dpert=D−(λmin⁡/2−ϵ)ImD_{\text{pert}} = D - (\lambda_{\min}/2 - \epsilon)I_m
  • Lösen Sie die entsprechende algebraische Riccati-Gleichung zur Initialisierung

Experimentelle Einrichtung

Testsysteme

  1. ACC-Benchmark-Problem: Kleinskaliges System (n=4,m=1n=4, m=1)
  2. CD-Player-Arm: Mittelskaliges System (n=120,m=2n=120, m=2)
  3. Hochgeschwindigkeits-Smartphone-Verbindungsleitung: Großskaliges System (n=800,m=4n=800, m=4)

Vergleichsmethoden

  • LMI: Standard-LMI-Methode basierend auf dem KYP-Lemma
  • LMI-TP: LMI-Methode mit Spurparametrisierung
  • Hamilton-Methode: Methode basierend auf Hamilton-Eigenwertstreuung

Bewertungsmetriken

  • H2\mathcal{H}_2-Fehler: ∥G−G^(⋅;C^)∥H2\|G - \hat{G}(\cdot; \hat{C})\|_{\mathcal{H}_2}
  • Rechenzeit und Iterationszahl
  • Erfolgsquote der Konvergenz zum globalen Optimum

Experimentelle Ergebnisse

Hauptergebnisse

Vergleich der Recheneffizienz

ModellMethodeIterationenGesamtzeit (s)Zeit pro Iteration (s)H2\mathcal{H}_2-Fehler
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⁻¹
CD-PlayerKLAP305,44×10⁻¹1,81×10⁻²1,06×10⁶
CD-PlayerLMI-TP1166,04×10²5,21×10⁰1,00×10⁶
SmartphoneKLAP22081,46×10²6,63×10⁻²8,32×10⁵

Wichtigste Erkenntnisse

  1. Recheneffizienz: KLAP ist 1-2 Größenordnungen schneller als traditionelle LMI-Methoden
  2. Globale Konvergenz: Ohne Vorwärtskopplungsterm sind alle lokalen Optima globale Optima
  3. Wirksamkeit der Neustartstrategien: Neustartstrategien können erfolgreich von nicht-globalen lokalen Optima erholen
  4. Großskalentauglichkeit: Funktioniert auch bei 800-dimensionalen Systemen effektiv

Fallstudien

ACC-Benchmark-Problem

  • Ohne Vorwärtskopplungsterm: Alle Initialisierungen konvergieren zum globalen Optimum
  • Mit Vorwärtskopplungsterm: 40% der zufälligen Initialisierungen konvergieren zu nicht-globalen lokalen Optima
  • Nach Verwendung der Neustartstrategien: Alle Initialisierungen konvergieren zum globalen Optimum

Hochgeschwindigkeits-Smartphone-Verbindungsleitung

  • Verbesserung des H2\mathcal{H}_2-Fehlers um etwa 31% gegenüber Referenzmethoden
  • Durch Diagonalisierungstransformation wird die Lösungszeit für eine einzelne Lyapunov-Gleichung von 550 ms auf 4 ms reduziert

Verwandte Arbeiten

Klassifikation von Passivierungsmethoden

  1. Methoden basierend auf dem KYP-Lemma: Erzeugen konvexe Optimierungsprobleme, aber hohe Rechenkosten
  2. Methoden basierend auf Hamilton-Spektrum: Mangelnde Konvergenzgarantien
  3. Methoden basierend auf Frequenzabtastung: Nur in bestimmten Frequenzbereichen wirksam

Vorteile dieses Artikels

  • Vermeidung großskaliger LMI-Lösung
  • Bereitstellung theoretischer Konvergenzgarantien
  • Anwendbarkeit auf großskalige Systeme
  • Klare Kriterien für globale Optimalität

Schlussfolgerungen und Diskussion

Hauptschlussfolgerungen

  1. Die KLAP-Methode transformiert erfolgreich das Optimierungsproblem mit Nebenbedingungen in ein uneingeschränktes Problem
  2. Garantiert globale Optimalität bei schiefsymmetrischen Vorwärtskopplungstermen
  3. Bietet effektive Mechanismen zur Erkennung lokaler Optimalität und zum Neustart
  4. Zeigt überlegene Recheneffizienz bei mehreren Benchmark-Tests

Einschränkungen

  1. Bei nicht-trivialen Vorwärtskopplungstermen können mehrere lokale Optima existieren
  2. Erfordert asymptotische Stabilitätsannahmen für das System
  3. Konzentriert sich derzeit hauptsächlich auf H2\mathcal{H}_2-Norm-Optimierung

Zukünftige Richtungen

  1. Erweiterung auf das beschränkte reelle Lemma zur Findung des nächsten kontraktiven Systems
  2. Anwendung auf parametrisierte Systeme und differential-algebraische Gleichungen
  3. Untersuchung von H∞\mathcal{H}_\infty-optimalen Passivierungsproblemen

Tiefgreifende Bewertung

Stärken

  1. Solide theoretische Beiträge: Vollständige theoretische Analyse einschließlich Existenz, Eindeutigkeit und globaler Optimalität
  2. Starke Methodische Innovation: Geschickte Nutzung der Niedrigrang-Zerlegung des KYP-Lemmas, Vermeidung von Rechnenengpässen traditioneller Methoden
  3. Hervorragende Praktikabilität: Algorithmus ist leicht zu implementieren und für großskalige Systeme geeignet
  4. Umfassende Experimente: Validierung der Methode auf Benchmark-Systemen verschiedener Größenordnungen

Schwächen

  1. Problem lokaler Optima: Bei allgemeinen Vorwärtskopplungstermen können lokale Optima auftreten
  2. Abhängigkeit von Initialisierung: Die Methodenleistung hängt in gewissem Maße von der Initialisierungsqualität ab
  3. Unvollständige theoretische Analyse: Die Analyse für den Fall D+DT⊁0D + D^T \not\succ 0 ist nicht ausreichend vollständig

Auswirkungen

  1. Akademischer Wert: Bietet neue theoretische Perspektiven und Lösungsmethoden für Passivierungsprobleme
  2. Praktischer Wert: Besonders geeignet für Passivierungsbedarf großskaliger Ingenieurssysteme
  3. Reproduzierbarkeit: Code und Daten sind öffentlich verfügbar, was Verifikation und Anwendung erleichtert

Anwendungsszenarien

  1. Passivierung großskaliger linearer Systeme
  2. Wiederherstellung der Passivität nach Modellreduktion
  3. Nachbearbeitung nach datengestützter Systemidentifikation
  4. Entwurf von Netzwerkverbindungssystemen

Literaturverzeichnis

Das Papier zitiert 58 relevante Literaturquellen, die hauptsächlich folgende Bereiche abdecken:

  • Theoretische Grundlagen dissipativer Systeme Willems, 1972
  • KYP-Lemma und positive Realität Anderson & Vongpanitlerd, 1973
  • Übersicht über Passivierungsmethoden Grivet-Talocia & Gustavsen, 2016
  • Numerische Optimierungsmethoden Boyd et al., 1989