function f=ispectre(Shat,Omega);

%                  f=ispectre(Shat,Omega);
%
% Calcule, etant donne un signal omega --> Shat (omega) donne comme un
% vecteur colonne, 
%           omega=-Omega:Omega/M:Omega-Omega/M, 
% les valeurs de l'antecedent approche de Shat par la prise de spectre 
% echantillonne aux points 
%              t= k pi/Omega, k=-M:1:M-1
% Si Omega = N*pi, on obtient la restitution sur 
% t = -1/2 :1/M : 1/2-1/M 
% du motif du signal 1-periodique ayant pour coefficients de Fourier
% entre -M et M-1 le vecteur Shat. 

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;
