实验报告——常微分方程及数值解法.docVIP

实验报告——常微分方程及数值解法.doc

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

liwenhua00 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档