[论文解读] 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

一句话总结

本文提出了一种切比雪夫谱配置方法，用于求解具有严格后验误差估计的线性周期性常微分方程（ODEs）和时滞微分方程（DDEs）。该方法引入了一项新颖的希尔伯特空间上算子的后验特征值扰动定理（Theorem 11），从而能够对单人迹算子的谱半径提供验证性界，实现了对具有整数时滞的DDEs的数学严格稳定性认证。

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.

研究动机与目标

• 开发一种可靠且经过验证的数值方法，用于计算具有整数时滞的线性周期性DDEs的解及其稳定性特性。
• 利用谱配置方法，为DDEs初值问题的解提供后验误差估计。
• 推导单人迹算子的真实特征值与其矩阵近似特征值之间距离的后验界。
• 将该方法推广至DDEs系统，并在时滞阻尼马修方程上展示其有效性。
• 解决谱方法中多项式求值、$L^\infty$ 范数估计以及条件数控制等方面的技术挑战。

提出的方法

• 在高斯-洛巴托点上使用切比雪夫谱配置对DDEs进行离散化，以高阶精度逼近基本解。
• 应用后验误差估计技术，以界定初值问题中残差的$L^\infty$-范数。
• 通过在子区间上使用谱配置求解一系列ODEs，并在配置点处插值解，构造单人迹矩阵近似。
• 将Bauer-Fike定理由矩阵推广至希尔伯特空间上的算子（定理11），以界定单人迹算子特征值扰动的界。
• 利用特征向量矩阵条件数的后验界，控制离散近似中特征值的精度。
• 在时滞马修方程上验证该方法，结果表明误差界随$N$增加呈指数衰减。

实验结果

研究问题

• RQ1 如何为具有整数时滞的线性DDEs的切比雪夫谱配置解，严格推导出后验误差估计？
• RQ2 单人迹算子的特征值与其矩阵近似特征值之间的关系是什么？该关系如何实现后验界？
• RQ3 能否利用谱配置方法，以数学上严格且可计算的误差界验证单人迹算子的谱性质？
• RQ4 单人迹矩阵近似特征值的精度和条件数如何随配置点数$N$变化？
• RQ5 该方法在多大程度上可推广至非整数时滞、中立型DDEs或泛函微分方程？

主要发现

• DDE初值问题的后验误差界随配置点数$N$的增加呈指数衰减，数值实验已验证此现象。
• 对于参数集$(a,b)=(-1.1,1)$，单人迹矩阵的最大特征值为$0.9369$，其与真实特征值距离的后验界为$r=0.0434$，由于$|\mu_1| + r < 1$，因此证明了系统的稳定性。
• 计算特征值周围误差圆盘的半径随$N$大致呈指数衰减，当$N \gtrapprox 220$时，浮点精度限制约为$10^{-5}$。
• 该方法通过逐像素计算单人迹矩阵特征值及其验证性误差界，成功生成了方程(1)的稳定性图。
• 特征值扰动定理（定理14）为单人迹算子的真实特征值与矩阵近似最近计算特征值之间的距离，提供了严格且可计算的界。
• 该方法可推广至DDEs系统，如在时滞阻尼马修方程上的应用所示，其收敛性和验证性特性保持一致。