初值问题的Euler方法和梯形法.doc

学 生 实 验 报 告 实验课程名称 偏微分方程 数统学院 学 院 数统 年级 2013 专业班 信计 02 学 生 姓 名 学 号 开 课 时 间 2015 至 2016 学年第 2 学期 总 成 绩 教师签名 数学与统计学院制 开课学院、实验室： 数统学院 2016年 月 日 实验项目 名 称 初值问题的Euler方法和梯形法 实验项目类型 验证 演示 综合 设计 其他 指导教师 曾芳 成 绩 是 一．实验目的 通过该实验，要求学生掌握初值问题的欧拉法和梯形法并能通过计算机语言编程实现两种算法。如下的初值问题： 该问题有解析解欧拉法求解，取步长将种步长的结果时刻的解析结果的误差列表显示 用梯形法求解该问题，取步长将种步长的结果时刻的解析结果的误差列表显示 在同一种方法下，请说明种网格大小的计算结果更加精确，并说明理由。相同的网格大小下，比较上述两种算法的计算结果，那种算法的结果要好一些，并说明理由。 function [x,y]=euler(fun,x0,xfinal,y0,n) if nargin5,n=50; end h=(xfinal-x0)/n; x(1)=x0;y(1)=y0; for i=1:n x(i+1)=x(i)+h; y(i+1)=y(i)+h*feval(fun,x(i),y(i)); end x=x; y=y; x1=0:0.2:1 y1=exp(2*x1) plot(x,y,x1,y1) function f=doty(x,y); f=2*y [x,y]=euler(doty,0,1,1,10) 梯形法： function [x,y]=tixing(fun,x0,xfinal,y0,n) if nargin5,n=50; end h=(xfinal-x0)/n; x(1)=x0;y(1)=y0; for i=1:n x(i+1)=x(i)+h; y(i+1)=y(i)+h*feval(fun,x(i),y(i)) y(i+1)=y(i)+h*(feval(fun,x(i),y(i))+feval(fun,x(i+1),y(i+1)))/2; end x=x; y=y; x1=0:0.2:1 y1=exp(2*x1) plot(x,y,x1,y1) 四．实验环境（所用软件、硬件等）及实验数据文件 Matlab 五．实验结果及实例分析 欧拉法输出分析： H=0.2 计算值 解析值 误差 1.0000 1.0000 0 1.4000 1.4918 -0.0918 1.9600 2.2255 -0.2655 2.7440 3.3201 -0.5761 3.8416 4.9530 -1.1114 5.3782 7.3891 -2.0108 H=0.1 1.0000 1.0000 0 1.2000 1.2214 -0.0214 1.4400 1.4918 -0.0518 1.7280 1.8221 -0.0941 2.0736 2.2255 -0.1519 2.4883 2.7183 -0.2300 2.9860 3.3201 -0.3341 3.5832 4.0552 -0.4720 4.2998 4.9530 -0.6532 5.1598 6.0496 -0.8899 6.1917 7.3891 -1.1973 H=0.05 1.0000 1.0000 0 1.1000 1.1052 -0.0052 1.2100 1.2214 -0.0114 1.3310 1.3499 -0.0189 1.4641 1.4918 -0.0277 1.6105 1.6487 -0.0382 1.7716 1.8221 -0.0506 1.9487 2.0138 -0.0650 2.1436 2.2255 -0.0820 2.3579 2.4596 -0.1017 2.5937 2.7183 -