Chapitre 1 : Rappels sur les filtres numériques (RIF et RII)

1.1 Transformée en Z
1- Définitions :
En traitement du signal il arrive souvent de ne devoir considérer que certaines valeurs d’un
signal temporel, par exemple à nT, n = 0; 1; 2 …. , où T est un nombre positif fixé, appelé
période d’échantillonnage.
Autrement dit, en partant d’un signal continu f (t), on considère un signal discret représenté
par (f(nT)) avec n ϵN. La transformée en Z permet d’associer une fonction d’une variable
complexe à ce type de signal discret, et d’employer certaines techniques de l’analyse
complexe pour étudier des propriétés du signal. La transformée en Z est également employée
en filtrage numérique ainsi que pour la résolution d’équations récurrentes.
Alors la transformée en z d’une séquence x(n) est définie comme la série X(z) calculée comme
suit : X(z)= où z ϵ₵ est une variable complexe et n ₵ N. On appelle encore
cette équation la transformée directe, car c’est la relation qui permet d’obtenir X(z) à partir de
x(n).
Cette transformée est aussi qualifiée de bilatérale par opposition à unilatérale. La transformée
en z unilatérale pour une séquence causale est définie par Xu(z) calculée comme suit :
Xu(z)=
L’opération inverse porte le nom de transformation inverse. Ainsi l'analyse d'un système
discret se fera en général au moyen de la transformée en Z.
On note que z étant une variable complexe et la fonction X(z) possède un domaine de
convergence (en z) ou région de convergence (RdC) qui est en général un anneau centré sur
l'origine et compris entre deux rayons R1 et R2, où R1 < |z | < R2
La transformée en z doit toujours indiquer sa région de convergence, puisque c'est une série
de puissance infinie.
– RdC ne contient pas de pôle de X(z). Elle correspond, en général, pour les signaux causaux à
l’extérieur d’un cercle |z|> a
– Pour les signaux physiques (qui ont une durée d’existence finie) :
RdC= ensemble du plan complexe avec l’exclusion possible de z=0 ou z=

Exemple : Déterminer les transformée en z des séquences finies suivantes :

X1(n) = { 1 , 2 , 5 , 7 , 0 , 1 }, x2(n) = { 1 , 2 , 5 , 7 , 0 , 1 } où pour x2 la valeur 5 signifie
x2(0) = 5
Alors : X1(z) = 1 + 2z-1 + 5z-2+ 7z-3+ z-5 R.O.C : tout le plan Z sauf z = 0
X2(z) = z2+ 2z+ 5 + 7z-1+ z-3 R.O.C : tout le plan Z sauf z = 0 et z = ∞

2- Propriétés de la transformée en Z :
a- Linéarité :
Soient α et β deux nombres complexes et x et y deux suites de nombres complexes.
Alors : Z( αx + βy) = αZ(x) + βZ(y)
b- Signal retardé
Soient x une suite numérique, k € N et y la suite de terme général
Soit y(t - n.Te) nulle pour t < n.Te, Zy(t - n.Te z-n .Zy(t)
On retiendra qu’un retard de Te se traduit par une multiplication par z-1.

c- Théorèmes de la valeur initiale et de la valeur finale


Y(0)= lim y(n ) = lim Y(z) et y(∞)= lim y(n) = lim 1-z-1 .Y(z) 
n0 z n∞ z
 d- Produit de convolution
La transformée en z du produit de convolution de deux signaux discrets g(n)
et u(n) est égale au produit des transformées en z de chacun des signaux.
Z {g (n) * u (n)} = Z ( = G(z) · U (z)

Ex : Calculer la transformée en z du produit de convolution des séquences suivantes :

x(n) = 2anu(n) et y(n) = (n-1)
X(z) = 2 / (1-az-1) Y(z)= z -1
X(z) . Y(z) = 2z-1 /(1-az-1)

2-Exemples de transformées :
2.1. Suite de Dirac :
X(z)
Comme on a, x(n) = 1 à n = 0 , on en déduit que :
X(z)= 1.z-0 alors X(z)=1

2.2. Suite échelon unité

La transformée en z de la suite échelon unité déterminée par :
u(n) = 1 pour tout n > 0.
On a :

Or on sait que :

On en déduit donc que :

Table de transformées en z :
1.2 Les filtres numériques
1.2.1 Classification des filtres numériques
Un filtre numérique peut être classé selon sa réponse :
Filtres numériques réponse impulsionnelle finie : la durée de sa réponse impulsionnelle est
finie : les filtres RIF ont leur réponse impulsionnelle à support fini
i.e. h(n) = 0 pour n<0 et n>N
Filtres numériques réponse impulsionnelle Infinie : les filtres RII ont leur réponse
impulsionnelle à support infini
i.e. h(n) ≠ 0 pour toute n
– le type de représentation temporelle
récursifs : la sortie y(n) dépend de l’entrée courante, des entrées précédentes et
des sorties précédentes
non récursifs : la sortie y(n) ne dépend que de l’entrée courante et des entrées
précédentes.
Réponse impulsionnelle
On sait, en analogique, que la réponse impulsionnelle d’un filtre le définit complètement. Il en
est de même pour les filtres numériques.
Notons toutefois une différence importante : l’impulsion analogique δ(t) n’est pas réalisable
physiquement, alors que la suite numérique impulsionnelle δ(n) est facilement réalisable :
δ(n) = 1 si n =0
δ(n) = 0 si n ≠0
On désignera par h(n) la réponse impulsionnelle du filtre.

1.2.2 Représentation d'un filtre numérique

Un filtre numérique est un algorithme de calcul qui fait correspondre à une suite d'échantillons
différences Y(n) s'écrit :

N M
Y(n)= ∑ bk .x(n-k) – ∑ ak .y(n-k)
k=0 k=1
Appelée équation aux différences linéaire à coefficients constants d’ordre N.
En utilisant la transformée en Z :
N M
Y(z)= ∑ bk. z-k.X(z) – ∑ ak. z-k .Y(z)
k=0 k=1

N M
Alors H(z) = Y(z) ∕ X(z) = ∑ bk. z ∕ (1 + ∑ ak. z-k )
-k
k=0 k=1

Si ak = 0, le filtre s'appelle non récursif ou RIF (Réponse Impulsionnelle Finie) ;

Si ak ≠ 0, le filtre s'appelle récursif ou RII (Réponse Impulsionnelle Infnie) ;
Si les pôles de H(z) sont à l'intérieur du cercle unité, le filtre est stable ;
Si la réponse impulsionnelle est causale h(n) = 0; pour tout n < 0, le filtre est causal.

Représentation par pôles et zéros

La fonction de transfert :

Les pôles de H(z) sont les valeurs de z pour lesquelles H(z) tend vers l'infini.
Les zéros de H(z) sont les valeurs de z pour lesquelles H(z) est nulle.
Les zéros et les pôles complexes de H(z) sont du type a ± jb.

1.2.3 Filtre à réponse impulsionnelle finie RIF

Les filtres RIF constituent une classe importante des filtres numériques car, ils sont
caractérisés par une phase linéaire dans la bande passante. Ils sont définis par l’équation
N
Y(n)= ∑ bk .x(n-k)
k=0
L’ordre du filtre numérique est donné par le degré de la récursivité de l’équation aux
différences finies associée : N.
La fonction de transfert en z est donnée par :
N
H(z) = ∑ bk. z-k
K=0
1.2.4 Propriétés des filtres RIF

-Les RIF sont toujours stables (système tout zéros).

- Une plus grande facilité d’implantation dans un système numérique de traitement.
- Une phase qui peut être exactement linéaire, par conséquent un temps de propagation de
groupe constant et une absence de distorsion harmonique dans le signal.
- A sélectivité équivalente, ils sont toujours plus coûteux (en temps de calcul) que leur
équivalent RII.

-Stabilité du filtre

Pour qu’un filtre numérique soit stable, il faut que ses pôles soient à l’intérieur du cercle
unité dans le plan complexe. Autrement dit, il faut que le cercle unité appartienne à la région
de convergence.
.
Où: z1, z2,…, zn sont les zéros, p1, p2,…, pn sont les pôles ou les valeurs qui annulent le
dénominateur.

Rque : Les filtres à réponse impulsionnelle finie sont toujours stables car ils n'admettent pas
de pôles.

1.2.5 Réalisation des filtres RIF et RII :

Les filtres numériques peuvent être réalisés à l'aide de trois éléments ou opérations de base,
soit l'élément gain, l'élément de sommation et le retard unitaire ; ces éléments sont suffisants
pour réaliser tous les filtres numériques linéaires possibles. La réalisation présentée dans la
figure ci-dessous est une réalisation directe de type filtre RII.

Figure 1.1. Réalisation directe de type filtre à réponse impulsionnelle infinie RII.

1.2.6 Analyse des filtres numériques :

Un filtre numérique linéaire peut être décrit par :
- une équation aux différences
- un produit de convolution
- sa fonction de transfert
- sa stabilité
- sa réponse impulsionnelle ou indicielle
- sa réponse fréquentielle et la nature du filtre
- son diagramme des pôles et des zéros
L’outil mathématique exploité pour faciliter son analyse est la transformée en Z.

Exemple :
Ex1. Soit le filtre numérique d´écrit par:
y(k) = 0.6x(k) + 0.1x(k − 1) + 0.6x(k − 2) + 1.2y(k − 1) − 0.4y(k − 2)

(a) Déterminer la fonction de transfert H(z).

(b) Préciser le type de filtre : RIF ou IIR. Justifier.

Ex 2. Soit le filtre numérique d´écrit par:

H(z) = 1 + 5z−1 − 3z−2 + 2.5z−3 + 4z−5

(a) D´déterminer l’´equation aux différences du filtre.

(b) Déterminer la réponse impulsionnelle.
(c) Préciser le type de filtre : RIF ou IIR. Justifier.

Solutions Exos filtres numériques :

Ex1 : y(k) = 0.6x(k) + 0.1x(k − 1) + 0.6x(k − 2) + 1.2y(k − 1) − 0.4y(k − 2)

Transformée en Z de l’équation :

Y(z) = 0.6X(z) + 0.1 z-1 X(z) + 0.6 z-2 X(z) + 1.2 z-1 Y(z) – 0.4z-2Y(z)

Y(z) [1 - 1.2 z-1 + 0.4z-2] = X(z) [0.6 + 0.1 z-1+0.6z-2]

La fonction de transfert est :

H(z) = Y(z)/X(z) = (0.6 + 0.1 z-1+0.6z-2) / (1 - 1.2 z-1 + 0.4z-2)

C’est un filtre à réponse impulsionnelle infinie : RII car H(z) possède des pôles.

Ex2: H(z) = 1 + 5z−1 − 3z−2 + 2.5z−3 + 4z−5

a- L’équation aux différences du filtre :

H(z)=Y(z)/X(z) = 1 + 5z−1 − 3z−2 + 2.5z−3 + 4z−5

Y(z)= X(z).[ 1+ 5z-1 – 3z-2 + 2.5 z-3 + 4z-5]

Y(k)= x(k) + 5x(k-1) -3x(k-2) +2.5x(k-3) +4x(k-5)

b- La réponse impulsionnelle :

h(k) = (k) + 5 (k-1) - 3(k-2) +2.5 (k-3) + 4(k-5)

avec (k)= 1 pour k=0 et (k)=0 pour k ≠0

c- Type de filtre : filtre RIF car il ne présente pas de boucle de retour dans sa structure
ou pas de pôles.
2- Synthèse des filtres numériques :

2.1 -Méthodes générales de synthèse des filtres numériques.

Les méthodes de synthèse les plus utilisées s’appuient sur les propriétés (bien connues !...) des
filtres analogiques.
On part d’un filtre analogique « prototype » de fonction de transfert Ganl (p) ayant
les propriétés souhaitées, puis on cherche le filtre numérique de fonction de transfert HNu (z)
ayant sensiblement les mêmes propriétés :

- même réponse impulsionnelle

- même réponse indicielle

- même réponse harmonique sur l’intervalle (fe/2 ; fe/2) … th de Shannon

2.1.1 Synthèse de filtres R.I.F.

Le calcul des coefficients des filtres FIR repose sur l'utilisation de la transformée de
Fourier numérique et des fenêtres spectrales (rectangulaire, Hamming, Hanning, ....).
Une suite d'échantillon h(n) séparés de Te=1/Fe admet pour transformée de Fourier
numérique H(jw) :

Un filtre RIF est définie par une fonction de transfert :

P
H(z) = ∑ a k .z − k
k=0

soit en régime harmonique :

P P
H( j.ω) = ∑ a k . ( e j.ω.TE )
−k
= ∑ a k .e − k. j.ω.TE
k=0 k=0

Il arrive souvent que ces coefficients aient des propriétés de symétrie :

a k = a P − k ou a k = − a P − k
Si on regroupe deux à deux les termes :
a k .e − k. j.ω.TE + a P − k .e − (P − k ). j.ω.TE = a k .  e− k. j.ω.TE + e− (P − k ). j.ω.TE 
P
− j. .ω.TE   p − k . j.ω.TE − . − k  j.ω.TE 
p 

= a k .e 2
.  e 2  + e 2  
 
 p   − j. P2 .ω.TE
= 2.a k .cos  − k  .ω.TE  .e
 2
 



∈\

 p   P
− j. .ω.TE
(si a k = − a P − k , on obtient 2.j.a k .sin  − k  .ω.TE  .e 2
)
 2  

On aura alors :
P
P
− j. .ω.TE 2
 p  
H( j.ω) = e 2
. ∑ 2.a k .cos  − k  .ω.TE 
k=0
 2
  

S
D’où :
P
Arg [ H( j.ω) ] = − .TE .ω ( + π ) (le terme + π intervient si S < 0)
2
 -2- Méthodes de synthèse

On cherche à réaliser un filtre de fonction de transfert :

P
H Num (z) = ∑ a k .z − k
k=0

ou, ce qui est équivalent, défini par la récurrence :

P
y(n) = ∑ a k .x(n − k)
k=0

Or, on a vu que pour un filtre RIF : a k = h(k ) , où h(n) est la réponse impulsionnelle
du filtre. Le problème se ramène donc à la détermination de la réponse impulsionnelle h(n) du
filtre numérique.
Deux cas peuvent se présenter :
• on sait calculer la réponse impulsionnelle g(t) du filtre analogique
prototype (technique de l’échantillonnage temporel)
• on ne sait pas calculer g(t), mais on connaît l’allure de la réponse en
fréquence du filtre analogique prototype. G anal ( j.ω) est connu, ou au moins
un gabarit (technique de l’échantillonnage en fréquence).
On travaillera sur un exemple pour mettre en œuvre ces deux techniques : on cherche à
F
réaliser sous forme numérique un filtre passe-bas idéal de fréquence de coupure FC = E
10
/\
G anal

I
f
- FE/2 -FC FC FE/2

-a) Echantillonnage temporel (invariance impulsionnelle)

Soit g(t) la réponse impulsionnelle du filtre analogique :

+∞

∫ δ(t).e
− j.2 π.f .t
δ(t) 
→ ∆ (f ) =
Fourier
.dt = 1
−∞

1 si − FC < f < FC

Filtrage
→ G(f ) = G anal (f ).∆ (f ) = 
 0 si f ∉ [ − Fc ; FC ]
+∞

∫ G(f ).e
j.2 π.f .t

→ g(t) =
Fourier Inverse
.df
−∞
+ FC

∫e
j.2 π.f .t
= .df
− FC

e j.2 π.FC .t − e − j.2 π.FC .t

=
j.2.π.t
sin ( 2π.FC .t )
= 2.FC .
2π.FC .t
Cette réponse impulsionnelle est échantillonnée pour fournir la réponse impulsionnelle
du filtre numérique :
sin ( 2.π.FC .n.TE )
h(n) = TE .g(n.TE ) = 2.FC .TE .
2.π.FC .n.TE
 FE
Pour FC = , on obtient numériquement :
10
 π
sin  n. 
5
h(n) = 0, 2. 
π
n.
5
Deux problèmes surgissent :
1. la réponse impulsionnelle est infinie. Pour la rendre finie, on tronque cette
réponse (fenêtre rectangulaire) ; voir TP : on ne prend que 17 échantillons,
de n = - 8 à n = + 8.
2. le filtre n’est pas causal : il répond avant l’impulsion. Il n’est donc pas
réalisable physiquement. Pour le rendre causal (et donc réalisable), on
translate cette réponse de 8 échantillons : h(-8) devient h(0), h(-7) devient
h(1) … h(8) devient h(16).
L’ensemble des opérations de synthèse est résumé page suivante.
Valeurs des coefficients :
k 0 1 2 3 4 5 6 7 8
ak -0.037841 -0.043247 -0.031183 0 0.046774 0.10091 0.15137 0.1871 0.2
k 9 10 11 12 13 14 15 16
ak
11 0.1871 I 0.15137 I 0.10091 I 0.046774 I 0 I -0.031183 I -0.043247 I -0.037841 I I
Dans le tableau on remarque qu'on retrouve la propriété h(k) = h(16 - k) et
l’argumenvariera de façon linéaire traduisant un retard de 8.

Synthèse d’un filtre passe-bas « idéal » avec une fréquence (FC = FE/10) :

C'est l'échantillonnage temporel:

sin ( 2π.FC .t )
Filtre analogique de réponse impulsionnelle g(t) evec g(t) = 2.FC .
2π.FC .t
g(t)

Echantillonnage de la réponseimpulsionnelle avec un pas Te donne la courbe :

h(n) = TE .g(n.TE )

h(n)

n
Fenêtre rectangulaire w(t) ou w(n) qui tronque la réponse impulsionnelle.(limitation du nombre
d’échantillons)

w
1

Réponse impulsionnelle ci-dessous est limitée à unnombre fini d’échantillons par la multiplication de:

h(n).w(n)

h(n)

Décalage de la réponse impulsionnelle vers la droite pour rendre le filtre causal :

(filtre causal, c’est à dire qu'il est réalisable physiquement)

h(n)

n
 -b) Echantillonnage en fréquence

on a la réponse en fréquence du filtre analogique prototype.

On décide d’échantillonner en fréquence G anal ( j.ω) . On définit donc la suite :

 ω 
H(k) = G anal  j.k. E 
 N
où N est le nombre d’échantillons. En reprenant le même exemple que précédemment
N = 17, H(0) = H(−1) = H(1) = 1 et H(2) = H(−2) = ... = H(8) = H(−8) = 0
On peut en déduire, par transformation de Fourier discrète (TFD), les valeurs de la
réponse impulsionnelle :
2.π
1 N/2 j. .k.n
h(n) = . ∑ H(k).e N
N k = − N/2
Dans notre cas :
1  2.π
− j. .n j. .n 
2.π
1   2.π  
h(n) = . 1 + e 17
+ e 17  = . 1 + 2.cos  .n   avec − 8 ≤ n ≤ 8
17 17   17  
 
Enfin, pour rendre ce filtre réalisable physiquement, on translate cette réponse
impulsionnelle de 8 échantillons. L’ensemble de ces opérations est résumé page suivante.
La figure ci-dessous indique la réponse en fréquence de chaque filtre.

1
Filtre idéal analogique
Echantillonnage temporel
Echantillonnage en fréquence

0.8
|G Num(j.ω )|

0.6

0.4

Echantillons H(k)

0.2

0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
f/FE
exemple de Synthèse d’un filtre passe-bas « idéal » (FC = FE/10) : un échantillonnage en fréquence

G anal
Réponse en fréquence du filtre
analogique prototype
G anal ( j.ω)

f
H(k)
Echantillonnage en fréquence :
 ω 
H(k) = G anal  j.k. E 
 N

h(n) Transformée de Fourier Discrète :

2.π
1 N/2 j. .k.n
h(n) = . ∑ H(k).e N
N k = − N/2

h(n) Décalage
(filtre réalisable physiquement)

2.1.2 Synthèse des filtres R.I.I.

1- Les différentes méthodes de synthèse:
 Il existe plusieurs méthodes de synthèse des filtres RII, suivant le critère qu’on
s’impose :
a- la réponse temporelle imposée est:
- une réponse impulsionnelle alors la méthode est appelée méthode
d’invariance impulsionnelle
- réponse indicielle alors la méthode est appelée méthode d’invariance
indicielle
b- simulation numérique d’une équation différentielle on a la méthode d’Euler et
la méthode des trapèzes
c- la réponse en fréquence imposée alors la méthode est appelée transformation
bilinéaire
Dans ces méthodes, on cherche à réaliser un filtre numérique à partir d’un
1
filtre analogique prototype de fonction de transfert G anal (p) = , avec τ = TE
1 + τ.p

2- Invariance impulsionnelle
Dans cette méthode on procède de la manière siuvante:
1 1 − τt
G anal (p) = 
Transformée de Laplace
→ g(t) = .e .γ (t)
1 + τ.p inverse
τ
n
T  − E
T

Echantillonnage
→ h(n) = TE .g(n.TE ) = E .  e τ  .u(n)
τ  
T 1

Transformée en Z
→ H(z) = E .
τ −
TE
1 − e τ .z − 1

1
Remarque : Ganal(p) possède un pôle (valeur de p rendant Ganal(p) infini) : p = −
τ
TE
−
H(z) possède un pôle : z = e τ
= eTE .p , où p est le pôle de Ganal(p).

3- Invariance indicielle
1 1 1
G anal (p) = ⇒ Y(p) = G anal (p).Γ(p) = .
1 + τ.p p 1 + τ.p
Dans cette méthode, on fait la transformée de Laplace inverse de Y(p) qui donne y(t) puis
on échantillonne y(t) qui va donnée y(n) (numérique) et enfin la TR en Z de y(n).

 − 
t
 Laplace
→ y(t) =  1 − e  .γ (t)
τ

 
 −
n.TE


Echantillonnage
→ y(n) =  1 − e τ  .u(n)
 
1 1 1
 
Transformée en Z
→ Y(z) = − = H(z).U(z) = H(z).
1 − z−1 T
− E 1 − z−1
1− e τ
.z −1
On obtient donc l’expression de H(z) :
 T
 −1
1 − e τ
 .z
1 − z−1
H(z) = (1 − z − 1 ) .Y(z) = 1 − TE = TE

− −
1 − e τ .z −1 1 e τ
.z −1