四阶龙格库塔法的基本思想_四阶龙格—库塔法的原理及其应用

● ● 高 教 视 野 龙榕一库锩法 摩理履冀麈厕 ◎冯建强 孙诗一 (扬州大学数学科学学院，江苏 扬州 225000) 【摘要】微分方程的数值模拟在工程应用 中有很重要的 意义．而在诸多的数值方法中，龙格一库塔法是使用最广 泛 、最有效 的数值 方法之一．龙格一 库塔 法的理论 基础是 泰 勒级数方法，龙格一库塔法吸收 了泰勒公式方法中的高精 度同时摒弃了泰勒级数方法中求高阶导数的弊端，利用复 合 函数 的思想合 理而 巧妙地 回避 了求 高 阶导数 这 一难 点， 使数值格式显得非常对称和紧凑．四阶龙格一库塔法是求 解微分方程的非常有用的工具，尽管其计算公式非常简单， 但该数值格式的基本原理却是非常深刻的． 首先，本文利用数学分析、数值分析的知识理解欧拉方 法及泰勒级数 方法 的局 限性 ，从 而理论 分析 导 出龙格一 库 塔 四阶格 式(四阶格 式的完整 导 出是任何 一本 参考 书上 所 没有的)．然后，探讨 了四阶龙格一库塔法在捷联惯性导航 中的应用．最后，对四阶龙格一库塔法的收敛性与稳定性进 行 了讨论． 【关键词】四阶龙格一库塔法；收敛性；稳定性；微 分 方程 【基金项目】该课题由扬州大学教改项 目YZUJX2016— 4A和扬州大学2016科创项 目资助． 一 、 概 述 对于一个复杂的函数，要用泰勒展开式求它的各阶导 数总是令人望而止步的，龙格一库塔法利用了泰勒级数法 的思想 ，但避免了对原来的函数进行解析求导的过程 ．所 以，本质上龙格一库塔法是建立在泰勒级数方法基础上的， 它摒弃了泰勒级数方法求导数的弊端，利用复合函数的思 想实 现求 导． 对于二阶、三阶的龙格一库塔法的推导，一些数值方法 的参考书都给出了具体的推导过程．但对于四阶龙格一库 塔法的推导，参考书上并没有指明，本文主要研究四阶龙格一 库塔法的推导以及四阶龙格一库塔法的具体应用实例． 二、四 阶龙格一库塔法 的理论推导 r 一 ．tJ、 设解初值问题{dt一八 ’的计算公式形如： 【 n(0)=Ⅱ。 kl 4-C2k2 4- 。3 uj+hb2lk1)， hb3Ikl+ uj+ hb4l l+ k3 4-C4k4)， ^632k2)， hb42 2+hb43k3)． 在 a2=b2l，a3=b31+bm a4=b4l+b42+b43的条件下 选择常数 Cl，c2，c3，c4，a2，a3，a4，b2l，b3l，b b4l，b42，b43使上 列的数值方法的截断误差的阶为 0(h )． 为此设 u(t)是方程 ：At,It)经过点( ， )的准确 解，将 u(t)在此点展开成幂级数 ： u( t)= ( )+hu ( )+ u”(0)+六^ u (々)+ “H ( )+O(h )． ● · · ● 司‘计算得到 H (t)=厂’11,”(t)= + ， t、={ j }{ {{ fjLu． “¨’(t)=／．Ⅲ+ +2[( + )，+ + ．， ]+ { f J {t { j、+ ’{ 、{t’{ {u l ’ j { f、 l ： f．舡 0 f + t 0j 0 s{0 ’{ u’ lt 1． 将 Ⅱ (t)，M”(t)，u (t)，u¨ (t)代入到 u(t )可得 u(tj+1)=u(tj)+hf+÷^ ( + )+六矿 + 0’ l +l{u 去 m } + l1jl l 0j+{ t 4-／J3． 对于原式 中 k。，k ，k，，k 在 (￡ ，uj)点 的幂 级数展开 式 ， 并注意应用公式和 已经写 出的 k．的表达式 ，可 以导 出 k。=f， 后