ADI와 Rational Krylov 방법의 동등성 및 의사 H₂ 최적 시프트

초록

본 논문은 Sylvester 방정식(특히 Lyapunov 방정식)의 저차원 근사 해를 구하기 위해 널리 사용되는 ADI 반복법과 Rational Krylov 투영법이 특정 시프트 선택, 즉 의사 H₂ 최적 시프트일 경우 완전히 동등함을 증명한다. 이러한 시프트는 Lyapunov 방정식에서 잔차가 Krylov 부분공간에 대해 직교성을 갖게 하며, 실험을 통해 거의 최적에 가까운 저랭크 근사 해를 제공함을 확인한다.

상세 분석

논문은 먼저 Sylvester 방정식 (AX+XB+Y=0) 를 다루며, 여기서 (Y=bc^{}) 와 같이 랭크‑1 형태를 가정한다. 기존에 대규모 시스템에서 사용되는 두 주요 방법인 ADI(Alternating Direction Implicit)와 Rational Krylov Projection Method(RKPM)의 성능은 시프트 선택에 크게 좌우된다는 점을 강조한다. 저자는 “의사 H₂ 최적 시프트(pseudo H₂‑optimal shifts)”라는 개념을 도입한다. 이는 전통적인 H₂ 최적 모델 축소에서 요구되는 인터폴레이션 조건 (H(-\lambda_i)=H_r(-\lambda_i)) 를 만족하지만, 파생 조건 (H’(-\lambda_i)=H’_r(-\lambda_i)) 은 포기한다. 이러한 시프트는 ( \lambda(Q^{}AQ) = -{\sigma_1,\dots,\sigma_r}) 혹은 ( \lambda(U^{*}BU) = -{\sigma_1,\dots,\sigma_r}) 를 만족하도록 선택된다.

핵심 정리(Theorem 2)는 다음을 보인다. 동일한 시프트 집합 (\sigma) 로 구성된 Krylov 부분공간의 정규직교 기저 (Q_r, U_r) 를 사용해 투영된 Sylvester 방정식 \