- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
数分析第九章
第一,共69。
第九章
常微分方程的数值解
一、Euler方法
二、Runge-Kutta方法
三、单步法的收敛性和稳定性
四、线性多步法
第二,共69。
很多科学技术和工程问题常用常微分方程的形
式建立数学模型.
但是对于绝大多数的微分方程问题,很难或者根本
不可能得到它的解析解.
本章重点考察一阶方程的初值问题
的数值解法,就是寻求解y(x)在一系列离散点
处的近似值的方法.
相邻两个节点间的距离称为步长.
第三,共69。
一、Euler方法
1欧拉公式
由初值条件表示积分曲线从
出发,并在处的切线斜率为
因此可以设想积分曲线在x=x0附近可以用切
线近似的代替曲线.
切线方程为
当x=x时,代入有
1
这样得到y(x)的近似值y的方法.
11
第四,共69。
重复上述方法,当x=x2时
依次可以计算出x,x,…处的近似值y,y,…
3434
由此得到Euler公式:
由于用折线近似代替方程的解析解,所以Euler方
法也称为Euler折线法.
例用Euler法计算初值问题的解在x=0.3时的近似
值,取步长h=0.1.
第五,共69。
解:
Euler公式的截断误差
局部截断误差:一步Euler公式产生的误差;
总体截断误差:Euler公式的累积总误差;
第六,共69。
在假设y=y(x),即第i步计算是精确的前
定义nn
提下,考虑的截断误差R=y(x)y称为局
nn+1n+1
部截断误差.
p+1
若某算法的局部截断误差为O(h),则称
定义
该算法有p阶精度.
欧拉法的局部截断误差:
所以欧拉法具有1阶精度.
第七,共69。
Lipschitiz条件:若存在正数L,使得对一切
x,y,y有
12
则称f(x,y)满足Lipschitiz条件.
欧拉法的总体截断误差:
设那么
为局部截断误差,所以
第八,共69。
第九,共69。
特别当n=m-1时,有
总体误差与h是同阶的.
上式还说明,当时,有即
也就是说ym收敛到方程的准确解
第十,共69。
后退
文档评论(0)