[論文レビュー] A small variation of the Taylor Method and periodic solutions of the 3-body problem

Óscar Perdomo|arXiv (Cornell University)|Jul 4, 2015

Spacecraft Dynamics and Control参考文献 4被引用数 3

ひとこと要約

本稿では、正確な解が不要なまま、全誤差および丸め誤差を追跡する修正されたテイラー法を提案し、力学系における厳密な数値的証明を可能にする。この手法を用いて、質量が100, 100, 200で周期T ≈ 96.241の新しい周期的解が三体問題に存在することを、正当化された数値積分と区間演算を用いて確認した。

ABSTRACT

In this paper we define a small variation of the Taylor method and a formula for the global error of this new numerical method that allows us to keep track of the round-off error and does not require previous knowledge of the exact solution. As an application we provide a rigorous proof of the construction/existence of a periodic solution of the three body problem. Some images of this periodic motion can be seen at https://www.youtube.com/watch?v=fSmQyeKcj5k

研究の動機と目的

• 力学系における常微分方程式を数値的に安定かつ厳密に解くための手法を開発すること。特に、正確な解が不明な場合を想定する。
• 特に三体問題のようなカオス的系において、長時間の数値積分における丸め誤差の追跡の課題に対処すること。
• 数値的に観測されたが未証明の解に基づき、三体問題における新しい周期的解の存在を正当化された数値的証明で示すこと。
• 記号的有理数演算に起因する計算コストの指数的増加を回避し、高精度な数値積分にコンピュータ代数システムを活用できるようにすること。
• 数値的に観測された周期的解が、新しい誤差境界式を用いることで、実際に数学的に周期的であることを厳密に証明できることを示すこと。

提案手法

• 各ステップで数学的演算を用いて高精度近似を計算し、その後に幅Hの有理数区間で結果をバウンディングする、テイラー法の微小な変種を提案する。
• 区間演算を用いて各反復の結果をユーザーが定義した幅Hの有理数区間内でバウンディングし、記号的正確な算術を必要とせずに誤差制御を実現する。
• 理論的貢献の中心は定理2.3であり、解の境界の事前知識を必要とせず、ステップサイズhと丸め許容誤差Hの関数として全誤差バウンディングを提供する。
• コンピュータ代数システムが有理数および平方根を含む式に対してタイトなバウンディングを計算できる能力を活用し、高精度な区間演算を可能にする。
• テイラー級数積分と正当化された数値的手法を組み合わせ、区間包摂を用いて切り捨て誤差と丸め誤差の両方を追跡する。
• 三体問題にこの手法を適用し、誤差バウンディングを制御した運動方程式の積分により、時間Tにおける状態が初期状態の小さな区間内に戻ることを確認した。

実験結果

リサーチクエスチョン

• RQ1 正確な解が不要な修正テイラー法を設計し、全誤差および丸め誤差を厳密に追跡できるか？
• RQ2 この手法を用いて、三体問題における数値的に観測された周期的解が実際に数学的に周期的であることを検証できるか？
• RQ3 カオス的系における周期性の厳密な数値的妥当性を達成するための最小限の計算コストは何か？
• RQ4 区間演算をテイラー級数積分と効果的に組み合わせ、式サイズの指数的増加を防ぎながら精度を維持できるか？
• RQ5 数値積分に誤差バウンディングを適用するだけで、三体問題における周期的解の存在を証明できるか？

主な発見

• 質量が100, 100, 200で周期T ≈ 96.241の三体問題における新しい周期的解が、厳密に証明された。
• 解は1周期後に初期状態と10^-14以内に一致することが数値的に観測され、本手法を用いてこれが正確であると正当化された。
• 解の境界の事前知識を必要とせず、ステップサイズhと丸め許容誤差Hのみを用いて全誤差を適切にバウンディングできた。
• 従来の記号的有理数算術に伴う式サイズの指数的増加を回避し、長時間積分を実行可能にした。
• 中間結果を幅が制御された有理数区間でバウンディングすることで、標準的なコンピュータ代数システムを用いた正当化された積分が可能になった。
• 証明により、図1.1に示された数値的に観測された解が実際に周期的であることが確認され、天体力学における数値解の正当化において重要な一歩を踏み出した。