[论文解读] A Superexponentially Convergent Functional-Discrete Method for Solving the Cauchy Problem for Systems of Ordinary Differential Equations

Makarov Volodymyr, Denis Dragunov|arXiv (Cornell University)|Dec 30, 2010

Fractional Differential Equations Solutions参考文献 6被引用 3

一句话总结

本文提出了一种针对一阶常微分方程组（SODEs）柯西问题的超指数收敛功能离散（FD）方法。该方法通过基于阿多米安多项式和分段常数离散化网格的递归格式，将FD方法推广至SODEs，在系统非线性项和强迫项满足充分条件时，实现了超指数收敛。主要贡献在于给出了一个严格的收敛定理，并对误差衰减速率提供了定量估计。

ABSTRACT

In the paper a new numerical-analytical method for solving the Cauchy problem for systems of ordinary differential equations of special form is presented. The method is based on the idea of the FD-method for solving the operator equations of general form, which was proposed by V.L. Makarov. The sufficient conditions for the method converges with a superexponential convergence rate were obtained. We have generalized the known statement about the local properties of Adomian polynomials for scalar functions on the operator case. Using the numerical examples we make the comparison between the proposed method and the Adomian Decomposition Method.

研究动机与目标

• 开发一种用于求解一阶常微分方程组（SODEs）柯西问题的数值-解析方法，具备优异的收敛特性。
• 将此前仅应用于特征值与两点边值问题的功能离散（FD）方法，推广至SODEs的初值问题。
• 建立保证FD方法在SODEs上实现超指数收敛的充分条件。
• 将所提出的FD方法与阿多米安分解法（ADM）进行比较，尤其关注ADM不收敛的情形。

提出的方法

• 该方法将解近似为向量函数部分和的级数，每一项通过在分段网格上递归求解线性柯西问题获得。
• 初始问题（j=0）求解一个具有分段常数系数的线性ODE，其系数由非线性算子在前一网格点处的取值确定。
• 高阶项（j≥1）通过包含非线性算子阿多米安多项式的递推公式计算，该方法以微扰方式线性化系统。
• 离散化网格为均匀网格，步长为h，且在每个网格点处施加跳跃条件，以确保解在区间间的连续性。
• 该方法利用非线性算子N(t, u)的结构，并通过阿多米安多项式实现类似泰勒展开的方式，高效处理非线性项。
• 收敛性分析基于对非线性算子的雅可比矩阵和强迫项的有界性进行控制，从而导出保证误差超指数衰减的充分条件。

实验结果

研究问题

• RQ1 功能离散方法能否推广至含非线性项的一阶常微分方程组（SODEs），并保持其超指数收敛速率？
• RQ2 非线性算子和强迫项需满足何种充分条件，才能保证FD方法在SODEs上的超指数收敛？
• RQ3 在收敛行为方面，FD方法与阿多米安分解法（ADM）相比如何，特别是在ADM发散的情况下？
• RQ4 FD方法是否可应用于ADM因不稳定或发散而失效的问题？

主要发现

• 在满足涉及雅可比矩阵谱性质和强迫项有界性的充分条件下，SODEs的FD方法可实现超指数收敛。
• 即使在ADM发散的情况下，该方法仍能实现超指数收敛，数值实验在具有振荡解的非线性系统中已验证此结果。
• 对于测试问题，当阶数p=4时，FD方法在区间[0,6]上的绝对误差低于2×10⁻⁶，且随着阶数增加，误差呈指数下降。
• 理论误差界通过雅可比矩阵对称部分的最小特征值和强迫项的L∞-范数推导得出，得到α≈0.47和κ≈2.9。
• 数值结果证实，FD方法的收敛速率显著快于ADM，后者在同一问题上于[0,2]区间表现出发散。
• 通过FD方法高精度恢复了精确解u₁(t)=sin(t)，u₂(t)=cos(t)，验证了理论收敛性结论。