function f=ispectre(Shat,Omega);

%                  f=ispectre(Shat,Omega);
%
% Calcule, etant donné un signal omega --> Shat (omega) donné comme un
% vecteur colonne, 
%           omega=-Omega:Omega/M:Omega-Omega/M, 
% les valeurs de l'antécédent approché de Shat par la prise de spectre 
% échantillonné aux points 
%              t=k pi/Omega, k=-M:1:M-1

N=length(Shat);
M=fix(N/2);
w=-ones(N,1);
for k=1:2:N
  w(k)=-w(k);
end
u=N*ifft(Shat);
u=w.*u;
u1=u(1:M);
u2=u(M+1:N);
f=[u2;u1];
f=(Omega/(N*pi))*f;