N=length(x); % nombre d'échantillons du signal
n=[0:N-1]; % indices des échantillons
Fe=8000; % fréquence d'échantillonnage
t=n/Fe; % échelle temporelle 
figure(1); plot(t,x); legend('signal original')
Xf=fft(x); % calcul de la TFD du signal
f=n*Fe/N; % échelle des fréquences calculées par TFD~;
figure(2); plot(f,abs(Xf)); legend('spectre du signal original')