[논문 리뷰] Lyapunov and Sylvester equations: A quadrature framework

Christian Bertram, Heike Faßbender|arXiv (Cornell University)|2019. 03. 13.

Matrix Theory and Algorithms참고 문헌 31인용 수 4

한 줄 요약

이 논문은 대규모 리아푸노프 및 실베스터 방정식을 해결하기 위해 구조를 유지하는 룬게쿠타 방법을 관련된 상미분방정식 시스템에 적용하는 새로운 적분 프레임워크를 제안한다. 대각선에 음의 룬게쿠타(디아르티) 방법을 통해 저랭크 잔차 유지 기법을 활용함으로써, ADI 반복과 동치가 되는 해법을 얻어내어 복소수 산술을 위한 실수화 전략을 포함한 효율적인 저랭크 근사가 가능해진다.

ABSTRACT

This paper introduces a novel framework for the solution of (large-scale) Lyapunov and Sylvester equations derived from numerical integration methods. Suitable systems of ordinary differential equations are introduced. Low-rank approximations of their solutions are produced by Runge-Kutta methods. Appropriate Runge-Kutta methods are identified following the idea of geometric numerical integration to preserve a geometric property, namely a low rank residual. For both types of equations we prove the equivalence of one particular instance of the resulting algorithm to the well known ADI iteration. As the general approach suggested here leads to complex valued computation even for real problems, we present a general realification approach based on similarity transformation.

연구 동기 및 목표

• 대규모 리아푸노프 및 실베스터 방정식을 풀기 위한 새로운 수치적 프레임워크를 개발하기 위해 상미분방정식의 수치적 적분을 활용한다.
• 기하학적 수치적 통합 원리를 통해 반복적 해법 단계에서 잔차의 저랭크 구조를 유지한다.
• 제안된 룬게쿠타 기반 알고리즘과 고전적인 ADI 반복 간의 동치성을 리아푸노프 및 실베스터 방정정에 대해 확립한다.
• 실수 문제에서 발생하는 복소수 연산을 다루기 위해 유사변환 기반 실수화 기법을 제공한다.

제안 방법

• 리아푸노프 방정식의 해는 시간에 따라 변하는 그람형의 적분으로 재정의되며, 이를 상미분방정식 시스템의 해로 해석한다.
• 해의 저랭크 근사를 생성하기 위해 상미분방정식 시스템에 룬게쿠타 방법을 적용하며, 안정성 함수는 부타흐 표기법에서 유도된다.
• 초기 잔차의 저랭크 성질을 유지하기 위해 안정성 함수가 특정 함수 항등식을 만족하도록 제약을 둔다.
• 실베스터 방정식의 경우 유사한 상미분방정식 시스템을 유도하고, 동일한 룬게쿠타 접근법을 적용함으로써 저랭크 반복적 방법을 도출한다.
• 이 프레임워크는 디아르티 방법이 저랭크 잔차를 유지함을 규명하여 ADI 반복과의 동치성을 이끌어낸다.
• 유사변환을 통한 실수화 전략을 도입하여 복소수 룬게쿠타 반복을 실수 산술로 변환함으로써 실수 문제에 적용 가능하게 한다.

실험 결과

연구 질문

• RQ1 상미분방정식 통합 기반의 적분 프레임워크를 통해 리아푸노프 및 실베스터 방정식에 대한 저랭크 반복 해법을 유도할 수 있는가?
• RQ2 어떤 유형의 룬게쿠타 방법이 해법 과정에서 잔차의 저랭크 구조를 유지하는가?
• RQ3 유도된 알고리즘이 리아푸노프 및 실베스터 방정식 모두에 대해 잘 알려진 ADI 반복과 동치인가?
• RQ4 이 방법에서 발생하는 복소수 연산은 실수 산술로 어떻게 변환할 수 있는가?
• RQ5 룬게쿠타 방법의 안정성 함수와 유리적 크릴로프 부분공간 프레임워크 사이의 관계는 무엇인가?

주요 결과

• 디아르티 방법을 사용할 경우, 제안된 리아푸노프 방정식에 대한 룬게쿠타 기반 알고리즘은 수학적으로 ADI 반복과 동치이다.
• 구조적으로 반복 과정 전반에 걸쳐 잔차의 저랭크 구조를 유지함으로써 계산 효율성을 보장한다.
• 실베스터 방정식의 경우 동일한 프레임워크를 통해 ADI 실베스터 변형과 동치인 알고리즘이 도출되었으며, 이는 알고리즘 4에서 확인된다.
• 저랭크 잔차를 유지하기 위해 룬게쿠타 방법의 안정성 함수는 특정 함수 항등식을 만족해야 하며, 이는 유효한 방법의 집합을 특성화한다.
• 유사변환을 통한 실수화 접근법은 정확도나 효율성을 잃지 않고 복소수 산술 방법을 실수 문제에 적용할 수 있게 한다.
• 이 프레임워크는 기하학적 수치적 통합과 유리적 크릴로프 및 모멘트 매칭 접근법을 연결하며, ADI 유사 방법의 통합적 유도를 제공한다.