function reponse_amplitude(b,a,fe);

% reponse_amplitude(b,a,fe);
% calcule la réponse en amplitude d'un filtre
% entrées:
% b,a vecteurs de coefficients du filtre
% si filtre RIF, mettre g,1
% fe fréquence d'échantillonnage. Si e paramètre est présent,
% les fréquences vraies en Hz apparaissent dans le graphe. sinon,
% ce sont les fréquences normalisées

h = freqz(b,a,2000);
if nargin==2,
    plot(linspace(0,0.5,2000),abs(h))
else
    plot(linspace(0,0.5*fe,2000),abs(h))
end
title('réponse en amplitude')
