一种改进型刚性常微分方程数值解法——T-R预估-校正法

摘要

Treanor显式方法是一种求解刚性常微分方程的数值积分法,用它来解刚性方程时仍然存在一些问题。本文对Treanor方法的稳定性进行了分析和讨论,以Treanor方法为基础,采用Richardson外推技巧,构造出一种改进型方法—T-R预估-校正法。分析和计算都表明,与Treanor方法相比,改进的方法具有更大的稳定域,且提高了数值解的精度。它更适合于刚性方程的求解,是数字仿真中一种较好的数值解法。