FILTRES A REPONSE IMPULSIONNELLE FINIE 27
• Propriétés générales:
– toujours stables
– définis par les éléments de leur réponse
impulsionnelle
– causals si leur réponse impulsionnelle g(k) est
nulle pour k < 0. S ’ils ne le sont pas, il suffit
de décaler cette réponse.
– Peuvent être à phase linéaire
Signal Processing Laboratory
Swiss Federal Institute of Technology, Lausanne
FILTRES A PHASE LINEAIRE (1) 28
• Si la réponse impulsionnelle g(k) est non nulle
entre k0 et k0 + L - 1, la réponse en fréquence du
filtre est:
k0 + L −1
− j2πfk
G( f ) = ∑ g (k ) e
k = k0
• En séparant partie réelle et imaginaire:
k0 + L −1 k0 + L −1
G ( f ) = ∑ g (k )cos(2 fk ) − j ∑ g (k )sin (2 fk )
k = k0 k = k0
Signal Processing Laboratory
Swiss Federal Institute of Technology, Lausanne
FILTRES A PHASE LINEAIRE (2) 29
• Pour obtenir une réponse en phase identiquement
nulle, il faut que la partie imaginaire soit
identiquement nulle. Ceci implique:
g (k ) = g (− k ) pour k0 ≤ k ≤ k0 + L − 1
• Cette relation lie k et L. On obtient finalement:
L −1
g (k ) = g (− k ) pour k ≤
2
Signal Processing Laboratory
Swiss Federal Institute of Technology, Lausanne
FILTRES A PHASE LINEAIRE (3) 30
• Cette réponse n’est donc pas causale. On peut la
rendre causal en la décalant de (L - 1)/2 pas vers
les valeurs positives de l’axe des k. La réponse du
filtre après décalage est:
G ' ( f ) = G ( f )exp[− jf (L − 1)]
• La réponse en phase est donc:
θ g ( f ) = − f (L − 1)
Signal Processing Laboratory
Swiss Federal Institute of Technology, Lausanne
FILTRES A PHASE LINEAIRE (4) 31
• Et la relation entre éléments de la réponse devient:
g (k ) = g (L − 1 − k )
g(k) g(k)
L impair L pair
k k
Signal Processing Laboratory
Swiss Federal Institute of Technology, Lausanne
SYNTHESE PAR SERIE DE FOURIER (1) 32
• On sait que:
12
g (k ) = ∫ G ( f )exp( j 2πfk ) df
-1 2
• Ceci exprime tout simplement le fait que les g(k)
sont les coefficients de la série de Fourier de la
fonction périodique G(f) de période 1.
Signal Processing Laboratory
Swiss Federal Institute of Technology, Lausanne
SYNTHESE PAR SERIE DE FOURIER (2) 33
• Si le filtre à réaliser est d ’un type standard (passe-bas,
…) présenté auparavant, ou si la forme de G(f) n ’est pas
trop compliquée, on peut obtenir facilement la forme
analytique de g(k).
• Par contre, si G(f) (ou ses dérivées) présente des
discontinuités, g(k) sera de longueur infinie. Tronquer
cette réponse entraîne l’apparition d’oscillations dans la
réponse en fréquence modifiée (phénomène de Gibbs).
• Il convient alors de multiplier la réponse impulsionnelle
tronquée par une fonction fenêtre.
Signal Processing Laboratory
Swiss Federal Institute of Technology, Lausanne
SYNTHESE PAR SERIE DE FOURIER (3) 34
• Exemple: filtre passe-bande, fréquences de
coupure 0.15 et 0.35
L = 31
L = 61
L = 61
+ Hamming
Signal Processing Laboratory
Swiss Federal Institute of Technology, Lausanne
APPROXIMATION NUMERIQUE (1) 35
• Il n ’est parfois pas possible d ’obtenir une
expression analytique pour g(k).
• On approxime la valeur de l ’intégrale par la
somme des surfaces de rectangles étroits
juxtaposés. Si N est le nombre de rectangles:
N 2−1
gˆ (k ) =
1 n nk
∑ G exp j2
N n=− N 2 N N
qui n ’est autre que la TFD inverse de G(n).
Signal Processing Laboratory
Swiss Federal Institute of Technology, Lausanne
APPROXIMATION NUMERIQUE (2) 36
• On a la relation:
+∞
ĝ (k ) = ∑ g (iN + k )
i =−∞
• Pour N suffisamment grand par rapport à la durée
présumée L de la réponse impulsionnelle:
gˆ (k ) ≈ g (k )
Signal Processing Laboratory
Swiss Federal Institute of Technology, Lausanne
ECHANTILLONNAGE EN FREQUENCE (1) 37
• Dénotons par Gd(f) la réponse fréquentielle
désirée, représentée par L échantillons Gd(n)
prélevés périodiquement avec une période 1/L sur
[-1/2, 1/2].
• Ces échantillons caractérisent complètement les L
échantillons de la réponse impulsionnelle:
1 L 2−1 2
g (k ) = ∑ Gd (n )exp j nk k0 ≤ k ≤ k0 + L − 1
L n=− L 2 L
Signal Processing Laboratory
Swiss Federal Institute of Technology, Lausanne
ECHANTILLONNAGE EN FREQUENCE (2) 38
• On sait qu ’il est possible d ’obtenir la transformée
en z par interpolation de ces échantillons Gd(n).
Avec k0 = 0:
− L L 2−1
1 − z 1
G (z ) = ∑ Gd (n )
L n=− L 2 1 − exp( j2 n L )z −1
• En évaluant G(z) sur le cercle unité, z = exp(j2πf),
on obtient la réponse fréquentielle G(f).
Signal Processing Laboratory
Swiss Federal Institute of Technology, Lausanne
ECHANTILLONNAGE EN FREQUENCE (3) 39
• On obtient:
− jf (1− L ) L 2−1
e sin (fL ) − j n L
G( f ) = ∑ Gd (n ) e
L n=− L 2 sin[( f − n L )]
• qui est une interpolation du type:
L 2−1
n
G ( f ) = ∑ Gd (n )ψ f −
n=− L 2 L
Signal Processing Laboratory
Swiss Federal Institute of Technology, Lausanne
ECHANTILLONNAGE EN FREQUENCE (4) 40
• Avec:
n sin (fL )
ψ f − = exp{j[( f − n L ) − fL ]}
L L sin[( f − n L )]
Parties réelle
et imaginaire
de ψ(f-n/L)
pour -1/2 < f < 1/2
L = 8, n = 1, n/L = 0.125
Signal Processing Laboratory
Swiss Federal Institute of Technology, Lausanne
ECHANTILLONNAGE EN FREQUENCE (5) 41
• Cette fonction d ’interpolation a les propriétés
suivantes:
n k
ψ f − = 0 pour f = avec k ≠ n
L L
n
lim ψ f − = 1
n L
f→
L
Signal Processing Laboratory
Swiss Federal Institute of Technology, Lausanne
ECHANTILLONNAGE EN FREQUENCE (6) 42
• Exemple: conception d ’un filtre passe-haut,
fréquence de coupure fc = 0.3125, L = 16
Réponse désirée
Réponse
impulsionnelle
non causale
Signal Processing Laboratory
Swiss Federal Institute of Technology, Lausanne
ECHANTILLONNAGE EN FREQUENCE (7) 43
Réponse fréquentielle
obtenue f
Réponse fréquentielle
obtenue (dB) f
Signal Processing Laboratory
Swiss Federal Institute of Technology, Lausanne
ECHANTILLONNAGE EN FREQUENCE (8) 44
• Même conception, mais avec L = 64
Réponse fréquentielle
obtenue
Réponse fréquentielle
obtenue (dB)
Signal Processing Laboratory
Swiss Federal Institute of Technology, Lausanne
ECHANTILLONNAGE EN FREQUENCE (9) 45
• Même conception, mais transition adoucie
Réponse fréquentielle
obtenue
Réponse fréquentielle
obtenue (dB)
Signal Processing Laboratory
Swiss Federal Institute of Technology, Lausanne