Rossler系统代码实现

1、Rossler系统的数学表现形式为

2、实现代码

clear,clc
tspn = [0 100];
y0 = [1,1,1];
a = 0.2; b= 0.2; c = 5.7;
[t,x]=ode45(@Rossler_ode,tspn,y0,[],a,b,c); 
%ode45(函数句柄，起止时间，初始值，flag（置空即可），传参)
plot3(x(:,1),x(:,2),x(:,3))

function dydt = Rossler_ode(t,y,a,b,c)
dxdt=-(y(2)+y(3));
dydt=y(1)+a*y(2);
dzdt=b+y(1)*y(3)-c*y(3);
dydt=[dxdt;dydt;dzdt];
end

3、结果