数值分析实验报告第三次.pdfVIP

下载本文档

6
0
约5.21千字
约 6页
2017-08-12 发布于重庆
举报
版权申诉

数值分析实验报告第三次.pdf

1、本文档共6页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
8、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

数值分析实验报告第三次.pdf

华中科技大学本科实验报告课程名称：数值分析姓名：姜福鑫学号：U201310044 专业班级：应数1302 指导老师：黄乘明实验题目：常微分方程的数值解法日期：2015-5-7 实验成绩： 1. 实验目的通过上机实验，理解利用计算机迭代求解微分方程数值解的整个过程，加深对所学计算方法的理论及算法特点的理解。 2. 实验题目 2.1 h=0.02, matlab 取编制改进欧拉公式的程序，求下列常微分方程初值问题的数值解： 0.9 = ，0≤x≤1， y  12x y(0)=1. 2.2 h=0.1, - matlab 取编制四阶经典龙格库塔公式的程序，求下列常微分方程初值问题的数值解： 2 = ,0 x 1 y  ≤ ≤ ， yx y(0)=1. 3. 程序文本 3.1function[x,y]=meauler(dyfun,xspan,y0,h) x=xspan(1):h:xspan(2);y(1)=y0; forn=1:(length(x)-1) k1=feval(dyfun,x(n),y(n)); y(n+1)=y(n)+h*k1; k2=feval(dyfun,x(n+1),y(n+1)); y(n+1)=y(n)+h*(k1+k2)/2; end x=x; y=y; 3.2function[x,y]=marunge4(dyfun,xspan,y0,h) x=xspan(1):h:xspan(2);y(1)=y0; forn=1:(length(x)-1) k1=feval(dyfun,x(n),y(n)); k2=feval(dyfun,x(n)+h/2,y(n)+h/2*k1); k3=feval(dyfun,x(n)+h/2,y(n)+h/2*k2); k4=feval(dyfun,x(n+1),y(n)+h*k3); y(n+1)=y(n)+h*(k1+2*k2+2*k3+k4)/6; end x=x; y=y; 4. 运行结果与分析 4.1 dyfun=inline(-0.9/(1+2*x),x,y); [x,y]=meauler(dyfun,[0,1],1,0.02); [x;y] 0 1 0.02 0.98235 0.04 0.96536 0.06 0.94899 0.08 0.9332 0.1 0.91794 0.12 0.90318 0.14 0.88889 0.16 0.87504 0.18 0.8616 0.2 0.84856 0.22 0.83588 0.24 0.82355