[논문 리뷰] A Structure Preserving Krylov Subspace Method for Large Scale Differential Riccati Equations

Antti Koskela, Hermann Mena|arXiv (Cornell University)|2017. 05. 21.

Matrix Theory and Algorithms인용 수 3

한 줄 요약

이 논문은 대규모 대칭 미분 리카티 방정식을 해결하기 위해 구조를 유지하는 켈러르 부분공간 방법을 제안한다. 이 방법은 행렬 A와 초기 조건 X₀, 외부 항 Q의 저랭크 인자들로 생성된 켈러르 부분공간에 문제를 투영함으로써 작동한다. 이 방법은 정확한 해 흐름의 양성과 단조성을 유지하며, 수치 실험을 통해 초선형 수렴성을 입증하고, 효율적인 사후 오차 추정 및 랭크 잘라내기 전략을 제공한다.

ABSTRACT

We consider a Krylov subspace approximation method for the symmetric differential Riccati equation $\dot{X} = AX + XA^T + Q - XSX$, $X(0)=X_0$. The method we consider is based on projecting the large scale equation onto a Krylov subspace spanned by the matrix $A$ and the low rank factors of $X_0$ and $Q$. We prove that the method is structure preserving in the sense that it preserves two important properties of the exact flow, namely the positivity of the exact flow, and also the property of monotonicity. We also provide a theoretical a priori error analysis which shows a superlinear convergence of the method. This behavior is illustrated in the numerical experiments. Moreover, we derive an efficient a posteriori error estimate as well as discuss multiple time stepping combined with a cut of the rank of the numerical solution.

연구 동기 및 목표

• 제어 및 모델 순서 축소에서 발생하는 대규모 대칭 미분 리카티 방정식을 위한 수치적으로 효율적이고 구조를 유지하는 방법을 개발하기 위해.
• 수치적 해가 정확한 흐름의 핵심 기하적 성질, 특히 양성과 단조성을 유지하도록 보장하기 위해.
• 수렴 분석과 실용적 구현을 위한 사전 및 사후 오차 추정을 제공하기 위해.
• 적응형 시간 스텝 및 저랭크 근사 기법을 통해 효율적인 시간 적분을 가능하게 하기 위해.

제안 방법

• 이 방법은 대규모 미분 리카티 방정식을 A와 X₀, Q의 저랭크 인자들로 생성된 켈러르 부분공간에 투영한다.
• 갈레르킨 투영을 사용하여 원래 방정식을 근사하는 저차원 동적 시스템을 유도한다.
• 저차원 시스템은 원래 방정식의 구조를 유지하므로 수치적 해가 정확한 흐름의 양성과 단조성을 물려받는다.
• 사후 오차 추정기를 유도하여 적응형 시간 스텝을 안내하고 해의 정확도를 보장한다.
• 계산 비용을 제어하고 저랭크 구조를 유지하기 위해 주기적으로 랭크 잘라내기를 적용한다.
• 이 방법은 적응형 시간 스텝과 저랭크 갱신을 결합하여 정확도와 효율성을 균형 잡는다.

실험 결과

연구 질문

• RQ1 대규모 미분 리카티 방정식에서 정확한 해 흐름의 양성과 단조성을 유지할 수 있는 켈러르 부분공간 방법을 설계할 수 있는가?
• RQ2 제안된 구조를 유지하는 켈러르 방법의 수렴 거동은 어떠한가? 초선형 수렴을 달성할 수 있는가?
• RQ3 저랭크 미분 리카티 방정식의 맥락에서, 효율적인 사후 오차 추정기를 어떻게 구성할 수 있는가?
• RQ4 랭크 잘라내기가 수치적 해의 정확도와 안정성에 미치는 영향은 어떠한가?
• RQ5 적응형 시간 스텝은 어떻게 저랭크 갱신과 효과적으로 조합하여 효율성과 정확도를 유지할 수 있는가?

주요 결과

• 제안된 방법은 정확한 해 흐름의 양성과 단조성을 유지하여 수치 근사의 기하학적 충실도를 보장한다.
• 이론적 사전 분석과 수치 실험을 통해 초선형 수렴성을 입증하였다.
• 효율적인 사후 오차 추정기를 도출하여 신뢰할 수 있는 적응형 시간 스텝과 오차 제어를 가능하게 하였다.
• 수치 실험은 적응형 시간 스텝과 저랭크 잘라내기의 조합이 정확도를 유지하면서 계산 비용을 줄이는 데 효과적임을 보였다.
• 랭크 감소에도 불구하고 높은 정확도를 유지하여 장시간 시뮬레이션에서의 강건성을 보였다.