10/11/2022
= , =
-fe/2 fe/2 f
-π π ω
Boularbah SOUISSI 1
Calcul de la transformée en z (TZ)
Exp1. Représenter la séquence suivante par sa transformée en z :
Exp2.
une série géométrique infinie: RDC
α=0.5
la région de convergence (ROC) est:
Boularbah SOUISSI 2
1
10/11/2022
3. Condition d'existence et la région de convergence (ROC)
Application:
Signal droitier: x1(n)=αnU(n), H1(z)=1/(1− z-1) et ROC: <
Signal Gaucher: x2(n)=βnU(-n-1), H2(z)=1/(1− β z-1) et ROC: < β
Signal bilateral: x3(n)=αnU(n)+nU(-n-1), H3(z)=1/(1− z-1)+ 1/(1− z-1) et ROC: < <
x5(n)=αnU(n), x6(n)=-βnU(-n-1), x7(n)=αnU(n)-nU(-n-1),
ROC d'un signal droitier ROC d'un signal gaucher ROC d'un signal bilatéral
ROC: < ROC: < β ROC: < <
H2(z)=1/(1- βz-1) H3(z)=1/(1- z-1)+ 1/(1- z-1)
H1(z)=1/(1- z-1)
Boularbah SOUISSI 3
Boularbah SOUISSI 4
2
10/11/2022
Fonction de transfert
Système discret linéaire
invariant dans le temps (DLI)
N M N M
y (n) ak y (n k ) bk x(n k ) Y ( z ) ak z k Y ( z ) bk z k X ( z )
k 1 k 0 k 1 k 0
M
z zk
z z1 z z 2 .........z z M b z M N .
M
k
Y ( z ) b0 b1z 1 ......... bM z M b z k
M N k 1
H ( z) k 0
H ( z ) b0 z . 0
X ( z ) 1 a1z 1 .......... ..aN z N N
1 ak z k z p1 z p2 .......... z pN N
z p
k 1 k 1
k
où zi sont les zéros, et pi sont les pôles.
1
0.8
Exemple: y(n)=x(n)+0.5y(n-1) 0.6
0.4
Imaginary Part
0.2
H(z)=z/(z-0,5) h(n)=(0.5)nu(n) 0
o x
-0.2
-0.4
0,5
-0.6
-0.8
-1
-1 -0.5 0 0.5 1
Real Part
Boularbah SOUISSI
Zeros: z1=0 et Poles: p1=0.5
5
Structure des filtres numériques récursifs
filtres à Réponse Impulsionnelle Finie (RIF) et Infinie (RII)
y(n)=αy(n-1)+x(n) y(n)=αx(n-1)+x(n)
Y z z 1Y z X z Y z .z 1 X z X z
H z 1/(1 z ) 1
H z Y ( z ) / X ( z ) .z 1 1
h[n] u[ n]
n h[ n] [ n 1] [ n]
Eléments de base z-1
+
Additionneur
retardateur multiplicateur.
Structure des filtres
RIF RII y(n)
x(n)
+ y(n) + y(n)=αy(n-1)+x(n)
x(n)
Z-1 Z-1
α
y(n)=αx(n-1)+x(n)
Filtres non récursifs de 1er ordre Filtres récursifs de 1er ordre
filtres RIF (Réponse Impulsionnelle finie). filtres RII (Réponse Impulsionnelle Infinie).
Boularbah SOUISSI 6
3
10/11/2022
Structure des filtres numériques récursifs
filtres à Réponse Impulsionnelle Finie (RIF) et Infinie (RII)
1. Filtres récursifs (RII): M
k
N M
Y ( z) b z
k 0
k
b0 b1 z 1 .........bM z M
y (n) ak y ( n k ) bk x( n k ) H ( z) N
X ( z) 1 a1 z 1 ............aN z N
k 1 k 0 1 ak z k
k 1
M
z z k
H ( z ) b0 z M N .
z z1 z z 2 ......... z z M b z M N .
k 1
Ce filtre récursif est associé z p1 z p2 .......... z pN 0 N
z p
k
exclusivement aux filtres RII k 1
(Réponse Impulsionnelle Infinie).
2. Filtres non récursifs (RIF):
M Y ( z) M
y (n ) bk x ( n k ) H ( z) bk z k
k 0 X ( z ) k 0
Ce filtre non récursif est associé exclusivement aux filtres
RIF (Réponse Impulsionnelle finie).
Boularbah SOUISSI 7
Stabilité et causalité d’un système et ROC
Diagramme de pôles et zéros
a) système instable et causal b) causal et stable c) stable et anti-causal
(Signal droitier) (Signal droitier) (Signal gaucher)
Stabilité :
•le cercle unité (z=1) est dans le domaine de convergence (ROC)
•Tous les pôles sont à l'intérieur de cercle unité
Causalité: si un système est causal, (h[n]=0, n<0), alors le ROC est à l'extérieur d'un cercle qui est l'extérieur
du pole le plus élevé.
Boularbah SOUISSI 8
4
10/11/2022
Transformée inverse
La définition formelle de la transformée en z inverse est :
Il existe des méthodes plus faciles pour faire la transformée inverse.
1. Forme polynomiale:
2.Longue division:
Boularbah SOUISSI 9
Transformée inverse
3. Expansion en fractions partielles
En multipliant par z, et en faisant la transformée inverse, on obtient :
4. Méthode des résidus: , =
= ,
1. Pôles simples:
= − ] : =
− 0.5
= 0.5 ( )
2. Pôles multiples:
1
= , = − ]
( − 1)!
: = , = ( )
Boularbah SOUISSI 10
5
10/11/2022
Réponse en fréquence
On peut obtenir la réponse en fréquence d'un système discret si on remplace z = ej2πF dans la fonction de
transfert H(z). C'est la DTFT de la réponse impulsionnelle h[n].
Exemple
Tracer la réponse en fréquence de
La première chose à faire est de remplacer z = ej2πF:
Et on trouve l'amplitude
0.5 sin 2
et la phase =− ( )
1 − cos 2
0.8
0.6
0.4
Imaginary Part
0.2
-0.2
-0.4
-0.6
-0.8
-1
-1 -0.5 0 0.5 1
Real P art
Boularbah SOUISSI 11
Description du système par ses pôles et ses zéros
Interprétation graphique d’une fonction de transfert dans le plan complexe z
La réponse d’un zéro ou d’un pôle (premier ordre)
Exemples (Filtre RIF) : = + −1
( ) +
= = +1 = ⇒ é = = 1+
( )
+ y(n)
x(n)
Z-1
α
M
=− . :
w=0 ==> = .
0.5
et w=π ==> = /
réponse fréquentielle
FIR Pass Haut
Boularbah SOUISSI
6
10/11/2022
La réponse d’un zéro ou d’un pôle (premier ordre)
Exemples (Filtre RIF) : = + −1
( ) +
= = +1 = ⇒ é = = 1+
( )
+ y(n)
x(n)
Z-1
α
= . :
w=0 ===> = /
et w= π ===> = .
% Filtre RIF Pass Haut: H(z)=1+0.5z^-1
b=[1 0.5];
a=[1 0];
% figure(1)
% freqz(b,a)
tf(b,a,-1)
[H,w]=freqz(b,a);
FIR Pass Bas
figure(3),plot(w,abs(H));grid Réponse fréquentielle
xlabel('w');ylabel('Abs(H)'),
title('Réponse fréquentielle: Filtre RIF
Passe Bas')
figure(4),zplane(b,a)
title('H(z)=1+0.5Z^-1')
Boularbah SOUISSI 13
Synthèse des filtres à réponse impulsionnelle finie (RIF)
1. Propriétés des
RIF à Phase linéaire
Filtre FIR
H(z)
H z=e = ( ) ( )
=− H z = h[n]z = b z
Un système à réponse de phase linéaire ne connaît pas
d'effet de distorsion de phase.
Symétrie
h0=h4, h1=h3
H z = h[n]z
= h0 + h1z + h2z + h3z + h4z
= h0 1 + z + h1(z + z )+h2 z
= h0z z + z + h1z (z + z )+h2z
=z {h0 z + z + h1 (z + z )+h2} Delete 2
H z=e =e h0(e +e + h1 (e
+ e )+h2 )
( )
symétrie impaire(5) et paire (4): (h[n]=h[N-n], n=0:N) =e (2h0cos 2ω + 2h1 2ω + h2) = ( )
et antisymétrie impaire (5) et paire (4): (h[n]=-h[N-n], n=0:N)
k est le centre de symétrie Boularbah SOUISSI = ⇨Phase linéaire 14
7
10/11/2022
Synthèse des filtres à réponse impulsionnelle finie (RIF)
2. Conception par la méthode des fenêtres
Faire la synthèse d’un filtre c’est rechercher la fonction de transfert correspondant à la spécification
d’un gabarit imposé. Plusieurs méthodes sont envisageables H(z)
Coefficients de filtre RIF: H z =∑ z ?
H z = [ ]z = z ⟹ = [ ]
Calcul des coefficients de h[n] par la méthode des fenêtres
le calcul des coefficients revient à trouver la réponse impulsionnelle du filtre. On part de l’amplitude
de la réponse en fréquence d’un filtre idéal et on prend sa transformée de Fourier inverse.
1 ( )
h n = H ω e dω Filtre ideal (rectangulaire): H ω = 1 h n =
2π
Boularbah SOUISSI 15
Réponse impulsionnelle idéale [n] des filtres.
Boularbah SOUISSI 16
8
10/11/2022
Effet de fenêtrage/Utilisation des fenêtres non rectangulaire :
Si on limite hd[n], cela revient à le multiplier par une fenêtre rectangulaire de même largeur, modifiant ainsi hd[n]. –
L’effet est la convolution de H(ω) par la transformée de Fourier de la fenêtre rectangulaire, ce qui affecte aussi la
réponse en fréquence du filtre. On peut minimiser l’effet en utilisant des fenêtres non rectangulaire (Tableau).
Note: l'ordre du filtre, N, dépendra du type de fenêtre et de la largeur de la bande de transition.
Fonction de Window W[n]
Boularbah SOUISSI 17
Synthèse des filtres RII
Ces filtres sont généralement obtenus par transposition dans le domaine numérique, d’une fonction de transfert
analogique. Cette correspondance confère aux filtres RII
• une facilité de conception et des possibilités de substitution directe avec leurs homologues analogiques.
• En terme de charge de calcul, ils sont souvent plus économiques que les filtres RIF car la récursivité réduit le
nombre de coefficients. En contrepartie les filtres RII présente une sensibilité importante vis-à-vis de la
quantification des coefficients ; des risques d’instabilité sont à prévoir pour les filtres à coefficient de qualité
élevé.
Chapitre 2:Filtres analogiques Filtres Numeriques RII
fonction de transfert analogique fonction de transfert numérique
( )
Méthodologie de conception
1. Définir le gabarit du filtre numérique
2. Convertir ce gabarit en un gabarit correspondant au filtre analogique
3. Faire la synthèse du filtre analogique (Butterworth, …) Ha(s)
4. Transformer le filtre analogique Ha(s) en un filtre numérique équivalent Hd(z) par deux métodes.
Méthodes
• Transformation bilinéaire
• Conservation de la réponse impulsionnelle du filtre analogique ("numérisation") (IIR by impulse invariance)
Boularbah SOUISSI 18
9
10/11/2022
Synthèse des filtres RII : Transformation bilinéaire
Chapitre 2:Filtres analogiques Filtres Numeriques RII
∗( − )
fonction de transfert analogique = fonction de transfert numérique
( + )
( )
= ( ) ∗( )
( )
Correspondances fréquentielles
2 2
= tan ou = tan ( ) Indice indique
2 2
d:digital
=2 =2 a:analogique
Exercisse
Concevoir un filtre pass-bas (FPB) avec une fréquence de coupure =0.2π et Te=1 en utilisant la
transformation bilinéaire appliquée à:
=
+
0.2π 0,35
= tan = 2 ∗ tan =0,65 Donc, =
+ 0,35
0,245(1 + )
= ( ) ) =
∗( SOUISSI
Boularbah
( )
1 − 0,509 19
10
