81欧接法.pptVIP

步长与误差 离散误差 设y=y(t)是初值问题的唯一解, 是离散近似解集 →从tk到tk+1步的计算误差 注 : ,每步忽略 那么在末端进行M步之后,误差积累到 欧拉方法的精度 设y=y(t)是初值问题 的解,若 是欧拉方法计算的近似值序列,则 注: 区间末端的误差称为最终全局误差. 实现 通过计算 求[a,b]上的初值问题 的近似解. Matlab Code f=inline(‘x’); a=0;b=10; Ya=1;M=20; h=(b-a)/M; T=zeros(1,M+1); Y=zeros(1,M+1); T=a:h:b; Y(1)=ya; for k=1:M Y(k+1)=Y(k)+h*f(T(k),Y(k)); end E=[T’ Y’]; Plot(E) 欧拉公式及其改进 tn tn+1 欧拉公式 欧拉公式的改进办法 后退的欧拉公式