function [omega,f]=spectreenergie(A,B,p)

%                [omega,f]=spectreenergie(A,B,p); 
%
% Etant donné un filtre ARMA de paramètres
%              A = [a0, a1,...,aN]      B= [b0,...,bM]
% (ces vecteurs sont déclarés en ligne) et un entier p tel que 
%                  2^p > sup(M,N), 
% cette routine calcule les valeurs du spectre d'énergie du filtre aux
% points 
%              omega = (2*pi/2^p)*[0:2^p-1] 


Den=fft(A,2^p);
Num=fft(B,2^p); 
omega=(2*pi/2^p)*[0:2^p-1];
f=(abs (Num./Den)).^2; 