[論文レビュー] Gradual approximation of the domain of attraction by gradual extension of the "embryo" of the transformed optimal Lyapunov function

Éva Kaslik, A.M. Balint+1 more|arXiv (Cornell University)|Dec 20, 2002

Quantum chaos and dynamical systems参考文献 6被引用数 5

ひとこと要約

本稿では、非線形解析的常微分方程式系における漸近安定な平衡点の吸引域（DA）を段階的に近似する、新規の反復的手法を提示する。この手法は、系列展開とコーシー＝アダマールの収束条件を用いて得られる変換された最適なリャプノフ関数の自然な解析接続性を活用し、関数の連続的な「胚」を計算することで、逐次的に吸引域の推定値を拡大する。その結果、ベンチマーク用多項式系において、従来手法を上回る精度を得た。

ABSTRACT

In this paper an autonomous analytical system of ordinary differential equations is considered. For an asymptotically stable steady state x0 of the system a gradual approximation of the domain of attraction DA is presented in the case when the matrix of the linearized system in x0 is diagonalizable. This technique is based on the gradual extension of the "embryo" of an analytic function of several complex variables. The analytic function is the transformed of a Lyapunov function whose natural domain of analyticity is the DA and which satisfies a linear non-homogeneous partial differential equation. The equation permits to establish an "embryo" of the transformed function and a first approximation of DA. The "embryo" is used for the determination of a new "embryo" and a new part of the DA. In this way, computing new "embryos" and new domains, the DA is gradually approximated. Numerical examples are given for polynomial systems. For systems considered recently in the literature the results are compared with those obtained with other methods.

研究の動機と目的

• 非線形解析的常微分方程式系における漸近安定な平衡点の吸引域（DA）を体系的に近似する手法を開発すること。
• 標準的なリャプノフ関数に基づくDA推定法の限界を克服すること。これらはしばしば過剰に保守的または不完全な近似をもたらす。
• 変換された最適なリャプノフ関数の自然な解析接続性を、DAの正確な特徴づけに活用すること。
• 計算的に実行可能で、各ステップで関数の「胚」を精緻化することで、段階的にDA推定値を拡大する反復的手順を提供すること。
• 既知のDAを有する多項式系に対して本手法を検証し、文献に登録された既存手法とその性能を比較すること。

提案手法

• 平衡点におけるヤコビ行列を対角化する同型写像 S を用いて元の系を変換し、新たな系 g = S⁻¹∘f∘S を得る。
• 最適なリャプノフ関数 V を W = V∘S に変換し、線形非同次偏微分方程式 ⟨∇W, g⟩ = -‖Sz‖² を満たし、W(0) = 0 を満たす。
• 定理3を用いて、次数2から始めて再帰的に W のテイラー級数の係数を計算し、高次項の係数は低次項に依存する。
• W の級数の収束領域 D⁰ はコーシー＝アダマールの公式により決定され、初期のDA推定値は DA⁰ = S(D⁰) として得られる。
• 新たな「胚」は、級数計算を高次項へと拡張することで生成され、収束領域 Dᵏ を逐次的に計算し、拡大された DA推定値 DAᵏ = S(Dᵏ) を得る。
• このプロセスは段階的に繰り返され、各反復で高次項の収束性を分析することで、関数の解析的構造を段階的に精緻化し、DA近似に新たな領域を追加する。

実験結果

リサーチクエスチョン

• RQ1 漸近安定な平衡点の吸引域は、変換された最適なリャプノフ関数の解析接続の反復的精緻化を通じて体系的に近似可能か？
• RQ2 線形非同次偏微分方程式の解の自然な解析接続性は、正確な吸引域を特徴づけるためにどのように利用可能か？
• RQ3 変換されたリャプノフ関数の「胚」を反復的に拡張することで、標準的手法に比べてDA推定の精度がどの程度向上するか？
• RQ4 多項式系（DAが既知または未知）に対して、この手法の計算的実行可能性と収束性はどのように評価できるか？
• RQ5 既存のDA推定技術と比較して、本手法はカバー範囲と精度の両面でどの程度優れているか？

主な発見

• 例1では、3ステップで真の円形吸引域と非常に近いDA推定値が得られ、境界が正確に捉えられた。
• 例2（3次元系）では、2ステップでDA推定値が著しく拡大し、真の双曲面境界をよく再現した。
• 例3（[10]からの2次元系）では、4ステップでの推定が[10]で報告された推定より広範囲をカバーし、優れたカバー範囲を示した。
• 例4（3次元系）では、2ステップ目で従来手法より広範なDA推定値が得られ、1ステップ目でも有意義な拡大が確認された。
• 表1から明らかなように、高次近似に伴い計算コストが著しく増加した。例2の2ステップ目では14.7時間かかっており、精度と計算時間のトレードオフが明確に示された。
• 反復的手順により、初期の収束領域を超えてDA推定値が段階的に拡大された。これは、関数の解析的構造を段階的に精緻化することで、真の吸引域を徐々に近似可能であることを確認した。