matlab第二章__常微分方程的数值解法.pptVIP

微分方程的数值解法 四阶龙格—库塔法 （The Fourth-Order Runge－Kutta Method） 常微分方程（Ordinary differential equations, ODE) 初值问题---给出初始值 边值问题---给出边界条件 一.解ODE的基本机理: 3. 根据式(2.2)编写计算导数的M函数文件-ODE文件 例题1：著名的Van der Pol方程 四阶Runge-Kutta公式 三. Runge-Kutta 法解Van der Pol 方程的Matlab 程序结构主程序：RK_vanderpol.m 子程序：RK_sub.m（函数文件） 解法2：采用Runge_Kutta法编程计算 主程序：RK_vanderpol.m t0=0; tN=20; y0=[0.25; 0]; h=0.001; t = t0 : h : tN; N = length (t); j = 1; for i = 1 : N t1 = t0 + h; K1 = RK_sub(t0, y0); K2 = RK_sub(t0 + h/2, y0 + h*K1/2); K3 = RK_sub(t0 + h/2, y0 + h*K2/2); K4 = RK_sub(t0 + h