Traitement Numérique du Signal
EXERCICES CORRIGÉS
Transformée de Fourier Discrète
Février 2022
A.ALLICHE
alliche_univebejaia@gmail,com
LAT 312 L3 Télécomunications
EXERCICE 1:
Trouvez la TFD de la séquence suivante,
x(n) = { ¼, ¼, ¼ }
Solution :
Produit
1. On détermine de Convolution
la longueur de la séquence, N
N = 3, k = 0,1,2
2. On utilise l’expression de la TFD X(k)
N −1
X(k) = ∑ x(n) e-j2πnk/N ,
n =0
k = 0, 1, 2
X(0) = ¼ + ¼ + ¼ = ¾
X(1) = ¼ + ¼e-j2π/3 + ¼e-j4π/3
= ¼ + ¼ [cos (2π/3) – jsin(2π/3) + ¼ [cos (4π/3) –
jsin(4π/3)
= ¼ + ¼ [-0.5 – j0.866] + ¼ [ -0.5 + j0.866]
= ¼ + ¼ [-1] = 0
LAT 312 L3 Télécomunications alliche_univebejaia@gmail,com
X(2) = ¼ + ¼ e-j4π/3 + ¼ e-j8π/3
= ¼ + ¼ [cos (4π/3) – jsin(4π/3)] +
¼ [cos (8π/3) – jsin(8π/3)]
= ¼Produit
+ ¼ [-0.5 +j0.866] + ¼[-0.5 –
de Convolution
j0.866]
=0
D’où,
X(k) = { ¾, 0, 0}
LAT 312 L3 Télécomunications alliche_univebejaia@gmail,com
EXERCICE 2: x(n) = {1,1,2,2,3,3}
Solution :
Produit
N −1
de Convolution
X(k) = ∑ x ( n) e
n =0
-j2πnk/N , k = 0,1,2,3,4,5
X(0) = 12, X(1) = -1.5 + j2.598
X(2) = -1.5 + j0.866, X(3) = 0
X(4) = -1.5 – j0.866, X(5) = -1.5 – j2.598
,
X(k) = {12, -1.5 + j2.598, -1.5 + j0.866, 0, -1.5 –
j0.866,
-1.5 – j2.598}
LAT 312 L3 Télécomunications alliche_univebejaia@gmail,com
EXERCICE 3: Calculez le Produit de convolution
x1(n) = {2, 1, 2, 1} & x2(n) = {1, 2, 3, 4}
Solution :
N= 4
(i) X1(0) = 6, X1(1) = 0, X1(2) = 2, X2(3) = 0
X1(k) = {6,0,2,0}
(ii) X2(0) = 10, X2(1) = -2+j2, X2(2) = -2, X2(3) = -2-j2
X2(k) = {10,-2+j2,-2,-2-j2}
3. Le produit deconvolution :
X3(k) = X1(k) X2(k) = {60, 0, -4, 0}
LAT 312 L3 Télécomunications alliche_univebejaia@gmail,com
TFDI: Transformée de Fourier Inverse
TFDI est la transformée de Fourier discrète inverse.
La séquence de longueur finie peut être obtenue à partir de
la transformée de Fourier discrète en effectuant TFDI.
Produit de Convolution
Définie par :
N −1
x(n) = 1/N ∑ X (k ) ej2πnk/N,
k =0
n = 0,1, …, N-1
LAT 312 L3 Télécomunications alliche_univebejaia@gmail,com
EXEMPLE 4:
Déterminer la TFDI de la séquence,
X(k) = {1, 2, 3, 4}
SOLUTION :
Produit de Convolution
N −1
x(n) = 1/4 ∑ X (k )
k =0
ej2πnk/4,
x(0) = ¼(1 + 2 + 3 + 4) = 5/2
x(1) = -0.5 – j0.5, x(2) = -0.5
x(3) = -0.5 + j0.5
x(n) = {2.5, -0.5-j0.5, -0.5, -0.5+j0.5}
LAT 312 L3 Télécomunications alliche_univebejaia@gmail,com
EXEMPLE 5:
Determiner la TFDI de la séquence,
X3(k) = {60, 0, -4, 0}
Produit de Convolution
SOLUTION :
x3(n) = 1/4 ∑ X (k ) ej2πnk/4,
3
k =0
x3(0) = 14, x3(1) = 16, x3(2) = 14, x3(3) = 16
x3(k) = {14, 16, 14, 16}
LAT 312 L3 Télécomunications alliche_univebejaia@gmail,com
EXEMPLE 6: x[n] = anU [n] a <1
SOLUTION :
Produit de Convolution
∞ − j 2 kπ
n ∞ − j 2 kπ
n
=x[k ] ∑
= x(nT )e
n = −∞
N
∑ U s ( n )e
a n
n = −∞
N
∞ n ∞ 2 kπ
− j 2 kπ −j
x[k ] = ∑ a ne N = ∑ (ae N )n
n=0 n=0
1 k = 0 , 1 , 2 , 3 , …………, N-1
x[ k ] =
2 kπ
−j
1 − ae N
LAT 312 L3 Télécomunications alliche_univebejaia@gmail,com
EXEMPLE 7: x[n] = anU [n+3] a <1
SOLUTION :
Produit
∞
de
− j 2kπ
n Convolution
∞ − j 2kπ
n
x[k ] = ∑ x(nT )e N =
∑ a nU s (n + 3)e N
n = −∞ n = −∞
∞ n ∞ 2 kπ −1 2 kπ
− j 2 kπ −j −j
x[k ] = ∑ e
a n N = ∑ (ae N )n +
∑ (ae N )n
n = −3 n=0 n = −3
a −3e j 6 kπ
x[k ] = 2 kπ
k = 0 , 1 , 2 , 3 , …………, N-1
−j
1 − ae N
LAT 312 L3 Télécomunications alliche_univebejaia@gmail,com
EXEMPLE 8 Transformée de Fourier Rapide
TFD de 2 échantillons “N=2”
1 1
X (0) = ∑ x(n)W2 = ∑ x(n) = x(0) + x(1)
1 n0
X (k ) = ∑ x(nProduit
)W2 nk k = 0,1
de Convolution n=0 n=0
n=0
1 1
X (1) = ∑ x ( n )W 2
n1
= ∑ x ( n )W 2
n
= x ( 0 )W 2 0 + x (1)W 2 1
n=0 n=0
= x ( 0 ) + x (1)W 2 (1 / 2 ) 2 = x ( 0 ) + x (1)( − 1) = x ( 0 ) − x (1)
X [0] = x[0] + x[1]
X [1] = x[0] − x[1]
Graphe de Fluence de FFT 2 points
LAT 312 L3 Télécomunications alliche_univebejaia@gmail,com
EXEMPLE 09: Transformée de Fourier Rapide
TFD de 4 échantillons “N=4”
3
X (k ) = ∑Produit
x(n)W4 nk k = 0Convolution
de ,1,2,3,
n=0
3 3
X (0) = ∑ x(n)W4 n0
= ∑ x(n) = x(0) + x(1) + x(2) + x(3)
n=0 n=0
3
X (1) = ∑ x(n)W4 = x(0)W4 + x(1)W4 + x(2)W4 + x(3)W4
n 0 1 2 3
n =0
= x(0) − jx(1) − x(2) + jx(3)
LAT 312 L3 Télécomunications alliche_univebejaia@gmail,com
EXEMPLE 09: Transformée de Fourier Rapide
TFD de 4 échantillons “N=4”
3
X (2) = ∑ x(n)W4 = x(0)W4 + x(1)W4 + x(2)W4 + x(3)W4
2n 0 2 4 6
n =0
Produit de Convolution
= x(0) + x(1)(−1) + x(2)(1) + x(3)W4
2
= x(0) − x(1) + x(2) − x(3)
3
X (3) = ∑ x(n)W4 = x(0)W4 + x(1)W4 + x(2)W4 + x(3)W4
3n 0 3 6 9
n =0
= x(0) + x(1)W4 + x(2)(1)W4 + x(3)W4
3 2 1
X (0) = [ x(0) + x(2)] + [ x(1) + x(3)]
= x(0) + jx(1) + (−1) x(2) + (− j ) x(3) X (1) = [ x(0) − x(2)] + (− j )[ x(1) − x(3)]
= x(0) + jx(1) − x(2) − jx(3) X (2) = [ x(0) + x(2)] − [ x(1) + x(3)]
X (3) = [ x(0) − x(2)] + j[ x(1) − x(3)]
LAT 312 L3 Télécomunications alliche_univebejaia@gmail,com
EXEMPLE 10: Transformée de Fourier Rapide
TFD de 8 échantillons “N=8”
ω4= eiπk ω2= eiπk/2 ω= eiπk/4
x0 x0 Produit
x +x
0 4
de Convolution
x +x +x +x x +x +x +x +x +x +x +x
0 4 2 6 0 4 2 6 1 5 3 7
F2
x1 x4 x0–x4 x0-x4+i(x2-x6) x0-x4+i(x2-x6)+eiπ/4 (x1-x5+i(x3-x7))
F4
x2 x2 x2+x6 x0+x4-(x2+x6) x0+x4-x2-x6+i(x1+x5-x3-x7)
x3 x6 x2–x6 x0-x4-i(x2-x6) x0-x4-i(x2-x6)+ eiπ3/4(x1-x5-i(x3-x7))
F8
x4 x1 x1+x5 x1+x5+x3+x7 x0+x4+x2+x6-(x1+x5+x3+x7)
x5 x5 x1–x5 x1-x5+i(x3-x7) x0-x4+i(x2-x6)-eiπ/4 (x1-x5+i(x3-x7))
x6 x3 x3+x7 x1+x5-(x3+x7) x0+x4-x2-x6-i(x1+x5-x3-x7)
x7 x7 x3–x7 x1-x5-i(x3-x7) x0-x4-i(x2-x6)- eiπ3/4(x1-x5-i(x3-x7))
Etage =1 Etage =2 Etage =3
LAT 312 L3 Télécomunications alliche_univebejaia@gmail,com
Graphe de Fluence de FFT 8 points
Produit de Convolution
LAT 312 L3 Télécomunications alliche_univebejaia@gmail,com
Graphe de Fluence de FFT 8 points
Produit de Convolution
LAT 312 L3 Télécomunications alliche_univebejaia@gmail,com
Graphe de Fluence de FFT/IFFT 4 points
Produit de Convolution
LAT 312 L3 Télécomunications alliche_univebejaia@gmail,com
Graphe de Fluence de FFT 8 points {x[n]} = {4, - 3, 2, 0, -1, - 2, 3, 1}
Produit de Convolution
{X[k]} = {4, 5 + j + j 2 , - 2 + 6j, 5 - j + j 2, 12, 5 + j - j 2 ,−2 + 6j, 5 - j - j 2}
LAT 312 L3 Télécomunications alliche_univebejaia@gmail,com
TFD Approche Matricielle
N −1 2π kn
X (k )
= ∑ x ( n )W
n =0
=
nk
N k 0,1,..., N − 1 W nk = e
N
−j
N
Produit X=Wx
de Convolution
X ( 0) W N
0
W N
0
W 0
N ... WN0 x ( 0)
0 ( N −1)
X (1 ) N
W W N
1
W N
2
... WN x ()
1
X= X ( 2 ) W= WN0 W 2
W 4
...
2( N −1) x= x 2
WN ( )
N N
... ... ... ... ... ... ...
X ( N − 1) x ( N − 1)
WN0 WN( N −1) WN2( N −1) ( N −1)
2
... WN
LAT 312 L3 Télécomunications alliche_univebejaia@gmail,com
TFD Approche Matricielle
EXEMPLE 11 : x(n)={1,1,0,0}
Produit de Convolution
1 1 1 1 1 2
1 j −1 − j 1 1 + j
X=Wx= =
1 −1 1 −1 0 0
1 − j −1 j 0 1 − j
LAT 312 L3 Télécomunications alliche_univebejaia@gmail,com