Henon映射分叉图

Henon映射的方程为:

 

分叉图是观察最终值范围的一个重要方式。

下面我们通过matlab实现henon映射的分叉图。

clearvars
clc
close all
a = 0:0.001:1.4;
b = 0.3;
N1 = 5000;
N2 = 100;
x = ones(N1+N2,length(a));
y = ones(N1+N2,length(a));
for j = 1:N1+N2-1
    x(j+1,:) = 1+y(j,:)-a.*x(j,:).^2;
    y(j+1,:) = b*x(j,:);
end
f = x(N1+1:end,:);
g = y(N1+1:end,:);
figure
plot(a,f,'b')
xlabel('a')
ylabel('x')
figure
plot(a,g,'b')
xlabel('a')
ylabel('y')
figure
for i=1:100
    subplot(10,10,i);
    plot(a,f(i,:),'b.')
    xlabel('a')
    ylabel('x')
end
figure
for i=1:100
    subplot(10,10,i);
    plot(a,f(i,:),'b.')
    xlabel('a')
    ylabel('y')
end

首先观察x