微分方程的数值解法 四阶龙格 库塔法 (The Fourth-Order Runge Kutta Method )常微分方程(Ordinary differential equations, ODE) n 初值问题- 给出初始值 n 边值问题- 给出边界条件 与初值常微分方程解算有关的指令 ode23 ode45 ode113 ode23t ode15s ode23s ode23tb一. 解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(si