刚性方程的数值求解法(二)

当关联系统矩阵A具有非常分散的特征值时,就会出现刚性性质。在这种情况下,要想找一个稳定的数值求解微分方程的方法是非常困难的。迄今尚无一种方法能合理地求解刚性方程。本文提出了一种特殊的方法,即仅对矩阵做一些乘方运算就能完成大量步数的数值计算。此法将能用足够微小的步长去保证计算方法的稳定与精度,因此能够解决线性刚性方程的求解问题。