Lorzen混沌方程求解-ode45

ode45函数

ode45是matlab中的求解微分方程数值解的函数，使用语法为：
[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需要在函数中编写好，以lorzen混沌方程为例

function dydt = lorenz_diff(t,y)
dydt = [ -10*(y(1)-y(2));
          30*y(1)-y(1)*y(3)-y(2);
          -8/3*y(3)+y(1)*y(2)];

调用ode45函数进行求解

clear
clc
y0 = [0,2,9];          
[t,y] = ode45(@lorenz_diff,[0,200],y0); 
plot(y(:,1),y(:,3),'-');
grid on;

lorzen系统相位图