Transforme de Fourier discrte et et transforme de

Transformée de Fourier discrète et et transformée de Fourier rapide Thomas LAMOTTE 1

La transformée de Fourier discrète l Signal analogique – La transformée de Fourier X( ) d ’un signal analogique x(t) est donnée par : l t représente le temps et f les fréquences – C ’est une opération de projection de x(t) sur l ’exponentielle. l Signal discrétisé et périodisé – Considérons une suite finie de N échantillons d’un signal discrétisé et périodisé. – On définit sa transformée de Fourier discrète (TFD) comme la suite Transformée de Fourier discrète et transformée de Fourier rapide 2

La transformée de Fourier discrète – De même, on définit la transformée de Fourier inverse (ITFD) par: l Interprétation vectorielle : – Les éléments de la suite peuvent être vus comme les composantes d ’un vecteur x dans un espace à N dimensions. X est alors une combinaison linéaire de N vecteurs de base wk où les composantes de chaque vecteur wk sont donnés par la suite – Exemple : l Considérons une TFD sur 16 valeurs, on peut tracer les parties réelles des composantes des cinq premier vecteurs de base wk Transformée de Fourier discrète et transformée de Fourier rapide 3

La transformée de Fourier discrète – La décomposition peut alors être exprimée sous la forme matricielle : l X = Wx. T où x. T désigne la transposée de x Transformée de Fourier discrète et transformée de Fourier rapide 4

La transformée de Fourier rapide l Complexité d ’une TFD – Pour la TFD, il y a l l N² multiplications complexes N(N-1) additions complexes – Les multiplications complexes ont une durée d ’exécution beaucoup plus longue que les additions. l Algorithmes de transformée de Fourier rapide (TFR) ou Fast Fourier Transform (FFT) – Dans les algorithme de transformée de Fourier rapide, le nombre d ’opération est considérablement réduit. – Il en existe plusieurs : l l TFR avec entrelacement temporel TFR avec entralacement temporel TFR en base 4 TFR en base double – Le plus connu est l’algorithme de Cooley. Tukey. Transformée de Fourier discrète et transformée de Fourier rapide 5

La transformée de Fourier rapide l Algorithme de Cooley-Tukey – La TFD peut s’écrire en séparant indices pairs et impairs: – On aboutit à deux transformées de longueur N/2. Xp correspond à la transformée des indices d’échantillons pairs et Xi à celle des indices impairs. – Il est possible de subdiviser encore Xp en Xpp et xpi en séparant les indices pairs et impairs et de même pour les indices impairs xi en xip et xii. Il est possible de réitérer plusieurs fois cette méthode. Transformée de Fourier discrète et transformée de Fourier rapide 6

La transformée de Fourier rapide l Exemple : – Prenons 8 échantillons, qui ont les valeurs successives suivantes: x 0, x 1, x 2, x 3, x 4, x 5 x 6, x 7, x 8. – La TFD se présente ainsi : Transformée de Fourier discrète et transformée de Fourier rapide 7

La transformée de Fourier rapide – Prenons la ligne X 1, en séparant les échantillons paires/impairs, puis en factorisant les échantillons impairs, on obtient: – Les termes des couples x 0 -x 1, x 2 -x 3, x 4 -x 5, x 6 -x 7 sont identiques à un facteur près, donc on peut encore subdiviser : – On peut montrer facilement, l a) l b) – D’où l’écriture suivante : Transformée de Fourier discrète et transformée de Fourier rapide 8

……. La transformée de Fourier rapide – Le calcul de certains terme revient plusieurs fois, on peut donc diminuer le nombre d’opérations à réaliser à l’aide d’ opérations « butterfly » – l Opérations « Butterfly » ou papillon – Le calcul « Butterfly » est le suivant : a b Transformée de Fourier discrète et transformée de Fourier rapide 9

La transformée de Fourier rapide l Application à la FFT x 0 X 0 x 4 X 1 x 2 X 2 x 6 X 3 x 1 X 4 x 5 X 5 x 3 X 6 x 7 X 7 c Signal l c Transformée de Fourier Remarque : – Ce schéma de calcul peut être implanté dans des DSP – Nombre de calculs N/2 log 2 N Transformée de Fourier discrète et transformée de Fourier rapide 10

La transformée de Fourier rapide l Incrément en « reverse carry » – Les couples d’échantillons doivent être choisis selon un ordre particulier. Cette incrémentation est appelée « reverse carry » (retenue inverse). – L’incrémentation conssite à additionner N/2 à l’indice puis à reporter la retenue à droite plutôt qu’à gauche. – Exemple: N= 8 N/2=4 soit 100 l l En retenue « non inverse » : 100 + 100 = 1000 = 8 En retenue « inverse » : 100+100 = 010 = 2 – Le 1 passe de la gauche vers la droite Exemple pour 8 échantillons – On peut aussi arranger les valeurs de fréquence selon l’incrémentation « reverse carry » : Transformée de Fourier discrète et transformée de Fourier rapide 11

La transformée de Fourier rapide x 0 X 0 x 1 X 4 x 2 X 2 x 3 X 6 x 4 X 1 x 5 X 5 x 6 X 3 x 7 X 7 c Signal c Transformée de Fourier discrète et transformée de Fourier rapide 12