matlab ode45用法

ode45表示采用四阶-五阶Runge-Kutta算法，它用4阶方法提供候选解，5阶方法控制误差，是一种自适应步长（变步长）的常微分方程数值解法，其整体截断误差为(Δx)^5。解决的是Nonstiff(非刚性)常微分方程。

ode45语法：

[T,Y] = ode45(odefun,tspan,y0)

[T,Y] = ode45(odefun,tspan,y0,options)

[T,Y,TE,YE,IE] = ode45(odefun,tspan,y0,options)

sol = ode45(odefun,[t0tf],y0...)

[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)

在设置了事件参数后的对应输出

TE 事件发生时间

YE 事件发生时之答案

IE 事件函数消失时之指针i

sol =ode45(odefun,[t0 tf],y0...)

sol 结构体输出结果

扩展资料：

如何在function里使用ode45输出值

(1) 主程式 (test.m)

边界值为 Y(1/1.5)=alpha=0 Y(1)=beta=0

用 shooting method 去解二阶 ode 的边界值问题，

解 ode 使用的指令为 ode45

(2)Function (funtest1.m)

解4 条first-order initial value problems

但a 的值是要从判断解出来的值运算後，是否有大於 1 来设定

H=0.25;

m=1.2;

si=((Y/x)^2-Y*Y'/x+(Y')^2)^0.5

if si>1

a=(si.^m-1)/(H*si)

elseif si<=1

a=0

end

参考资料：百度百科-ode45