ode45是最常用的求解非刚性微分方程(组)的指令,它采用变步长四、五阶runge-kutta法,求解精度比ode23、ode15要高。但ode45的计算量比较大,而ode23计算量小.
要求调用ode45函数求解
(1)建立函数文件funt.m:function up=funt(t,u) up=u-(2*t/u);(2)求微分方程:clc;clear; t0=0;tf=1;u0=1; [t,u]=ode45('funt',[t0,tf],u0);
[t,x]=ode45('myfun',[0,600],x0,[]);提示两个错误:错误使用 odearguments .
t>=0即可用ode45,它的物理意义类似时间,但也不局限于时间。你的问题比较模糊,建议将相关的指令图表以及微分解析再复习一遍