[논문 리뷰] Lagrange approximation of transfer operators associated with holomorphic data

Oscar F. Bandtlow, Julia Slipantschuk|arXiv (Cornell University)|2020. 04. 07.

Quantum chaos and dynamical systems인용 수 6

한 줄 요약

이 논문은 해석적 데이터를 갖는 전이 연산자에 대해 라그랑주 보간 기반 스펙트럼 근사 방법을 제안하며, 고유값과 고유함수의 지수 수렴을 달성한다. 이 방법은 조건점과 다항식 보간을 사용하여 행렬 근사식을 구성하며, 수렴 속도는 기저 역동계의 복소수 수축 비율에 의해 결정된다.

ABSTRACT

We show that spectral data of transfer operators given by holomorphic data can be approximated using an effective numerical scheme based on Lagrange interpolation. In particular, we show that for one-dimensional systems satisfying certain complex contraction properties, spectral data of the approximants converge exponentially to the spectral data of the transfer operator with the exponential rate determined by the respective complex contraction ratios of the underlying systems. We demonstrate the effectiveness of this scheme by numerically computing eigenvalues of transfer operators arising from interval and circle maps, as well as Lyapunov exponents of (positive) random matrix products and iterated function systems, based on examples taken from the literature.

연구 동기 및 목표

• 해석적 역동계에서 유도된 전이 연산자에 대한 수치적으로 효율적인 스펙트럼 근사 방법을 개발하기 위해.
• 기존의 프로젝션 기반 방법(예: 욘의 방법)이 지닌 비최적 수렴 속도를 극복하기 위해.
• 조건점에서의 라그랑주 보간을 활용해 고유값과 고유함수의 지수 수렴을 확립하기 위해.
• 리아푸노프 지수와 정적 측도를 포함한 스펙트럼 데이터를 계산하기 위한 실용적이고 실행 가능한 알고리즘을 제공하기 위해.
• 구간 사상, 원주 사상, 랜덤 행렬 곱, 반복 기하 체계에서의 수치 예제를 통해 방법의 효과성을 입증하기 위해.

제안 방법

• 이 방법은 라그랑주 보간을 활용해 $n$개의 조건점에서의 다항식 보간을 정의하는 프로젝션 연산자 $\Pi_n$를 정의한다.
• 전이 연산자 $\mathcal{L}$는 $\mathcal{L}_n = \Pi_n \mathcal{L} \Pi_n$로 근사되며, 조건점 $x_i$에서 $\mathcal{L}e_j$의 평가를 통해 $n \times n$ 행렬로 표현된다.
• 조건점은 체비셰프 다항식의 근 또는 타원형 경로 위에 선택되어 스펙트럼 수렴을 보장한다.
• 이 방법은 갈레르킨 또는 유한 섹션 방법에서 요구되는 고비용의 적분과 도함수 계산을 피하기 위해, 조건점 또는 스펙트럼 갈레르킨 접근 방식을 구현한다.
• 이 방법은 사상 성분의 해석적 연장성을 활용해 타원 영역 $E_{\gamma,R}$를 정의함으로써 복소수 수축 분석이 가능해진다.
• 근사식 $\mathcal{L}_n$의 고유함수 기능을 활용해 불변 측도에 대한 적분과 리아푸노프 지수를 계산한다.

실험 결과

연구 질문

• RQ1 라그랑주 보간이 해석적 데이터를 갖는 전이 연산자의 스펙트럼 근사에 대해 지수 수렴을 제공할 수 있는가?
• RQ2 라그랑주-체비셰프 조건점 방법의 수렴 속도는 전통적 방법(예: 욘의 유한 섹션 방법)과 비교해 어떻게 되는가?
• RQ3 복소수 수축 비율이 고유값과 고유함수의 지수 수렴 속도를 결정하는 데 어떤 역할을 하는가?
• RQ4 이 방법은 해석적 생성자를 갖는 반복 기하 체계에 대해 리아푸노프 지수와 정적 측도를 효과적으로 계산할 수 있는가?
• RQ5 이 방법은 해석적 연장이 가능한 구간 및 원주 사상의 전이 연산자에 대해 어느 정도 적용 가능한가?

주요 결과

• 근사식의 스펙트럼 데이터는 기저 시스템의 복소수 수축 비율 $r/R$에 의해 지배되는 지수 속도로 전체 전이 연산자의 스펙트럼 데이터로 수렴한다.
• 구간 사상 예제에서 $\Phi_1(x) = \frac{1}{6}\sin(\pi x/4) + \frac{1}{4}$, $\Phi_2(x) = \frac{1}{3}\sin(\pi x/4) + \frac{2}{3}$ 일 때, $n=100$에서 리아푸노프 지수를 89자리 정밀도로 계산하였으며, 주기 18 궤도로부터 이전에 얻은 결과와 일치한다.
• 최적 수축 비율은 수치적으로 $r/R \approx 0.4138$로 도출되었으며, 이는 리아푸노프 지수 근사에 대해 $O((0.4138)^n)$의 지수 수렴 속도를 의미한다.
• 원주 사상 예제에서는 수축 비율 $r/R \approx 0.53$로 지수 수렴을 달성하였으며, 이는 이론적 예측과 일치한다.
• 해석성 가정 하에 기존의 욘 유형 방법(수렴 속도 $O(\log n / n)$)에 비해 지수 수렴을 달성함으로써 이 방법이 우월함을 입증하였다.
• 이 방법은 수치적 적분과 고차 도함수 계산을 피하며, 체비셰프-가우스 quadrature를 활용해 행렬 평가를 수행할 수 있어 실용적으로 효율적이다.