《建模方法教学资料》微分方程数值解法.pptVIP

第九节 常微分方程的数值解法 一阶常微分方程的初值问题： 节点：x1x2 xn 步长 为常数 一 欧拉方法（折线法） yi+1=yi+hf(xi,yi)(i=0,1, , n 1) 优点：计算简单。 缺点：精度不高。 二 改进的欧拉方法 三 龙格—库塔法(Runge-Kutta) 欧拉公式可改写为 （它每一步计算f(xi,yi) 一次，截断误差为O(h2)) 改进的欧拉公式可改写为 （它每一步计算f(x,y)两次，截断误差为O(h3)) 标准四阶龙格—库塔公式 （每一步计算f(x,y)四次，截断误差为O(h5)) 例 分别用改进的欧拉格式和四阶龙格—库塔格式解初值问题（取步长h=0.2）： 表7—4 节点 改进欧拉法 四阶龙格—库塔法 准确解 0 1 1 1 0.2 1.186667 1.183229 1.183216 0.4 1.348312 1.3