微分方程的数值解法matlab(四阶龙格—库塔法).pptxVIP

微分方程的数值解法四阶龙格—库塔法（The Fourth-Order Runge－Kutta Method）常微分方程（Ordinary differential equations, ODE)初值问题---给出初始值边值问题---给出边界条件与初值常微分方程解算有关的指令ode23 ode45 ode113 ode23tode15s ode23sode23tb一.解ODE的基本机理:1.列出微分方程(2.1)初始条件令2. 把高阶方程转换成一阶微分方程组 (2.2)(2.3)例：著名的Van der Pol方程令 降为一阶初始条件3. 根据式(2.2)编写计算导数的M函数文件-ODE文件把t,Y作为输入宗量，把 作为输出宗量 %M function file name: dYdt.m function Yd = f (t, Y) Yd = f (t,Y) 的展开式例Van der Pol方程 %M function file name: dYdt.m function Yd = f (t, Y) Yd=zeros(size(Y));[t, Y]=solver (‘ODE函数文件名’, t0,tN, Y0, tol);ode454. 使编写好的ODE函数文件和初值 供微分方程解算指令（solver）调用说明：t0：初始时刻；tN：终点时刻Y0：初值；tol：计算精度Solver解算指令的使用格式输出宗量形式例题1：著名的Van der Pol方程解法1：采用ODE命令% 主程序 （程序名：VanderPol _ex1.m） t0 = 0;tN = 20;tol = 1e-6; Y0 = [0.25; 0.0]; [t, Y]=ode45 (‘dYdt’, t0,tN, Y0, tol); subplot (121), plot (t, Y) subplot (122), plot (Y( :, 1), Y( :, 2))Van der Pol方程% 子程序 （程序名： dYdt.m ） function Ydot = dYdt (t, Y)Ydot=[Y(2);-Y(2)*(Y(1)^2-1)-Y(1)];或写为function Ydot = dYdt (t, Y)Ydot=zeros(size(Y));Ydot(1)=Y(2);Ydot(2)=-Y(2)*(Y(1).^2-1)-Y(1)];各种solver 解算指令的特点解法指令解题类型特 点适合场合ode45非刚性采用4、5阶Runge－Kutta法大多数场合的首选算法ode23非刚性采用Adams算法较低精度（10－3）场合ode113非刚性多步法；采用Adams算法；高低精度均可（10－3～10－6）ode45计算时间太长时取代ode45ode23t适度刚性采用梯形法则算法适度刚性ode15s刚性多步法；采用2阶Rosenbrock算式，精度中等当ode45失败时使用；或存在质量矩阵时ode23s刚性一步法；采用2阶Rosenbrock算式，低精度低精度时，比ode15s有效；或存在质量矩阵时ode23tb刚性采用梯形法则－反向数值微分两阶段算法，低精度低精度时，比ode15s有效；或存在质量矩阵时二. 四 阶 Runge-Kutta 法对 I=[a,b]作分割步长初值问题的数值解法分为两大类计算 的近似值 时只用到 ，是自开始方法 单步法-Runge-Kutta 方法多步法-Admas方法Runge-Kutta法是常微分方程的一种经典解法MATLAB 对应命令：ode45四阶Runge-Kutta公式四 阶 Runge-Kutta 法计算流程图初始条件： ； 积分步长： 迭代次数：for i = 1 : N Next i开始子程序计算输出结果EndPlot三. Runge-Kutta 法解Van der Pol 方程的Matlab 程序结构主程序：RK_vanderpol.m 子程序：RK_sub.m（函数文件） 解法2：采用Runge_Kutta法编程计算主程序：RK_vanderpol.mt0=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_