1 pièce(s) jointe(s)
Calcul et tracé de la FFT
Bonjour,
Je vous propose un nouvel élément à utiliser : Calcul et tracé de la FFT
Voici la fonction plot_fft qui permet de tracé en toute simplicité la FFT (Fast Fourier Transform) d'un signal sur des fréquences positives.
Le tracé se décompose en deux subplots :
- celui du dessus qui a des ordonnées linéaires et des fréquences normalisées en abscisses ;
- celui du dessous avec des ordonnées en dB et des fréquences « normales ».
Les deux objets Axes ont un datacursor personnalisé affichant les valeurs de l'amplitude et des deux types de fréquence.
La fonction vous donne la possibilité de faire du zero padding ainsi que d'appliquer une fenêtre.
Vous pouvez récupérer en sortie le couple fréquences/FFT de la partie positive qui est affichée.
Exemple de code utilisant la fonction plot_fft :
Code:
1
2
3
4
5
6
7
8
| N = 1000;
fs = 1e4;
t = (0:N-1)/fs;
fo = 300;
f1 = 1.2e3;
x = 1.2 * sin(2*pi*fo*t) + 3 * sin(2*pi*f1*t);
[X f] = plot_fft(x, fs, 2, @hann); |
Qu'en pensez-vous ?
