[論文レビュー] Computed Chaos or Numerical Errors

Lun-Shin Yao|arXiv (Cornell University)|Jun 21, 2005

Scientific Research and Discoveries参考文献 15被引用数 4

ひとこと要約

この論文は、有限時間刻みを用いた非線形微分方程式の数値的解法において、真の解を表すのではなく、数値誤差の拡大によって生じる解が得られる可能性があると主張している。ローレンツ方程式を事例として用い、有限時間刻みが安定・不安定多様体によって形成される仮想的分離界面を越えて軌道がジャンプすることを示し、解の一意性が破れ、誤差が爆発的に増幅されることで、計算されたカオスが真のカオス的力学の信頼できる近似であるかどうかに疑問を呈する。

ABSTRACT

Discrete numerical methods with finite time-steps represent a practical technique to solve initial-value problems involving nonlinear differential equations. These methods seem particularly useful to the study of chaos since no analytical chaotic solution is currently available. Using the well-known Lorenz equations as an example, it is demonstrated that numerically computed results and their associated statistical properties are time-step dependent. There are two reasons for this behavior. First, chaotic differential equations are unstable so that any small error is amplified exponentially near an unstable manifold. The more serious and lesser-known reason is that stable and unstable manifolds of singular points associated with differential equations can form virtual separatrices. The existence of a virtual separatrix presents the possibility of a computed trajectory actually jumping through it due to the finite time-steps of discrete numerical methods. Such behavior violates the uniqueness theory of differential equations and amplifies the numerical errors explosively. These reasons imply that, even if computed results are bounded, their independence on time-step should be established before accepting them as useful numerical approximations to the true solution of the differential equations. However, due to these exponential and explosive amplifications of numerical errors, no computed chaotic solutions of differential equations independent of integration-time step have been found. Thus, reports of computed non-periodic solutions of chaotic differential equations are simply consequences of unstably amplified truncation errors, and are not approximate solutions of the associated differential equations.

研究の動機と目的

• 非線形微分方程式の数値的計算によるカオス的解が、真の解の真正な近似であるのか、それとも数値誤差に起因する人工的結果であるのかを調査すること。
• 離散的数値解法における有限時間刻みの役割と、カオス系における軌道の安定性に与える影響を検討すること。
• 特異点の多様体によって形成される仮想的分離界面（仮想的分離界面）の形成とその影響を分析すること。
• カオス系における数値誤差が指数的かつ爆発的に増幅され、計算結果の信頼性を損なう仕組みを示すこと。
• 時間刻みに依存しない計算されたカオス的解が、まだ発見されていないことから、それらが真の力学的現象ではなく、数値的アーティファクトであると主張すること。

提案手法

• カオス系の代表例としてローレンツ方程式を用い、異なる時間刻みでの数値的解法を検討する。
• 不安定多様体および特異点付近の数値的軌道の挙動を分析し、誤差増幅のメカニズムを同定する。
• 平衡点の安定・不安定多様体が形成する仮想的分離界面（物理的境界ではないが、有限時間刻みの離散化によって越えられる）を特定する。
• 異なる時間刻みを用いた数値積分を実施し、計算された軌道および統計的性質が刻み幅に依存するかを検証する。
• 初期値問題に対する解の一意性の破れが、仮想的分離界面を越える軌道のジャンプによって生じ、常微分方程式の理論的性質に矛盾することを評価する。
• 計算された軌道の統計的分析を用い、時間刻み依存性を示し、真の解への収束が見られないことを評価する。

実験結果

リサーチクエスチョン

• RQ1 数値積分における有限時間刻みが、カオス的微分方程式の解をどの程度歪めているのか。
• RQ2 数値的軌道は、特異点の安定・不安定多様体によって形成される仮想的分離界面を越えることができるのか。その結果は何か。
• RQ3 計算されたカオス的解が有界であるにもかかわらず、なぜ依然として積分時間刻みに依存しているのか。
• RQ4 カオス系における数値誤差がどのように爆発的に増幅され、計算されたカオスの妥当性にどのような含みを持つのか。
• RQ5 カオス系において報告される非周期的解は、元の微分方程式の真の解であるのか、それとも切り捨て誤差に起因するアーティファクトにすぎないのか。

主な発見

• ローレンツ方程式の数値的カオス的解は時間刻みに依存しており、真の解への収束的近似ではないことが示された。
• 特異点の安定・不安定多様体によって形成される仮想的分離界面の存在が、有限時間刻みによる離散化の結果、数値的軌道がそれらを越えてジャンプすることを可能にしている。
• このようなジャンプは、初期値問題に対する一意性定理に反するため、数値的解に根本的な不整合をもたらす。
• カオス系における数値誤差は不安定多様体付近で指数的に増幅され、仮想的分離界面を越えると爆発的に増幅される。
• 時間刻みに依存しない計算されたカオス的解は、まだ発見されておらず、すべての解が真の解ではなく、数値的アーティファクトである可能性が高い。
• したがって、カオス系における非周期的解の報告は、不安定に増幅された切り捨て誤差に起因する可能性が高く、真のカオス的力学とは無関係である。