原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
PAGE
4 -
实验报告填写说明
（实验项目名称、实验项目类型必须与实验教学大纲保持一致）
1．实验环境：
实验用的硬件、软件环境。
2．实验目的：
根据实验教学大纲，写出实验的要求和目的。
3．实验原理：
简要说明本实验项目所涉及的理论知识。
4．实验步骤：
这是实验报告极其重要的容。对于验证性验，要写清楚操作方法，需要经过哪几个步骤来实现其操作。对于设计性和综合性实验，还应写出设计思路和设计方法。对于创新性实验，还应注明其创新点。
5．实验结论：
根据实验过程中得到的结果，做出结论。
6．实验总结：
本次实验的收获、体会和建议。
7．指导教师评语及成绩：
指导教师依据学生的实际报告内容，给出本次实验报告的评价和成绩。
1 实验环境：
MATLAB R2014a
2 实验目的：
（1）掌握Matlab基本知识，能够编写简单程序
?(2)?熟练掌握求解常微分方程初值问题的常用数值解法
3实验原理：
显式Euler法：
隐式Euler法：
改进Euler法：
四阶Adams显式公式：
四阶Runge-Kutta公式：
4实验步骤：
显式Euler法
①在Matlab中创建M文件。
②取步长h=0.1，变量t根据步长h从0取到1，初值y（1）=1。
③进行显式Euler迭代：
for i=1:n-1;
y(i+1)=y(i)+h*(t(i)+1-y(i))
end。
④运行M文件，得到解。
隐式Euler法：
①进行隐式Euler迭代（在迭代过程中加上对y1（i+1）的赋值）
for i=1:n-1;
y1(i+1)=y(i)+h*(t(i)+1-y(i));
y(i+1)=y(i)+h*(t(i+1)+1-y1(i+1))
end
②运行M文件，得到解。
改进Euler法：
①进行改进Euler迭代在迭代过程中加上对y1（i+1）的赋值）
for i=1:n-1;
y1(i+1)=y(i)+h*(t(i)+1-y(i));
y(i+1)=y(i)+h/2*(t(i)+1-y(i)+t(i+1)+1-y1(i+1))
end
②运行M文件，得到解。
四阶Runge-Kutta法：
①
对进行赋值，其中
②根据迭代公式y(i+1)=y(i)+h/6*(K1+2*K2+2*K3+K4)以及①中K1,K2,K3,K4的值求解问题。运行M文件，得到解。
四阶Adams显式公式：
①根据四阶Runge-Kutta法中所求得的前四个迭代解分别赋值到y（1），y（2），y（3），
y（4）。
②根据迭代公式
y(i+1)=(1/24)*(18.5*y(i)+5.058*y(i-1)-3.7*y(i-2)+0.9*y(i-3)+0.24*i+3.24)进行迭代，不同的是需要迭代次数需减少，因为在①中已经为前四个解进行了赋值。
③运行M文件，得到解。
5实验结论
显式Euler法
隐式Euler法
改进Euler法
四阶RK法
四阶Adams显式
1.0000
1.0000
1.0000
1.0000
1.0000
1.0000
1.0100
1.0050
1.0048
1.0048
1.0100
1.0281
1.0190
1.0187
1.0187
1.0290
1.0536
1.0412
1.0408
1.0408
1.0561
1.0857
1.0708
1.0703
1.0746
1.0905
1.1240
1.1071
1.1065
1.1133
1.1314
1.1679
1.1494
1.1488
1.1574
1.1783
1.2168
1.1972
1.1966
1.2051
1.2305
1.2703
1.2500
1.2493
1.2565
1.2874
1.3279
1.3072
1.3066
1.3109
1.3487
1.3894
1.3685
1.3679
1.3679
表1 迭代结果表格
6实验总结（收获、体会和建议）
初步掌握了运用Matlab编写简单程序求解常微分方程初值问题的方法，巩固并熟悉了求解常微分方程初值问题的常用数值解法。
7指导教师评语及成绩：
评 语
评语等级
优
良
中
及格
不及格
实验方案设计的合理程度
实验结论的记录情况
实验总结情况
实验报告是否按时完成,字迹清楚,文字叙述流畅,层次清晰
成 绩：
指导教师签名：
批阅日期：
附录1：源 程 序
%显式Euler法
clear
clc
h=0.1;
t=0:h:1;
n=length(t);
y(1)=1
for i=1:n-1;
y(i+1)=y(i)+h*(t
文档评论（0）