第六章节常微分方程数值解法幻灯片.pptVIP

第六章 常微分方程初值问题数值解法 已知 求 y(xi) (i=1,2,…,n)的近似值yi。 作为入门，我们只讨论n=1的情况。即 已知 求 y(xi) (i=1,2,…,n)的近似值yi。 求常微分方程初值问题数值解法的途径 假定 建立数值方法(求近似值方法)常 用的途径有： （1）用差商代替导数 （2）使用数值积分的方法 （3）使用泰勒公式 第一节 欧拉法（尤拉法） 一、欧拉法及其截断误差 欧拉法就是用差商代替导数求常微分方程初值问题的方法。它的几何意义是用一条初始点重合的折线来近似表示曲线。 y (x1,y1) (x2,y2) (x0,y0) y=y(x) 0 x0 x1 x2 x3 x 其公式为： 假定yi=y(xi),则 y(xi+1)=y(xi)+hf(xi,yi)+o(h2),所以欧拉公式的局部截断误差可表示为o(h2)。 用欧拉方法求初值问题举例 例 取步长h=0.1用欧拉方法求初值问题： y=x+y y(0)=1 的数值解。 解 有 f(x,y)=x+y , x0=0 , y0=1 y1=y0+hf(0,1) =1+0.1(0+1) =1.10000 y2=y1+hf(0.1,y1) =1.22000 …… 二、改进的欧拉法及预测-校正公式 1. 改进的欧拉法 改进的欧拉法即用公式 计算 y(xi) (i=1,2,…,n)的近似值yi的方法。 由于该公式是一个隐式，实际计算时，可将欧拉法和该方法联合使用，构成如下的计算公式 进行 yi+1 的迭代计算，直到 并把 2. 预测-校正公式 在步长h较小的情况下，用欧拉公式算出 后，通常只进行一次迭代就能满足精度要求。只对yi+1 进行一次迭代计算的做法称为预测-校正法，即 可以证明，预测-校正法的局部截断误差为o(h3)。 所以，y1=1.24 y2(p)=y1+hf(x1,y1)=1.24+0.2(0.2+1.24)=1.528 y2(c)=y1+ [f(x1,y1)+f(x2,y2(p))]=1.24+0.2/2(0.2+1.24+0.4+1.528)=1.577 所以，y2=1.577 类似地, 可计算y3等。 改进的欧拉法求初值问题举例 例 取步长h=0.2,用改进的欧拉法(预测-校正公式)求初值问题 y=x+y y(0)=1 的数值解。 解 f(x,y)=x+y , x0=0 , y0=1, 用欧拉公式计算y1的初始值y1(p) y1(p)=y0+hf(x0,y0)=1+0.2(0+1)=1.2 然后计算y1(c) y1(c) =y0+ [f(x0,y0)+f