TRAITEMENT NUMÉRIQUE DU SIGNAL :

TD 3 : TZ, RIF, RII

Enseignant : BEN MESSAOUD Mohamed Anouar
Maître de Conférences à FST, Master SSD, ENIT
LAB LSITI, ENIT

1
 2

Exercice 2
Soit le filtre numérique dont l’équation récurrente est donnée par :
y[n]=0.75 x[n]-2x[n-1]+x[n-2]
1) Calculer la réponse impulsionnelle h(n) de ce filtre.
yk   xk   hk   Y z   X z .H z 
H z  est la transformée en z de la réponse à l’impulsion hk 
H z  est appelée la fonction de transfert du système linéaire et invariant.

h[n]=0.75d(n)-2d[n-1]+d[n-2]

2) De quel type de filtre s’agit-il ?

C’est un filtre RIF

3) Calculer la réponse harmonique de ce filtre.

 j 2f
En remplacent z par e dans H(z) , on obtient la réponse harmonique H(f).

z = e jwTe avec Te = 0,5  z = e (jw/2)

On a H(z) = 0,75-2z-1+z-2
 H(f) = 0,75 -2 e -(jw/2) +e -(jw)
 3

Exercice 2
4) a) On associe en série à ce filtre un autre filtre numérique de façon que le comportement de
ce système soit équivalent au filtre analogique suivant :
H(p) = [(p-1)(p+1)]/ [(p+0.5)(p+3)]
Calculer la fonction de transfert du filtre recherché sachant que la période d’échantillonnage T =
0.5 s.
y(k)= [x(k)-x(k-1)]/T Y(z) =[X(z)-z-1X(z)]/T Y(z) = X(z)[1-z-1]/T
p = (1-z-1)/Te  p = 2 (1-z-1)
H(z) = [4(1-z-1)2)-1]/[(2-2z-1+0.5) (2-2z-1+3)]

= [4(1-2z-1+z-2)-1]/[(2.5-2z-1)(5-2z-1)]

= [(3-8z-1+4z-2)]/[(25/2-15z-1+4z-2)]

= [(3-8z-1+4z-2)]/[(25/2-15z-1+4z-2)]

= [(0.75-2z-1+z-2)]/[1/(25/8-(15/4)z-1+z-2)]

On a H(z) = H1(z) . H2(z)

D’où H2(z) = 1/[(25/8-(15/4)z-1+z-2)]

 4

Exercice 2

 a yk  n    b xk  m
N M

On a : n m alors
n 0 m0

a z Y z   b z  m X z 
N M
n
n m
M
m
Y z 
n 0 m0
 bm z
H z   H z   m 0
X z 
Comme alors
N
n
 an z
n 0
Par factorisation du numérateur et du dénominateur, il est possible d’écrire la
fonction de transfert en fonction des pôles et des zéros:
M
 1  zm z
1

H z   H 0 m 1

 1  p  d’où H2(z) = (8/25)/[(1-(6/5)z-1+(8/25)z-2)]

N
1
nz
n 0
4) b) En déduire son type, son ordre et ses coefficients.

C’est un filtre RII d’ordre 2 a pour coefficients

b0 = 8/25 ; a1 = 6/5 ; a2 = -8/25

 5

Exercice 3
Soit un filtre numérique de réponse impulsionnelle h1(k) =( ½)ku(k) où u(k) est l’échelon unité.
1) Que peut-on dire sur la stabilité de ce filtre ?
Ce filtre est stable car h(k) converge

2) Calculer la fonction de transfert en H1(z)

H1(z) = z/(z-1/2)

3) On considère le filtre numérique avec une fonction de transfert en H2(z) = 1/[4-z], que peut-on
dire de la stabilité de ce filtre ? Calculer sa réponse impulsionnelle.
Les deux filtres sont mis en cascade pour obtenir un filtre H équivalent.
a) Quelle est sa fonction de transfert H(z) ?
H2(z) = 1/(4-z)

H2(z) = -[(z z-1)]/[(z-4)]

 TZ-1(H2(z) ) = h2(k) = (-1/4)*(4)k

 TZ-1 h2(k) = (-4)k-1u(k)

H(z) = -z/(z-1/2)(z-4)

 H(z) = a/(z-1/2) + b/(z-4)

 6

Exercice 3
 H(z) = a/(z-1/2) + b/(z-4)

a = (-1/z-4)z=1/2 = -1/-3.5 = 2/7

b = (-1/z-1/2)z=4 = -1/3.5 = -2/7

D’où H(z) = (2/7). (z/z-1/2) - (2/7). (z/z-4)

b) Calculer sa réponse impulsionnelle.

h(k) = [(2/7) (1/2)k-(2/7) (4)k]u(k)

 7

Exercice 5
Soit le système décrit par :
y[n]=0.83x[n]+0.166y[n-1]]

1) Calculer sa fonction de transfert.

C’est un filtre RII

y[n]=0.83x[n]+0.166y[n-1]  y[n] - 0.166y[n-1]=0.83x[n]

H(z)=0.83/1-0.166z-1

D’où H(z)=0.83z/z-0.166

2) Retrouver le filtre analogique équivalent (T = 1s).

On a p = (1-z-1)/T  p = 1-z-1  z-1 = 1 - p

On a H(z)=0.83/1-0.166z-1  H(p) = 0.83/(1-0.166(1-p))

H(p) = 0.83/(1-0.166+0.166p)

H(p) = 0.83/(0.834+0.166p)
 8

Exercice 5
La réponse indicielle unitaire s(n) d’un système discret vaut 2(1/3)nu(n), où u(n) est l’échelon
unité.
Calculer la réponse impulsionnelle de ce système.

s(n) = 2(1/3)nu(n)S(z) = 2*(z/z-1) = 2z/(z-1/3)

E(z) = z/z-1

H(z) = S(z)/E(z) = [2z/(z-1/3)]*[z-1/z] = [2(z-1)]/(z-1/3)

TZ-1 H(z) = [2z]/[z-1/3]-[[z-1*[2(z)]/(z-1/3)]

h(n) = 2(1/3)n-2(1/3)n-1 = 2(1/3)n-2(1/3)n (1/3)-1

h(n) = 2(1/3)n -6(1/3)n

h(n) = -4(1/3)n
 9

Exercice 4
Il faut rendre ce filtre causal
On pose k = n+1  n = k-1
y(k)= y(k-1) – 0.5 y(k-2) +0.5 x(k-1)
TZ: Y(z)[1-z-1+0.5z-2] = 0.5 z-1X(z)
H(z) = Y(z)/X(z) = 0.5z-1/(1-z-1+0.5z-2)=0.5z/(z2-z+0.5)
x(n)= d(n)  X(z) =1
S(z) = e-aTsin(wT) z/(z2-2 e-aTcos(wT)z +e-2aT)
Avec e-aT sin(wT)=0.5 ; 2 e-aTcos(wT)=1; e-2aT=0.5 alors e-aTcos(wT)=0.5
tg(wT) =1 ; e-2aT=0.5
wT=pi/4 ; e-aT=0.25
Donc h(n)=(1/4)nsin(n*pi/4)
Réponse impulsionnelle d’entrée X(z)=z/z-1
Y(z)= [0.5z/(z2-z+0.5)]*[z-1/z]
Y(z)=0.5/(z-1)(z2-z+0.5)
Delta = j2
z’=(1-j)/2
z’’=(1+j)/2
Y(z)=0.5/(z-1)(z-z’)(z-z’’)on a 3 pôles simples