matlab用ode45求解微分方程组
matlab用ode45求解微分方程组?
ode45是最常用的求解非刚性微分方程(组)的指令,它采用变步长四、五阶Runge-Kutta法,求解精度比ode23、ode15要高。但ode45的计算量比较大,而ode23计算量小,且误差大
调用格式:[t,y]=ode45(odefun,tspan,y0)
odefun 用以表示f(t,y)的函数句柄或inline函数,t是标量,y是标量或向量
tspan 如果是二维向量[t0,tf],表示自变量初值t0和终值tf;如果是高维向量[t0,t1,t2,。。。,tn],则表示输出结点列向量
y0 表示初值向量
t 表示结点列向量(t0,t1,t2,。。。,tn)^T
y 表示数值解矩阵,每一列对应y的一个分量。
实例说明:
例1 y=y-2t/y,y(0)=1,0