[논문 리뷰] Methods for verified stabilizing solutions to continuous-time algebraic Riccati equations

Tayyebe Haqiri, Federico Poloni|arXiv (Cornell University)|2015. 09. 07.

Numerical Methods and Algorithms참고 문헌 29인용 수 10

한 줄 요약

이 논문은 개선된 크라체프스키 및 고정점 방법을 활용하여 연속시간 대수적 린카티 방정식(CAREs)의 안정화 해를 검증 가능한 계산 방법을 제시한다. 기저 변환, 간격 산술, 그리고 새로운 기울기 행렬 공식화를 통해 O(n³) 복잡도를 갖는 검증 가능한 포함 구간을 도출하며, 기존 최첨단 알고리즘보다 여러 벤치마크에서 성능을 뛰어넘고, 비대각화 가능한 닫힌 루프 행렬의 경우에도 신뢰성 유지한다.

ABSTRACT

We describe a procedure based on the Krawczyk method to compute a verified enclosure for the stabilizing solution of a continuous-time algebraic Riccati equation $A^*X+XA+Q=XGX$ building on the work of [B.~Hashemi, \emph{SCAN} 2012] and adding several modifications to the Krawczyk procedure. We show that after these improvements the Krawczyk method reaches results comparable with the current state-of-the-art algorithm [Miyajima, \emph{Jpn. J. Ind. Appl. Math} 2015], and surpasses it in some examples. Moreover, we introduce a new direct method for verification which has a cubic complexity in term of the dimension of $X$, employing a fixed-point formulation of the equation inspired by the ADI procedure. The resulting methods are tested on a number of standard benchmark examples.

연구 동기 및 목표

• 연속시간 대수적 린카티 방정식(CAREs)의 유일한 안정화 해를 포함하는 엄밀한 간격 행렬 포함 구간을 계산하는 검증 가능한 수치적 방법을 개발하는 것.
• 계산 복잡도를 감소시키고 간격 산술에서의 감싸임 효과를 최소화하여 CAREs에 대한 검증 가능한 계산의 효율성과 신뢰성을 향상시키는 것.
• 기존 방법의 한계—특히 비대각화 가능하거나 악조건인 닫힌 루프 행렬에서의 실패—를 해결하기 위해 기저의 대각화 가능성에 의존하지 않는 새로운 고정점 공식화를 도입하는 것.
• 계산된 해 포함 구간이 안정성 특성(행렬 A - GX의 허르비츠 안정성)을 검증할 수 있도록 보장하는 것.
• 표준 CARE 테스트 문제 세트에 걸쳐 제안된 방법들을 최첨단 알고리즘과 비교하기 위한 벤치마크 수행 및 성능 비교

제안 방법

• 해밀토니안 행렬의 불변부공간 구조를 통한 기저 변환을 활용하여 크라체프스키 방법을 CAREs에 적용함으로써 간격 산술 연산을 감소시키고 감싸임 효과를 최소화한다.
• 비검증 Riccati 해법에 영감을 얻은 기저 변환을 적용하여 CARE를 재구성함으로써 안정화 해의 노름이 유한해지도록 하여 더 좁은 간격 포함 구간을 가능하게 한다.
• 기울기 행렬의 표준 간격 자료 평가 대신 대수적으로 유도된 표현식을 도입하여 더 작고 정확도가 높은 간격 포함 구간을 얻는다.
• ADI 반복 프레임워크를 기반으로 한 새로운 고정점 공식화를 도입하여, 닫힌 루프 행렬의 대각화 가능성에 의존하지 않고도 반복마다 O(n³) 복잡도를 갖는다.
• 모든 고유값이 열린 왼쪽 반평면에 위치함을 확인하여 계산된 간격 해의 안정성 특성을 검증하는 절차를 구현한다.
• 간격 산술 및 검증된 선형대수(예: 검증된 행렬 분해)를 사용하여 고유값 계산을 포함한 모든 연산이 엄밀히 경계됨을 보장한다.

실험 결과

연구 질문

• RQ1 기저 변환과 개선된 기울기 행렬 추정을 통한 크라체프스키 방법의 개선은 경쟁 가능한 성능과 더 좁은 포함 구간을 갖는 검증 가능한 CARE 해를 달성할 수 있는가?
• RQ2 불변부공간 변환을 통한 CARE 재구성은 검증 계산에서 더 나은 조건 문제와 향상된 간격 포함 구간을 가져오는가?
• RQ3 ADI 절차를 기반으로 한 고정점 방법은 닫힌 루프 행렬이 열화되거나 거의 열화된 경우 크라체프스키 기반 방법의 대안으로 신뢰할 수 있는가?
• RQ4 제안된 방법들은 [25]의 최첨단 알고리즘과 비교하여 계산 비용, 신뢰성, 포함 구간 품질 측면에서 어떻게 성능을 내는가?
• RQ5 검증된 방법의 확장성은 행렬 차원 n에 대해 어떻게 되며, 실질적으로 O(n³) 복잡도를 유지하는가?

주요 결과

• 개선된 크라체프스키 방법은 [25]의 최첨단 알고리즘과 유사한 성능을 보이며, 여러 벤치마크 예제에서 더 좁은 해 포함 구간을 제공한다.
• 새로운 고정점 방법은 기저의 대각화가 불가능하거나 악조건인 경우에도 성공적으로 처리하며, 크라체프스키 기반 방법이 실패할 수 있는 상황에서 유용하다.
• 테스트된 네 가지 방법(H, M, K, F) 모두 행렬 차원 n에 대해 약 O(n³)로 스케일링됨을 실험으로 확인하였다. (크기 10에서 1000까지의 문제 세트 기준)
• 실제로 가장 빠른 방법은 [25] 기반의 Method M이지만, Method K(개선된 크라체프스키)는 더 높은 신뢰성과 낮은 실패 빈도를 보이며, 특히 악조건 문제에서 유리하다.
• Method F는 가장 높은 실패율을 보이지만, 닫힌 루프 행렬이 대각화 불가능한 경우에만 특별히 효과적이므로, 그 특수한 용도에 유용하다.
• 가장 큰 테스트 케이스(n = 1000)에서 유일하게 해와 그 안정성 특성을 모두 검증한 것은 Method K로, 극단적인 경우에서도 뛰어난 강건성을 입증하였다.