推荐一个学习混沌的网站

这上面有大量的关于混沌和分形matlab程序，大家可以上去看看。

http://luobo.yculblog.com/archive.56999.html

Henon系统

function newx=Henon(x);
% Henon方程[差分方程]
% newx=Henon(x)
%
% 方程如下:
%   x(k+1)=-p*x(k)*x(k)+y(k)+1
%   y(k+1)=q*x(k)
%
% Example:
%   Y=Iterative(1000,[1;0]);
%   plot(Y(:,1),Y(:,2),'.');
%
% Author's email: ustb03-07@yahoo.com.cn
%
p=1.4;
q=.3;
newx=zeros(2,1);

 

newx(1)=-p*x(1)*x(1)+x(2)+1;
newx(2)=q*x(1);

 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

 

function Y=Iterative(n,X);
% 方程迭代
% 参数如下:
%   Y=Iterative(n,X);
%   n-迭代次数
%   X-函数的初值
%
% Example:
%   Y=Iterative(1000,[1;0]);
%
% Author's email: ustb03-07@yahoo.com.cn
%
for i=1:n
   
% 对Henon方程迭代
  X=feval(@Henon,X);
  
  Y(i,:)=X;
end

Rossler系统

function dx=Rossler(t,x);
% Rossler方程[不显含时间t的自治系统]
% 方程如下:
% dx/dt=-y-z
% dy/dt=x+a*y
% dz/dt=b+z*(x-c)
%
% Example:
% [T,Y]=ode45('Rossler',[0,400],[1;1;1]);
% plot3(Y(:,1),Y(:,2),Y(:,3));
%
% Author's email: ustb03-07@yahoo.com.cn
%
a=.15;
b=.2;
c=10;
dx=zeros(3,1);

dx(1)=-x(2)-x(3);
dx(2)=x(1)+a*x(2);
dx(3)=b+x(3)*(x(1)-c);

Chen系统

function dx=Chen(t,x);
% Chen方程[不显含时间t的自治系统]
% 方程如下:
% dx/dt=a*(y-x)
% dy/dt=(c-a)*x+c*y-x*z
% dz/dt=x*y-b*z
%
% Example:
% [T,Y]=ode45('Chen',[0,50],[1;1;1]);
% plot3(Y(:,1),Y(:,2),Y(:,3));
%
% Author's email: ustb03-07@yahoo.com.cn
%
a=35;
b=3;
c=28;
dx=zeros(3,1);

dx(1)=a*(x(2)-x(1));
dx(2)=(c-a)*x(1)+c*x(2)-x(1)*x(3);
dx(3)=x(1)*x(2)-b*x(3);