[bs] Ikeda Map

Referensi : wikipedia

Matlab code

clc

clear all

n=5000;

mu=0.9;

x0=0.1;

y0=0.1;

% Initialize

t=0.4-6/(1+x0^2+y0^2);

x(1,1)=1+mu*(x0*cos(t)-y0*sin(t));

y(1,1)=mu*(x0*sin(t)+y0*cos(t));

% Simulate

for i=2:n

    t=0.4-6/(1+x(i-1,1)^2+y(i-1,1)^2);

    x(i,1)=1+mu*(x(i-1,1)*cos(t)-y(i-1,1)*sin(t));

    y(i,1)=mu*(x(i-1,1)*sin(t)+y(i-1,1)*cos(t));

end

plot(x,y,'b.','MarkerSize',4)

axis equal

axis([1.1*min(x) 1.1*max(x) 1.1*min(y) 1.1*max(y)])

xlabel('X','FontSize',12);ylabel('Y','FontSize',12)

title('Ikeda Map','FontSize',12)