常微分方程数值解法欧拉格式与梯形格式的比较.doc_第1页
常微分方程数值解法欧拉格式与梯形格式的比较.doc_第2页
常微分方程数值解法欧拉格式与梯形格式的比较.doc_第3页
常微分方程数值解法欧拉格式与梯形格式的比较.doc_第4页
全文预览已结束

VIP免费下载

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

年级00000班号0000学号0000专业000姓名00000实验名称000000实验类型设计型综合型创新型实验目的或要求常微分的数值解法用欧拉格式与梯形格式进行比较。实验原理（算法流程图或者含注释的源代码）已知准确解是首先建立准确解的函数文件function f=fun(x)f=1-exp(-5*x.2);欧拉格式程序如下：h=0.1;x=0:0.1:1;y=zeros(1,11);y(1)=0;for i=1:10 y(i+1)=y(i)+h*(10*x(i)*(1-y(i);endy1=fun(x);plot(x,y,r*,x,y1,*-)梯形格式程序如下：function f=funx(x,y)f=10*x.*(1-y);h=0.1;x=0:0.1:1;y=zeros(1,11);y(1)=0;for n=1:10 x(n+1)=x(n)+h; z0=y(n)+h*feval(funx,x(n),y(n); for k=1:3 z1=y(n)+h/2*(feval(funx,x(n),y(n)+feval(funx,x(n+1),z0); if abs(z1-z0)1e-3 break; end z0=z1; end y(n+1)=z1;endy1=fun(x);plot(x,y,r*,x,y1,*-)实验结果分析及心得体会程序运行结果如下：一、欧拉格式-. 是真实值，* 是数值解。从图形来看还是比较接近真实解的。这是梯形格式的数值解与真实解的图形。-. 是真实值，* 是数值解从图形看出，此题用梯形格式的数值解是非常接近真实值的。心得体会：欧拉格式的程序较简单，对于梯形格式的数值解法，也是逐步显式的过程，用Euler方法提供迭代处置，而得到迭代公式，从而建立程序的编写过程，这是关键。从
人人文库> 全部分类> 应用文书 > 合同范本

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论