常微分方程的数值解法分析.docVIP

第八章 常微分方程的数值解法 一．内容要点 考虑一阶常微分方程初值问题: 微分方程的数值解：设微分方程的解y(x)的存在区间是[a,b]，在[a,b]内取一系列节点a= x0 x1… xn =b,其中hk=xk+1-xk ；(一般采用等距节点，h=(b-a)/n称为步长)。在每个节点xk求解函数y(x)的近似值:yk≈y(xk)，这样y0 , y1 ,...,yn称为微分方程的数值解。 用数值方法，求得f(xk)的近似值yk，再用插值或拟合方法就求得y(x)的近似函数。 (一)常微分方程处置问题解得存在唯一性定理 对于常微分方程初值问题:如果： (1) 在的矩形内是一个二元连续函数。 (2) 对于y满足利普希茨条件，即 则在上方程的解存在且唯一，这里C=min((A-x0),x0+B/L),L是利普希茨常数。 定义：任何一个一步方法可以写为，其中称为算法的增量函数。 收敛性定理：若一步方法满足： (1)是p解的. (2) 增量函数对于y满足利普希茨条件. (3) 初始值y0是精确的。则kh=x-x0,也就是有 (一)、主要算法 1．局部截断误差 局部截断误差：当y(xk)是精确解时，由y(xk)按照数值方法计算出来的的误差y(xk+1)- 称为局部截断误差。 注意：yk+1和的区别。因而局部截断误差与误差ek+1=y(xk+1) -yk+1不同。 如果局部截断误差是O(