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 : Approximation de faible rang basée sur le lemme KYP pour la passivation optimale en H2\mathcal{H}_2

Informations fondamentales

  • Identifiant de l'article : 2501.05178
  • Titre : KLAP: KYP lemma based low-rank approximation for H2\mathcal{H}_2-optimal passivation
  • Auteurs : Jonas Nicodemus, Matthias Voigt, Serkan Gugercin, Benjamin Unger
  • Classification : math.OC (Optimisation et contrôle mathématiques)
  • Date de publication : 14 octobre 2025 (prépublication arXiv)
  • Lien de l'article : https://arxiv.org/abs/2501.05178

Résumé

Cet article propose une nouvelle méthode de passivation appelée KLAP, basée sur le lemme de Kalman-Yakubovich-Popov (KYP) et les équations de Lur'e associées pour la passivation de systèmes linéaires invariants dans le temps. Dans ce cadre, le problème de passivation correspond à la recherche d'une perturbation d'un système non-passif donné pour le rendre passif, tout en minimisant la distance H2\mathcal{H}_2 ou H2\mathcal{H}_2 pondérée en fréquence entre le système non-passif original et le système passif résultant. L'étude montre que ce problème peut être formulé comme un problème d'optimisation sans contrainte dont la fonction objectif est efficacement différentiable même dans les contextes de grande dimension. Il est démontré que tout minimiseur du problème sans contrainte produit le même système passif, et qu'en l'absence de terme de précommande, chaque minimiseur local est également un minimiseur global.

Contexte et motivation de la recherche

Importance du problème

  1. Besoins de modélisation des systèmes physiques : De nombreux domaines physiques tels que les circuits électriques, les systèmes électriques, les systèmes mécaniques et l'élasticité poreuse nécessitent des garanties de passivité pour obtenir des résultats de simulation physiquement significatifs
  2. Exigences d'interconnexion en réseau : Les systèmes passifs servant de blocs de construction pour les grands modèles de réseau, l'interconnexion préservant la puissance de systèmes passifs produit un système globalement passif
  3. Défis de modélisation pratique : Bien que les processus physiques eux-mêmes soient passifs, les modèles obtenus par des méthodes de réduction de modèle non structurées ou par identification de systèmes pilotée par les données sont souvent non-passifs

Limitations des méthodes existantes

Les méthodes de passivation existantes se divisent principalement en trois catégories :

  1. Méthodes LMI basées sur le lemme KYP : Le coût de calcul augmente rapidement avec la taille du système en raison de la nécessité d'existence de la matrice de Lyapunov
  2. Méthodes basées sur les caractéristiques spectrales de la matrice hamiltonienne : Manquent de garanties de convergence et peuvent nécessiter plusieurs itérations
  3. Méthodes basées sur l'échantillonnage discret en fréquence : Ne peuvent garantir la passivité que sur des plages de fréquence spécifiques

Motivation de la recherche

Cet article vise à développer une méthode de passivation efficace capable de :

  • Traiter les systèmes de grande dimension
  • Fournir des garanties de convergence
  • Trouver la solution optimale au sens de la norme H2\mathcal{H}_2

Contributions principales

  1. Paramétrisation explicite : En exploitant l'existence de solutions de rang minimal de l'inégalité KYP, une paramétrisation explicite de tout système passif avec nmnm variables de décision est obtenue
  2. Reformulation en optimisation sans contrainte : Le problème d'optimisation convexe avec contraintes est reformulé en problème d'optimisation non-convexe sans contrainte, établissant la résolubilité, l'unicité et les méthodes de calcul du gradient
  3. Théorie de l'optimalité globale : Il est prouvé que dans le cas d'un terme de précommande antisymétrique (D+DT=0D + D^T = 0), tout minimiseur local est également un minimiseur global
  4. Détection des minimiseurs locaux : Un nouveau critère est fourni pour vérifier si un minimiseur local est un minimiseur global en utilisant les solutions extrémales de l'inégalité KYP
  5. Stratégies algorithmiques pratiques : Des stratégies d'initialisation et de redémarrage basées sur la modification du terme de précommande sont proposées

Détails de la méthode

Définition de la tâche

Étant donné un système dynamique linéaire invariant dans le temps : Σ:{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'objectif est de trouver le système modifié : Σ^(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}

tel que Σ^(C^)\hat{\Sigma}(\hat{C}) soit passif et minimise la distance H2\mathcal{H}_2 avec le système original.

Fondements théoriques principaux

Lemme KYP et paramétrisation de la passivité

Selon le lemme KYP, un système est passif si et seulement s'il existe des matrices L∈Rn×mL \in \mathbb{R}^{n \times m} et M∈Rm×mM \in \mathbb{R}^{m \times m} telles que : C=BTL−1(−LLT)+MLTC = B^T\mathcal{L}^{-1}(-LL^T) + ML^TD+DT=MMTD + D^T = MM^T

où L\mathcal{L} est l'opérateur de Lyapunov : L(X)=ATX+XA\mathcal{L}(X) = A^TX + XA.

Fonction objectif et gradient

La fonction objectif peut être exprimée comme : 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))

où PP est la grammienne de commandabilité. Le gradient est : ∇J(L)=2XL−2P(CT−C^(L)T)M\nabla J(L) = 2XL - 2P(C^T - \hat{C}(L)^T)M

Architecture de l'algorithme

Flux de l'algorithme KLAP

  1. Initialisation : Obtenir L0L_0 initial en utilisant l'algorithme 1
  2. Optimisation : Résoudre le problème sans contrainte en utilisant L-BFGS
  3. Détection de la globalité : Vérifier les valeurs propres de Y∗=A−B(D+DT)−1M(L∗)TY^* = A - B(D+D^T)^{-1}M(L^*)^T
  4. Stratégie de redémarrage : Si un optimum local est détecté, exécuter une étape de gradient et redémarrer

Stratégie d'initialisation

Rendre le système passif en perturbant le terme de précommande DD :

  • Calculer λmin⁡=min⁡ωλmin⁡(Φ(iω))\lambda_{\min} = \min_\omega \lambda_{\min}(\Phi(i\omega))
  • Définir Dpert=D−(λmin⁡/2−ϵ)ImD_{\text{pert}} = D - (\lambda_{\min}/2 - \epsilon)I_m
  • Résoudre l'équation algébrique de Riccati correspondante pour l'initialisation

Configuration expérimentale

Systèmes de test

  1. Problème de référence ACC : Système de petite dimension (n=4,m=1n=4, m=1)
  2. Bras de lecteur CD : Système de dimension moyenne (n=120,m=2n=120, m=2)
  3. Liaison d'interconnexion de smartphone haute vitesse : Système de grande dimension (n=800,m=4n=800, m=4)

Méthodes de comparaison

  • LMI : Méthode LMI standard basée sur le lemme KYP
  • LMI-TP : Méthode LMI avec paramétrisation de trace
  • Méthode hamiltonienne : Méthode basée sur la perturbation des valeurs propres hamiltoniennes

Indicateurs d'évaluation

  • Erreur H2\mathcal{H}_2 : \|G - \hat{G}(\cdot; \hat{C})\|_{\mathcal{H}_2
  • Temps de calcul et nombre d'itérations
  • Taux de succès de convergence vers l'optimum global

Résultats expérimentaux

Résultats principaux

Comparaison de l'efficacité de calcul

ModèleMéthodeItérationsTemps total (s)Temps par itération (s)Erreur 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⁻¹
Lecteur CDKLAP305,44×10⁻¹1,81×10⁻²1,06×10⁶
Lecteur CDLMI-TP1166,04×10²5,21×10⁰1,00×10⁶
SmartphoneKLAP22081,46×10²6,63×10⁻²8,32×10⁵

Résultats clés

  1. Efficacité de calcul : KLAP est 1 à 2 ordres de grandeur plus rapide que les méthodes LMI traditionnelles
  2. Convergence globale : En l'absence de terme de précommande, tous les optima locaux sont des optima globaux
  3. Efficacité de la stratégie de redémarrage : La stratégie de redémarrage peut récupérer avec succès les optima locaux non-globaux
  4. Applicabilité à grande dimension : Fonctionne efficacement sur des systèmes de dimension 800

Analyses de cas

Problème de référence ACC

  • Sans terme de précommande : Toutes les initialisations convergent vers l'optimum global
  • Avec terme de précommande : 40 % des initialisations aléatoires convergent vers un optimum local non-global
  • Après utilisation de la stratégie de redémarrage : Toutes les initialisations convergent vers l'optimum global

Liaison d'interconnexion de smartphone haute vitesse

  • Amélioration de l'erreur H2\mathcal{H}_2 d'environ 31 % par rapport à la méthode de référence
  • Grâce à la transformation diagonale, le temps de résolution d'une seule équation de Lyapunov est réduit de 550 ms à 4 ms

Travaux connexes

Classification des méthodes de passivation

  1. Méthodes basées sur le lemme KYP : Produisent des problèmes d'optimisation convexe mais avec un coût de calcul élevé
  2. Méthodes basées sur le spectre hamiltonien : Manquent de garanties de convergence
  3. Méthodes basées sur l'échantillonnage en fréquence : Efficaces uniquement sur des plages de fréquence spécifiques

Avantages de cet article

  • Évite la résolution de LMI de grande dimension
  • Fournit des garanties théoriques de convergence
  • Applicable aux systèmes de grande dimension
  • Dispose de critères clairs d'optimalité globale

Conclusions et discussion

Conclusions principales

  1. La méthode KLAP transforme avec succès un problème d'optimisation avec contraintes en problème sans contrainte
  2. Garantit l'optimalité globale dans le cas d'un terme de précommande antisymétrique
  3. Fournit des mécanismes efficaces de détection des optima locaux et de redémarrage
  4. Démontre une efficacité de calcul supérieure sur plusieurs tests de référence

Limitations

  1. Pour les termes de précommande non triviaux, plusieurs optima locaux peuvent exister
  2. Nécessite que le système satisfasse l'hypothèse de stabilité asymptotique
  3. Actuellement orienté principalement vers l'optimisation en norme H2\mathcal{H}_2

Directions futures

  1. Extension au lemme de réalité bornée pour trouver le système contractif le plus proche
  2. Application aux systèmes paramétrés et équations différentielles-algébriques
  3. Étude du problème de passivation optimale en H∞\mathcal{H}_\infty

Évaluation approfondie

Points forts

  1. Contributions théoriques solides : Fournit une analyse théorique complète, incluant l'existence, l'unicité et l'optimalité globale
  2. Forte innovativité de la méthode : Exploite intelligemment la décomposition de faible rang du lemme KYP, évitant les goulots d'étranglement de calcul des méthodes traditionnelles
  3. Praticité remarquable : L'algorithme est facile à implémenter et applicable aux systèmes de grande dimension
  4. Expériences suffisantes : La méthode est validée sur des systèmes de référence de différentes dimensions

Insuffisances

  1. Problème des optima locaux : Pour les termes de précommande généraux, la méthode peut toujours converger vers des optima locaux
  2. Dépendance à l'initialisation : Les performances de la méthode dépendent dans une certaine mesure de la qualité de l'initialisation
  3. Limitations de l'analyse théorique : L'analyse pour le cas D+DT⊁0D + D^T \not\succ 0 n'est pas suffisamment complète

Impact

  1. Valeur académique : Fournit une nouvelle perspective théorique et une méthode de résolution pour le problème de passivation
  2. Valeur pratique : Particulièrement adaptée aux besoins de passivation des systèmes d'ingénierie de grande dimension
  3. Reproductibilité : Le code et les données sont publiquement disponibles, facilitant la vérification et l'application

Scénarios d'application

  1. Passivation de systèmes linéaires de grande dimension
  2. Récupération de la passivité après réduction de modèle
  3. Post-traitement après identification de système pilotée par les données
  4. Conception de systèmes d'interconnexion en réseau

Références

L'article cite 58 références connexes, couvrant principalement :

  • Fondements de la théorie des systèmes dissipatifs Willems, 1972
  • Lemme KYP et théorie de la réalité positive Anderson & Vongpanitlerd, 1973
  • Synthèse des méthodes de passivation Grivet-Talocia & Gustavsen, 2016
  • Méthodes d'optimisation numérique Boyd et al., 1989