cos ( 2 * np. pi / T1 * t) + np. sin ( 2 * np. pi / T2 * t) # affichage du signal plt. plot ( t, signal) # calcul de la transformee de Fourier et des frequences fourier = np. fft ( signal) n = signal. size freq = np. fftfreq ( n, d = dt) # affichage de la transformee de Fourier plt. plot ( freq, fourier. real, label = "real") plt. imag, label = "imag") plt. legend () Fonction fftshift ¶ >>> n = 8 >>> dt = 0. 1 >>> freq = np. fftfreq ( n, d = dt) >>> freq array([ 0., 1. 25, 2. 5, 3. 75, -5., -3. 75, -2. 5, -1. 25]) >>> f = np. fftshift ( freq) >>> f array([-5., -3. 25, 0., 1. Transformée de fourier python programming. 75]) >>> inv_f = np. ifftshift ( f) >>> inv_f Lorsqu'on désire calculer la transformée de Fourier d'une fonction \(x(t)\) à l'aide d'un ordinateur, ce dernier ne travaille que sur des valeurs discrètes, on est amené à: discrétiser la fonction temporelle, tronquer la fonction temporelle, discrétiser la fonction fréquentielle.
append ( f, f [ 0]) # calcul d'une valeur supplementaire z = np. append ( X, X [ 0]) Exemple avec translation ¶ x = np. exp ( - alpha * ( t - 1) ** 2) ( Source code)
Exemples simples ¶ Visualisation de la partie réelle et imaginaire de la transformée ¶ import numpy as np import as plt n = 20 # definition de a a = np. zeros ( n) a [ 1] = 1 # visualisation de a # on ajoute a droite la valeur de gauche pour la periodicite plt. subplot ( 311) plt. plot ( np. append ( a, a [ 0])) # calcul de A A = np. fft. fft ( a) # visualisation de A B = np. append ( A, A [ 0]) plt. subplot ( 312) plt. real ( B)) plt. ylabel ( "partie reelle") plt. Transformée de Fourier. subplot ( 313) plt. imag ( B)) plt. ylabel ( "partie imaginaire") plt. show () ( Source code) Visualisation des valeurs complexes avec une échelle colorée ¶ Pour plus d'informations sur cette technique de visualisation, voir Visualisation d'une fonction à valeurs complexes avec PyLab. plt. subplot ( 211) # calcul de k k = np. arange ( n) # visualisation de A - Attention au changement de variable plt. subplot ( 212) x = np. append ( k, k [ - 1] + k [ 1] - k [ 0]) # calcul d'une valeur supplementaire z = np. append ( A, A [ 0]) X = np.