Henon映射的数学表示为:
clear,clc
%% when a =0.5,b=0.05
a = 0.5;
b = 0.05;
x(1) = 0;
y(1) = 0;
for i = 2:2000
x(i) = y(i-1) - a*x(i-1)^2 +1;
y(i) = b.*x(i-1);
end
subplot(2,2,1)
scatter(x,y)
txt = {'a = 0.5','b= 0.05'};
text(0.1,0.04,txt)
set(gca,'xlim',[min(x),max(x)])
xlabel('n');
ylabel('x(n)');
%% when a=1,b=0.05
a = 1;
b = 0.05;
x(1) = 0;
y(1) = 0;
for i = 2:2000
x(i) = y(i-1) - a*x(i-1)^2 +1;
y(i) = b.*x(i-1);
end
subplot(2,2,2)
scatter(x,y)
txt = {'a = 1','b= 0.05'};
text(0,0.04,txt)
set(gca,'xlim',[min(x),max(x)])
xlabel('n');
ylabel('x(n)');
%% when a=1.2,b=0.1
a = 1.2;
b = 0.1;
x(1) = 0;
y(1) = 0;
for i = 2:2000
x(i) = y(i-1) - a*x(i-1)^2 +1;
y(i) = b.*x(i-1);
end
subplot(2,2,3)
scatter(x,y)
txt = {'a = 1.2','b= 0.1'};
text(-0.2,0.08,txt)
set(gca,'xlim',[min(x),max(x)])
xlabel('n');
ylabel('x(n)');
%% when a=1.4,b=0.3
a = 1.4;
b = 0.3;
x(1) = 0;
y(1) = 0;
for i = 2:2000
x(i) = y(i-1) - a*x(i-1)^2 +1;
y(i) = b*x(i-1);
end
subplot(2,2,4)
scatter(x,y)
txt = {'a = 1.4','b= 0.3'};
text(-1,0.2,txt)
set(gca,'xlim',[min(x),max(x)])
xlabel('n');
ylabel('x(n)');
结果: