第7章节常微分方程数值解法.pptVIP

略去高阶项，再做变换即可得到的形式 u?=λu. 对于m个方程的方程组, 可线性化为y?=Ay, 这里A为m×m雅可比矩阵(?fi/?yj)，若A有m个特征值λ1,λ2,?,λm，其中λi可能是复数，所以，为了使模型方程结果能推广到方程组，方程(4.8)中λ为复数. 为保证微分方程本身的稳定性，还应假定Re(λ)0. 　　下面先研究欧拉方法的稳定性. 模型方程y?=λy的欧拉公式为 设在节点yn上有一扰动值εn，它的传播使节点值yn+1产生大小为的扰动值εn+1，假设用yn*=yn+εn，按欧拉公式得出yn+1*=yn+1+εn+1的计算过程不再有新的误差，则扰动值满足 可见扰动值满足原来的差分方程(4.9). 这样，如果差分方程的解是不增长的，即有 则它就是稳定的. 这一论断对于下面将要研究的其它方法同样适用. 　　显然，为要保证差分方程(4.9)的解是不增长的，只要选取h充分小，使 　　在μ=hλ的复平面上，这是以(-1,0)为圆心，1为半径的单位圆. 称为欧拉法的绝对稳定域，一般情形可由下面定义. 定义6 单步法(4.1)用于解模型方程y?=λy，若得到的解yn+1=E(hλ)yn，满足|E(hλ)|1，则称方法(4.1)是绝对稳定的. 在μ=hλ的平面上, 使|E(hλ)|1的变量围成的区域，称为绝对稳定区域，它与实轴的交称为绝对稳定区间.