[논문 리뷰] Gradual approximation of the domain of attraction by gradual extension of the "embryo" of the transformed optimal Lyapunov function

Éva Kaslik, A.M. Balint외 1명|arXiv (Cornell University)|2002. 12. 20.

Quantum chaos and dynamical systems참고 문헌 6인용 수 5

한 줄 요약

이 논문은 자동적 분석 ODE 시스템에서 점점 더 정확한 근사치를 제공하는 점점 확장되는 영역을 제공하는 새로운 반복적 방법을 제시한다. 이는 변환된 최적 라플라스 함수의 자연스러운 해석적 영역을 이용하여, 급격히 수렴하는 급수 전개와 카우치-하다르 수렴 조건을 통해 유도된 해석적 성질을 활용한다. 이 방법은 반복적으로 함수의 '배태기'를 계산함으로써 점차적으로 영역을 확장하며, 기존의 방법보다 벤치마크 다항식 시스템에서 더 정확한 근사치를 제공한다.

ABSTRACT

In this paper an autonomous analytical system of ordinary differential equations is considered. For an asymptotically stable steady state x0 of the system a gradual approximation of the domain of attraction DA is presented in the case when the matrix of the linearized system in x0 is diagonalizable. This technique is based on the gradual extension of the "embryo" of an analytic function of several complex variables. The analytic function is the transformed of a Lyapunov function whose natural domain of analyticity is the DA and which satisfies a linear non-homogeneous partial differential equation. The equation permits to establish an "embryo" of the transformed function and a first approximation of DA. The "embryo" is used for the determination of a new "embryo" and a new part of the DA. In this way, computing new "embryos" and new domains, the DA is gradually approximated. Numerical examples are given for polynomial systems. For systems considered recently in the literature the results are compared with those obtained with other methods.

연구 동기 및 목표

• 자기적 분석 ODE 시스템에서 점점 더 정확한 근사치를 제공하는 점점 확장되는 영역을 제공하는 새로운 반복적 방법을 개발하는 것.
• 일반적인 라플라스 기반 DA 추정 기법이 종종 보수적이거나 불완전한 근사치를 제공하는 데서 비롯되는 한계를 극복하는 것.
• 변환된 최적 라플라스 함수의 자연스러운 해석적 영역을 DA의 정확한 특성화로 활용하는 것.
• 계산적으로 실현 가능한 반복 절차를 제공하여 각 단계에서 함수의 '배태기'를 정밀화함으로써 점차적으로 DA 추정을 확장하는 것.
• 기존 문헌의 접근 방식과 비교하여 알려진 DA를 가진 다항식 시스템에서 이 방법의 성능을 검증하는 것.

제안 방법

• 균형점에서의 자코비안을 대각화하는 동형사상 S를 통해 원래 시스템을 변환하여 새로운 시스템 g = S⁻¹∘f∘S를 도출한다.
• 최적 라플라스 함수 V는 W = V∘S로 변환되며, 이는 선형 비동차 PDE ⟨∇W, g⟩ = -‖Sz‖²를 만족하고, W(0) = 0을 충족한다.
• W의 테일러 급수 계수는 제2차부터 시작하여 정리 3를 이용해 재귀적으로 계산되며, 고차수 계수는 저차수 계수에 의존한다.
• W의 급수 수렴 영역 D⁰은 카우치-하다르 공식을 통해 결정되며, 초기 DA 추정치는 DA⁰ = S(D⁰)로 얻어진다.
• 새로운 '배태기'는 고차수 항까지의 급수 계산을 확장함으로써 생성되며, 반복적으로 새로운 수렴 영역 Dᵏ를 계산하여 확장된 DA 추정치 DAᵏ = S(Dᵏ)를 도출한다.
• 이 과정은 단계별로 반복되며, 각 반복에서 점차 고차수 항의 수렴 특성을 분석함으로써 DA 근사치에 새로운 영역을 추가한다.

실험 결과

연구 질문

• RQ1 점점 더 정확한 근사치를 제공하는 점점 확장되는 영역을 제공하는 새로운 반복적 방법을 개발하는 것.
• RQ2 선형 비동차 PDE의 해에 대한 자연스러운 해석적 영역이 정확한 DA를 특성화하는 데 어떻게 활용될 수 있는가?
• RQ3 변환된 라플라스 함수의 '배태기'를 반복적으로 확장함으로써 표준 기법 대비 DA 추정의 정확도가 얼마나 향상되는가?
• RQ4 이 방법은 알려진 또는 알려지지 않은 DA를 가진 다항식 시스템에서 계산적으로 실현 가능하고 수렴 행동이 어떻게 되는가?
• RQ5 기존의 DA 추정 기법과 비교하여 이 방법은 커버리지와 정밀도 측면에서 어떻게 성능을 발휘하는가?

주요 결과

• 예제 1에서, 이 방법은 세 번의 단계를 거쳐 진짜 원형 영역의 주변을 정확하게 포착한 DA 추정치를 달성하였다.
• 3차원 시스템인 예제 2에서는 두 번의 단계를 거친 후 DA 추정치가 크게 확장되었으며, 진짜 쌍곡면 경계에 매우 가까운 근사치를 제공하였다.
• 2차원 시스템인 예제 3([10]에서 유래)에서는 4단계 추정치가 [10]에서 보고된 추정치보다 더 넓은 영역을 커버하여 우수한 커버리지가 입증되었다.
• 3차원 시스템인 예제 4에서는 이전의 접근 방식보다 두 번째 단계에서 더 광범위한 DA 추정치를 생성하였으며, 첫 번째 단계의 추정치 역시 의미 있는 확장을 보였다.
• 이 방법의 계산 비용은 고차수 근사치에 따라 크게 증가하였으며, 표 1에서 볼 수 있듯이 예제 2의 두 번째 단계는 14.7시간이 소요되어 정확도와 계산 시간 사이의 상충 관계가 있음을 시사한다.
• 반복 절차는 초기 수렴 영역을 초월하여 DA 추정치를 점차적으로 확장하는 데 성공하였으며, 이는 함수의 해석적 구조를 점차 정밀화함으로써 진짜 영역을 점차적으로 근사할 수 있음을 확인한다.