[論文レビュー] Estimation of initial conditions from a scalar time series

Anil Maybhate, R. E. Amritkar|ArXiv.org|Feb 17, 2000

Time Series Analysis and Forecasting参考文献 1被引用数 3

ひとこと要約

本稿では、1つのスカラー時系列から多次元力学系の初期状態をニュートン・ラフソン法を応用した新規手法により推定する。時間発展を組み込んだ変更された反復スキームを用いることで、条件付きリャプノフ指数が正であっても高速かつ2次収束を達成し、スカラー信号から得られる同一初期状態を用いた力学系の即時同期を可能にする。

ABSTRACT

We introduce a method to estimate the initial conditions of a mutivariable dynamical system from a scalar signal. The method is based on a modified multidimensional Newton-Raphson method which includes the time evolution of the system. The method can estimate initial conditions of periodic and chaotic systems and the required length of scalar signal is very small. Also, the method works even when the conditional Lyapunov exponent is positive. An important application of our method is that synchronization of two chaotic signals using a scalar signal becomes trivial and instantaneous.

研究の動機と目的

• 1つの状態変数のスカラー時系列から多次元力学系の全初期状態を再構成する課題に対処すること。
• 既存手法がしばしば失敗する、条件付きリャプノフ指数が正である状態でも信頼性高く動作する手法を開発すること。
• スカラー信号からの初期状態の正確な推定により、同一初期状態を用いた力学系の即時かつ容易な同期を可能にすること。
• 標準的な埋め込み遅延時間よりも短い長さの時系列でも、初期状態の回復が正確に可能であることを示すこと。

提案手法

• 本手法は、時間ステップの逐次更新を通じてシステムの時間発展を組み込んだ、変更された多次元ニュートン・ラフソン法を用いる。
• 試行軌道 y^n と真の軌道 x^n の差を w^n = y^n - x^n と定義し、w^n = 0 を求める。
• 差のベクトルの時間発展は行列更新式 W^n = A^{n-1} A^{n-2} ... A^0 W^0 でモデル化され、ここで A^n = I + Δt J^n であり、J^n は y^n におけるシステムのベクトル場のヤコビ行列である。
• 時間発展された誤差ベクトルに基づいて、システムのダイナミクスの勾配を用いて初期推定値 y^0 を反復的に補正する。
• 収束は誤差ノルム e_i = ||x^0 - y^0||_i で評価され、収束速度は e_{i,n+1} ≈ α e_{i,n}^2 におけるパrameter α で定量化される。
• 本手法は、ローレンツ系、チャア回路、ローレンツ系、ディスクダイナモ、三波プラズマ結合系、4次元位相変換回路など、周期的・カオス的・間欠的ダイナミクスを示す多様な系に適用され、多様なダイナミクス行動において一貫した性能を示した。

実験結果

リサーチクエスチョン

• RQ1 多次元力学系の全初期状態ベクトルは、1つの成分のスカラー時系列から正確に再構成可能か？
• RQ2 条件付きリャプノフ指数が正であっても、本手法は有効に機能するか？（通常は同期を困難にする条件である。）
• RQ3 信頼性のある初期状態推定に必要な時系列長はどの程度短くできるか？特に標準的埋め込み手法と比較して。
• RQ4 収束速度および耐障害性は、ダイナミクスの状態（周期的・カオス的・間欠的）にどの程度依存するか？

主な発見

• 本手法は2次収束を達成し、ローレンツ系の x₁ 時系列からの初期状態推定において、e₁ について α ≈ 1.97、e₃ について α ≈ 1.95 となる。
• x₂ 時系列を用いた場合、条件付きリャプノフ指数が正であっても、収束速度は2次（e₁ について α ≈ 1.97、e₃ について α ≈ 1.95）を維持する。
• x₃ 時系列を用いた場合、2次未塔の収束（e₁ について α ≈ 1.28、e₂ について α ≈ 1.30）となるが、依然として信頼性高く収束する。
• チャア回路の限界円形軌道においては、収束が2次未塔（e₂ について α ≈ 1.29、e₃ について α ≈ 1.25）となるが、初期状態の推定に成功している。
• 本手法は、ローレンツ系、ディスクダイナモ、三波プラズマ結合系、4次元位相変換回路など、多様な系において、周期的・カオス的・間欠的ダイナミクスのすべての状態で有効に機能する。
• 条件付きリャプノフ指数が正であっても、スカラー信号から推定した初期状態をレプリカの初期状態に設定することで、同一のカオス的系の即時同期が可能になる。