KLAP: KYP Lemma Based Low-Rank Approximation for $\mathcal{H}_2$-Optimal Passivation

Basic Information

• Paper ID: 2501.05178
• Title: KLAP: KYP lemma based low-rank approximation for $\mathcal{H}_2$-optimal passivation
• Authors: Jonas Nicodemus, Matthias Voigt, Serkan Gugercin, Benjamin Unger
• Classification: math.OC (Mathematical Optimization and Control)
• Publication Date: October 14, 2025 (arXiv preprint)
• Paper Link: https://arxiv.org/abs/2501.05178

Abstract

This paper proposes a novel passivation enforcement method called KLAP for linear time-invariant systems based on the Kalman-Yakubovich-Popov (KYP) lemma and associated Lur'e equations. Within this framework, the passivation problem corresponds to finding a perturbation to a given non-passive system that renders it passive while minimizing the $\mathcal{H}_2$ or frequency-weighted $\mathcal{H}_2$ distance between the original non-passive system and the resulting passive system. The study demonstrates that this problem can be formulated as an unconstrained optimization problem whose objective function can be efficiently differentiated even in large-scale settings. It is proven that any minimizer of the unconstrained problem yields the same passive system, and in the absence of feedthrough terms, every local minimizer is also a global minimizer.

Research Background and Motivation

Problem Significance

1. Physical System Modeling Requirements: Systems in numerous physical domains including circuits, power systems, mechanical systems, and poroelasticity require passivity guarantees to obtain physically meaningful simulation results
2. Network Interconnection Requirements: Passive systems serve as building blocks for large-scale network models, and power-preserving interconnections of passive systems yield overall passive systems
3. Practical Modeling Challenges: Although physical processes are inherently passive, models obtained through unstructured model reduction methods or data-driven system identification techniques are often non-passive

Limitations of Existing Methods

Existing passivation methods fall into three main categories:

1. LMI Methods Based on KYP Lemma: Computational cost grows rapidly with system size due to the requirement for Lyapunov matrix existence
2. Methods Based on Hamiltonian Matrix Spectral Characteristics: Lack convergence guarantees and may require multiple iterations
3. Discrete Frequency-Based Methods: Can only guarantee passivity within specific frequency ranges

Research Motivation

This work aims to develop an efficient passivation method capable of:

• Handling large-scale systems
• Providing convergence guarantees
• Finding optimal solutions in the $\mathcal{H}_2$ norm sense

Core Contributions

1. Explicit Parametrization: Leveraging the existence of rank-minimizing solutions to KYP inequalities, an explicit parametrization of any passive system with $nm$ decision variables is obtained
2. Unconstrained Optimization Reformulation: Reformulates the convex constrained optimization problem as a non-convex unconstrained optimization problem, establishing solvability, uniqueness, and gradient computation methods
3. Global Optimality Theory: Proves that for skew-symmetric feedthrough terms ($D + D^T = 0$), any local minimizer is also a global minimizer
4. Local Optimality Detection: Provides new criteria using KYP inequality extremal solutions to check whether a local minimizer is a global minimizer
5. Practical Algorithm Strategies: Proposes initialization and restart strategies based on feedthrough term modification

Methodology Details

Problem Formulation

Given a linear time-invariant dynamic system: $\Sigma : \begin{cases} \dot{x}(t) = Ax(t) + Bu(t) \\ y(t) = Cx(t) + Du(t) \end{cases}$

The objective is to find a modified system: $\hat{\Sigma}(\hat{C}) : \begin{cases} \dot{x}(t) = Ax(t) + Bu(t) \\ y(t) = \hat{C}x(t) + Du(t) \end{cases}$

such that $\hat{\Sigma}(\hat{C})$ is passive and the $\mathcal{H}_2$ distance to the original system is minimized.

Core Theoretical Foundation

KYP Lemma and Passivity Parametrization

Based on the KYP lemma, a system is passive if and only if there exist matrices $L \in \mathbb{R}^{n \times m}$ and $M \in \mathbb{R}^{m \times m}$ such that: $C = B^T\mathcal{L}^{-1}(-LL^T) + ML^T$$D + D^T = MM^T$

where $\mathcal{L}$ is the Lyapunov operator: $\mathcal{L}(X) = A^TX + XA$.

Objective Function and Gradient

The objective function can be expressed as: $J(L) = \text{tr}((C - \hat{C}(L))P(C^T - \hat{C}(L)^T))$

where $P$ is the controllability Gramian. The gradient is: $\nabla J(L) = 2XL - 2P(C^T - \hat{C}(L)^T)M$

Algorithm Architecture

KLAP Algorithm Flow

1. Initialization: Obtain initial $L_0$ using Algorithm 1
2. Optimization: Solve the unconstrained problem using L-BFGS
3. Global Optimality Detection: Check eigenvalues of $Y^* = A - B(D+D^T)^{-1}M(L^*)^T$
4. Restart Strategy: If local optimality is detected, execute a gradient step and restart

Initialization Strategy

Render the system passive by perturbing the feedthrough term $D$:

• Compute $\lambda_{\min} = \min_\omega \lambda_{\min}(\Phi(i\omega))$
• Set $D_{\text{pert}} = D - (\lambda_{\min}/2 - \epsilon)I_m$
• Solve the corresponding algebraic Riccati equation to obtain initialization

Experimental Setup

Test Systems

1. ACC Benchmark Problem: Small-scale system ($n=4, m=1$)
2. CD Player Arm: Medium-scale system ($n=120, m=2$)
3. High-Speed Smartphone Interconnect Link: Large-scale system ($n=800, m=4$)

Comparison Methods

• LMI: Standard LMI method based on KYP lemma
• LMI-TP: LMI method with trace parametrization
• Hamiltonian Method: Method based on Hamiltonian eigenvalue perturbation

Evaluation Metrics

• $\mathcal{H}_2$ error: $\|G - \hat{G}(\cdot; \hat{C})\|_{\mathcal{H}_2}$
• Computation time and iteration count
• Success rate of convergence to global optimum

Experimental Results

Main Results

Computational Efficiency Comparison

Key Findings

1. Computational Efficiency: KLAP is 1-2 orders of magnitude faster than traditional LMI methods
2. Global Convergence: In the absence of feedthrough terms, all local optima are global optima
3. Restart Strategy Effectiveness: The restart strategy successfully recovers from non-global local optima
4. Large-Scale Applicability: Remains effective on systems of dimension 800

Case Studies

ACC Benchmark Problem

• Without feedthrough terms: All initializations converge to global optimum
• With feedthrough terms: 40% of random initializations converge to non-global local optima
• Using restart strategy: All initializations converge to global optimum

High-Speed Smartphone Interconnect Link

• Compared to reference methods, $\mathcal{H}_2$ error improvement of approximately 31%
• Through diagonalization transformation, single Lyapunov equation solution time reduced from 550ms to 4ms

Related Work

Classification of Passivation Methods

1. KYP Lemma-Based Methods: Produce convex optimization problems but with high computational cost
2. Hamiltonian Spectrum-Based Methods: Lack convergence guarantees
3. Frequency Sampling-Based Methods: Only effective within specific frequency ranges

Advantages of This Work

• Avoids large-scale LMI solving
• Provides theoretical convergence guarantees
• Applicable to large-scale systems
• Offers explicit global optimality criteria

Conclusions and Discussion

Main Conclusions

1. The KLAP method successfully transforms the constrained optimization problem into an unconstrained one
2. Guarantees global optimality for skew-symmetric feedthrough terms
3. Provides effective local optimality detection and restart mechanisms
4. Demonstrates superior computational efficiency on multiple benchmark tests

Limitations

1. For non-trivial feedthrough terms, multiple local optima may exist
2. Requires the system to satisfy asymptotic stability assumptions
3. Currently focused primarily on $\mathcal{H}_2$ norm optimization

Future Directions

1. Extension to bounded real lemma for finding nearest contractive systems
2. Application to parametric systems and differential-algebraic equations
3. Investigation of $\mathcal{H}_\infty$-optimal passivation problems

In-Depth Evaluation

Strengths

1. Solid Theoretical Contributions: Provides comprehensive theoretical analysis including existence, uniqueness, and global optimality
2. Strong Method Innovation: Cleverly exploits low-rank decomposition of KYP inequalities, avoiding computational bottlenecks of traditional methods
3. Outstanding Practicality: Algorithm is easy to implement and applicable to large-scale systems
4. Comprehensive Experiments: Validates method effectiveness on benchmark systems of varying scales

Weaknesses

1. Local Optimality Issue: For general feedthrough terms, may still be trapped in local optima
2. Initialization Dependency: Method performance depends to some extent on initialization quality
3. Incomplete Theoretical Analysis: Analysis for the case $D + D^T \not\succ 0$ is incomplete

Impact

1. Academic Value: Provides new theoretical perspective and solution method for passivation problems
2. Practical Value: Particularly suitable for large-scale engineering systems requiring passivation
3. Reproducibility: Code and data are publicly available, facilitating verification and application

Applicable Scenarios

1. Passivation of large-scale linear systems
2. Recovery of passivity after model reduction
3. Post-processing of data-driven system identification
4. Design of network interconnected systems

References

The paper cites 58 related references, primarily covering:

• Dissipative systems theory foundations Willems, 1972
• KYP lemma and positive realness theory Anderson & Vongpanitlerd, 1973
• Passivation methods survey Grivet-Talocia & Gustavsen, 2016
• Numerical optimization methods Boyd et al., 1989