[論文レビュー] Approximation of solutions of DDEs under nonstandard assumptions via Euler scheme

Natalia Czyżewska, Paweł Morkisz+1 more|arXiv (Cornell University)|Jun 7, 2021

Numerical methods for differential equations参考文献 22被引用数 2

ひとこと要約

本稿では、右辺関数の正則性が非標準的である場合の遅れ微分方程式（DDEs）に対するオイラースキームの厳密な誤差解析を提示する。具体的には、グローバルに線形成長、グローバル片側リプシッツ連続性、局所ホルダー連続性という条件を想定する。主な貢献は、ホルダー指数に明示的に依存する点誤差の理論的上界を導出したことである。この上界は、非リプシッツ非線形性を有する多次元DDEに対する数値実験でも検証されている。

ABSTRACT

We deal with approximation of solutions of delay differential equations (DDEs) via the classical Euler algorithm. We investigate the pointwise error of the Euler scheme under nonstandard assumptions imposed on the right-hand side function $f$. Namely, we assume that $f$ is globally of at most linear growth, satisfies globally one-side Lipschitz condition but it is only locally H\"older continuous. We provide a detailed error analysis of the Euler algorithm under such nonstandard regularity conditions. Moreover, we report results of numerical experiments.

研究の動機と目的

• 右辺関数 f がグローバルリプシッツでないが、より弱い正則性条件を満たす場合の古典的オイラースキームの収束性を分析すること。
• グローバル線形成長、片側リプシッツ条件、局所ホルダー連続性という非標準的仮定の下で、オイラー法の厳密な点誤差境界を確立すること。
• 通常はグローバルリプシッツ連続性を仮定する既存の誤差解析フレームワークを、物理的・工学的モデルに一般的に見られるより現実的で滑らかでない設定へと拡張すること。
• 非リプシッツ非線形性を有する代表的なDDEに対して、理論的結果を数値実験で検証すること。
• 類似した弱い正則性仮定の下で、高次精度数値スキームへの解析の拡張に基盤を提供すること。

提案手法

• 時間区間 [jτ, (j+1)τ] を順次処理する際、均等なステップサイズ h = τ/N を用いてオイラースキームを繰り返し適用し、直前の区間の解を次回の初期値として用いる。
• 理論的解析は、グローバルリプシッツ連続性に代えてホルダー連続性と片側リプシッツ条件を考慮した、グロンウォール型の帰納的誤差不等式に依拠する。
• 重要なステップとして、同様の弱い正則性仮定の下でODE解の挙動を解析する補助的補題の証明を行う。特に、ホルダー指数 β_i に対して、局所誤差が h^α および h^β_i の形で抑えられることを示す。
• 局所切捨て誤差の見積もりと誤差のグローバルな伝播を組み合わせることで誤差境界を導出する。片側リプシッツ条件を用いて誤差の増大を制御する。
• 本解析の新規性は、関数 f のホルダー指数 β1, ..., βp に依存する誤差の依存関係を明示的に追跡している点である。
• 非リプシッツ右辺関数を有する3つのテストDDEに対して数値実験を実施し、理論的収束率とホルダー指数が誤差に与える影響を確認した。

実験結果

リサーチクエスチョン

• RQ1 右辺関数 f がグローバルリプシッツでないが、局所ホルダー連続かつ片側リプシッツ条件を満たす場合、DDEに対するオイラースキームの収束速度はどのようになるか？
• RQ2 これらの非標準的仮定下で、オイラースキームの点誤差は関数 f のホルダー指数にどのように依存するか？
• RQ3 グローバルリプシッツ連続性が成り立たない状況下でも、線形成長と片側リプシッツ条件に依拠して、オイラー法の厳密な誤差境界を確立できるか？
• RQ4 数値実験は理論的誤差見積もり、特にホルダー指数への依存性をどの程度確認できるか？
• RQ5 提案された解析フレームワークは、同様の正則性条件を満たすDDEに対して高次精度数値スキームへと拡張可能か？

主な発見

• 本稿では、右辺関数 f のホルダー指数 β1, ..., βp に明示的に依存するオイラースキームの点誤差に対する理論的上界を確立した。
• 誤差境界は O(h^α + ∑h^β_i) の形をとり、α は時間方向のホルダー指数、β_i は空間方向のホルダー指数である。これは、ホルダー正則性が低下するにつれて収束速度が劣化することを示している。
• 片側リプシッツ定数 H+ がゼロの場合、誤差境界は ∆ + C(1 + ∥ξ∥)(1 + ∆)(b−a)(h^α + ∑h^β_i) の形を取り、有利な条件下で h に関して線形収束を示す。
• H+ > 0 の場合、誤差境界には指数関数的要因 e^{H+(b−a)} と H+ および問題パラメータに依存する乗法的定数が含まれ、片側リプシッツ定数の影響が反映されている。
• 数値実験により理論的収束率が確認され、ホルダー指数が小さくなるほど誤差が増大する傾向が示され、理論的正則性依存性が妥当であることが検証された。
• 複数のテストケースにわたり結果が安定しており、非リプシッツ非線形性を有する多次元DDEに対しても適用可能であることが示された。これは、材料科学や相転移モデルにおける現実的応用への適応可能性を示している。