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