数值分析报告.docVIP

《常微分方程数值解》实验报告 实验序号：　　1 实验项目名称：欧拉方法 学　　号 201141210230 姓　　名 谢冬冬 专业、班 11级信计二班 实验地点 K7-405 指导教师 汪玉霞 时 间 2014-3-11 实验目的及要求 （1）熟悉运用显式计算格式 向前Euler方法： 中点公式： 改进Euler方法： 或 求解常微分方程和方程组的初值问题 的过程。 （2）熟悉运用隐式计算格式 向后Euler方法： 求解一阶线性常微分方程初值问题，并将计算结果与上两种方法的结果对照。 （3）对李普希兹常数巨大的微分方程，用向前Euler方法在不同步长下计算，对方法的数值稳定性进行分析，得出向前Euler方法的绝对稳定区间的有限性。 二、实验设备要求与实验原理 1、支持Intel Pentium Ⅲ及其以上CPU，内存256MB以上、硬盘1GB以上容量的微机；软件配有Windows98/2000/XP操作系统及MATLAB软件； 2、实验过程中，务必将实验所得数据列表，分析数据，得出实验结论。（建议同时网上提交电子版实验报告：780668212@） 实验内容与步骤 1. 分别用向前Euler、向后Euler、改进Euler方法求解 并与精确解在各节点处的值对照，计算绝对误差，分析各算法的性能。 2．分别用向前Euler方法、向后Euler方法求解. 步长分别为，计算处的数值解，并与函数值比较。 实验程序、结果与数据处理 （一）题目一 向前Euler方法： function [x,y]=xiedongdong1(a,b,h,y0) %qianouler n=(b-a)/h; x=zeros(1,n+1);y=zeros(1,n+1);x(1)=a;y(1)=y0; for i=1:n x(i+1)=x(i)+h; y(i+1)=y(i)+h*youxianf(x(i),y(i)); end dsolve(Dy=x-y+1,y(0)=1,x) ans = x+exp(-x) jqz=1+exp(-1) jqz = 1.3679 format long jqz=1+exp(-1) jqz = 1.367879441171442 [x,y]=xiedongdong1(0,1,0.1,1) Columns 8 through 11 1.178296900000000 1.230467210000000 1.287420489000000 1.348678440100000 wucha=jqz-y(11) wucha = 0.019201001071443 [x,y]=xiedongdong1(0,1,0.01,1) Columns 97 through 101 1.341047118104550 1.347236646923505 1.353464280454270 1.359729637649727 1.366032341273230 wucha2=jqz-y(101) wucha2 = 0.001847099898213 [x,y]=xiedongdong1(0,1,1/1000,1) Columns 997 through 1001 1.365169890791096 1.365800720900305 1.366431920179404 1.367063488259225 1.367695424770966 wucha3=jqz-y(1001) wucha3 = 1.840164004764944e-004 列表比较不同步长所得结果 步长 数值解 精确解 误差（取绝对值） 0.1 1.348678440100000 1.367879441171442 0.019201001071443 0.05 1.358485922408542 1.367879441171442 0.009393518762900 0.01 1.366032341273230 1.367879441171442 0.001847099898213 改进Euler方法： function [x,y]=xiedongdong2(a,b,h,y0) n=(b-a)/h; x=zeros(1,n+1);y=zeros(1,n+1);x(1)=a;y(1)=y0; for i=1:n x(i+1)=x(i)+h; k1=youxianf(x(i),y(i)); k2=youxianf(x(i+1),y(i)+h*k1); y(i+