第5讲(实践课续)格龙-库塔方法.docVIP

第四讲（续） 3.3 单步法的收敛性与绝对稳定性.3.1 单步法的收敛性定义.1 设y(x)是初值问题(1.2.1)的精确解，是单步法(3.2.2)在处产生的近似解，若 则称方法(3.2.2)产生的数值解收敛于. 实际上，定义中是一固定点，当h→0时n→∞，n不是固定的.因显然方法收敛，则在固定点处的整体误差，当p≥1时. 下面定理给出方法(3.2.2)收敛的条件定理.1　设初值问题(.2.1)的单步法(3.2.2)是p阶方法(p≥1)，且函数对y满足Lipschitz条件，即存在常数L＞0使对，均有　　　　 则方法(3.2.2)收敛，且. 定理证明略..3.2 绝对稳定性用单步法(3.2.2)求数值解，由于原始数据及计算过程舍入误差影响，实际得到的不是而是，其中是误差，再计算下一步得到 以Euler法为例，若令,则　　　　(.3.1) 如果，则从计算到误差不增长，它是稳定的但如果条件不满足就不稳定例.4 y′=-100y,y(0)=1,精确解为，用Euler法求解得 若取h=0.025，则,当，而，显然计算是不稳定的.如果用后退Euler法(3.1.5)解此例，仍取h=0.025,则　　　　　　,即 显然当，计算是稳定的.　　由此看到稳定性与方法有关，也与有关，在此例中.在研究方法的稳定性时，通常不必对一般的f(x,y)进行讨论，而只针对模型方程　　　　　　　　　　　　　　(