蔡氏电路与logistic映射

蔡氏电路:

function ydot = DyDt( t,y)

ydot=[9*(y(2)-y(1)+0.68*y(1)-0.5*(-1.27+0.68)*(abs(y(1)+1)-abs(y(1)-1)));...
      y(1)-y(2)+y(3);...
      -14.87*y(2)];
 
end  


tspan=[0,100];
options=odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-4]);
y0=[0.025;-0.022;0.8];
[tt,yy]=ode45(@DyDt,tspan,y0);
%plot3(yy(:,1),yy(:,2),yy(:,3))
subplot(1,3,1);plot(yy(:,1),yy(:,2));title('\fontname{隶书}x-y相平面图');
subplot(1,3,2);plot(yy(:,1),yy(:,3));title('\fontname{隶书}x-z相平面图');
subplot(1,3,3);plot(yy(:,2),yy(:,3));title('\fontname{隶书}y-z相平面图');

tspan&#