[論文レビュー] Chebyshev collocation for linear, periodic ordinary and delay differential equations: a posteriori estimates

Ed Bueler|arXiv (Cornell University)|Sep 23, 2004

Numerical methods for differential equationsMathematics参考文献 28被引用数 18

ひとこと要約

本稿では、整数遅れを伴う線形周期的常微分方程式（ODE）および遅れ微分方程式（DDE）を解くためのチェビシェフコロケーション法を提示し、厳密な事後誤差推定値を提供する。ヒルバート空間上の作用素に対する新規な事後固有値摂動定理を導入することで、モノドロミー作用素のスペクトル半径に対する検証可能な上限が得られ、整数遅れを伴うDDEの数学的に厳密な安定性認証が可能となる。

ABSTRACT

We present a Chebyshev collocation method for linear ODE and DDE problems. We first give a posteriori estimates for the accuracy of the approximate solution of a scalar ODE initial value problem. Examples of the success of the estimate are given. For linear, periodic DDEs with integer delays we define and discuss the monodromy operator U as our main goal is reliable estimation of the stability of such DDEs. We prove a theorem which gives a posteriori estimates for eigenvalues of U, our main result. This result is based on a generalization to operators on Hilbert spaces of the Bauer-Fike theorem for (matrix) eigenvalue perturbation problems. We generalize these results to systems of DDEs. A delayed, damped Mathieu equation example is given. The computation of good bounds on ODE fundamental solutions is an important technical issue; an a posteriori method for such bounds is given. Certain technical issues are also addressed, namely the evaluation of polynomials and the estimation of L^\infty norms of analytic functions. Generalization to the non-integer delays case is also considered.

研究の動機と目的

• 整数遅れを伴う線形周期的DDEの解と安定性特性を信頼性があり検証可能な数値的手法で計算すること。
• スペクトルコロケーションを用いてDDEの初期値問題の解に対する事後誤差推定値を提供すること。
• モノドロミー作用素の固有値とその行列近似の固有値との間の距離に対する事後誤差境界を導出すること。
• この手法をDDEの連立系に拡張し、遅れ付き減衰型マチエ方程式への応用を示すこと。
• スペクトル法における多項式評価、$L^\infty$ノルム推定、条件数制御の技術的課題に対処すること。

提案手法

• ガウス＝ロバッチョフ点上でチェビシェフコロケーションを用いてDDEを離散化し、高次の精度で基本解を近似する。
• 初期値問題における残差の$L^\infty$ノルムをバインドするための事後誤差推定技術を適用する。
• スペクトルコロケーションを用いて部分区間上で逐次的ODEを解き、コロケーション点における解を補間することでモノドロミー行列近似を構築する。
• ヒルバート空間上の作用素に対するバウアー＝フィックの定理の一般化（定理11）を用いて、モノドロミー作用素の固有値摂動をバインドする。
• 固有ベクトル行列の条件数に対する事後誤差境界を用いて、離散近似における固有値の精度を制御する。
• 遅れ付きマチエ方程式に対して手法を検証し、$N$の増加に伴い誤差境界が指数関数的に収束することを示した。

実験結果

リサーチクエスチョン

• RQ1 整数遅れを伴う線形DDEのチェビシェフコロケーション解に対して、どのように厳密な事後誤差推定値を導出できるか？
• RQ2 モノドロミー作用素の固有値とその行列近似の固有値との間の関係は何か？ そして、これを事後にバインドすることは可能か？
• RQ3 スペクトルコロケーションを用いて、モノドロミー作用素のスペクトル的性質を数学的に厳密な誤差境界で検証できるか？
• RQ4 コロケーション点数$N$の増加に伴い、モノドロミー行列近似の精度と条件数はどのように変化するか？
• RQ5 この手法は、非整数遅れ、中立型DDE、または関数微分方程式へまで一般化可能か？

主な発見

• DDE初期値問題に対する事後誤差境界は、数値実験によりコロケーション点数$N$の増加に伴い指数関数的に減少することが示された。
• パrameterセット$(a,b)=(-1.1,1)$に対して、モノドロミー行列の最大固有値は$0.9369$であり、真の固有値との距離に対する事後誤差境界は$r=0.0434$であった。$|\mu_1| + r < 1$を満たすため、安定性が証明された。
• 計算された固有値の周囲の誤差円の半径は、$N$の増加に伴い概ね指数関数的に減少し、$N \gtrapprox 220$で浮動小数点精度が約$10^{-5}$程度に制限される。
• 本手法により、方程式(1)の安定性図が、ピクセル単位でモノドロミー行列固有値とその検証可能な誤差境界を計算することで正しく生成された。
• 固有値摂動定理（定理14）により、モノドロミー作用素の真の固有値とその行列近似の最も近い計算固有値との距離に対する、数学的に厳密で計算可能な境界が得られた。
• 遅れ付き減衰型マチエ方程式への応用例から、連立DDEへへの一般化が可能であることが示され、同様の収束性と検証性を示した。