日志

Henon映射的分叉图的Matlab实现

已有 2538 次阅读2008-9-26 11:00 |个人分类:分叉|

% Author: Thomas Lee
% E-mail: lixf1979@126.com
% Corresponding: School of Mathematics, Physics and Software Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China

%Plotting the bifurcation of the 2D Henon Map

b=0.3;
niter=2000;
x=0;
y=0;
A=[];
X=[];
%for
a=0:0.001:1.4;
for i=1:niter
     xprev=x;
     yprev=y;
     x=a-xprev.*xprev+b*yprev;
     y=xprev;
    % xprev=x;
    if i>1960
%hold on
A=[A;a];
X=[X;x];
end
end
%end
plot(A,X,'k.','markersize',1);
%set(h,{'MarkerSize'},{1});
title('Henon Bifurcation, a=0-1.4, b=0.3');