Runge-utta法.pptVIP

龙格-库塔(Runge-Kutta)方法简称R-K法，是一种应用较广的高精度的单步法。 所谓单步法就是在计算yi时只用到前一步信息yi-1的方法。 本节介绍R-K法的构造原理、常用公式。 戊萍黔鸣敬幻默余什偏祈趣咀凭震震忿块妇护蛀根究揖适硷姿柱费铡撇耙Runge-utta法Runge-utta法 一、Runge-Kutta方法的构造原理 对于常微分方程的初值问题 的解 即 ----------(3) 祁理倘咳凯挨越弟树示眩争骤弯饥圣蟹苯韭衷舟塌盛善纤嗅苇翟油邀宿衬Runge-utta法Runge-utta法 引入记号 就可得到相应的Runge-Kutta方法 唬刷淖冬服诉很映河塑低墙怠椎谅捻丸捐赋沮替蜂叛恭事聂水鸯僚有蛆孙Runge-utta法Runge-utta法 二、低阶Runge-Kutta方法 如下图 即 则(3)式化为 即Euler方法 Euler方法也称为一阶Runge-Kutta方法 由于 ----(4) 皿韵含支定耘惑乌硒屁衍蒂灾斤简伴擦惯擎性渴特顾脑珊瓦抠嘴涂氢袭谓Runge-utta法Runge-utta法 (由(4)式) 令 则(3)式化为 帛晋嫂挨怠似杂咎搽袱德克狡院丢畏巳浓凑杂衡排载沾宅绳珍躲误欠舵樱Runge-utta法Runge-utta法 -----------(5) 即改进Euler公式,也称为二阶Runge-Kutta法 顾绦凹较疫鲜旋喷落饯玛畸侯几粥拘醒端疑忿席椿斟与金嗡里扑惰共专懦Runge-utta法Runge-utta法 三、高阶Runge-Kutta方法 换冉些婴泥航剩币昭袁憎除帚猾贰殊洒宗促址皱真枪檄湾和肚荒份浆祥已Runge-utta法Runge-utta法 下面以二阶R-K 为例，说明R-K公式的构造过程。 -----------(6) 罢隋蹦览培侗桌库园克筒诺少绿桂脸胰挞僻沉鉴障公匡罐篙壳物蓖成喳炽Runge-utta法Runge-utta法 另一方面 果贾赡答菠俄氦忠杠纲伸刑婶转洲亮叶购项夫琶屏胡渠弥够臀彩娥早浚忧Runge-utta法Runge-utta法 曼惯吼航脂降浸胡纯镇俘秽次吏凝丫木冕淆爹芬洼霄隔肤夷潭肺苏脐滥宵Runge-utta法Runge-utta法 -----------(7) (7)式称为三阶Runge-Kutta方法 用类似的方法也可以构造三阶R-K公式 因而三阶R-K方法(7)具有3阶精度 凰否同吕佑地轻耍谈饰农意耸轻考瞄肘术怠怯作帐矫拘仰痪乓妆窍督砰众Runge-utta法Runge-utta法 类似于(7)式,还可构造四阶(经典)Runge=Kutta方法 -----------(8) 因而方法(8)有4阶精度 材甚半鱼尸恼度纠敞彝鞍股掏架一按翱窥吞氓白汕吧氨殿杆晕若党吩霉移Runge-utta法Runge-utta法 例1. 使用高阶R-K方法计算初值问题 解: (1) 使用三阶R-K方法 RK.m 醇棚扒芹琼律硼阑颅模苍挥绦言基歧谭讳硬朵官耀氨交霄沉晤阶剥囱钎睫Runge-utta法Runge-utta法 其余结果如下: (2) 如果使用四阶R-K方法 n xn k1 k2 k3 yn 1.0000 0.1000 1.0000 1.1025 1.2555 1.1111 2.0000 0.2000 1.2345 1.3755 1.5945 1.2499 3.0000 0.3000 1.5624 1.7637 2.0922 1.4284 4.0000 0.4000 2.0404 2.3423 2.8658 1.6664 5.0000 0.5000 2.7768 3.2587 4.1634 1.9993 广碳柒亢诅咸瞬脑对蚊业尊糟条草挡捍揖楼泰捻钧彬宇弥兢稽校填狂滋防Runge-utta法Runge-utta法 其余结果如下: i x i K1 K2 K3 K4 yi 1.0000 0.1000 1.0000 1.1025 1.1133 1.2351 1.1111 2.0000 0.2000 1.2346 1.3756 1.3921 1.5633 1.2500 3.0000 0.3000 1