- 1、本文档共7页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
实验报告——常微分方程及数值解法
实验报告
实验项目名称 常微分方程的数值解法
实 验 室 数学实验室
所属课程名称 微分方程数值解
实 验 类 型 上机实验
实 验 日 期 2013年3月11日
班 级 10信息与计算科学
学 号 2010119421
姓 名 叶达伟
成 绩
实验概述: 【实验目的及要求】
运用不同的数值解法来求解具体问题,并通过具体实例来分析比较各种常微分方程的数值解法的精度,为以后求解一般的常微分方程起到借鉴意义。
【实验原理】
各种常微分方程的数值解法的原理,包括Euler法,改进Euler法,梯形法,Runge-Kutta方法,线性多步方法等。
【实验环境】(使用的软硬件)
Matlab软件
实验内容: 【实验方案设计】
我们分别运用Euler法,改进Euler法,RK方法和Adams隐式方法对同一问题进行求解,将数值解和解析解画在同一图像中,比较数值解的精度大小,得出结论。
【实验过程】(实验步骤、记录、数据、分析)
我们首先来回顾一下原题:
对于给定初值问题:
1. 求出其解析解并用Matlab画出其图形;
2. 采用Euler法取步长为0.5和0.25数值求解,并将结果画在同一幅图中,比较两者精度;
3. 采用改进Euler法求解,步长取为0.5;
4. 采用四级Runge-Kutta法求解,步长取为0.5;
5. 采用Adams四阶隐格式计算,初值可由四级Runge-Kutta格式确定。
步骤一,求出解析解并用Matlab画出其图形;
图一 初值问题的解析解的图像
步骤二,采用Euler法取步长为0.5和0.25数值求解(2.16),并将结果画在同一幅图中,比较两者精度;
采用Euler法取步长为0.5和0.25数值求解
图二 Euler方法的计算结果与解析解的比较
从图像中不难看出,采用Euler法取步长为0.5和0.25数值求解
图三 不同步长的Euler法的计算结果与解析解的绝对误差的比较
从图像中我们不难看出,步长为0.25的Euler法比步长为0.5的Euler法的精度更高。
步骤三,采用改进Euler法求解,步长取为0.5;
图四 改进Euler法的计算结果与解析解的比较
步骤四,采用四级Runge-Kutta法求解,步长取为0.5;
四级Runge-Kutta
图五 四级RK方法的计算结果与解析解的比较
步骤五,采用Adams四阶隐格式计算,初值可由四级Runge-Kutta格式确定。
Adams四阶隐格式计算
图六 Adams四阶隐格式计算
由于采用Adams四阶隐格式计算
为了将Adams四阶隐格式计算
图七 Adams四阶隐格式计算
最后,我们总结一下各种方法计算出来的数值解与解析解之间的绝对误差,详见下表:
Euler法
改进Euler法
四级RK方法
四阶Adams隐格式方法
0.0
0.5
1.0
1.5
2.0
2.5
3.0
3.5
4.0
4.5
5.0
表一 步长为0.5时各种数值解法的计算结果与解析解之间的绝对误差
由于本题中步长h取为取值0.5,步长取值较大,计算结果与精确解之间的误差较大,难以判断各种数值方法的计算精度。因此,我在本题的基础上进行了拓展,将步长取为0.05,进一步研究各种方法的计算精度。
从先前的例子中可以看出,误差主要是在t取值接近5的时候出现,所以在步长h=0.05时重点研究的计算结果与解析解之间的绝对误差,以比较各种方法的计算精度。绝对误差情况详见下表:
Euler法
改进Euler法
四级RK方法
四阶Adams隐格式方法
4.50
4.55
4.60
4.65
4.70
4.75
4.80
4.85
4.90
4.95
5.00
表二 步长为0.05时各种数值解法的计算结果与解析解之间的绝对误差
通过步长的减小,我们得出了更加精确的数值解,同时也可以比较各种算法的精度大小。不难发现,对于这个问题,四级RK方法的计算精度最大,其次是四阶Adams隐格式方法和改进Euler法,精度最差的是Euler法。(这不代表哪种方法最优,因为我们只计算了一个实例)
【结论】(结果)
通过一个具体实例,我们比较了各种数值计算方法的精度大小,对于这个问题,四级RK方法的计算精度最大,其次是四阶Adams隐格式方法和改进Euler法,精度最差的是Euler法。
【小结】
本次实验,将不同的常微风方程的数值计算方法在计算机上进行实现,使我们的理论知识掌握得更加扎实
文档评论(0)