常微分方程数值解法欧拉格式与梯形格式的比较

常微分方程数值解法欧拉格式与梯形格式的比较

ID:35240935

大小:55.50 KB

页数:4页

时间:2019-03-22

常微分方程数值解法欧拉格式与梯形格式的比较_第1页
常微分方程数值解法欧拉格式与梯形格式的比较_第2页
常微分方程数值解法欧拉格式与梯形格式的比较_第3页
常微分方程数值解法欧拉格式与梯形格式的比较_第4页
资源描述:

《常微分方程数值解法欧拉格式与梯形格式的比较》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库

1、年级00000班号0000学号0000专业000姓名00000实验名称000000实验类型设计型综合型创新型实验目的或要求常微分的数值解法用欧拉格式与梯形格式进行比较。实验原理(算法流程图或者含注释的源代码)已知准确解是首先建立准确解的函数文件functionf=fun(x)f=1-exp(-5*x.^2);欧拉格式程序如下:h=0.1;x=0:0.1:1;y=zeros(1,11);y(1)=0;fori=1:10y(i+1)=y(i)+h*(10*x(i)*(1-y(i)));endy1=fun(x);plot(x,y,'r*',x,y1,'*-')梯形格式程序如下:functionf=f

2、unx(x,y)f=10*x.*(1-y);h=0.1;x=0:0.1:1;y=zeros(1,11);y(1)=0;4forn=1:10x(n+1)=x(n)+h;z0=y(n)+h*feval(@funx,x(n),y(n));fork=1:3z1=y(n)+h/2*(feval(@funx,x(n),y(n))+feval(@funx,x(n+1),z0));ifabs(z1-z0)<1e-3break;endz0=z1;endy(n+1)=z1;endy1=fun(x);plot(x,y,'r*',x,y1,'*-')4实验结果分析及心得体会程序运行结果如下:一、欧拉格式-.是真实值,

3、*是数值解。从图形来看还是比较接近真实解的。4这是梯形格式的数值解与真实解的图形。-.是真实值,*是数值解从图形看出,此题用梯形格式的数值解是非常接近真实值的。心得体会:欧拉格式的程序较简单,对于梯形格式的数值解法,也是逐步显式的过程,用Euler方法提供迭代处置,而得到迭代公式,从而建立程序的编写过程,这是关键。从图形来看,梯形格式显然比欧拉格式更接近真实值。成绩评定教师签名:年月日4

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。