function f=hilbertanalog(s,T);

%              f=hilbertanalog(s,T);
%
%Calcule, étant donné un signal analogique S=S(t), t=-T:T/M:T-T/M, donné
%comme un vecteur colonne, la transform\'ee de Hilbert du signal S, soit
%f(t)=CaH[S](t), t=-T:T/M:T-T/M. La sortie est donnée (comme l'entrée) sous
%forme de vecteur colonne. 

s=s'; 
N=length(s);
M=fix(N/2);
w=-ones(1,N);
for k=1:2:N
  w(k)=-w(k);
end
u=fft(s);
u=w.*u;
u1=2*u(2:M);
spectre=[zeros(1,M),u(1), u1];
spectre=(2*T/N)*spectre;
f=imag(ispectre(spectre',M*pi/T)'); 