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、结果