第三章数值积分和常微分方程数值解法.ppt

3-2-2-1 Runge-Kutta法——问题的提出 例： 基元反应： k2 已知：[A]0=[B]0=1mol·L-1, [C]0=[D]0=[E]0=0 1mol·L-1, k1=1.0min-1， k2=0.25min-1， k3=0.5min-1， 求各组分浓度随时间的变化（0-20min）. 基本思想： 求平均斜率 差商： 由微分中值定理，存在 使 得 由 (xi,xi+1)上平均斜率 3-2-2-2 Runge-Kutta法——方法原理 二阶RK公式：（改进Euler） 改进Euler： 二点斜率 通式： 满足： 3-2-2-2 Runge-Kutta法——方法原理 ——变形Euler 三阶RK公式：（三点斜率） 或 3-2-2-2 Runge-Kutta法——方法原理 四阶RK公式：（四点斜率） 或 3-2-2-2 Runge-Kutta法——方法原理 常微分方程组四阶RK公式： 变步长：判据： （方程） （方程组） 3-2-2-2 Runge-Kutta法——方法原理 RK4(F,X,Y,N,H) H2=H/2,H6=H/6 DO I=1，N-1 XI=X(I), YI=Y(I) X12=XI+H2 , X11=XI+H X(I+1)=X11 RK1=F(XI,YI) RK2=F(X12,YI+H2*RK1) RK3=F(X12, Y