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 based low-rank approximation for $\mathcal{H}_2$ -optimal passivation

基本信息

论文ID: 2501.05178
标题: KLAP: KYP lemma based low-rank approximation for $\mathcal{H}_2$ -optimal passivation
作者: Jonas Nicodemus, Matthias Voigt, Serkan Gugercin, Benjamin Unger
分类: math.OC (数学优化与控制)
发表时间: 2025年10月14日 (arXiv预印本)
论文链接: https://arxiv.org/abs/2501.05178

摘要

本文提出了一种名为KLAP的新型被动性强制方法，用于基于Kalman-Yakubovich-Popov (KYP)引理和相关Lur'e方程的线性时不变系统被动化。该框架中的被动化问题对应于寻找对给定非被动系统的扰动，使系统变为被动，同时最小化原始非被动系统与结果被动系统之间的 $\mathcal{H}_2$ 或频率加权 $\mathcal{H}_2$ 距离。研究表明，该问题可以表述为无约束优化问题，其目标函数即使在大规模设置下也能高效微分。证明了无约束问题的任何最小化器都产生相同的被动系统，并且在没有前馈项的情况下，每个局部最小化器也是全局最小化器。

研究背景与动机

问题的重要性

物理系统建模需求: 电路系统、电力系统、机械系统和多孔弹性等众多物理领域的系统需要被动性保证，以获得物理上有意义的仿真结果
网络互连要求: 被动系统作为大型网络模型的构建块，保功率互连的被动系统能产生整体被动系统
实际建模挑战: 尽管物理过程本身是被动的，但通过非结构化模型降阶方法或数据驱动系统辨识技术得到的模型往往是非被动的

现有方法的局限性

现有被动化方法主要分为三类：

基于KYP引理的LMI方法: 计算成本随系统规模增长很快，由于需要Lyapunov矩阵的存在
基于Hamiltonian矩阵谱特征的方法: 缺乏收敛保证，可能需要多次迭代
基于离散频率的方法: 只能在特定频率范围内保证被动性

研究动机

本文旨在开发一种高效的被动化方法，能够：

处理大规模系统
提供收敛保证
在 $\mathcal{H}_2$ 范数意义下找到最优解

核心贡献

显式参数化: 利用KYP不等式的秩最小化解的存在性，获得了具有 $nm$ 个决策变量的任何被动系统的显式参数化
无约束优化重构: 将凸约束优化问题重构为非凸无约束优化问题，建立了可解性、唯一性和梯度计算方法
全局最优性理论: 证明了在斜对称前馈项情况下( $D + D^T = 0$ )，任何局部最小化器也是全局最小化器
局部最优检测: 提供了使用KYP不等式极值解检查局部最小化器是否为全局最小化器的新准则
实用算法策略: 提出了基于前馈项修改的初始化策略和重启策略

方法详解

任务定义

给定线性时不变动态系统： $\Sigma : \begin{cases} \dot{x}(t) = Ax(t) + Bu(t) \\ y(t) = Cx(t) + Du(t) \end{cases}$

目标是找到修改后的系统： $\hat{\Sigma}(\hat{C}) : \begin{cases} \dot{x}(t) = Ax(t) + Bu(t) \\ y(t) = \hat{C}x(t) + Du(t) \end{cases}$

使得 $\hat{\Sigma}(\hat{C})$ 是被动的，并最小化与原系统的 $\mathcal{H}_2$ 距离。

初始化: 使用算法1获得初始 $L_0$
优化: 使用L-BFGS求解无约束问题
全局性检测: 检查 $Y^* = A - B(D+D^T)^{-1}M(L^*)^T$ 的特征值
重启策略: 如果检测到局部最优，执行梯度步并重启

初始化策略

通过扰动前馈项 $D$ 使系统变为被动：

计算 $\lambda_{\min} = \min_\omega \lambda_{\min}(\Phi(i\omega))$
设置 $D_{\text{pert}} = D - (\lambda_{\min}/2 - \epsilon)I_m$
求解相应的代数Riccati方程获得初始化

实验设置

测试系统

ACC基准问题: 小规模系统( $n=4, m=1$ )
CD播放器摆臂: 中等规模系统( $n=120, m=2$ )
高速智能手机互连链路: 大规模系统( $n=800, m=4$ )

对比方法

LMI: 标准基于KYP引理的LMI方法
LMI-TP: 带迹参数化的LMI方法
Hamiltonian方法: 基于Hamiltonian特征值扰动的方法

评价指标

$\mathcal{H}_2$ 误差： $\|G - \hat{G}(\cdot; \hat{C})\|_{\mathcal{H}_2}$
计算时间和迭代次数
收敛到全局最优的成功率

模型	方法	迭代次数	总时间(s)	每次迭代时间(s)	$\mathcal{H}_2$ 误差
ACC	KLAP	12	2.29×10⁻⁴	1.91×10⁻⁵	8.71×10⁻¹
ACC	LMI	13	4.61×10⁻³	3.54×10⁻⁴	8.71×10⁻¹
ACC	LMI-TP	11	3.59×10⁻²	3.26×10⁻³	8.71×10⁻¹
CD播放器	KLAP	30	5.44×10⁻¹	1.81×10⁻²	1.06×10⁶
CD播放器	LMI-TP	116	6.04×10²	5.21×10⁰	1.00×10⁶
智能手机	KLAP	2208	1.46×10²	6.63×10⁻²	8.32×10⁵