初值问题地Euler方法和梯形法.doc
《初值问题地Euler方法和梯形法.doc》由会员分享,可在线阅读,更多相关《初值问题地Euler方法和梯形法.doc(13页珍藏版)》请在三一文库上搜索。
1、学生实验报告实验课程名称偏微分方程数值解开课实验室 数统学院学院数统年级2013 专业班信计02学生姓名学号开课时 间 2015 至2016 学年第 2 学期总成绩教师签名数学与统计学院制开课学院、实验室:数统学院实验时间:2016年 月曰实验项目名称初值问题的Euler方法和梯形法实验项目类型验证演示综合设计其他指导教师曾芳成绩是一. 实验目的通过该实验,要求学生掌握求解初值问题的欧拉法和梯形法,并能通过计算机语言编程实现这两种算法。二. 实验内容考虑如下的初值问题:du一 2u,t 0,1dt u 01该问题有解析解u t = e21。1. 用欧拉法求解该问题,取步长h= 0.2,0.1,
2、0.05 ,将3种步长的计算结果(tn - nh,n 1,2,.,1/ h时刻的计算结果),解析结果和相应的绝对误差列表显示。2. 用梯形法求解该问题,取步长h- 0.2,0.1,0.05 ,将3种步长的计算结果(tn - nh,n 1,2,.,1/ h时刻的计算结果),解析结果和相应的绝对误差列表显示。3. 在同一种方法下,请说明哪种网格大小的计算结果更加精确,并说明理由。在相同的网格大小下,比较上述两种算法的计算结果,那种算法的结果要好一些,并说明理由。三实验原理、方法(算法)、步骤欧拉法的迭代格式及误差估计:Un+1 =U n+hf(t n ,Un). IUn-U( tn ) =0(h)
3、 欧拉法:fun cti on x,y=euler(fu n, xO,xfi nal,y 0,n)if n arg in<5,n=50;endh=(xfi nal-x0)/n;x(1)=x0;y(1)=y0;for i=1: nx(i+1)=x(i)+h;y(i+1)=y(i)+h*feval(fu n, x(i),y(i);endx=x'y=y'x1=0:0.2:1y1=exp(2*x1)plot(x,y,x1,y1)fun cti on f=doty(x,y);f=2*yx,y=euler('doty',0,1,1,10)梯形法:fun cti on
4、x,y=tix in g(fu n, x0,xfi nal,y 0,n)if n arg in<5,n=50;endh=(xfi nal-x0)/n;x(1)=x0;y(1)=y0;for i=1: nx(i+1)=x(i)+h;y(i+1)=y(i)+h*feval(fu n,x(i),y(i)y(i+1)=y(i)+h*(feval(fu n, x(i),y(i)+feval(fu n,x(i+1),y(i+1)/2;endx=x'y=y'x1=0:0.2:1y1二exp(2*x1)plot(x,y,x1,y1)四实验环境(所用软件、硬件等)及实验数据文件Matlab
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 初值问题 Euler 方法 梯形
链接地址:https://www.31doc.com/p-12472657.html