clc
close all
clear all

phi=20 ;%en degrés
rau=0.9;
N=100;
noise=normrnd(0,1,1,N);

a1=-2*cosd(phi)*rau;
a2=rau^2;
%filterAR = filt(1,[1,a1,a2]); % Le filtre AR
traj=zeros(1,N);
for k =3:N
    traj(k)=noise(k)-a1*traj(k-1)-a2*(traj(k-2));
end

%filterAR;
plot(traj) % Affichage de trajectoire
hold on
figure
plot(noise)