原创力文档- 1、本文档共10页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
数值计算方法52
第五章 常微分方程数值解
§ 5.2 Runge-Kutta法
k 1,2,3,
§ 5.2 Runge-Kutta法
考虑改进Euler法
yk yk 1 hf (xk 1 , yk 1 )
h
y y [f (x , y ) f (x , y )]
k k 1 k 1 k 1 k k
2
h
如果将其改成 y k y k 1 (K1 K2 )
2
K1 f (xk 1 , yk 1 )
(1)
K f (x , y hK )
2 k k 1 1
y y (x )
0 0
2
1
改进Euler法是由梯形公式和Euler公式复合而成
梯形公式具有2阶精度
同样可以证明,改进Euler法也具有2阶精度
形如(1)式的求解公式称为二阶Runge-Kutta法
对于Simpson求解公式：
h h h
y y [f (x , y ) 4f (x , y (x )) f (x , y )]
k k 1 k 1 k 1 k 1 k 1 k k
6 2 2
这是隐式多步法 选取适当的显化方法,可得类似(1)
的高阶Runge-Kutta方法
以下使用中值定理进行推导 3
为了同学们课后复习的方便, 以下的内容将k写成n.
一、Runge-Kutta方法的导出
对于常微分方程的初值问题
y f (x , y ) a x b
y (a) y 0
的解y y (x), 在区间[x , x ]上使用微分中值定理, 有
n1 n
(x , x )
y (x ) y (x ) y ( )(x x ) n1 n1 n
n n1 n1 n n1
即 y (x ) y (x ) hy ( )
文档评论（0）