我要加入 登录
声振论坛 返回首页

liliangbiao的个人空间 http://home.vibunion.com/?69372 [收藏] [复制] [分享] [RSS]

日志

Henon映射的分叉图的Matlab实现

已有 2609 次阅读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');

发表评论 评论 (2 个评论)

回复 kpgkpg 2008-11-18 20:22
if i>1960
这句代码什么意思呢?
回复 liliangbiao 2008-11-19 09:04
if i> 1960

facelist doodle 涂鸦板

您需要登录后才可以评论 登录 | 我要加入

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-2 21:31 , Processed in 0.036318 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

返回顶部