微分方程的数值解法n 四阶龙格 库塔法n (The Fourth-Order Runge Kutta Method )常微分方程(Ordinary differential equations, ODE)n 初值问题- 给出初始值n 边值问题- 给出边界条件与初值常微分方程解算有关的指令ode23 ode45 ode113 ode23tode15s ode23sode23tb一. 解ODE 的基本机理:2. 把高阶方程转换成一阶微分方程组1. 列出微分方程初始条件令(2.1) (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);4. 使编写好的O