Discussion :

[Missay]

Bonjour,

En fait j'ai un sinusoïde de fréquence f0= 100 Hz et fréquence échantillonnage
Fe=1000
Normalement je devrais observer mes deux raies en f1=-100 Hz et f2=100 Hz

mais en programmant, j'ai deux raies une à f1=-400 Hz et f2=400 Hz ce qui n'es pas NORMAL

voila mon code :

Code :

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

1
2
3
4
5
6
7
8
f0  = 100;
fs  = 1000; 
t   = 0:1/fs:1;
x   = sin(2*pi*f0*t);
y   = abs(fft(x)); 
f   = linspace(-fs/2,fs/2,length(x));
%f=0:fs/length(x):fs-(fs/length(x))
plot(f,y);

Un petit coup de pouce SVP

[Invité]

Bonjour,

La sortie de la fonction FFT correspond aux fréquences allant de 0 à fs, et non -fs/2 à fs/2. L'intervalle [fs/2 fs] étant identique à [-fs/2 0] (périodicité de la FFT), on utilise la fonction fftshift par dessus pour avoir l'ordre que tu recherches

Code :

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

y   = abs( fftshift(fft(x)) );

[Missay]

Merci beaucoup