龙格-库塔方法基本原理.ppt
54页1、在连续系统的仿真中 主要的计算工作是求解一阶微分方程y f x y y x0 y0解析法只能用来求解一些特殊类型的方程 实际仿真问题中归结出来的微分方程主要靠数值解法 由于实际运算只能完成有限项或有限步运算 因此要将有些需用极限或无穷过程进行的运算有限化 对无穷过程进行截断 这样产生的误差成为截断误差 根据实际情况建立的数学模型往往难以求解 通常需要通过近似替代 将所求解的数学模型简化为易求解的数值计算问题后再进行求解 数学模型的理论解与数值计算问题的精确解之间的误差称为截断误差 这是计算方法本身带来的误差 所以也成为方法误差 2020 4 3 3 得到高精度方法的一个直接想法是利用Taylor展开假设式y f x y a x b 中的f x y 充分光滑 将y xi 1 在xi点作Taylor展开 若取右端不同的有限项作为y xi 1 的近似值 就可得到计算y xi 1 的各种不同截断误差的数值公式 例如 取前两项可得到 3 2龙格 库塔方法 2020 4 3 4 其中 P阶泰勒方法 若取前三项 可得到截断误差为O h3 的公式 类似地 若取前P 1项作为y xi 1 的近似值 便得
2、到 2020 4 3 5 显然p 1时 yi 1 yi hf xi yi 它即为我们熟悉的Euler方法 当p 2时 要利用泰勒方法就需要计算f x y 的高阶微商 这个计算量是很大的 尤其当f x y 较复杂时 其高阶导数会很复杂 因此 利用泰勒公式构造高阶公式是不实用的 但是泰勒级数展开法的基本思想是许多数值方法的基础 R K方法不是直接使用Taylor级数 而是利用它的思想 2020 4 3 6 龙格 库塔 R K 法的基本思想 Euler公式可改写成 则yi 1的表达式与y xi 1 的Taylor展开式的前两项完全相同 即局部截断误差为O h2 Runge Kutta方法是一种高精度的单步法 简称R K法 2020 4 3 7 同理 改进Euler公式可改写成 上述两组公式在形式上共同点 都是用f x y 在某些点上值的线性组合得出y xi 1 的近似值yi 1 且增加计算的次数f x y 的次数 可提高截断误差的阶 如欧拉法 每步计算一次f x y 的值 为一阶方法 改进欧拉法需计算两次f x y 的值 为二阶方法 局部截断误差为O h3 2020 4 3 8 于是可考虑用
3、函数f x y 在若干点上的函数值的线性组合来构造近似公式 构造时要求近似公式在 xi yi 处的Taylor展开式与解y x 在xi处的Taylor展开式的前面几项重合 从而使近似公式达到所需要的阶数 既避免求高阶导数 又提高了计算方法精度的阶数 或者说 在 xi xi 1 这一步内多计算几个点的斜率值 然后将其进行加权平均作为平均斜率 则可构造出更高精度的计算格式 这就是龙格 库塔 Runge Kutta 法的基本思想 一般龙格 库塔方法的形式为 2020 4 3 9 其中ai bij ci为待定参数 要求上式yi 1在点 xi yi 处作Tailor展开 通过相同项的系数确定参数 称为P阶龙格 库塔方法 10 Runge Kutta方法的推导思想 对于常微分方程的初值问题 的解y y x 在区间 xi xi 1 上使用微分中值定理 有 即 2020 4 3 11 引入记号 就可得到相应的Runge Kutta方法 2020 4 3 12 如下图 即 则上式化为 即Euler方法 Euler方法也称为一阶Runge Kutta方法 2020 4 3 二阶龙格 库塔法在 xi xi
4、1 上取两点xi和xi a2 xi a2h 以该两点处的斜率值K1和K2的加权平均 或称为线性组合 来求取平均斜率k 的近似值K 即 式中 K1为xi点处的切线斜率值K1 hf xi yi hy xi K2为xi a2h点处的切线斜率值 比照改进的欧拉法 将xi a2视为xi 1 即可得 2020 4 3 13 确定系数c1 c2 a2 b21 可得到有2阶精度的算法格式 2020 4 3 14 因此 将y xi 1 在x xi处进行Taylor展开 将在x xi处进行Taylor展开 2020 4 3 15 K1 hf xi yi 2020 4 3 16 这里有4个未知数 3个方程 存在无穷多个解 所有满足上式的格式统称为2阶龙格 库塔格式 令 对应项的系数相等 得到 2020 4 3 17 注意到 就是二阶龙格 库塔公式 也就是改进的欧拉法 因此 凡满足条件式有一簇形如上式的计算格式 这些格式统称为二阶龙格 库塔格式 因此改进的欧拉格式是众多的二阶龙格 库塔法中的一种特殊格式 若取 就是另一种形式的二阶龙格 库塔公式 2020 4 3 18 此计算公式称为变形的二阶龙格 库塔法 式
5、中为区间的中点 也称中点公式 Q 为获得更高的精度 应该如何进一步推广 2020 4 3 19 二级R K方法是显式单步式 每前进一步需要计算两个函数值 由上面的讨论可知 适当选择四个参数c1 c2 a2 b21 可使每步计算两次函数值的二阶R K方法达到二阶精度 能否在计算函数值次数不变的情况下 通过选择不同的参数值 使得二阶R K方法的精度再提高呢 答案是否定的 无论四个参数怎样选择 都不能使公式的局部截断误差提高到三阶 这说明每一步计算两个函数值的二阶R K方法最高阶为二阶 若要获得更高阶得数值方法 就必须增加计算函数值的次数 三阶龙格 库塔法 2020 4 3 20 为进一步提高精度 在区间 xi xi 1 上除两点xi和xi a2 xi a2h 以外 再增加一点xi a3 xi a3h 用这三点处的斜率值K1 K2和K3的加权平均得出平均斜率K 的近似值K 这时计算格式具有形式 2020 4 3 21 同理推导二阶公式 将y xi 1 和yi 1在x xi处进行Taylor展开 使局部截断误差达到O h4 使对应项的系数相等 得到系数方程组 参数的选择不唯一 从而构成一类不同
《龙格-库塔方法基本原理.ppt》由会员自***分享,可在线阅读,更多相关《龙格-库塔方法基本原理.ppt》请在金锄头文库上搜索。
药店营业员培训-终端拦截技巧.ppt
穿越A4纸三种图解.ppt
虚拟仿真实训软件的制作技术.ppt
社会调查的历史发展.ppt
铰链四杆机构的基本性质课件.ppt
黑暗中的碰撞-盲人体验馆设计环艺毕业设计.ppt
老僧与蝎子材料作文讲评及例文.ppt
青岛版六年级科学上册《2生物的启示》课件.ppt
项羽之死优秀课件课件.ppt
菩萨蛮-书江西造口壁课件.ppt
鸟的天堂PPT课件.ppt
辩证唯物主义知识结构图.ppt
镇静催眠药-护理药物学.ppt
菩萨蛮优秀课件(温庭筠)PPT.ppt
蒲松龄《狼》(复习总结课).ppt
近代史人物介绍-张灵甫.ppt
韦庄《菩萨蛮》(其二)ppt课件课件.ppt
荷史诗中的战争.ppt
粉磨对水泥质量的影响.ppt
财政学-社会保障制度.ppt
2023-04-02 50页
2023-04-02 53页
2022-08-17 23页
2022-08-17 30页
2022-08-17 16页
2022-08-17 18页
2022-08-17 17页
2022-08-17 33页
2022-08-15 64页
2022-08-15 71页