Discussion :

[chneu87]

Bonjour,

Je réalise une FFT sur un signal et je désire récupérer l'amplitude et la valeur (fréq) du Peak (le max).

Code :

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

1
2
3
4
5
6
7
8
9
 
L = length(Signal);
Y = fft(Signal);
FT = abs(Y/L);
DispFT = 2*FT(1:L/2+1);
freq = Sample_Freq*(0:L/2)/L;
 
[Amplitude, IFreq] = max(DispFT);
Frequence = freq(IFreq);

Ce code fonctionne (Frequency = 200Hz) mais je voudrai maintenant avoir le Peak entre 100 Hz et 500 Hz (faire abstraction des autres peak max si < à 100Hz ou > à 500Hz)

J'ai tenté ça :

Code :

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

1
2
3
4
5
6
7
 
 
FTCalcul = DispFT(100:500);
FreqCalcul = Freq(100:500);
 
[Amplitude, IFreq] = max(FTCalcul);
Frequency = FreqCalcul(IFreq);

Je n'obtient pas avec ce code le même résultat qu'avant (Frequency = 186.4), alors que je test avec une sinusoïde de 200Hz (entre 100 et 500Hz) je devrai avoir le même résultat non ?

Merci d'avance.

[chneu87]

Un petit exemple :

Dans ce cas là (image ci-dessous), ma fonction me retourne le peak max (200Hz) mais j'aimerai maintenant qu'il me retourne le peak max entre 400Hz et 800Hz (exemple 600Hz). Je n'arrive pas à réaliser cela.



Merci d'avance.

Chneu

[chneu87]

J'ai trouvé ma réponse, je la poste si besoin :

Code :

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

1
2
3
4
5
6
7
8
9
 
IFreqMin = find(abs(freq-400) < 0.08, 1, 'first'); % récupération de l'index pour la fréq min à 400Hz
IFreqMax = find(abs(freq-800) < 0.08, 1, 'last'); % récupération de l'index pour la fréq max à 800Hz
 
FTCalcul = DispFT(IFreqMin:IFreqMax);
FreqCalcul = Freq(IFreqMin:IFreqMax);
 
[Amplitude, IFreq] = max(FTCalcul);
Frequency = FreqCalcul(IFreq);

Bonne journée.

Chneu