Transformée De Fourier Python
Thursday, 1 August 2024ylabel ( r "Amplitude $X(f)$") plt. title ( "Transformée de Fourier") plt. subplot ( 2, 1, 2) plt. xlim ( - 2, 2) # Limite autour de la fréquence du signal plt. title ( "Transformée de Fourier autour de la fréquence du signal") plt. tight_layout () Mise en forme des résultats ¶ La mise en forme des résultats consiste à ne garder que les fréquences positives et à calculer la valeur absolue de l'amplitude pour obtenir l'amplitude du spectre pour des fréquences positives. L'amplitude est ensuite normalisée par rapport à la définition de la fonction fft. # On prend la valeur absolue de l'amplitude uniquement pour les fréquences positives X_abs = np. abs ( X [: N // 2]) # Normalisation de l'amplitude X_norm = X_abs * 2. 0 / N # On garde uniquement les fréquences positives freq_pos = freq [: N // 2] plt. Transformée de fourier python online. plot ( freq_pos, X_norm, label = "Amplitude absolue") plt. xlim ( 0, 10) # On réduit la plage des fréquences à la zone utile plt. ylabel ( r "Amplitude $|X(f)|$") Cas d'un fichier audio ¶ On va prendre le fichier audio suivant Cri Wilhelm au format wav et on va réaliser la FFT de ce signal.
Transformée De Fourier Python Pour
array ([ x, x]) y0 = np. zeros ( len ( x)) y = np. abs ( z) Y = np. array ([ y0, y]) Z = np. array ([ z, z]) C = np. angle ( Z) plt. plot ( x, y, 'k') plt. pcolormesh ( X, Y, C, shading = "gouraud", cmap = plt. cm. hsv, vmin =- np. pi, vmax = np. Transformée de fourier python program. pi) plt. colorbar () Exemple avec cosinus ¶ m = np. arange ( n) a = np. cos ( m * 2 * np. pi / n) Exemple avec sinus ¶ Exemple avec cosinus sans prise en compte de la période dans l'affichage plt. plot ( a) plt. real ( A)) Fonction fftfreq ¶ renvoie les fréquences du signal calculé dans la DFT. Le tableau freq renvoyé contient les fréquences discrètes en nombre de cycles par pas de temps. Par exemple si le pas de temps est en secondes, alors les fréquences seront données en cycles/seconde. Si le signal contient n pas de temps et que le pas de temps vaut d: freq = [0, 1, …, n/2-1, -n/2, …, -1] / (d*n) si n est pair freq = [0, 1, …, (n-1)/2, -(n-1)/2, …, -1] / (d*n) si n est impair # definition du signal dt = 0. 1 T1 = 2 T2 = 5 t = np. arange ( 0, T1 * T2, dt) signal = 2 * np.
Pour remédier à ce problème, on remplace la fenêtre rectangulaire par une fenêtre dont le spectre présente des lobes secondaires plus faibles, par exemple la fenêtre de Hamming: def hamming(t): return 0. 54+0. 46*(2**t/T) def signalHamming(t): return signal(t)*hamming(t) tracerSpectre(signalHamming, T, fe) On obtient ainsi une réduction de la largeur des raies, qui nous rapproche du spectre discret d'un signal périodique.