Program2
% ---------------------------------------
% correlation des echantillons en x de la fonction de
Henon
% a=1.4, b=0.3, x0=y0=0
% ---------------------------------------
clear;
n = input('nombre d''iterations = ');
x=zeros(1,n+1);
y=zeros(1,n+1);
x(1)=0; y(1)=0; % initialisation
des parametres.
a=1.4; b=0.3;
% --------------------
% remplissage des vecteurs x et y
% --------------------
for i=1:n
y(i+1)=b*x(i);
x(i+1)=1+y(i)-a*(x(i)^2);
end
%------------------
% affichage des echantillons
%------------------
figure1=plot(x(1:length(x)-1), x(2:length(x)), 'r.');
% ici, affiche x(i+1) en fonction de x(i).
% figure1=plot(x(1:length(x)-2), x(3:length(x)), 'r.');
--> pour afficher x(i+2) en fonction de x(i).
grid;
title('Correlation des echantillons x(i+1) et x(i)');
xlabel('x(i)');
ylabel('x(i+1)');
zoom;
% permet de "zoomer et dezoomer" sur une partie du graphe avec la souris.