23/03/2023

Transformée de Fourier Discrète

(DFT)

Chapitre 3

1
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

Rappel
 Transformée de Fourier:
 Signal analogique:
La transformée de Fourier du signal x(t ) :


 x(t ) e
 j t
 X ( )  dt


où   2 f est la fréquence radiale (en radian/seconde).

• Le module X ()  spectre du signal x(t )

• La phase  ()

1
 X ( ) e
j t
 x(t )  d
2 

2
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

1
 23/03/2023

Transformée de Fourier à temps discret (TFTD)

 Transformée de Fourier:
 Signal discret:

La transformée de Fourier à temps discret (TFTD) (ou DTFT)

d’une séquence x(n) :


 X ( w)   x ( n)
n 
e  j wn

X ( w) fonction continue et périodique de 2


1
 x ( n) 
2  X (w)

e j wn dw

3
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

Transformée de Fourier à temps discret (TFTD)

 Transformée de Fourier:

 X(f )   x ( n) e
n 
 2 j nf

X ( f ) fonction continue et périodique de 1

1
2
 x ( n)  
1 2
X ( f ) e 2 j nf df

Ou bien

 X (e jw )   x(n) e  j wn
n 

X (e ) : est périodique avec w  2

jw

 X (e jw )  X re (e jw )  jX im (e jw )  X (e jw ) e j ( w) ,  ( w)  arg X (e jw ) .
4
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

2
 23/03/2023

Transformée de Fourier à temps discret (TFTD)

 Exemple 1:
On considère une séquence de durée finie :

 x(n)  1, 2,3, 4,5


 X (e jw )  e jw  2  3e jw  4e2 jw  5e3 jw

 Exemple 2:

On considère une séquence causale:

 x(n)  (0.5) n u (n)

1
 X (e jw ) 
1  0.5 e jw

5
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

Transformée de Fourier à temps discret (TFTD)

 Exemple 3:
On considère une fonction porte rectangle:
1 si n  0,..., N  1
 x(n)  rect N (n)  u (n)  u (n  N )  
0 sin on

N 1
 X ( f )   e 2 j nf  1  ...  e 2 j ( N 1) f
n 0

N pour f  0

 1  e 2 j Nf
 j ( N 1) f sin( N  f )
pour f  0
 1  e2 j f  e sin( f )


6
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

3
 23/03/2023

Transformée de de Fourier Discrète (TFD)

 Définition
• Considérons une séquence finie: x(n)   x(0),..., x( N  1)
sa TFTD s ' écrit :

N 1
 X ( f )   x(n) e 2 j nf , où f   0,1
n 0

• Soit N le nombre de valeurs de f pour calculer la DFT.

7
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

Transformée de de Fourier Discrète (TFD)

• Soit N le nombre de valeurs de f pour calculer la DFT.

• On prend N valeurs uniformément réparties sur l’intervalle  0,1 :

k
f  , avec k  0,..., N  1
N

Alors :
N 1
k
 X( )   x(n) e 2 j nk / N
N n 0

8
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

4
 23/03/2023

Transformée de Fourier discret (TFD)

 Définition
 On appelle TFD de la suite x(n), n  0,..., N  1 , la suite X (k ),

k  0,..., N  1 , définie par:

N 1
 X (k )   x(n) WNn k , où : WN  e 2 j / N , WNN  1
n 0

 La transformée inverse:

N 1
1
 x ( n) 
N
 X (k ) W
k 0
n k
N , k  0,...N  1. WNN  1

9
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

Transformée de Fourier discrète (TFD)

 Exemples: TFD
 On considère une séquence de longueur N:
1, n  0
 x ( n)  
0, ailleurs
N 1
 X (k )   x(n) WNn k  1
n 0

 On considère une séquence de longueur N:

1, n  m
 y ( n)  
0, ailleurs

 Y (k ) WNkm
10
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

5
 23/03/2023

Transformée de Fourier discrète (TFD)

 Présentation matricielle: TFD

 La TFD définie para:

N 1
 X (k )   x(n) WNn k , où : WN  e 2 j / N , WNN  1
n 0

 X   X (0) X (1) ... X ( N  1) 

t

 x   x(0) x(1) ... x( N  1) 

t

 X  DN x, DN : est la matrice TFD N X N

11
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

Transformée de Fourier discrète (TFD)

 Présentation matricielle: TFD
• Forme matricielle:  X  DN x

1 1 1 ....... 1 
 1 2 N 1 
1 W W N N....... W N 
1 WN2 WN4 ....... WN2( N 1) 
 
. 
 DN 
.


. 
 
. 
. 
 N 1 2( N 1) ( N 1)( N 1)


1 WN WN ....... WN 


12
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

6
 23/03/2023

Transformée de Fourier discrète (TFD)

 Présentation matricielle: TFD inverse (IDFT)
 La TFD inverse définie para:
N 1
1
 x ( n) 
N
 X (k ) W
n 0
n k
N , n  0,...N  1.
• Forme matricielle:
 x(0)   X (0) 
 x(1)   X (1) 
   
1 1 1 ... 1  .  
1 .

 1 2  ( N 1)      DN  
1 W W N N ... W N  .  . 
1 WN2 WN4 ... WN2( N 1)  .  . 
1     
1
 DN  .   x( N  1)   X ( N  1) 
N 
. 
. 
  ( N 1) 2( N 1)  ( N 1)( N 1) 
1 WN WN ... WN 
13
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

Transformée de Fourier discrète (TFD)

 Propriétés de TFD
• Linéarité:  Si x(n) 
TFD
 X (k ) et y (n) 
TFD
 Y (k )
Alors,  a x(n)  b y (n) 
TFD
 a X (k )  bY (k )
où a, b 
• Décalage temporel:
k
2 j n0
 x((n  n0 ) mod N ) 
TFD
 X (k ) e N

• Décalage fréquentiel: k0
2 j n
 x ( n) e N

TFD
 X ((k  k0 ) mod N )

• Formule de Parseval:
N 1 N 1
1
  x ( n)   X (k )
2 2

n 0 N n 0
14
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

7
 23/03/2023

Transformée de Fourier discrète (TFD)

 Réduction Modulo-N
• La TFD est évaluée sur une séquence x(n) de longueur finie: 0  n  N 1
Si la séquence x(n) est de longueur infinie ou de longueur L  N, on
divise la séquence en bloques de longueur N.
par exemple, si L  4N :

N N N N
x(n)  x0 x1 x2 x3

 x(n)  x0  x1  x2  x3

15
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

Transformée de Fourier discrète (TFD)

• On exprime les sous bloques en terme de la séquence x(n),0  n  L 1:
pour m  0, 1, ....
xm (n)  x(mN  n), n  0, 1, ,...., N  1

Où le mième bloque occupe l’intervalle:  mN , ( m  1) N 

 x(n)  x0  x1  x2  x3  ...
 x(n)  x( N  n)  x(2 N  n)  ...

ou bien, x(n)   x(mN  n), n  0, 1, ...., N  1
m 0

• De manière plus générale:


 x ( n)   x(mN  n),
m
  n  

 on note: x(n) est périodiquede N .

16
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

8
 23/03/2023

Transformée de Fourier discrète (TFD)

 Exemple: x(n)  1, 2,  2,3, 4,  2,  1,1 T

Trouver x(n) pour N  4 et N  3

 pour N  4 :
 1   4  5 
 2   2  0 
x ( n)         
 2   1   3
     
 3  1   4 
 pour N  3 :
 1  3   1 3 
x(n)   2    4   1    7 
 2   2  0   4 
 On note: on ajoute un zéro au troisième sous bloque pour N=3
17
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

Transformée de Fourier discrète (TFD)

 Convolution circulaire
• La convolution circulaire est notée par:
 z(n)  x(n)  y(n)

Exemple: x(n)  1, 2,3,1 , y ( n)  1, 0,1,1

z (n)  x(n)  y (n)  6, 6,5, 4

• la suite Z(k) = X(k) Y(k) a TFD inverse:

1
 z ( n)  X (k )  Y (k )
N

18
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

9
 23/03/2023

Transformée de Fourier discrète (TFD)

 Exercice

Soit x(n)  1,1, 2, 2,3,3 et h(n)  2,1, 2

Calculer:
Produit de convolution y(n) et en déduire le produit de convolution
circulaire y(n) mod3

19
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

Transformée de Fourier Rapide (FFT)

 FFT: Fast Fourier Transform

• Publiée pour la première fois en 1965 par Cooly et Tucky, est

une technique de calcul rapide de la TFD.
• La TFD de x(n) :
N 1
 X (k )   x(n) WNn k , k  0,1,....,N  1
n 0

avec : WN  e2 j / N

 Propriétés de WN :
 WNk n  WNk ( n N )  WN( k  N ) n
k n N
 WN 2
 WNk n
 WN2 k n  WNn k
2
20
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

10
 23/03/2023

Transformée de Fourier Rapide (FFT)

 Algorithme de la FFT
• On considère que la longueur de la séquence x(n) est une
puissance de 2: N  2 B  B  log ( N )
2

• La TFD de x(n) peut s’écrire comme la somme d’un terme lié

aux indices de rang pair et d’un terme lié aux indices de rang
impair:

N 1 N 1
2 2
 X (k )   x(2n) W
n 0
N
k (2 n )
  x(2n  1) W
n 0
N
k (2 n 1)

21
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

Transformée de Fourier Rapide (FFT)

 Algorithme de la FFT
On pose :
 g (n)  x(2n),
N
 h(n)  x(2n  1), n  0, 1, . . ., 1
2
N 1
2
 G (k )  W
n 0
kn
N
2
g (n)
N 1
2 N
 H (k )  W
n 0
kn
N
2
h(n), k  0, 1, . . .,
2
1

 Les deux termes (pair et impair) de X(k) peuvent être

exprimés en termes de G(k) et H(k):
22
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

11
 23/03/2023

Transformée de Fourier Rapide (FFT)

 Algorithme de la FFT

 On note que WN et WN
2
sont reliés par:
WN  e  j 2 /( N /2)  e  j 4 / N  WN2
2

On obtient :N 21 N 1
2
 X (k )   g ( n) W
n 0
N
kn
2
W k
N  h( n) W
n 0
N
kn
2

 X (k )  G (k )  WNk H (k ), k  0,1,...., N  1
 On écrit ces N équations comme deux groupes de N/2
équations:
 X (k )  G (k )  WNk H (k )
(kN )
 X (k  N )  G (k  N )  WN 2
H (k  N ), k  0,1,...., N  1
2 2 2 23
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

Transformée de Fourier Rapide (FFT)

 Algorithme de la FFT

• La périodicité de DFT: G (k  N )  G (k ) et H (k  N )  H (k )
2 2

2  j
Et la propriété:
N N
WN 2
 (e N
) 2
 e j    1

D’où l’algorithme de la FFT:

 X (k )  G (k )  WNk H (k ),
 X (k  N )  G (k )  WNk H (k ).
2
k  0,1,...., N  1
2

24
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

12
 23/03/2023

TFD &FFT
 Exercices (voir TD-Série N°4 + exercice supplémentaire )

25
Université Mohammed V de Rabat, FSR, SMP6 - Parcours 1, 2022-2023, Pr. Lh. Masmoudi

13