clear
close all
%dynamique des fréquences
freq=[-4:0.01:4];

%paramètres du système
%en continu 
aprime = -(2*sqrt(2)+1);
bprime=-4;
epsilon=0.1;

%analyse en temps continu
s= j*2*pi*freq;
Hs=ft_pilote(aprime, bprime, epsilon, s);

%Affichage
n1=figure;
hold on
plot(freq, abs(Hs), 'r-')
grid on
xlabel('Fréquence')
ylabel('Module')
legend('|Hs|')

disp('fin de l''analyse en temps continu')
pause

%analyse en discret
%indiquer une période d'echantillonnage
Te1=1/4
%indiquer l'équation de passage du plan des "s" au plan des "z" 
z=exp(2*pi*j*freq*Te1);
s=(1-(1./z))/Te1;

Hz1=ft_pilote(aprime, bprime, epsilon, s);

figure(n1)
hold on
plot(freq, abs(Hz1),'g-')
legend('|Hs|', '|Hz1|')

%calcul de l'ecart à 1hz
indice = 401+100;
freq(indice)
abs(Hs(indice))
abs(Hz1(indice))

%on change de période d'echantillonage
Te2=0.0192
%indiquer l'équation de passage du plan des "s" au plan des "z" 
z=exp(2*pi*j*freq*Te2);
s=(1-(1./z))/Te2;

Hz2=ft_pilote(aprime, bprime, epsilon, s);

figure(n1)
plot(freq, abs(Hz2),'b-')
legend('|Hs|', '|Hz1|', '|Hz2|')

%calcul de l'ecart à 1hz
abs(Hz2(indice))

erreur1 = 100*abs(abs(Hs(indice))-abs(Hz1(indice)))/abs(Hs(indice))
erreur2 = 100*abs(abs(Hs(indice))-abs(Hz2(indice)))/abs(Hs(indice))