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 보조정리 기반 $\mathcal{H}_2$ -최적 수동성 강제를 위한 저차 근사

기본 정보

논문 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년 1월 14일 (arXiv 사전인쇄본)
논문 링크: https://arxiv.org/abs/2501.05178

초록

본 논문은 Kalman-Yakubovich-Popov (KYP) 보조정리 및 관련 Lur'e 방정식을 기반으로 선형 시불변 시스템의 수동성 강제를 위한 KLAP이라는 새로운 방법을 제안한다. 본 프레임워크의 수동성 강제 문제는 주어진 비수동 시스템에 대한 섭동을 찾아 시스템을 수동으로 만들면서 동시에 원래의 비수동 시스템과 결과 수동 시스템 간의 $\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⁵