chua.m

实现了三种混沌系统：蔡氏混沌系统

%蔡氏电路
function xdot=chua(num,chua_x)
apha=10;beta=14.87;a=-1.27;b=-0.68;%标准参数
%apha=9.4664;beta=14.0205;a=-1.4057;b=-0.5279;%禹思敏<多涡卷混沌吸引自研究的新结果>归一化后的参数
%apha=21.3,beta=23.9,a=-1.1,b=-0.7;

xdot=[-apha*chua_x(1)+apha*chua_x(2)-apha*[b*chua_x(1)+(a-b)*(abs(chua_x(1)+0.5)-abs(chua_x(1)-0.5))*0.5];chua_x(1)-chua_x(2)+chua_x(3);-chua_x(2)*beta];

%function xdot=chua(t,x);

%xdot=[(x(2)-x(1))*0.119782*(10^6)-(0.00032*x(1)+0.5*(-0.532)*(abs(x(1)+0.2)-abs(x(1)-0.2)))*0.197628*(10^9)
%    (x(1)-x(2))*0.012653*(10^6)-0.020876*(10^9)*x(3)
%    -x(2)*0.107526*(10^3)];


%蔡氏方程的相图
%clear all;
% option=odeset('AbsTol',1e-6,'Reltol',1e-4); 
%[t,y]=ode45(@chua,[0,100],[0.000296 -0.000635 -0.000654],option);
%plot3(y(:,1),y(:,2),y(:,3));
%title('chua方程相图');
%xlabel('x');ylabel('y');zlabel('z');grid;