- 1、本文档共9页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
10
Runge-Kutta Methods–II
Absolute Stability
10.1 Absolute Stability of RK Methods
The notion of absolute stability developed in Chapter 6 for LMMs is equally
relevant to RK methods. Applying an RK method to the linear ODE x (t) =
λx(t) with (λ) 0, absolute stability requires that xn → 0 as n → ∞.
We begin with the most general two-stage, second-order RK method dis-
cussed in Section 9.4.
Example 10.1
Investigate the absolute stability of the RK method with Butcher array
0 0
a a 0 a = 1/(2θ).
1 − θ θ
Applying this method to the solution of the ODE x (t) = λx(t) and writing
D.F. Griffiths, D.J. Higham, Numerical Methods for Ordinary Differential Equations,
Springer Undergraduate Mathematics Series, DOI 10.1007/978-0-85729-148-6_10,
© Springer-Verlag London Limited 2010
136 10. Runge-Kutta Methods–II Absolute Stability
1
h = λh, we find that
hk1 = hf (tn , xn ) = hλxn = hxn ,
hk = hf (t + ah, x + ahk ),
2 n n 1
= h(x + ahk ) = h(1 + ah)x .
n 1 n
Hence,
xn+1 = xn + (1 − θ)hk1 + θhk2
= 1 + h(1 + θah) xn
and, since aθ = 1 , it follows that
2
您可能关注的文档
- pachelbel-canon-in-d-major外文电子书籍.pdf
- Ph_D defense (12)外文电子书籍.pdf
- Reputation and Public Administration外文电子书籍.pdf
- The major transitions in evolution外文电子书籍.pdf
- The mental representation of causal conditional reasoning外文电子书籍.pdf
- U1P1 Warming up1外文电子书籍.ppt
- Vivid Dreams Unlocking Shadow Memories外文电子书籍.pdf
- 按揭专员月工作总结外文电子书籍.doc
- 设备点检管理制度(宝钢).docx
- 天津医科大学天然药物化学课程试卷.docx
文档评论(0)