[論文レビュー] A computable realization of Ruelle's formula for linear response of statistics in chaotic systems

Nisha Chandramoorthy, Qiqi Wang|arXiv (Cornell University)|Feb 10, 2020

Mathematical Dynamics and Fractals参考文献 24被引用数 16

ひとこと要約

本稿では、力学系のカオス的挙動に対して、Ruelleの線形応答公式の計算可能な再定式化である「空間分割感度（S3）」を提案する。S3は感度を安定部分と不安定部分に分離し、安定部分については接線解の平均化、不安定部分についてはエルゴード的・双曲的構造を用いた正則化を施す。これにより、$N$ 個の位相点を用いた ${\cal O}(1/\sqrt{N})$ の収束速度を達成し、一様双曲的吸引子に対して、有限差分法に比べて顕著に高い精度と効率を発揮する。

ABSTRACT

We present a computable reformulation of Ruelle's linear response formula for chaotic systems. The new formula, called Space-Split Sensitivity or S3, achieves an error convergence of the order ${\cal O}(1/\sqrt{N})$ using $N$ phase points. The reformulation is based on splitting the overall sensitivity into that to stable and unstable components of the perturbation. The unstable contribution to the sensitivity is regularized using ergodic properties and the hyperbolic structure of the dynamics. Numerical examples of uniformly hyperbolic attractors are used to validate the S3 formula against a naïve finite-difference calculation; sensitivities match closely, with far fewer sample points required by S3.

研究の動機と目的

• カオス的力学系における統計的感度を計算する課題に取り組むこと。従来の接線／随伴法は不安定性のため失敗する。
• シャドーイング法や遷移作用素の摂動法のような既存手法の代替手段として、収束問題や高コストに悩まされる計算可能性の高い代替手法を開発すること。
• カオス的力学の双曲的構造とエルゴード的性質を活用して、Ruelleの線形応答公式を数値計算に適した形に再定式化すること。
• 位相点 $N$ 個を用いて ${\cal O}(1/\sqrt{N})$ の収束速度を達成し、有限差分法に比べてはるかに少ないサンプル数で正確な感度推定が可能になるようにすること。
• 一様双曲的吸引子に対して数値的検証を行い、有限差分法の結果と密接に一致することを示し、はるかに少ない軌道数で実現可能であることを示すこと。

提案手法

• 位相空間を不安定な葉に測度的分割することで、感度を安定・不安定多様体からの寄与に分解する。
• 安定寄与については、安定寄与が良好に振る舞い、サンプリングに適していることを利用し、軌道に沿った接線方程式解の有限サンプル平均を計算する。
• 不安定寄与は、不安定多様体上で部分積分を施し、SRB測度の分解と条件付き測度の絶対連続性を用いて正則化する。
• 再定式化された公式は、不安定な葉の商空間上の積分として表現され、被積分関数には観測関数の不安定方向への微分と条件付き測度の密度が含まれる。
• 相関の指数的減衰と中心極限定理を用いてサンプリング誤差を抑え、収束速度 ${\cal O}(1/\sqrt{N})$ を保証する。
• 最終的な公式は、$N$ 個の初期条件に対するモンテカルロサンプリングにより実装され、不安定部分は時間発展状態への有限和として計算され、不安定接線ベクトルの直接積分を回避する。

実験結果

リサーチクエスチョン

• RQ1 Ruelleの線形応答公式は、カオス的力学系に対して、数値的に安定かつ計算可能な形に再定式化可能か？
• RQ2 エルゴード的・双曲的構造を用いて、不安定成分の感度を正則化することで、接線法で見られる指数的発散を回避できるか？
• RQ3 安定・不安定寄与を分離する空間分割アプローチにより、$N$ 個のサンプルで ${\cal O}(1/\sqrt{N})$ の収束速度が達成可能か？
• RQ4 S3手法は、有限差分法に比べてはるかに少ない軌道数で正確な感度推定が可能か？
• RQ5 S3公式は、一様双曲的系において、基準となる有限差分計算と照合した場合に、ロバストで正確であるか？

主な発見

• S3手法は、中心極限定理と相関の指数的減衰を用いた解析により、不安定感度成分について ${\cal O}(1/\sqrt{N})$ の収束速度を達成することが証明された。
• 不安定寄与は、不安定多様体上で部分積分を施すことにより正則化され、不安定接線ベクトルを直接計算する必要がなくなる。
• 一様双曲的吸引子に対する数値的検証では、S3の感度が、有限差分計算結果と密接に一致しており、サンプル数がはるかに少ない状況下でも成立している。
• 本手法は、安定寄与と不安定寄与を明確に分離できており、安定部分は有限サンプル平均で計算され、不安定部分は時間発展状態への有限和で計算される。
• 不安定部分の誤差境界は $e_{N,M} \leq c^\mathrm{u}/\sqrt{N} + c_2 \gamma^M / (1 - \gamma)$ と示され、$M \to N$ のとき、支配的項は ${\cal O}(1/\sqrt{N})$ になる。
• 理論的枠組みは、SRB測度の分解と不安定な葉における条件付き測度の絶対連続性に依拠しており、再定式化における数学的厳密性を保証する。