Matlab中解常微分方程的ode.docVIP

Matlab中解常微分方程的ode45 [原创]Matlab中解常微分方程的ode 452010-11-03 11：44ode是专门用于解微分方程的功能函数，他有ode23,ode45,ode23s等等，采用的是Runge-Kutta算法。ode45表示采用四阶，五阶runge-kutta单步算法,截断误差为(Δx)^3。解决的是Nonstiff(非刚性)的常微分方程.是解决数值解问题的首选方法,若长时间没结果,应该就是刚性的,换用ode23来解.其他几个也是类似的用法 使用方法 [T,Y]=ode45(odefun,tspan,y0) odefun是函数句柄,可以是函数文件名,匿名函数句柄或内联函数名 tspan是区间 [t0 tf]或者一系列散点[t0,t1,.,tf] y0是初始值向量 T返回列向量的时间点 Y返回对应T的求解列向量 [T,Y]=ode45(odefun,tspan,y0,options) options是求解参数设置,可以用odeset在计算前设定误差,输出参数,事件等 [T,Y,TE,YE,IE]=ode45(odefun,tspan,y0,options) 每组(t,Y)之产生称为事件函数。每次均会检查是否函数等于零。并决定是否在零时终止运算。这可以在函数中之特性上设定。例如以events或@events产生一函数。[value,isterminal,d