Discussion :

[Newenda]

Bonjour,

J'ai un signal sismique (signal) enregistré à une fréquence d'échantillonnage (Fs) de 125 Hz et en calculant sa FFT de cette manière:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
 
NFFT = 2^nextpow2(numel(signal)); 
Y = fft(signal)/numel(signal);
YA=2*abs(Y(1:NFFT/2+1));
f = Fs/2*linspace(0,1,NFFT/2+1);

j'obtiens:




J'aimerais savoir d'où vient cette symétrie étrange à environ 37Hz dans la FFT ? il n'y a aucune raison d'en avoir le signal étant sismique (mélange de bruit de partout, aléatoire etc.)

Et par quelle miracle la FFT s'arrête à la fréquence de Nyquist ?


Merci d'avance.

[Matersss]

Salut,

Je ne suis pas un pro dans ce domaine mais je sais que quand un système est discrétisé à un certain pas et que l'on trace son diagramme de Bode, il est tracé jusqu'à la fréquence de Nyquist car après c'est symétrique. Je me demande s'il peut y avoir un rapport avec ton cas.

Matersss

[Newenda]

merci Materss pour ta réponse,

Oui tu as raison, cependant ici la symétrie commence bien avant la fréquence Fs/2.

Je ne vois aucune raison de traitement qui pourrait expliquer ceci.. peut être une coupure analogique du système d'acquisition non?

[Invité]

Bonjour,

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

Y = fft(signal, NFFT)/numel(signal);

[Newenda]

super, merci Winjerome.

[anita.crus]

Bonjour

j'ai un signal fréquentiel que j'ai obtenu en effectuant des essais, je voudrais le lisser afin de se débarrasser du bruit, mais je ne sais pas comment le faire. on m'a conseillé d'appliquer une fft au signal ensuite une ifft pour récupérer la zone des pics qui m'intéressent mais ça ne marche pas!