ISSAT

ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Chapitre 3. Analyse fréquentielle des systèmes

discrets LTI

Pr. Djemal Ridha

Departement d’Electronique
ISSAT Sousse -2021-2022

3.1. Les fonctions de corrélation des signaux discrets

3.2. Décomposition en séries de Fourier des signaux discrets (DTFS)
3.3. Exemples d’auto-corrélation des signaux discret avec Matlab
3.4. Analyse de Fourier des signaux discrets (DTFT)
3.5. Exemples d’application des systèmes discrets
3.6. Transformée de Fourier des systèmes discrets

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 1

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Décomposition en séries de Fourier des signaux discrets et périodiques

Les séries de Fourier des signaux discrets (DTFS)
Comme indiqué dans le chapitre 1 un signal discret est périodique de période N
si: x[n]=x[n+N]
La période fondamentale est le plus petit entier N positif pour lequel
x[n]=x[n+N] avec 0=2/N est dite la fréquence fondamentale.
 2 
 Par exemple: j
est périodique de période N.
e  N 
n

 L’ensemble des exponentielles complexes des signaux discrets qui sont

périodiques de période N sont définis par:

 2 
jk  n
jk0n ( )  e j
 k [ n]  e  e   N
, k  0,1,2,... (3.2)
 Ces signaux ont des fréquences qui sont des multiples de 0=2/N, sont
reliées harmoniquement. Tous les signaux ont des fréquences multiples de la
fréquences dite fondamentale [Link] signaux { k[n]} sont tous périodiques de
période N : 2 2
 [n]  e j ( k  N ) . n e jk ( ) n . j 2n   [n]
kN N e N k

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 2

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Décomposition en séries de Fourier des signaux discrets et périodiques

Les séries de Fourier des signaux discrets (DTFS)
 Il y a N composantes périodiques distinctes. Cette décomposition est une
conséquence du fait que les exponentielles complexes et discrètes et qui
différent par un multiple de 2 sont identiques. En effet:  0[n]=N[n],
 1[n]=  N+1[n] et en général:  k[n] =  k+N[n]

o Le terme des séries pour k=0 est appelée la composante DC

o Les deux termes pour k=1: es composantes fondamentales ou les
premières composantes harmoniques.
o Les deux termes pour k=2 les composantes des secondes
harmoniques.
o Les deux termes pour k=M les composantes des Mieme harmoniques.

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 3

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Décomposition en séries de Fourier des signaux discrets (DTFS)

Les séries de Fourier des signaux discrets (DTFS)
 Supposons que x[n] est un signal discret et périodique de période N. La
représentation en séries de Fourier d’un tel signal s’exprime par:
 2 
jk  n
x[n]   ak  k [n]   ak e  N 
(3.3)
k n k n

x[0]   ak
k n.

x[1]   ak e j 2k / N
k n
…

x[ N  1]   ak e j 2k ( N 1) / N
k n

 Ces équations représentent N équations linéaires pour le N coefficients ak

inconnus où k varie à travers N valeurs entières successives

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 4

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Décomposition en séries de Fourier des signaux discrets (DTFS)

Les séries de Fourier des signaux discrets (DTFS)
 La série de Fourier est une série finie avec N termes. Les coefficients ak
peuvent être calculées en utilisant l’expression suivante:
 2 
1 1  jk  n
 jk 0n (3.4)
ak   x [ n ]e   x[n]e  N 

N k n N k n

 Pour k allant de 0 à N-1: .

x[n]  a00[n]  a11[n]  .....  aN 1 N 1[n]
 Pour k allant de 1 à N:
x[n]  a11[n]  a22[n]  .....  aN  N [n]
 Les coefficients ak ont la propriété qui est:

ak  N  ak

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 5

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Exemple de manipulation des signaux discret avec Matlab

 Exemple1:
% Manipulation des signaux % Signal retardé (avec délai)
N=128; f1=150; f2=450; x1d=[zeros(1,20),x1(1:N-20)];
f3=1500; fs=8000; n=0:N-1; subplot(2,3,4);
x1=sin(2*pi*(f1/fs)*n); plot(n,x1d); grid;
x2=(1/3)*sin(2*pi*(f2/fs)*n); title('Delayed x(n), [x1(n-
x3=sin(2*pi*(f3/fs)*n); 20)]');
figure(1); % Addition de signaux
subplot(1,1,1); xadd=x1+x2;
subplot(2,3,1); subplot(2,3,5);
plot(n,x1); grid; plot(n,xadd);
title('x1=sin(2*pi(f1/fs)n)'); grid;
subplot(2,3,2); title('x1(n)+x2(n)');
plot(n,x2);grid; % Multiplication de signaux
title('x2=(1/3)sin(2*pi(f2/fs) xmult=x1.*x3;
n)'); subplot(2,3,6);
subplot(2,3,3); plot(xmult);
plot(n,x3); grid; grid;
title('x3=sin(2*pi(f3/fs)n)'); title('x1*x3');

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 6

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Exemple de manipulation des signaux discret avec Matlab

 Résultats de simulation sur Matlab:
x1=sin(2*pi(f1/fs)n) x2=(1/3)sin(2*pi(f2/fs)n) x3=sin(2*pi(f3/fs)n)
1 0.4 1

0.5 0.2 0.5

0 0 0

-0.5 -0.2 -0.5

-1 -0.4 -1
0 50 100 0 50 100 0 50 100

x(n)décalé, [x1(n-20)] Addition:x1(n)+x2(n) Multiplication: x1*x3

1 1 1

0.5 0.5 0.5

0 0 0

-0.5 -0.5 -0.5

-1 -1 -1
0 50 100 0 50 100 0 50 100

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 7

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Décomposition en séries de Fourier des signaux discrets (DTFS)

Exemple 1:
Déterminons les coefficients ak de la DTFS pour un signal périodique x[n] :
x[n]
2  2 
1  jk   n
n ak   x[n]e  N 
1

-5 -1 0 1 5 N nN 
N=6 ( k allant de 0 à 5)
 2  .
1  jk  n
ak   x[n]e  6 
N n N 
  jk     jk     
1  jk  3   jk   
   
  
 1  2e  3    2 e  3   e  3    1
1
 2e
6
  6    
 
1  1 2 
 1  4 cos k    cos k
6 3 6 3 3

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 8

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Décomposition en séries de Fourier des signaux discrets (DTFS)

Exemple 2: Soit: 6 2
x[n]  sin n  sin( 3n)
5 5
Trouver ak: 6
x[n] est périodique de période N = 5 car:
5 3
2 5
1  j 3 5  n  j 3 5  n 
 2   2 

x[n]  e e   2 
2j  N 1 jk  n
  x[n]   ak e  N 
1 .
2 j k  3 k 0

ak   a k  a k 5
 1
k  3
 2 j ak 1/2j

k
-1 0 1
-1/2j

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 9

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Décomposition en séries de Fourier des signaux discrets (DTFS)

Exemple 3: Trouver les coefficients de la série de Fourier ak du signal
discret et périodique x[n] comme le montre la figure suivante:

 2 
1  jk   n
ak   x[n]e  N 
N nN 
x[n]
2
1
2 3 n
N=6 -5 -1 0 1. 5

2 2  1  jk  n
0    ak   x[n]e 0
N 6 3 N nN 
1

  2e  jk 0 2  2e  jk 0  1  2e jk 0  2e jk 0 2
6

 1  4 cos k0  4 cos 2k0 
1
6

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 10

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Décomposition en séries de Fourier des signaux discrets (DTFS)

Exemple 4: Soit le signal périodique x[n] de période N=8:
k  3k 
 Déterminer x[n] si: ak  cos  sin  
4  4 
ak  
1 jk0
e  e  jk0 2 2 
2 0   

1
2j

e j 3k0  e  j 3k0  N 8 4

 2 
1
8

 4e jk0  4e  jk0 ak 
1
 x[n]e
 jk 0n

1
 x[n]e
 jk  n
 N 

N k n N k n
 4 je j 3k0  4 je  j 3k0 

4, n  1
 4 j , n  3

x[n]  
4 j , n  3
0, si non (dansune seule periode)

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 11

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Analyse de Fourier des signaux discrets (DTFT)

Calcul de la transformée de Fourier des signaux discrets DTFT

 Nous étudions la transformée de Fourier des signaux périodiques en

premier lieu. Ensuite, nous enchainons pour les signaux apériodiques.
 Nous verrons également les propriétés de la transformée de Fourier des
signaux discrets.
Soit x[n] un signal apériodique. On pourrait construire un signal périodique
x[n] qui est égale au signal apériodique. sur une période N1 n  N2
Plus la période augmente (allant jusqu'a l’infini), plus le signal périodique
devient assez proche du signal apériodique x[n].

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 12

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Analyse de Fourier des signaux discrets (DTFT)

Calcul de la transformée de Fourier des signaux périodiques et discrets
En prenant un signal discret x[n] apériodique, on pourrait construire uns
signal ~
x [n] périodique.

N 2 x[n]
x [n]   ak e jk N n étant un signal périodique
~
k  N

1 2 n
~  jk
ak   x [ n ] e N
n
-N1 N1
N n ( N )
~
x [ n]
1 N1 2
  x [n] e jk N n
N n N 1
1  2 n
  x [n] e jk 0n 0 
N n N N -N1 N1 N

~
Comme x[n]  x [n] -N1 n  N1 , il est plus commode choisir l’intervalle de
sommation sur (N) de -N1 à N1.

Puisque x[n] = 0 en dehors de l’intervalle [-N1,N1],

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 13

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Analyse de Fourier des signaux discrets (DTFT)

Transformée de Fourier des signaux périodiques discrets DTFT
Définissons la transformée de Fourier d’un signal discret x[n]:
    x[n]e

j  jn
Xe
n

ak 
1
N

X e jk 0  avec 0 
2
N
1 
et  0
N 2
Remarque: X(ej) est appelée la transformée de Fourier des signaux
discrets. Elle est périodique de période 2. Notons que Nak correspond
à des échantillons de la DTFT en une période. Ainsi, les coefficients ak
sont proportionnels aux échantillons de x[n]. Par abus de langage, X(ej)
est appelée X(j) ou X().
1 B
~
x [ n]   X (e jk0 ) e jk0n
k  N  N - k0 +
1
~
x [ n] 
2
 X (e jk0 ) e jk0n 0
k  N  B

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 14

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Analyse de Fourier des signaux discrets (DTFT)

Transformée de Fourier des signaux périodiques discrets DTFT
Plus N augmente, 0 diminue.
N, , 0 tend vers 0
x [n]  x[n] La sommation devient une intégrale.
~

TF
Xe 
j

 x[n]e jn
n  

TF inverse X (e j )e jnd
1
x[n]  
2 2
 Il y a une différence frappante entre la Transformée de Fourier et
celle en série de Fourier qui est la périodicité de 2 pour la fréquence
numérique . Soit:e j ( 2n)  e j .En se basant sur la périodicité de la TFD,
nous utilisons la transformée en série de Fourier pour des fonctions
périodiques pour exprimer la transformée de Fourier Inverse. En utilisant:
1 
x[n]   ak e jk0n
2 k  

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 15

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Exemples d’application des systèmes discrets

Transformée de Fourier des signaux périodiques discrets DTFT

Exemple 1: Trouver le DTFT de la séquence x[n]:

x[n]  a nu[n] a 1

X e  
 
j n  jn
a e   (ae  j ) n
n  n 0
1

1  ae  j

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 16

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Exemples d’application des systèmes discrets

Exemple 2:
Trouver la transformée de Fourier de la séquence:
n
x[n]  a a 1

j
 a e  j n
n
X (e )[n] 
n  
 1
n  j n
a e   a  n e  j n
n 0 n  

 ae    ae  
 n  n
j
 j
1
n 0 n 0
1 1
  j
 j
1
1  ae 1  ae
1  a2

1  2a cos  a 2
Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 17
 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Exemples d’application des systèmes discrets

Exemple 3: On considère le signal rectangulaire dans le domaine fréquentiel

On se propose de calculer signal discret x[n]
X(ej)


- -B B 

1 j
x[n]   X (e )e jn d
2 2
B
1 j n
  e d
2 B
sin Bn

n

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 18

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Exemples d’application des systèmes discrets

Cas des signaux périodiques
j0n
On considère la transformée de Fourier suivante d’un signal discret x[n]  e :

X (e j )   2 (  0  2k ) X (e j )
k 
2 …
…

0 -4 0 -2

0 0 0 +2 0 +4
Transformée de Fourier discrète de: x[n]  e j0n

Afin de vérifier la validité de cette expression, calculons le TF inverse de

X(ej), Soit: x[n]  1 X (e j )e jn d
 2 2

 2 (  0  2k )e jn d
1
 
2 2 k 

 e j (0  2r )n  e j0n

Remarque: Dans un intervalle de 2, on a une seule impulsion dans la
sommation.

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 19

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Exemples d’application des systèmes discrets

Application de la transformée de Fourier des signaux périodiques

On considère maintenant une séquence périodique x[n] avec une période N
et ayant la représentation suivante en série de Fourier, Soit :
x[n]   ak e jk(2 / N )n
k  N 
x[n] est périodique de période N et de fréquence 0=2/N
La TF d’un tel signal discret et périodique est donnée par:
  2k 
X (e j )   2 ak     
k   N 
En prenant l’intervalle pour k allant de 0, 1, …. N-1, on a:

x[n]  a0  a1e j (2 / N )n  a2e j 2(2 / N )n  .... a N 1e j ( N 1)(2 / N )n

x[n] est une combinaison linéaire des signaux avec 0= 0, 2/N, 4/N, … ,
(N-1)2/N

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 20

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Exemples d’application des systèmes discrets

Exemple 4: Soit x[n] une signal périodique discret: x[n]  cos(0 n)

Cherchons alors X(ej)
1 j0n 1  j0n
x[n]  e  e
2 2

X (e j )    ( (  0  2k )   (  0  2k ))
k 
Remarque: x[n] est périodique avec une fréquence du fondamentale

0=2/N , on a ainsi: x[n]   ak e jk0n X (e j )   2ak  (  k0 )
k  N  k 
j
Sur un intervalle -: X (e )   (  0 )   (  0 )
Pour , on a:
2 2
-: X (e j )   (  )   (  )
5 X (e j )
5

… …
-2 2

-2-0 -2+0 -0 0 0 2-0 2+0

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 21

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Exemples d’application des systèmes discrets

Application de la transformée de Fourier des signaux périodiques
Soit x[n] un signal périodique défini par:

x[n]    (n  kN )
k  
x[n] est périodique de période N et de fréquence 0=2/N
 2 
 jk n
1
ak   x[n]e  N 
N n N 
En prenant l’intervalle de sommation pour 0 k  N-1, on a: ak =1/N

2  2k 
X (e j )      
k  N  N 
… x[n]

-N 0 2N
… n
N
X (e j )

… 2
N

0
… 

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 22

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Exemples d’application des systèmes discrets

Exemple 2: On considère le signal x[n]   [n]


X  e    x[n]e  jn  1
 j 
  n  
Ceci veut dire que la représentation de la transformée de Fourier résulte
d’une contribution égale des différentes composantes des fréquences.

j
Exemple 3: Soit: X (e )   2  (  0  2k )
k  
1 j jn
Cherchons x[ n]  2 X (e )e d
2
1  jn
 2 k  2 (  0  2k )e d
2
 e j (0 2r ) n
 e j0n
0 n’est pas nécessairement un nombre rationnel. Ainsi, e j0n peut être
périodique ou non.

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 23

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Transformée de Fourier des systèmes discrets

TF des signaux périodiques
Rappelons qu’un system discret LTI s’exprime par:

y[n]   h[k ]x[n  k ]
k  

j0 n
Supposons que: x[n]  e

y[ n] 

 h[k ]e
j0 ( n  k )    nh[n]e jn
He j 

k  
 Réponse fréquentielle
j0n  j0k
e  h[k ]e
k  
 e j0n H (e j )

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 24

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Quelques propriétés
Décalage en temps et en fréquence
SF
TF
x[n]  X ( ) x[n]  ak
TF
 j n0
SF  jk n0
x[n  n0 ]  X ( )e x[n  n0 ]  ak e
SF
x[n]  X ( ) x[n]  ak
x[n]e j0n  X (  0 ) jk  n SF
x[n]e 0 0  ak  k
0
Sommation et différentiation
n
Sommation : x[n]   y[i]
i 
Differentiation : y[n]  x[n]  x[n  1]

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 25

 ISSAT
ISSAT Sousse–
Kairouan–
Département
Département
Génie
Génie
électronique
électrique 2020-2021
2021-2022

Quelques propriétés
Convolution
Des propriétés similaires pour la convolution des signaux non périodiques
s’expriment par:
TF
si x[n]  X (e j )
TF
et h[n]  H (e j )
TF
y[n]  x[n]  h[n]  Y (e j )  X (e j )  H (e j )
Si x1[n] et x2[n] sont deux signaux périodiques discrets ayant une
période commune N, alors on peut définir::
y[n]   x1[n]x2[n  k ]  x1[n]  x2[n]
kN 
Ayant: ak les coef. de la série de Fourier de x1[n]
bk les coef. de la série de Fourier de x2[n]
ck les coef. de la série de Fourier de y[n]
Dans ce cas: ck =N ak bk

Traitement Numérique du signal Pr. Ridha Jemal Chapitre 3– page: 26