梯形法则 matlab代码,matlab第二章常微分方程的数值解法

《matlab第二章常微分方程的数值解法》由会员分享，可在线阅读，更多相关《matlab第二章常微分方程的数值解法(36页珍藏版)》请在金锄头文库上搜索。

1、 微分方程的数值解法四阶龙格库塔法 (The Fourth-Order RungeKutta 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。

2、,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. 使编写好的ODE函数文件和初值 供微分 方程解算指令(solver)调用Solver解算指令的使用格式t, Y=solver (ODE函数文件名, t0, tN, Y0, tol);ode45输出宗量形式说明：t0：初始时刻；tN：终点时刻 Y0：初值； 。

3、tol：计算精度例题1：著名的Van der Pol方程% 主程序 (程序名：VanderPol _ex1.m)t0 = 0; tN = 20; tol = 1e-6; Y0 = 0.25; 0.0;t, Y=ode45 (dYdt, t