[論文レビュー] On the fractal dimension of the Duffing attractor

Mariusz Tarnopolski|arXiv (Cornell University)|May 29, 2013

Mathematical Dynamics and Fractals参考文献 9被引用数 5

ひとこと要約

本研究では、ボックスカウント法と相関次元法を用いて、数値的にDuffingアトラクターのフラクタル次元を推定し、サンプル点数の増加に伴い安定値への収束を示した。結果として、ボックスカウント次元は dC ≈ 1.43、相関次元は dG ≈ 1.38 となり、位相空間内に非整数のフラクタル構造を有する奇妙な自己相似アトラクターであることが示された。

ABSTRACT

The box counting dimension $d_C$ and the correlation dimension $d_G$ change with the number of numerically generated points forming the attractor. At a sufficiently large number of points the fractal dimension tends to a finite value. The obtained values are $d_C\approx 1.43$ and $d_G\approx 1.38$.

研究の動機と目的

• 数値的に生成された点数の増加に伴うDuffingアトラクターのフラクタル次元推定値の収束行動を調査すること。
• カオス的状態におけるDuffingアトラクターのボックスカウント次元（dC）および相関次元（dG）の安定値を特定すること。
• 非線形力学分野で広く用いられる2つの手法を用いて次元を推定することで、Duffingアトラクターのフラクタル性を検証すること。
• 有限な点集合に対して、相関次元アルゴリズムの数値的複雑性と信頼性をボックスカウント法と比較して評価すること。
• 位相空間内でのアトラクターの幾何的複雑性の定量的測定を提供し、それが奇妙なアトラクターとしての分類を支援すること。

提案手法

• Mathematicaのデフォルト手法を用いて、Duffing方程式（¨x + δẋ + βx + αx³ = γ cos(ωDt)）の数値的積分を実行し、軌道を生成した。
• 位置、速度、駆動の位相を状態変数として用いて、2階の常微分方程式を3次元の非線形自己同調系に変換した。
• 駆動周期ごとにサンプリングするPoincaré断面を適用し、アトラクターの構造を抽出した。
• MathematicaのBoxCount関数を用いてボックスカウント次元（dC）を計算し、εが小さい範囲で ln N(ε) と ln(1/ε) の関係に線形回帰を適用した。
• 相関次元（dG）は、相関和 C(r) = (1/N²)Σᵢ<ⱼ H(r − ||xi − xj||) を実装した並列Pythonプログラムにより算出。rは小さな値の範囲で変化させた。
• log C(r) と log r、および log N(ε) と log(1/ε) の線形回帰を実施し、傾きを推定することでフラクタル次元を算出し、誤差最小化を適用した。

実験結果

リサーチクエスチョン

• RQ1 サンプル点数の増加に伴い、Duffingアトラクターのボックスカウント次元（dC）および相関次元（dG）はどのように変化するか？
• RQ2 十分に大きな点集合に対して、フラクタル次元推定値が有限の値に収束するか。これは真のフラクタル構造を示唆するか？
• RQ3 与えられたカオス的パrameter設定下で、Duffingアトラクターのボックスカウント次元 dC の収束値は何か？
• RQ4 Duffingアトラクターの相関次元 dG の収束値は何か。また、dC と比較してどのように異なるか？
• RQ5 相関次元アルゴリズムの数値的複雑性は点数の増加に伴いどのように変化するか。大規模データセットに対しても計算上の実行可能性は保証されるか？

主な発見

• 点数の増加に伴い、ボックスカウント次元 dC は約 1.43 に収束し、1から2の間のフラクタル構造を示していることがわかった。
• 相関次元 dG は約 1.38 に収束し、非整数かつ自己相似的な幾何構造を有する奇妙なアトラクターの存在を裏付けた。
• N ≥ 1592点の場合、log C(r) と log r、および log N(ε) と log(1/ε) の線形回帰により、dG ≈ 1.357 および dC ≈ 1.330 の安定的かつ一貫性のある推定値が得られた。
• 1つの点を追加した場合、相関次元推定値の標準偏差は21%増加したが、次元値自体の変化はわずか4%にとどまり、点数の変化に対して感度が高いことが示された。
• 相関次元アルゴリズムの計算複雑性は O(n²) であり、実行時間に二次多項式をフィッティングすることで確認された。31,838点の場合、予測時間は12,943秒であり、実際の時間12,610秒と非常に近い値であった。
• 両手法とも有限で非整数の値に収束し、Duffingアトラクターのフラクタル性を確認した。これは、散逸的カオス系における奇妙なアトラクターとしての分類を支援するものである。