Exemples
Exemple #1:
// Création d'un signal t = 0:0.1:1000; x = 3*sin(t)+8*sin(3*t)+0.5*sin(5*t)+3*rand(t); // Calcul de sa fft y = fft(x, -1); // Inversion du résultat shift = fftshift(y); // Inversion du résultat inversé invShift = ifftshift(shift); // On vérifie qu'on a bien recréé le résultat de départ and(y == invShift) |  |  |
Example #2:
// Création d'une image 2D t = 0:0.1:30; x = 3*sin(t')*cos(2*t)+8*sin(3*t')*sin(5*t)+.. 0.5*sin(5*t')*sin(5*t)+3*rand(t')*rand(t); // Calcul de la fft y = fft(x,-1); // Inversion du résultat shift = fftshift(y); // Inversion du résultat inversé invShift = ifftshift(shift); // On vérifie qu'on a bien recréé le résultat de départ and(y == invShift) | | |