原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
微分方程数值解法实验告
湖南工程学院 微分方程数值解法实验报告
专业班级 姓名 学号 同组实验人员 信息与计算科学1101 王珍 09
实验日期 2014年5 月6日 第 1次实验 指导老师 杨继明 评分 实验名称 用欧拉方法解常微分方程的初值问题 实验目的 熟悉掌握常微分方程的初值问题的数值格式并程序实现 考虑下列常微分方程的初值问题：
采用欧拉方法求解初值问题
采用Matlab程序设计语言编程实现该问题的数值求解。
在Matlab命令窗口中输入命令：
E=euler(0,1,1,20)
运行程序得到实验结果为：
E =
0 1.00000000000000
0.05000000000000 0.75000000000000
0.10000000000000 0.56250000000000
0.15000000000000 0.42187500000000
0.20000000000000 0.31640625000000
0.25000000000000 0.23730468750000
0.30000000000000 0.17797851562500
0.35000000000000 0.13348388671875
0.40000000000000 0.10011291503906
0.45000000000000 0.07508468627930
0.50000000000000 0.05631351470947
0.55000000000000 0.04223513603210
附Matlab程序代码：
function E=euler(a,b,ya,M)
%Input - y=f is the function
% - a and b are the left and right endpoints
% - ya is the initial condition y(a)
% - M is the number of steps
%Output - E=[T Y] where T is the vector of abscissas and
% - Y is the vector of ordinates
%If f is defined as an M-file function use the @ notation
% call E=euler(@f,a,b,ya,M).
%If f is defined as an anonymous function use the
% call E=euler(f,a,b,ya,M).
h=(b-a)/M;
T=zeros(1,M+1);
Y=zeros(1,M+1);
T=a:h:b;
Y(1)=ya;
for j=1:M
Y(j+1)=Y(j)+h*f(T(j),Y(j));
end
E=[T Y];
function f=f(t,y)
f=-5*y;
指导教师评语：
签字：
年 月 日
实验日期 2014年5 月6日 第 2次实验 指导老师 杨继明 评分 实验名称 利用peYF格式计算近似解 实验目的 熟悉掌握常微分方程的初值问题的数值格式并程序实现 考虑下列常微分方程的初值问题：
其精确解为
采用Matlab程序设计语言编程实现该问题的数值求解。
在Matlab命令窗口中输入命令：
u = peYF(-2,0.02,11,0,1,5)
运行程序得到实验结果为：
u =
Columns 1 through 4
1.74530439815628 1.75031196341828 1.46872586082017 1.00810241079480
Columns 5 through 8
0.54438411523664 0.25469560184372 0.24968803658172 0.53127413917983
Columns 9 through 11
0.99189758920520 1.45561588476336
文档评论（0）