61-刚性问题认识v1.docVIP

1、 举例1 求解时出现方程组的各分量数量级差别很大，给数值求解带来困难，这类问题都称为刚性（Stiff）问题。 如某化学反应方程为： 方程右端系数矩阵为，特征值为，微分方程的精确解为 分析： 当t(∞时，u(1，v(t)(1，均为稳态解。 u与v中第二项为快变分量，大约计算到t=0.005秒时近似为0，称为时间常数；而为满变分量，其时间常数为，大约计算到秒时慢变分量接近0。 若使用4阶RK方法求解，根据绝对稳定要求，h-2.78/2000.5=0.00138，而从慢分量来看，需要20秒才能稳定，因此需要计算20/0.00138，约为14493步才能停止。 1.1 线性系统一般情形 给定方程： 假设A有m个特征根，以及相应的特征向量 令，则有：AQ=Qdiag(λ1, λ2,…λm)=QD； 令Q-1AQ=B，则A=QDQ-1； 令x(t)=Q-1y(t) 则成为： D为对角阵 由上可知：系统的解为： 为稳态解。 定义1：若线性系统中A的特征值满足：λj)0,j=1,…,m s