[論文レビュー] Lyapunov exponents in constrained and unconstrained ordinary differential equations

Michael D. Hartl|arXiv (Cornell University)|Mar 20, 2003

Mathematical Control Systems and Analysis参考文献 2被引用数 5

ひとこと要約

本稿では、制約付きおよび非制約付き常微分方程式系におけるリャプノフ指数の計算に耐障害性のある数値的手法を提示する。制約によって生じる課題に特に注目し、偏差ベクトル法とヤコビアンベース法を比較する。特に関数値分解（SVD）とグラム・シュミット直交化に焦点を当て、非スケーリングSVDが正確な結果をもたらす一方、スケーリングによって軸ラベルの曖昧さが生じ、系統的な誤差を引き起こすことを示している。

ABSTRACT

We discuss several numerical methods for calculating Lyapunov exponents (a quantitative measure of chaos) in systems of ordinary differential equations. We pay particular attention to constrained systems, and we introduce a variety of techniques to address the complications introduced by constraints. For all cases considered, we develop both deviation vector methods, which follow the time-evolution of the difference between two nearby trajectories, and Jacobian methods, which use the Jacobian matrix to determine the true local behavior of the system. We also assess the merits of the various methods, and discuss assorted subtleties and potential sources of error.

研究の動機と目的

• 常微分方程式系、特に制約を含む系におけるリャプノフ指数を信頼性高く計算するための数値技術の開発。
• 冗長な座標が独立した自由度を表さないといった、制約によって生じる複雑さの解消。
• 特にスケーリング条件下での、偏差ベクトル法とヤコビアンベース法の精度および安定性の比較。
• 相対論的回転する物体を含む複雑な力学系への応用を想定した、汎用的なC++実装の提供。
• SVDベースの手法におけるスケーリング時の軸ラベルの曖昧さを解消し、指数推定における系統的誤差を回避すること。

提案手法

• 初期条件の小さな球の楕円体の時間発展を用いて、主軸の成長率からすべてのn個のリャプノフ指数を計算する。
• ヤコビアン行列の特異値分解（SVD）を用いて、定期的な間隔で偏差ベクトルを直交化およびスケーリングする。
• 特にスケーリングが必要な状況において、SVDの代替として頑健な手法としてグラム・シュミット法を採用する。
• 元のODEおよびその変分方程式の数値積分を用いて、偏差ベクトルの時間発展を実装する。
• 主軸の順序を一貫して（s1 ≥ s2 ≥ ... ≥ sn）とることで、最大および最小の指数を正確に追跡する。
• スケーリングを避けることで、軸ラベルの誤りを回避し、指数推定の長期的整合性を保証する非スケーリングSVDを用いる。

実験結果

リサーチクエスチョン

• RQ1 座標数が真の自由度を上回る制約付き力学系において、リャプノフ指数をどのように正確に計算できるか？
• RQ2 ODEにおけるリャプノフ指数計算において、SVDとグラム・シュミット法の数値的利点と落とし穴は何か？
• RQ3 SVDベースの手法で偏差ベクトルをスケーリングすると、なぜ指数推定に系統的誤差が生じるのか？
• RQ4 さまざまな数値積分法および直交化スキームは、カオス的系において安定性および精度の観点でどのように比較できるか？
• RQ5 SVDのスケーリングを避けることで、軸ラベルの曖昧さを回避し、信頼性のあるリャプノフ指数推定が可能か？

主な発見

• SVD手法はスケーリングが適用されると、主軸のラベル付けの曖昧さにより失敗し、リャプノフ指数が系統的に過大評価される。
• 非スケーリングSVDは、グラム・シュミット法と良好に一致し、スケーリングを避けることでその正確性が裏付けられる。
• スケーリングにより、時間ステップ間で最も長い軸が誤ってラベル付けされ、成長率測定が歪められ、指数がバイアスを受ける。
• グラム・シュミット法は、一貫した直交化とラベル付けのおかげで、実用上より頑健であり、特にスケーリング条件下で優位である。
• テストケースとしてローレンツ系と強制減衰振動子を用い、非スケーリング条件下でSVDとグラム・シュミット法の結果が強く一致した。
• スピンするコンパクト物体がカー時空に存在するような複雑な系でも手法を検証し、実世界の相対論的力学に応用可能であることが確認された。