1、 9.3 Runge-Kutta法法 第九章第九章 常微分方程數值解常微分方程數值解 8.3 Runge-Kutta法法 1 (,) kkkk yyhf xy 111 (,)(,) 2 kkkkkk h yyf xyf xy 考慮改進Euler法 如果將其改成 00 ()yy x 1 (,) kk Kf xy 211 (,) kk Kf xyhK 112 () 2 kk h yyKK -(1) 改進Euler法是由梯形公式和Euler公式復合而成 梯形公式具有2階精度 形如(1)式的求解公式稱為二階二階Runge-Kutta法法 同樣可以證明,改進Euler法也具有2階精度 0 11 ( (

2、, ), ) kkkk yy y f xx ay xy b hyx a 對于微分方程初值問題 如已知需要求， 1 1 01() ( , ,) kk k k kkk h h x f xhh yy y yy xhy 則由微分中值定理 可得 1 ),(, (,),) kkkk kkkkk fxxxhhx fx y xx y yyyy 其中為在區間上的平均斜率。 為在點上的斜率( ) 基本思路基本思路 0 1 1 2 (,) (,)(,) ( ) , () kkkk kkk k k k k yy yyyh a f x xf x x yy x hO y 在中 僅用一個點處的斜率 來近似代替區間上的平均斜

3、率， 局部 歐拉公式 截斷誤差為。 1 1211 121 1 (,) (,) (,)(,) 1 1 () () 2 2 pkk k ckpk kk kp k k c k k f x Kf h x f yyy y yyyy yy y xKf xK KK hh h yy 在 改進的歐拉公式 中 11 111 3 1 (,)(,) (,) (,) ,() kkk k k k k kk k k k h h xf x xK f xK xx yyy y O y h y 則是用點處的斜率和由此點 處信息預估的點處的斜率 的算術平均值來近似代替 區間上的平均斜率，局部截斷誤差為。 所以如果在區間上多預估幾個點

4、的斜率值， 再將它們的線性組合作為平均斜率的近似值， 則就有可能構造出精度更高的計算格式。 推廣 1 (, ) kkkk yyhxy h ( , )f x y其中 是用在一些點上值的線性組合來構成 這種單步法稱為Runge-Kutta方法,簡記為簡記為R-K公式公式. ,Runge-Kutta.RRf若 是由 個 值線性組合構成 則稱線性方法級 三階三階Runge-KuttaRunge-Kutta 1 4 , () kk h xx O 若在區間上再增加一個新點， 即用三個點上的斜率進行加權平均作為平均斜率， 則可望得到截斷誤差為的計算公式， 1 123 123 1 21 32 1213231

5、, , (,) (,) (,) () 1,1 kkkkk k k k k kk k k xphqhxxxx KKK Kf x Kf xK Kf x phph qhqh h pq y y yK KyyKK 預估預估 其中為區間上的三個點； ；00； 是三個斜率的線性組合系數。 即 三階龍格三階龍格- -庫塔公式庫塔公式 4 1 1 2 -() pq hO 如果取中點和終點的斜率， 則可得到三階龍一種截斷誤差為格的庫塔公式 1 21 32 3 231 12 1 (,) (,) 22 (,) (4) 6 141 , 666 k k k k k k k k y y y y Kf x Kf xK Kf

6、xK KK hh hh h Ky 即 三階龍格三階龍格- -庫塔公式庫塔公式 1 21 32 13 1 1 23 (,) (,) 33 22 (,) 33 (3) 4 13 ,0, 44 k k k kk k k k Kf x Kf xK Kf xK hh h y h h K y yyK y 即 4 12 33 ()- pq hO 如果取任意兩點，如和終點的斜率， 則可得到另一種截斷誤差為的三階龍格 庫塔公式 四階龍格四階龍格- -庫塔公式庫塔公式 1 5 1 ,1 2 () - kk pq h xx O : 若在區間上仍取三個點 （，），但在中點處又校正， 則可望得到截斷誤差為 四階古典形式

7、的龍格 庫 的計 塔公式 算公式， 1 21 32 4 12341 3 / (,) (,) 22 (,) 22 (,) (22) 6 k k k k k k k k k k hh hh Kf x Kf xK Kf xK Kf x y y y y yy K KKKK hh h 多校即(RK4)正一次 構造一般的R級Runge-Kutta方法 1 22211 11 1 ,11 1122 (,) (,) () (,) () k k k kkR k k RkRRR RR R Kf x Kf xpqK RK Kf xpqKqK y h KKK y y yy h hhh h ,Taylor iiis p

8、q其中等均為待定的參數, 根據展開 并由期望的階數確定, 且一般不唯一. Runge-kutta方法的階與級的關系方法的階與級的關系 在Runge-kutta計算格式(RK)中.計算函數值 f 的次數 R 稱為級級, 級數與階數是不同的, 可以證明R級Runge-kutta公式 的 最高階數是 R . 通常所說的 R 級 m 階Runge-kutta公式 指要計算 R個f(x,y)的函數值, 且對應的計算公式是 m 階的. Butcher得出如下Runge-kutta方法的級數級數R與階數階數m的對應關系: 因此, 通常使用4級4階Runge-kutta公式(RK4). ()fR每步計算 的個數級數 2 3 4 5 6 7 R8 2 3 4 4 5 6 R-2 可達到的最高精度階數 應當注意，高階應當注意，高階R-K公式的推導是基于初公式的推導是基于初 值問題的解值問題的解y(x)的的