matlab ode45精辟说明

ode45
求非刚性微分方程的解，中阶方法。
语法 [自变量,应变量] = ode45(@（函数名）,步进,初值向量)

例：

分析：
y1,y2,y2是应变量,t是自变量，所以dy初始化zeros(3,1)确保是列向量
代码

function dy = odefun(t,y)
dy = zeros(3,1);    %列向量
dy(1) = y(2) * y(3);
dy(2) = -y(1) * y(3);
dy(3) = -0.51 * y(1) * y(2);
end

tspan = [0:0.1:12];% 时间跨度取0,0.1...,12
y0 = [0,1,1];% 初始值
[t,y] = ode45( @odefun, tspan, y0 );% 调用语句
% 绘图
plot(t,y(:,1),'-',t,y(:,2),'-.',t,y(:,3),'.')
legend('x','y','z')