시간 제한 H₂ 최적 근사 모델 차원 축소 알고리즘

논문은 먼저 모델 차원 축소(MOR)의 필요성을 강조한다. 현대의 물리 시스템은 수천 개의 상태 방정식으로 기술되며, 직접 시뮬레이션은 계산 비용이 prohibitive하다. 기존 MOR 기법으로는 Balanced Truncation(BT), Modal Truncation, Moment Matching, 그리고 H₂‑optimal IRKA 등이 있다. BT는 안정성 보장과 a priori error bound를 제공하지만, 대규모 Lyapunov 방정식 해결이 필요해 실용성이 떨어진다. Moment Matching은 Krylov 서브스페이스를 이용해 효율적으로 보간점을 맞추지만, 전통적인 H₂‑optimal 조건을 완전 만족하지는 않는다. IRKA는 Wilson 조건을 만족하도록 반복적으로 보간점을 업데이트하지만, 수렴 보장이 없고, 안정성도 보장되지 않는다. 시간‑제한 모델 차원 축소는 시스템의 특정 시간 구간(예: 0~t초) 내에서 정확도를 높이는 것이 목표이며, 이를 위해 Time‑Limited BT(TLBT)와 Time‑Limited IRKA가 제안되었다. 그러나 TLBT는 두 개의 대규모 Lyapunov 방정식을 필요로 하고, 안정성 보장이 약하며, TL‑IRKA는 큰 시스템에 적용하기 위해 pole‑residue 형태를 요구하는 등 실용적 한계가 있다. 본 논문은 이러한 한계를 극복하고자, 시간‑제한 H₂‑norm(‖·‖_{H₂,t})에 대한 1차 최적조건을 기반으로 한 새로운 알고리즘을 제시한다. 핵심 아이디어는 Pseudo‑Optimal Rational Krylov(PORK) 알고리즘을 시간‑제한 상황에 맞게 확장하고, 파라미터화 기법을 도입해 보간점·접선 방향을 사전에 선택한 뒤, 반복 없이 최적조건의 일부를 만족하도록 ROM을 구성하는 것이다. 구체적인 절차는 다음과 같다. 1. 보간점 σ_i (i=1…r)와 오른쪽 접선 방향 ĉ_i 를 선택한다. 2. 입력 Krylov 서브스페이스 V̂_r = span{(σ_i I−A)^{-1}B ĉ_i } 를 구축한다. 3. 출력 서브스페이스 Ŵ_r = V̂_r (V̂_r^T V̂_r)^{-1} 로 정의한다. 4. 투영을 통해 Â_r = Ŵ_r^T A V̂_r, B̂_r = Ŵ_r^T B, Ĉ_r = C V̂_r 를 얻는다. 5. B⊥ = B − V̂_r B̂_r, L̂_r = (B⊥^T B⊥)^{-1} B⊥^T (A V̂_r − V̂_r Â_r) 를 계산한다. 6. Ŝ_r = Â_r − B̂_r L̂_r 로 정의하고, Sylvester 방정식 A V̂_r + V̂_r(−Ŝ_r) + B(−L̂_r)=0 를 만족하도록 V̂_r 를 조정한다. 이때 σ_i는 Ŝ_r의 고유값이 되며, L̂_r는 각 보간점에서의 잔여를 조정한다. 파라미터 ξ를 ξ = −Q̂_s^{-1} L̂_r^T 로 설정하면, H₂‑norm 감소식 (21)을 만족하게 되며, 이는 전체 H₂‑optimal 조건 중 ‘pseudo‑optimal’ 부분을 만족한다는 의미이다. 시간‑제한 최적조건을 반영하기 위해, 위 파라미터화 과정에 추가적인 제약을 부여한다. 구체적으로, 시간‑제한 Gramians P_T, Q_T 를 근사하는 행렬 ˜P_T, ˜Q_T 를 Sylvester 형태로 정의하고, 이들에 대한 연립식 (11)–(12)를 만족하도록 ξ 를 미세 조정한다. 이렇게 하면 ROM이 시간‑제한 H₂‑norm ‖E‖_{H₂,t} 를 최소화하는 방향으로 동작한다. 알고리즘의 안정성 보장은 Â_r가 Hurwitz(모든 고유값이 좌반평면)임을 증명함으로써 확보한다. 이는 V̂_r, Ŵ_r 가 A의 안정성을 보존하도록 구성되었기 때문이다. 또한, 제안된 적응형 프레임워크는 CURE(Cumulative Reduction) 아이디어를 차용한다. 초기 보간점·접선 방향을 임의로 선택한 뒤, 매 단계마다 새로운 보간조건을 추가하고 기존 조건을 유지한다. 이 과정에서 오차 ‖E‖_{H₂,t} 가 단조 감소한다는 정리를 증명했으며, 이는 보간점 선택에 대한 민감성을 크게 완화한다. 수치 실험에서는 두 개의 표준 벤치마크(예: 1D 열전달 모델, 대규모 전력 시스템 모델)를 사용했다. 시스템 차원 n≈10⁴, 목표 차원 r≈50 로 설정했으며, 시간 구간 t=15 s 로 제한하였다. 결과는 다음과 같다. - 제안 알고리즘은 TLBT와 TL‑IRKA 대비 시간‑제한 H₂‑오차를 평균 2.5배 감소시켰다. - 모든 실험에서 ROM의 고유값이 좌반평면에 위치해 안정성을 보장했다. - 중요한 물리적 모드(예: 저주파 진동 모드)의 고유값과 잔여가 원본과 거의 일치했다. - 계산 시간은 Krylov 기반 사전 계산 단계와 투영 단계만 포함해 전체 TLBT 대비 30 % 이하였으며, 메모리 사용량도 O(n·r) 수준에 머물렀다. 결론적으로, 본 논문은 (i) 반복‑무료 구현, (ii) 대규모 시스템 적용 가능성, (iii) 시간‑제한 H₂‑optimal 조건의 부분 만족, (iv) 오차 단조 감소와 안정성 보장이라는 네 가지 핵심 장점을 동시에 제공하는 새로운 MOR 프레임워크를 제시한다. 이는 제한된 시간 구간에서의 정확도가 중요한 전력 시스템, 항공우주, 자동차 제어 등 다양한 분야에 직접적인 활용 가치를 제공한다.