- 1、本文档共8页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
初值问题的Euler方法和梯形法重点讲义
学 生 实 验 报 告
实验课程名称 偏微分方程 数统学院
学 院 数统 年级 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 -
您可能关注的文档
最近下载
- 2024.10政法干警锻造新时代政法铁军专题研讨班发言材料(5篇).docx VIP
- 医疗器械出库复核程序.docx
- 董责险-PPT_可编辑.ppt VIP
- 后勤岗位竞聘演讲稿PPT.pptx
- 历年华二自招考试数学试卷汇编(共5套,附答案).pdf
- 高州风土人情资料.ppt
- 食品加工技术专业及农产品加工类专业群建设项目.pdf
- 高一英语必修一单元精练Unit 3 Family Matters 重点单词变形词组短语句型(外研版2019).pdf VIP
- 英汉语言对比(华中科技大学)中国大学MOOC慕课 客观题题库答案.docx
- 《回弹法检测水泥基灌浆材料抗压强度技术规程》标准全文.docx VIP
文档评论(0)