数值分析：一阶常微分方程初值问题数值解法.pptx

第七章常微分方程数值解;近似方法有两类:

一类是近似解析方法.如级数解法,逐次逼近法等;另一类就是数值解法,它可以给出解在一些离散点上的近似值.;定理7.1设f(x,y)在D={(x,y)|a?x?b,y?R}

上有定义且连续,同时满足如下茨(Lipschitz)条件;求解之前,应检验是否满足Lipschitz条件.;要计算出数值解,实际上就是求对应于一系列已知节点

a=x0x1…xn=b处的函数值y0,y1,…,yn;7.2(Euler)方法;2.数值微分

由导数的定义知:当h充分小时有:;3数值积分

在[xn,xn+1]上对积分得;在公式(7.5)中分别用yn,yn+1近似代替y(xn),

y(xn+1)得:;例1（p173）：用法求解初值问题;7.2.2隐式公式及改进的方法;利用数值积分的梯形公式得;通常采用的算法:;改进的法也可以写成:;设;定义：假设;公式的局部截断误差分析

设yn=y(xn),则有

此时公式可写成;2.改进的公式的局部截断误差分析;(*)-(**)得到改进的公式的局部截断误差为;7.3-(Runge-Kutta)方法;其中:;以2阶方法为例;而y(xn+1)在xn处的展开式为;这是一个有三个方程四个未知数的方程组.