dodo

[H,f]=freqz(num,den,512,fe);
amp=abs(H);
phase=unwrap(angle(H));
% Tracé de la réponse fréquentielle en amplitude et phase
figure
subplot(2,1,1)
plot(f,amp),grid
title('Réponse fréquentielle en amplitude et phase pour [0;fe/2]'),
ylabel('Amplitude')
subplot(2,1,2)
plot(f,phase*180/(pi)), grid
ylabel('Phase en °'),
xlabel('fréquence (Hz)')
% Tracé du diagramme de Bode en amplitude et phase
figure
subplot(2,1,1)
plot(f,20*log10(amp)), grid
title('Diagramme de Bose en amplitude et phase pour[0;fe/2]'),
ylabel('Amplitude en dB')
subplot(2,1,2)
plot(f,phase*180/(2*pi)), grid
ylabel('Phase en °'),
xlabel('fréquence (Hz)')