공유 선형 해석을 이용한 다중 Lyapunov Sylvester Riccati 방정식 통합 ADI 프레임워크

초록

**
본 논문은 대규모 Lyapunov, Sylvester, Riccati 방정식들을 동시에 해결하기 위해 ADI(Alternating Direction Implicit) 방법의 선형 시스템 풀이를 공유하는 통합 프레임워크를 제안한다. 저차원 투영을 통한 Petrov‑Galerkin 해석과 보간 특성을 분석하고, 기존 방법과 달리 극점 배치는 작은 규모 연산에 국한시켜 선형 시스템 풀이만 두 번 수행하면 6개의 Lyapunov, 1개의 Sylvester, 10개의 Riccati 방정식을 동시에 얻을 수 있음을 보인다. 수치 실험을 통해 연산량 감소와 정확도 유지가 확인되었다.

**

상세 분석

**
이 논문은 ADI 방법이 실제로는 Petrov‑Galerkin 투영 알고리즘이며, 각 반복 단계에서 ADI 시프트의 거울 이미지에 대한 보간(interpolation)을 강제한다는 점을 명확히 증명한다. 기존 연구에서는 Lyapunov 방정식에만 적용된 이 보간‑극점 관계가 Sylvester와 Riccati 방정식에도 동일하게 존재하지만, 후자들의 경우 극점(pole) 배치가 보간점과 일치하지 않는다. 대신, 극점은 투영된 작은 차원의 방정식이 유일해지도록 선택되며, 이는 전적으로 작은 행렬 연산에 의해 수행된다. 이러한 차이를 인식함으로써 저자들은 ‘선형 시스템 풀이(shifted linear solves)’가 모든 방정식에 공통적으로 필요하고, 실제 비용의 대부분을 차지한다는 사실을 활용한다.

핵심 아이디어는 두 개의 시프트된 선형 시스템
\