蔡氏电路与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=[0,100];
options=odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-4]);
%y0=[0.025;-0.022;0.8];
y0=[0.025;-0.02;0.08];
[tt,yy]=ode45(@DyDt,tspan,y0);
plot3(yy(:,1),yy(:,2),yy(:,3))

view(-15.5,26);
xlabel('\fontsize{16}x')
ylabel('\fontsize{16}y')
zlabel('\fontsize{16}z')

losistic映射
clear;clc;
n=100;
a=3.83;
x=0.9*ones(n,1);
for i=2:n
x(i)=a*x(i-1)*(1-x(i-1));
end
plot(x)
xlabel('\fontsize{16}a=3.83'),ylabel('\fontsize{16}x')
title('\fontsize{16}x\fontname{隶书}的变化与a的关系')

%---------------feigenbaum图---------
f=@(a,x)a*x*(1-x);
x=0.9;
for a=1:0.01:4
for i=1:200;
x=f(a,x);
if i>100
plot(a,x,'k','linewidth',1);
hold on
end
end
end
xlabel('\fontsize{16}a'),ylabel('\fontsize{16}x')

hold off