Discussion :

[jericho]

Bonjour,
j'aimerai calculer la fréquence d'un signal de type sinusoidal. Pour cela j'utilise la FFT qui me permet de passer du domaine temporel au domaine fréquentiel.
Cependant, après avoir fait ma FFT, je trouve des resultats incohérents. Voici un exemple du code de la doc Matlab:

Fs = 1000; % Sampling frequency
T = 1/Fs; % Sample time
L = 1000; % Length of signal
t = (0:L-1)*T; % Time vector

y = 7*sin(2*pi*50*t) % je crée une sinusoide de 50hz et d'amplitude 7 volts

NFFT = 2^nextpow2(L); % Next power of 2 from length of y

Y = fft(y,NFFT)/L; % Calcul de la FFT

A ce stade Y est une matrice de 1024 points complexes représentant la FFT de y.
Pour trouver la fréquence j'ai mis dans une matrice le module de Y, trouvé le maximum, calculé la matrice avec la phase de Y, lu la valeur dont l'indice est le maximum pour la matrice précédente. Code:

module = abs(Y) % calcul du module (donc de l'amplitude?)
[valeur,Indice]=max(module) % valeur=3.2978 et I=52 (deja la je m'attendais à avoir valeur=~7)
phase = angle(Y) % calcul de la phase de Y (donc des fréquences?)
phase(52) % retourne -0.9572 !!!

Pouvez vous m'indiquer ou sont mes erreurs de raisonnement?
j'ai enfait 2 questions principales:
1)que représente le "Sampling frequency" ?
2)comment déduire la fréquence depuis la FFT?

Merci.

[jericho]

J'ai trouvé un début de réponse:
1)le Sample Time c'est le temps séparant 2 points de l'échantillon (et le Sampling Frequency son inverse = le nombre de points par seconde)
2)code pour trouver la fréquence:

y = 7*sin(2*pi*50*t); % je crée une sinusoide de 50hz et d'amplitude 7 volts
Y = fft(y,NFFT)/L; % Calcul de la FFT
amplitude = 2*abs(Y(1:NFFT));
[amplitudeMAX,FrequenceSignal] = max(amplitude)

% On obtient les bon résultats:
amplitudeMAX = 6.5814
FrequenceSignal = 52

Maintenant ce pose un autre problème, lorsque je mets une fréquence de 1MegaHz par exemple, j'obtiens à nouveau de faux résultats.
Quelqu'un aurait-il une idée d'où cela provient?

[jericho]

Après une bonne nuit de sommeil et après avoir rêvé de Fourier voici la solution! Dans mon cas j'ai besoin de calculer une fréquence entre 1Mhz et 10Mhz, j'ai donc mis le Sampling Frequency à 20 000 000.

%Parametres

Fs = 20000000; % Sampling frequency
T = 1/Fs; % Sample time
L = 1000; % Length of signal
t = (0:L-1)*T; % Time vector
NFFT = 2^nextpow2(L); % Next power of 2 from length of y

%Calculs

y = 7*sin(2*pi*3500000*t); % je crée une sinusoide de 3.5Mhz et d'amplitude 7 volts
Y = fft(y,NFFT)/L; % Calcul de la FFT
amplitude = 2*abs(Y(1:NFFT));
f = Fs/2*linspace(0,1,NFFT/2);
plot(f,2*abs(Y(1:NFFT/2)))
[amplitudeMAX,indiceFrequence] = max(amplitude);
frequence = f(indiceFrequence)

frequence =

3.5029e+006