function estim_DSP(x,fe,flim)

% function estim_DSP(x,fe,flim);
% affiche l'estimation de la DSP du signal x
% fe fréquence d'échantillonnage
% flim limite de fréquence pour la représentation

L = length(x);
[Px,f] = pwelch(x-mean(x),L/2,[ ],[ ],fe);
Ind = find(f<=flim);
plot(f(Ind),Px(Ind));
xlabel('Hz')



