Université Ibn Zohr A.
U : 2023-2024
Faculté des Sciences Appliquées
Ait Melloul Responsable : M. Alaoui
Solution de TD5 : Traitement du signal
Filière : SMP
Parcours Electronique / S6
Questions de cours :
1./ Donner les conditions de stabilité pour un filtre RII ; que deviennent ces conditions pour
un filtre RIF ?
2./ Déterminer la causalité du filtre numérique dont l’équation aux différences suivante :
y [ n ] = x [ n + 1] + x [ n ] + x [ n − 1] ,
- Que peut-on faire pour corriger le problème ?
3./ Un filtre RII possède la fonction de transfert en Z :
1
H ( z) =
1 − 0,8 z −5
3-1./ Tracer la location des pôles et des zéros du filtre dans le plan complexe.
3-1./ Que peut-on dire de la stabilité du filtre ?
4./ Soit un système linéaire discret invariant dans le temps (SLIT), causal, caractérisé par la
fonction de transfert suivante :
z +1
H (z) =
z + 0, 25
Le signal d’entrée x [ n ] est donné par :
x [ n ] = ( 0, 4 ) U [ n ]
n
4-1./ Déterminer la réponse y [ n ] du SLIT causal.
4-2./ Le système est stable ou instable, justifier votre réponse ?
Correction : Question de cours
1./ Donner les conditions de stabilité pour un filtre RII ; que deviennent ces conditions pour
un filtre RIF ?
Réponse 1./ :
La transformée en Z (TZ) d’une séquence numérique x [ n ] , représentant un signal discret
(numérique) ou la réponse d’un système est définie comme suit :
n =+∞
X [ Z ] = TZ { x [ n ]} = ∑ x [ n] Z −n
n =−∞
Ou, Z est un variable complexe Z ∈ ^ et n ∈ ] .
- l’ensemble des valeurs de Z pour lesquelles la TZ existe est désigné par le terme « région de
la convergence » (RDC) de X [ Z ] et est défini par la condition de l’équation suivante :
1
PH : Laaboubi Mostafa SMP/ S6: Parcours Electronique / 2023-2024
n =+∞
∑ x [ n] Z
n =−∞
−n
<∞
- Particulièrement la transformée en Z de la réponse impulsionnelle d’un système linéaire
invariant dans le temps (SLIT). C’est la fonction de transfert de ce système. Nous
présenterons aussi une méthode permettant de déterminer si un SLIT est stable, à partir des
propriétés de sa fonction de transfert. Tout système linéaire invariant dans le temps (SLIT)
peut être exprimé sous la forme d’une équation aux différences :
M N
y [ n ] = ∑ bk x [ n − k ] + ∑ ak y [ n − k ] ,
k =0 k =1
Où, x [ n ] est le signal d’entrée du système et y [ n ] son signal de sortie.
La transformée en Z aux deux membres de l’équation aux différences, nous obtenons alors :
M N
Y [ Z ] = ∑ bk X [ Z ]Z − k + ∑ ak Y [ Z ]Z − k ,
k =0 k =1
Le rapport entre Y [ Z ] , qui est la transformée en Z du signal de sortie y [ n ] , et X [ Z ] , qui est
celle du signal d’entrée x [ n ] , est donné par :
M
Y [Z ] ∑b Z k
−k
H [Z ] = = k =0
,
X [Z ] N
1 − ∑ ak Z −k
k =1
- H [ Z ] est la fonction de transfert du SLIT. Il est possible d’obtenir la fonction de transfert
d’un SLIT à partir de l’équation aux différences qui le caractérise et vice-versa.
- H [ Z ] est aussi la transformée en Z de la réponse impulsionnelle h [ n ] du SLIT. On
remarquera que puisque le signal de sortie y [ n ] de tout SLIT peut être obtenu en effectuant le
produit de convolution de son signal d’entrée x [ n ] par sa réponse impulsionnelle h [ n ] , on a :
Y [Z ]
y [ n ] = x [ n ] ∗ h [ n ] YZZ
ZZZ Z Y [ Z ] = H [ Z ]. X [ Z ] U H [ Z ] =
X
TZ
,
X [Z ]
Les racines du numérateur et du dénominateur d’une fonction de transfert telle que H [ Z ] sont
respectivement désignées comme étant les zéros et les pôles de cette dernière.
La fonction de transfert rationnelle :
N [Z ]
H [Z ] = ,
D [Z ]
SLIT est stable si et seulement si sa région de convergence de H [ Z ] contient le cercle
unité Z = 1 ⇒ Les pôles de H [ Z ] sont à l’intérieur du cercle unité.
NB : ordre du numérateur <= celui du dénominateur de H [ Z ]
Proposition : Un filtre linéaire est stable si et seulement si sa fonction de transfert converge
sur le cercle unité.
{Z ∈ ^ / Z = 1} ∈ RDC ,
- Si le filtre est causal de réponse impulsionnelle infinie RII, la stabilité est assurée si et
seulement si le rayon de convergence (RDC) est inférieur à 1. Les pôles sont à l’intérieur du
cercle unité (Fig.1).
2
PH : Laaboubi Mostafa SMP/ S6: Parcours Electronique / 2023-2024
Fig. 1 : Pôles situés à l’intérieur de Z = 1 .
- Si le filtre est anti-causal de réponse impulsionnelle infinie RII, la stabilité est assurée si et
seulement si le rayon de convergence est supérieur à 1. Les pôles sont à l’extérieur du cercle
unité (Fig.2).
Fig. 2 : Pôles situés à l’extérieur de Z = 1 .
- Si le filtre est un filtre de réponse impulsionnelle finie RIF, ce filtre est toujours stable si
l’entrée est bornée.
2./ Déterminer la causalité du filtre numérique dont l’équation aux différences suivante :
y [ n ] = x [ n + 1] + x [ n ] + x [ n − 1] ,
- Que peut-on faire pour corriger le problème ?
Réponse 2./ :
Non causal à cause de l’indice k+1. Cependant la nature LIT du filtre, permet de décaler tout
les indices de x de 1
3./ Un filtre RII possède la fonction de transfert en Z :
1
H ( z) =
1 − 0,8 z −5
3-1./ Tracer la location des pôles et des zéros du filtre dans le plan complexe.
3-1./ Que peut-on dire de la stabilité de ce filtre ?
3
PH : Laaboubi Mostafa SMP/ S6: Parcours Electronique / 2023-2024
Réponse 3./ :
1 1
H (z) = = −5 5
z ( z − 0,8 )
−5
1 − 0,8 z
z5
H ( z) = 5
( z − 0,8)
⇒ :
La fonction de transfert H [ z ] possède 5 zéros à z = 0 , et 5 pôles.
z 5 − 0,8 = 0 ,
z 5 = 0,8 ,
i
( 2 k +1) π
z = 5 0,8e 5
, k = 0,1,...., 4
Im( z )
z2
z1
z3 zéros
Re ( z )
z5
z4
Fig. 3 : Pôles et les zéros de H [ z ] situés sur le cercle unité Z = 1 .
⇒ Le système ou le filtre est stable puisque les pôles sont à l’intérieur du cercle unité Z = 1 .
4./ Soit un système linéaire discret invariant dans le temps (SLIT), causal, caractérisé par la
fonction de transfert suivante :
z +1
H (z) =
z + 0, 25
Le signal d’entrée x [ n ] est donné par :
x [ n ] = ( 0, 4 ) U [ n ]
n
4-1./ Déterminer la réponse y [ n ] du SLIT causal.
4-2./ Le système est stable ou instable, justifier votre réponse ?
Réponse 4./ :
4-1./ Le système est linéaire discret et invariant dans le temps (SLIT), causal, caractérisé par
la fonction de transfert H ( z ) ⇒ :
Y ( z ) = X ( z ) .H ( z ) ,
Or :
4
PH : Laaboubi Mostafa SMP/ S6: Parcours Electronique / 2023-2024
x [ n ] = ( 0, 4 ) U [ n ]
n
+∞
⇒ : TZ { x [ n ]} = X [ Z ] = ∑ ( 0, 4 ) .Z − n ,
n
n=0
1
X [Z ] = .
1 − 0, 4 Z −1
Alors :
⎛ z +1 ⎞ ⎛ 1 ⎞
Y ( z ) = X ( z ) .H ( z ) = ⎜ ⎟ .⎜ −1 ⎟
,
⎝ z + 0, 25 ⎠ ⎝ 1 − 0, 4 Z ⎠
z ( z + 1)
Y (z) = ,
( z + 0, 25 )( z − 0, 4 )
Le degré du numérateur de H [ z ] étant du même ordre que celui de son dénominateur,
considérons :
Y ( z) ( z + 1)
= ,
z ( z + 0, 25)( z − 0, 4 )
⇒ :
Y (z) a b
= + ,
z ( z − 0, 4 ) ( z + 0, 25 )
- Calcul de a :
a = ( z − 0, 4 ) .
( z + 1) = 2,15 .
( z + 0, 25)( z − 0, 4 ) z =0,4
- Calcul de b :
b = ( z + 0, 25 ) .
( z + 1) = −1,15 .
( z + 0, 25 )( z − 0, 4 ) z =−0,25
Donc :
Y (z) 2,15 1,15
= − ,
z ( z − 0, 4 ) ( z + 0, 25 )
et :
2,15 1,15
Y (z) = − ,
(1 − 0, 4 z −1 ) (1 + 0, 25 z −1 )
En utilisant la transformée en Z inverse, on obtient :
⎧⎪ 2,15 ⎫⎪
y [ n ] = TZ −1 {Y ( z )} = TZ −1 ⎨
1,15
− −1 ⎬
,
⎩⎪ (1 − 0, 4 z ) (1 + 0, 25 z ) ⎭⎪
−1
Donc, la réponse du système en sortie y [ n ] est :
y [ n] = ⎡ 2,15. ( 0, 4 ) − 1,15. ( −0, 25) ⎤ U [ n] .
n n
⎣ ⎦
z +1
4-2./ La fonction de transfert H [ z ] est : H ( z ) = ,
z + 0, 25
5
PH : Laaboubi Mostafa SMP/ S6: Parcours Electronique / 2023-2024
possède un zéro à z0 = −1 , et un pôle à z1 = −0, 25 , comme z1 ≤ 1 donc la stabilité est assurée
car le rayon de convergence (RDC) est inférieur à 1. Les pôles sont à l’intérieur du cercle
unité.
Exercice 1 :
Soit le filtre numérique défini par l’équation suivante:
1 1 1
y (n) = x(n) + x(n − 1) + x(n − 2)
4 2 4
1./ Déterminer la fonction de transfert H(z) du système.
2./ Déduire la réponse impulsionnelle du filtre numérique.
3./ Déterminer la réponse impulsionnelle fréquentielle H(f).
4./ Donner le module de la fonction de transfert H ( f ) , tracer le comportement fréquentiel du
filtre et déduire son type.
Correction : Exercice 1
Soit le filtre numérique défini par l’équation suivante:
1 1 1
y (n) = x(n) + x(n − 1) + x(n − 2)
4 2 4
1./ Déterminer la fonction de transfert H(z) du système.
Réponse 1./ :
- Le filtre est défini par l’équation suivante:
1 1 1
y (n) = x(n) + x(n − 1) + x(n − 2)
4 2 4
Donc, la transformée en Z est :
⎧1 1 1 ⎫
TZ { y (n)} = TZ ⎨ x( n) + x(n − 1) + x(n − 2) ⎬
⎩4 2 4 ⎭
1 1 1
= TZ { x(n)} + TZ { x(n − 1)} + TZ { x(n − 2)}
4 2 4
Alors,
1 1 1
Y ( Z ) = X ( Z ) + X ( Z ) Z −1 + X ( Z ) Z −2
4 2 4
Or, la fonction de transfert est :
Y (Z )
H (Z ) =
X (Z )
Donc :
1 1 1
H ( Z ) = + Z −1 + Z −2
4 2 4
= (1 + 2 Z −1 + Z −2 )
1
4
2./ Déduire la réponse impulsionnelle du filtre numérique.
Réponse 2./ :
6
PH : Laaboubi Mostafa SMP/ S6: Parcours Electronique / 2023-2024
- La réponse impulsionnelle h ( n ) :
D’après la définition de la transformée en Z :
+∞
H ( Z ) = ∑ h ( n )Z − n avec n∈]
−∞
1 1 −1 1 −2
H (Z ) = + Z + Z
4 2 4
Donc :
⎧1
⎪4 si n = 0 et n = 2
⎪
⎪1
h ( n) = ⎨ si n = 1
⎪2
⎪0 sin on
⎪
⎩
3./ Déterminer la réponse impulsionnelle fréquentielle H(f).
Réponse 3./ :
- Réponse impulsionnelle fréquentielle H ( f ) :
H ( f ) = H ( Z ) Z =ei 2 πfTe
1 1 − i 2 πfTe 1 i 4 πfTe
= + e + e
4 2 4
1
(
= 1 + 2e − i 2 πfTe + ei 4 πfTe
4
)
1
Avec Te = est la période d’échantillonnage.
fe
1
(
H ( f ) = 1 + 2e − i 2 πfTe + e − i 4 πfTe
4
)
1 ⎛ e −i 2 πfTe ⎞
= ⎜ −i 2 πfTe + 2e− i 2 πfTe + e− i 2 πf 2Te ⎟
4⎝e ⎠
e− i 2 πfTe ⎛ 1 − i 2 πfTe ⎞
= ⎜ − i 2 πfTe + 2 + e ⎟
4 ⎝e ⎠
− i 2 πfTe
=
e
4
(
2 + ei 2 πfTe + e− i 2 πfTe )
− i 2 πfTe
=
e
4
( 2 + 2 cos ( 2πfTe ) )
Finalement :
e − i 2 πfTe
H(f )=
2
(1 + cos ( 2πfTe ) )
7
PH : Laaboubi Mostafa SMP/ S6: Parcours Electronique / 2023-2024
4./ Donner le module de la fonction de transfert H ( f ) , tracer le comportement fréquentiel du
filtre [ − f e , f e ] et déduire son type.
Réponse 4./ :
- Module de H ( f ) :
e − i 2 πfTe
H(f) =
2
(1 + cos ( 2πfTe ) )
e− i 2 πfTe
=
2
(1 + cos ( 2πfTe ) )
1
2
(1 + cos ( 2πfTe ) )
=
- Comportement fréquentiel du filtre (Figure 7):
H(f)
fe fe f
− fe − fe
2 2
Fig. 4 : Comportement fréquentiel du filtre.
H ( f ) : est appelée le filtre cosinusöde surélevée elle conserve la composante à la fréquence
fe
f = 0 et élimine celle à f = (Figure 7), donc c’est un filtre basse bas de fréquence de
2
fe
coupure FC = .
2
Exercice 2 :
Soit un filtre numérique causal RII d’ordre 1, caractérisé par l’équation aux différences
suivante :
y [ n ] = x [ n ] + a. y [ n − 1] .
1./ Etudier la stabilité du système ?
Réponse 1./ :
8
PH : Laaboubi Mostafa SMP/ S6: Parcours Electronique / 2023-2024
L’analyse d’un système linéaire invariant dans le temps SLIT (filtre numérique) consiste à
extraire les principales caractéristiques de ce filtre, qui sont :
- La fonction de transfert H [ z ] .
- La réponse fréquentielle H ( f ) et la nature du filtre.
- La réponse impulsionnelle temporelle h [ n ] .
- La stabilité.
Considérons le filtre défini par l’équation aux déférences suivante :
y [ n ] = x [ n ] + a. y [ n − 1] .
- La fonction de transfert H [ z ] :
La fonction de transfert H [ z ] est obtenue en appliquant la transformée en Z (TZ) aux
deux membres de l’équation :
Y [ z ] = X [ z ] + a.z −1Y [ z ] ,
Y [ z] 1
=
X [ z] 1 − a.z −1
Soit :
Y [ z] 1
H [ z] = = .
X [ z] 1 − a.z −1
- Stabilité du filtre :
On peut étudier la stabilité du filtre à partir de la position des pôles de H [ z ] :
z
H [ z ] ==
z−a
Pôles de H [ z ] ,
z −a = 0,
Soit :
z0 = a .
On constate pour :
- z0 = a ≤ 1 , la valeur de la réponse impulsionnelle décroit vers zéro avec le temps
⇒ donc le système est stable dans ce cas.
- z0 = a > 1 , la valeur de la réponse impulsionnelle croit vers l’infini avec le temps
⇒ donc le système est instable dans ce cas.
2./ Déterminer la réponse impulsionnelle h [ n ] du système.
Réponse 2./ :
- La réponse impulsionnelle temporelle h [ n ] :
La réponse impulsionnelle temporelle peut être obtenue par de deux façons soit par :
- Le calcul direct, qui consiste à utiliser l’équation aux différences et à
remplacer l’entrée du système par l’impulsion de Dirac :
⎧1 si n = 0,
x [ n] = δ [ n] = ⎨
⎩0 sin on
Alors,
x [ n] = δ [ n] ⇒ y [ n] = h [ n] .
9
PH : Laaboubi Mostafa SMP/ S6: Parcours Electronique / 2023-2024
- La transformée en Z inverse ( TZ −1 ).
h [ n ] = TZ − { H [ z ]} =
1
∫ H [ z ] z n −1dz ,
2iπ cercle unité
h [ n] = ∑ résidus ( H [ z ] .z n −1 ) .
Pôles: H [ z ]
Détermination de La réponse impulsionnelle temporelle h [ n ] : par la méthode directe
- Le SLIT est causal : ⇒ ∀ n < 0 : y [ n] = 0 .
On a :
y [ n ] = x [ n ] + a. y [ n − 1] ,
Alors :
- Pour n = 0 ,
y [ 0] = x [ 0] + a. y [ −1] ,
x [ 0] = δ [ 0] ⇒ y [ 0] = h [ 0] ,
y [ −1] = 0
Donc :
h [ 0] = δ [ 0] + a. y [ −1] = δ [ 0] ,
h [ 0] = 1 .
- Pour n = 1 ,
y [1] = x [1] + a. y [ 0] ,
y [ 0] = h [ 0] ⇒ y [ 0] = 1 ,
x [1] = δ [1] = 0 ⇒ h [1] = a. y [ 0] .
Donc :
h [1] = a .
- Pour n = 2 ,
y [ 2] = x [ 2] + a. y [1] ,
y [1] = h [1] ⇒ y [1] = a ,
x [ 2] = δ [ 2] = 0 ⇒ h [ 2] = a. y [1] .
Donc :
h [ 2] = a 2 .
- Pour n = 3 ,
y [3] = x [3] + a. y [ 2] ,
y [ 2] = h [ 2] ⇒ y [ 2] = a 2 ,
x [3] = δ [3] = 0 ⇒ h [3] = a. y [ 2] .
Donc :
h [3] = a 3 .
Par récurrence, on trouve :
∀n≥0 h [ n] = a n ,
D’où :
10
PH : Laaboubi Mostafa SMP/ S6: Parcours Electronique / 2023-2024
∀ n∈] h [ n ] = a n .u [ n ] .
Où :
⎧a n si n ≥ 0
h [ n] = ⎨
⎩0 sin on
3./ Donner la fonction de transfert en fréquence.
Réponse 3./ :
- On peut calculer la réponse fréquentielle H ( f ) par la transformée de Fourier discrète (TFD)
de h [ n ] . Cependant, on l’obtient en posant z = e+ i 2 πfTe :
H ( f ) = H [ z] z = ei 2 πfTe
,
1
H(f)=
1 − a.e−2iπfTe
Ou encore,
i .arg ( H ( f ) )
H(f)= H(f )e .
Avec :
⎧ 1
⎪ H ( f ) = 1 + a 2 − 2a cos ( 2πfT )
2
⎪ e
⎨ .
⎪arg H ( f ) ) = Arctg ⎡ a sin ( 2 πfT ) ⎤
⎪ (
e
⎢ ⎥
⎩ ⎣ 1 − a cos ( 2 πfT e ) ⎦
Pour tracer H ( f ) , on s’intéresse à la bande de fréquence [ 0, Fe ] , car le signal est
2
1
échantillonné à Fe = . Donc :
Te
H ( f + nFe ) = H ( f )
2 2
pour n∈] .
- on définit le Gain statique G s par :
G s = H ( f = 0 ) = H ( Z = 1) ,
1
G s= .
1− a
Fig. 5 : Spectre en module de H ( f ) entre 0 et Fe .
11
PH : Laaboubi Mostafa SMP/ S6: Parcours Electronique / 2023-2024
D’après la représentation (spectre en module de H ( f ) ), on constate que ce filtre est un filtre
Fe
passe-bas de fréquence de coupure Fc = du premier ordre (ordre 1).
2
4./ Soit un système SLIT causal, caractérisé par la fonction de transfert suivante :
H ( z) =
( z + 4 )( z + 0, 25 )
( z + 0,5 ) ( z 2 + 0, 2 z + 0,5 )
- Le système est stable ou instable, justifier votre réponse ?
Réponse 4./ :
- Calcul les pôles de H [ z ]
Le dénominateur de la fonction de transfert H [ z ] est :
( z + 0,5) ( z 2 + 0, 2 z + 0,5) = 0 ,
Soit :
⎧ z + 0,5 = 0
⎪
⎨Où
⎪ z 2 + 0, 2 z + 0,5 = 0
⎩
Alors :
⎧ z1 = −0,5
⎪où
⎪⎪
⎨ z2 = −0,1 − 0, 7i
⎪où
⎪
⎪⎩ z3 = −0,1 + 0, 7i
Donc les pôles de H [ z ] sont : z1 , z2 et z3 .
Les modules respectifs sont :
z1 = 0,5 ,
et,
z2 = z3 = ( 0,1) + ( 0, 7 ) = 0,5 = 0, 707 .
2 2
Tous les pôles de H [ z ] se trouvent à l’intérieurs du cercle unité et le SLIT est causal
⇒ le systèmes est stable.
12
PH : Laaboubi Mostafa SMP/ S6: Parcours Electronique / 2023-2024