2025-11-12T23:07:10.120463

Symplectic model order reduction of port-Hamiltonian systems

Glas, Mamunuzzaman, Mu et al.

This work proposes a novel structure-preserving model order reduction (MOR) method for linear, time-invariant port-Hamiltonian (pH) systems. Our goal is to construct a reduced order pH system, which can still be interpreted in the physical domain of the full order model. By this we mean, that if an electrical circuit is the initial high-dimensional pH system, we want the reduced order model to be still interpretable as an electronic circuit. In the case of the well-known mass spring damper (MSD) system, there are MOR methods available, which already guarantee the preservation of this particular structure. Moreover, we show that our new structure-preserving MOR method, which is based on symplectic MOR methods, will recover the known second-order Arnoldi method in the case of MSD systems. However, for the example of an electrical circuit pH model (and more models of similar block structure), our method yields a novel model reduction method. We present numerical results on the aforementioned electronic circuit model, highlighting the advantages of the proposed method.

academic

포트-해밀턴 시스템의 심플렉틱 모델 차수 축소

기본 정보

논문 ID: 2203.07751
제목: Symplectic model order reduction of port-Hamiltonian systems
저자: Silke Glas, Mir Mamunuzzaman, Hongliang Mu, Hans Zwart
분류: math.OC (최적화 및 제어), math.CA (고전 해석)
발표 시간: 2022년 3월 (arXiv 사전인쇄본)
논문 링크: https://arxiv.org/abs/2203.07751

초록

본 논문은 선형 시불변 포트-해밀턴(pH) 시스템을 위한 새로운 구조 보존 모델 차수 축소(MOR) 방법을 제안한다. 연구 목표는 축소된 pH 시스템이 원래 전체 차수 모델의 물리적 영역에서 여전히 해석 가능하도록 하는 것이다. 구체적으로, 초기 고차원 pH 시스템이 전자 회로인 경우, 축소된 모델도 전자 회로로 해석 가능해야 한다. 질량-스프링-감쇠(MSD) 시스템의 경우 이러한 특정 구조를 보존하는 MOR 방법이 이미 존재한다. 본 논문은 심플렉틱 MOR 방법 기반의 새로운 구조 보존 MOR 방법이 MSD 시스템의 경우 알려진 2차 Arnoldi 방법을 복원할 수 있음을 보여준다. 그러나 전자 회로 pH 모델과 같은 유사한 블록 구조를 가진 모델의 경우, 이 방법은 새로운 모델 차수 축소 기법을 생성한다.

연구 배경 및 동기

문제 설명

포트-해밀턴 시스템은 고전 해밀턴 프레임워크의 일반화로, 고전 해밀턴 방법을 개방 물리 시스템 모델링의 강력한 프레임워크와 결합하기 때문에 공학 문제에서 매우 유용하다. 그러나 제어 설계에서 제어기는 보통 원래 시스템과 동일한 차수를 가지므로, 원래 시스템의 기본 동작을 포착할 수 있는 축소된 모델에 대한 필요성이 대두된다.

기존 방법의 한계

구조 손실 문제: 전통적인 MOR 방법을 직접 적용하면 pH 구조가 손실됨
물리적 해석성 부재: 기존 방법은 축소된 모델의 물리적 영역에서의 해석 가능성을 보장할 수 없음
블록 구조 파괴: 특수한 블록 구조를 가진 pH 시스템의 경우, 기존 방법은 이러한 구조를 유지할 수 없음

연구 동기

물리 모델을 pH 형식으로 재구성할 때, 행렬 J, R 및 H는 보통 특수한 블록 구조를 가진다. 본 논문의 목표는 차수 축소 과정에서 이러한 블록 구조를 유지하여 축소된 모델이 여전히 전체 차수 모델과 동일한 유형의 물리 시스템으로 해석 가능하도록 하는 것이다.

핵심 기여

심플렉틱-pH 방법 도입: 심플렉틱 MOR 기반의 새로운 pH 시스템 MOR 방법을 제안하여 블록 구조 형식(2)에 적용 가능하게 하고, 축소된 모델이 동일한 형식의 pH 시스템을 유지하도록 함
이론적 연결: MSD 시스템의 특수한 경우에서 제안된 방법이 2차 Arnoldi 방법(SO-Arnoldi)을 복원함을 증명
수치 검증: 전자 회로(LRCR 회로)에서 수치 결과를 제시하고 문헌의 다른 알려진 방법과 비교
안정성 이점: SO-Arnoldi 방법이 특정 경우에 불안정한 축소 모델을 생성하는 반면, 본 논문의 방법은 안정성을 유지함을 발견

방법 상세 설명

작업 정의

선형 시불변 pH 시스템이 주어졌을 때:

ẋ(t) = (J - R)Hx(t) + Bu(t)
y(t) = B*Hx(t)

여기서 x(t) ∈ C²ⁿ은 상태, u(t) ∈ Cᵖ은 입력, y(t) ∈ Cᵖ은 출력이다. 목표는 동일한 형식이지만 상태 공간 차원이 더 작은 시스템을 찾아서, 두 전달 함수가 미리 정해진 주파수 점에서 같도록 하는 것이다.

핵심 이론 프레임워크

심플렉틱 기하학 기초

정의 1 (심플렉틱 형식): V를 C 위의 짝수 차원 벡터 공간이라 하고, 심플렉틱 형식 Ω: V × V → C가 다음을 만족한다고 하자:

반대칭성: Ω(v₁, v₂) = -Ω(v₂, v₁)
비퇴화성: 모든 v₁ ∈ V에 대해 Ω(v₁, v₂) = 0 ⟺ v₂ = 0

정의 2 (심플렉틱 사상): 행렬 Q가 (J₂ₙ, J₂ₖ)-심플렉틱이라 함은:

Q*J₂ₙQ = J₂ₖ

무손실 경우의 심플렉틱 MOR

정리 4: 무손실 pH 시스템(R = 0)에 대해, Q ∈ C²ⁿˣ²ᵏ가 (J₂ₙ⁻¹, J₂ₖ⁻¹)-심플렉틱 행렬이고 다음을 만족한다면:

ran(B) ⊆ ran(Q)
(sₘI₂ₙ - J₂ₙH)⁻¹Bμₘ ∈ ran(Q)

축소된 pH 시스템은 주어진 주파수에서 원래 시스템과 같다.

손실이 있는 경우의 심플렉틱 MOR

정리 5: 손실이 있는 pH 시스템에 대해, 확장 행렬을 구성한다:

QE := [Q₁  0  ]
      [0   Q₂]

여기서 Q₂는 (J₂ₙ, J₂ₖ)-심플렉틱이고, Q₁*는 Q₂의 좌역이다. 적절한 조건 하에서, 축소된 시스템은 보간 조건을 만족한다.

블록 구조 보존 알고리즘

알고리즘 1 (심플렉틱-pH MOR):

보간 데이터에 따라 V₁ 구성
V₂ := HV₁ 계산
V₁과 V₂를 블록 구조에 따라 분할
공식(27)을 통해 Ψ₁, Ψ₂, Φ₁, Φ₂ 구성
Q₁과 Q₂ 구축
Petrov-Galerkin 투영 적용

구성 공식:

Ψ₂ = orth(V₂₂)
Ψ₁ = V₁₂(Ψ₂*V₁₂)⁻¹
Φ₂ = V₂₁(Ψ₂*Ĵₙ*V₂₁)⁻¹Ĵₖ*
Φ₁ = V₁₁(Φ₂*V₁₁)⁻¹

기술적 혁신점

심플렉틱 구조 보존: 심플렉틱 사상을 통해 해밀턴 시스템의 기본 기하학적 구조 보존
블록 대각 구조 유지: 특수한 Q₁과 Q₂ 구성으로 물리적 의미의 블록 구조 보존 보장
보간 조건: 축소된 시스템이 미리 정해진 주파수 점에서 원래 시스템과 정확히 일치하도록 보장
통일된 프레임워크: MSD 시스템의 SO-Arnoldi 방법을 특수한 경우로 포함

실험 설정

수치 예제: LRCR 회로

그림 1에 표시된 LRCR 회로를 고려하며, n=50개 단위 포함, 매개변수 설정:

커패시턴스: Cᵢ = 1 pF
인덕턴스: Lᵢ = 1 nH
저항: Rᵢ,₁ = 1 mΩ, Rᵢ,₂ = 1 kΩ

평가 지표

H₂ 오차:

H₂error = √(1/2π ∫ Trace((G(iω) - Gred(iω))*(G(iω) - Gred(iω)))dω)

H∞ 오차:

H∞error = max‖G(iω) - Gred(iω)‖₂

비교 방법

pH-Arnoldi: 정리 2 기반의 구조 보존 방법
SO-Arnoldi: 2차 Arnoldi 방법
심플렉틱-pH: 본 논문에서 제안한 방법

보간점 선택

구간 10⁻², 10² 내에서 로그 스케일상 균등하게 보간점 선택, 예를 들어 6개의 보간점: 10⁻²i, -10⁻²i, 10⁰i, -10⁰i, 10²i, -10²i

실험 결과

주요 결과

그림 3에서 다음을 관찰할 수 있다:

심플렉틱-pH vs pH-Arnoldi: 심플렉틱-pH가 대부분의 경우 더 낮은 H₂ 및 H∞ 오차 제공
심플렉틱-pH vs SO-Arnoldi:
- 축소 차원 <24일 때 성능 유사
- 축소 차원=32일 때 SO-Arnoldi 성능 급격히 저하
- 축소 차원 28과 36일 때 SO-Arnoldi가 불안정한 시스템 생성