- 1、本文档共56页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
7: ODE * 求此问题数值解的标准四阶Runge-Kutta公式为 7: ODE * 7.5.3 高阶常微分方程组初值问题的数值解法 对于高阶常微分方程组的初值问题, 考察下面的二阶常微分方程组初值问题: 令 将此问题转化为一阶常微分方程组初值问题 7: ODE * 等价的一阶常微分方程组初值问题 7: ODE * 本章小结 基本知识 建立初值问题数值解的计算格式的方法 数值微分方法、Taylor 展开法和数值积分方法 截断误差与方法的精度 Runge-Kutta方法 基本思想 一阶、二阶、四阶Runge-Kutta方法 一阶常微分方程组数值解法 高阶方程(组)的数值解法 * * * 当前的整体截断误差=上一步的整体截断误差+上一步整体截断误差引起的导函数f(x,y)的误差+当前的局部截断误差。 * 7: ODE * 7.1.2 截断误差与方法的精度 定义 7.1.1 称误差 en+1 = y ( xn+1 ) - yn+1 为数值方法在点 xn+1 的截断误差, 又称整体截断误差. 设 yk= y ( xk ) (k = 0,1,. . . , n),则 为数值方法在点 xn+1 的局部截断误差. 7: ODE * 整体截断误差 en+1 是在没有引进舍入误差的情况下, 纯粹因为不准确的计算格式造成的, 故又称为方法误差.它不仅与 x = xn+1 这一步的计算有关, 而且和 xn , xn-1 ,. . . , x1 这几步的计算都有关系. 局部截断误差是假设 xn 之前各数值解没有误差, 仅由 xn 到 xn+1 这一步计算由计算格式引起的误差. 7: ODE * 如Euler公式 在点 xn+1 的整体截断误差 en+1 = y (xn+1)- yn+1 局部截断误差 当前的整体截断误差=上一步的整体截断误差+上一步整体截断误差引起的导函数 f (x, y) 的误差+当前的局部截断误差。 7: ODE * 定义7.1.2 若某数值方法的局部截断误差为 则称该方法具有 P 阶精度, 或称其为 P 阶方法. 可以证明: Euler 方法的局部截断误差 其具有一阶精度. 梯形方法的局部截断误差 其具有二阶精度. 改进的 Euler 方法的局部截断误差 具有二阶精度. 7: ODE * 7.2 Runge-Kutta 方法 继续讨论前面的 Taylor 展开法。 设 y∈C2[a , b], 由 Taylor 公式有 由 故上式即为 略去余项, 并以 yn , yn+1 代替 y (xn) 及 y (xn+1), 得到Euler 公式. 7: ODE * 进一步假设 y∈Cp+1[a, b],由 Taylor 公式有 其中 由 故式(7.2.1)即为 7: ODE * 7.2 Runge-Kutta 方法 略去余项, 并以数值解 yn , yn+1 代替 (7.2.3) 中的 解析解y (xn) 及 y (xn+1), 可得到一个差分方程,即 其中余项可写成 注:这里 7: ODE * 在(7.2.3)中略去余项,用yn, yn+1代替y(xn)及y(xn+1) 其中 称 (7.2.4) 式为求解常微分方程初值问题数值解Taylor的格式 . 7: ODE * 由于局部截断误差 可知它是一个 p 阶方法。当p=1时,上式正是Euler 公式。但当 p 2 时,需要计算f (x, y(x) ) 的高阶导数,特别是对于复杂函数 f (x, y(x)) 的求导,这无疑是大大增加计算量,这是它最大的缺点。因此高阶的Taylor方法是不实用的。 德国数学家C.Runge及M.W.Kutta提出了一种改进策略,得到了至今还作为高精度单步法被广泛使用龙格-库塔法(Runge-Kutta method)。 7: ODE * 7.2.1 Runge-Kutta方法的基本思想 Runge-Kutta方法是利用 f 在某些点处函数值的线性组合替代(7.2.4)步长 h 后面括号中的因子来构造差分方程, 从而避免了高阶导数的计算, 这就是 Runge-Kutta 方法的基本思想. 用f 在某些点处函数值的线性组合替代这一部分 7: ODE * 其一般形式为: 其中 r 是上式中调用 f 的个数, r 称为级数,
您可能关注的文档
最近下载
- 2024年昆明滇池国家旅游度假区国有资产投资经营管理(集团)有限责任公司人员招聘笔试备考题库及答案解析.docx
- 【教案】胚胎工程的理论基础教学设计高二下学期生物人教版(2019)选择性必修3.docx VIP
- 2024年新教材七年级上册道德与法治1.2《规划初中生活》教案.docx
- 小学五年级体育《蹲踞式跳远:助跑与起跳》教学设计.doc
- 残疾儿童送教上门一生一案资料.doc
- 铅火法冶炼工(稿).doc
- 与时代同行--中国电视文艺晚会的发展趋势初探.pdf
- 中国-象棋基础知识(适合儿童教学用)课件108页.ppt
- 安川伺服驱动器使用说明书..docx
- 第九章 马利冷却塔NC8400用户手册.pdf
文档评论(0)