_关于单步法于多步法部分方法的分析.docVIP

摘要: 通过《数值分析》这门课程我们学习了求解常微分方程初值问题的方法分为单步法和多步法,单步法主要有欧拉方法。改进的欧拉方法和龙格法,多步法主要有Adams法等方法,关于一阶常微分方程初值问题求解用数值方法求的是的近似值.所以本文主要用改进的欧拉方法和Adams法分别作为单步法和多步法的为代表进行误差分析等比较. 方法介绍： 一．单步法 1、欧拉公式： 1、方法构造的思想：微分方程初值问题 再用前项差代替一阶导数,则 此时x0, y0均已知y(x0)= y0,则由此式可近似求出y(x1)的近似值y, 般的利用在处的微分方程可得:，于这个差分格式我们称为欧拉公式. 特点 (1).单步方法; (2)显示格式 (3)局部截断误差因而是一阶精度. 2.改进欧拉方法 对其从到进行定积分得: 将右端的定积分用梯形公式来进行近似计算: 二.多步法 前面讲的方法：欧拉方法、改进欧拉方法、龙格－库塔方法均是单步方法，即在每一步要计算时，只要前面一个值已知的条件下秒可以计算出了。可以自成系统进行直接计算，因为初始条件只有一个已知，由可以计算，，用前面已知计算出来的，这前面计算好的个值来计算，这样自由度的增加来提高格式的精度。 如: 1.Adams方法： 1、方法构造的思想：微分方程初值问题将微分方程在上积分，. 即若已知来计算， 简记，用的拉格朗日插值多项式代替f 截断部分，用等距步长，得到的方法就是显式Adams方法. 的局部截断误差是因而是四阶精度的. K步的Adams（阿达姆斯）公式：. 1步显式Adams公式:,即Euler公式. 2步显式Adams公式: 3 步显式Adams 公式：4 步显式Adams 公式(较常用公式之一)： 类似可推导隐式 Adams 公式： 1 步隐式Adams 公式： 2 步隐式Adams 公式：3 步隐式Adams 公式(较常用公式之一)： 预测-校正方法：用显式公式预测，用隐式公式校正。 4 阶Adams 公式为例 校正： 三．实例分析： 运用matlab分别用改进欧拉方法和Adams（三步法）法进行求解. 求得精确解为: 分别在窗口中运行: 方法1欧拉： x0=0;xn=1.2;y0=2;h=0.1; [y,x]=Eulergj(x0,xn,y0,h); n=length(x); fprintf( i x(i) y1(i)\n); for i=1:n fprintf(%2d %4.4f %4.4f\n,i,x(i),y(i)); end 方法2 Adams（三步法）x0=0;xn=1.2;y0=2;h=0.1; [y,x]=Adamswai3(x0,xn,y0,h); n=length(x); fprintf( i x(i) y2(i)\n); for i=1:n fprintf(%2d %4.4f %4.4f\n,i,x(i),y(i)); end Matlab运行结果： 1.改进欧拉方法： i x(i) y1(i) 1 0.0000 2.0000 2 0.1000 1.6328 3 0.2000 1.4388 4 0.3000 1.3200 5 0.4000 1.2409 6 0.5000 1.1852 7 0.6000 1.1445 8 0.7000 1.1140 9 0.8000 1.0907 10 0.9000 1.0726 11 1.0000 1.0584 12 1.1000 1.0472 13 1.2000 1.0382 2. Adams（三步法）法i x(i) y2(i) 1 0.0000 2.0000 2 0.1000 1.5934 3 0.2000 1.4076 4 0.3000 1.2198 5 0.4000 1.1877 6 0.5000 1.1160 7 0.6000 1.1038 8 0.7000 1.0738 9 0.8000 1.0630 10 0.9000 1.0484 11 1.0000 1.0400 12 1.1000 1.0319 13 1.2000 1.0261 把数据放到excel中处理得: 改进的欧拉方法 Adams（三步法）法 精确解 误差 误差 0 2.0000 2.0000 0 2.0000 0 0.1 1.6097 1.6328 -0.023143 1.5934 0.016227 0.2 1.4181 1.4388 -0.020732 1.4076 0.010525 0.3 1.3037 1.3200 -0.016363 1.2198 0.083848 0.4 1.2281 1.2409 -0.0