% 产生 Henon 序列
% 使用平台 - Matlab6.5
% 作者：陆振波，海军工程大学
% 欢迎同行来信交流与合作，更多文章与程序下载请访问我的个人主页
% 电子邮件：luzhenbo@sina.com
% 个人主页：luzhenbo.88uu.com.cn
% modified by superdont:http://blog.youkuaiyun.com/superdont
%
Henon映射x(n
+
1
)
=
1
-
a
*
x(n)
^
2
+
y(n);y(n
+
1
)
=
b
*
x(n)
clc
clear
closeall
a
=
1.4
;
b
=
0.3
;
x0
=
0
;
y0
=
0
;
N1
=
2000
;
%
前面的迭代点数
N2
=
8000
;
%
后面的迭代点数
z
=
zeros(N1
+
N2,
2
);
for
i
=
1
:N1
+
N2
x
=
1
-
a
*
x0
^
2
+
y0;
y
=
b
*
x0;
x0
=
x;
y0
=
y;
z(i,
1
)
=
x;
z(i,
2
)
=
y;
end
z
=
z(N1
+
1
:end,:);
X
=
z(:,
1
);
Y
=
z(:,
2
);
plot(X,Y,
'
r.
'
,
'
MarkerSize
'
,
1
)
xlabel(
'
z
'
);ylabel(
'
y
'
)
title(
'
Henonattractor
'
)
% 使用平台 - Matlab6.5
% 作者：陆振波，海军工程大学
% 欢迎同行来信交流与合作，更多文章与程序下载请访问我的个人主页
% 电子邮件：luzhenbo@sina.com
% 个人主页：luzhenbo.88uu.com.cn
% modified by superdont:http://blog.youkuaiyun.com/superdont
%
Henon映射x(n
+
1
)
=
1
-
a
*
x(n)
^
2
+
y(n);y(n
+
1
)
=
b
*
x(n)clcclearclosealla
=
1.4
;b
=
0.3
;x0
=
0
;y0
=
0
;N1
=
2000
;
%
前面的迭代点数N2
=
8000
;
%
后面的迭代点数z
=
zeros(N1
+
N2,
2
);
for
i
=
1
:N1
+
N2x
=
1
-
a
*
x0
^
2
+
y0;y
=
b
*
x0;x0
=
x;y0
=
y;z(i,
1
)
=
x;z(i,
2
)
=
y;endz
=
z(N1
+
1
:end,:);X
=
z(:,
1
);Y
=
z(:,
2
);plot(X,Y,
'
r.
'
,
'
MarkerSize
'
,
1
)xlabel(
'
z
'
);ylabel(
'
y
'
)title(
'
Henonattractor
'
)
扫一扫
7095
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
