Synthèse des Filtres à Réponse Impulsionnelle Infinie (RII)
2.1 Méthodologie générale de synthèse des filtres RII
Introduction
La synthèse d'un filtre numérique RII consiste à déterminer les coefficients d'une fonction de transfert
numérique H(z) à partir de spécifications fréquentielles données. Cette méthode s'appuie sur la
conception de filtres analogiques bien établis, puis leur transformation vers le domaine numérique.
a) Étapes successives de la synthèse d'un filtre RII
1. Spécification du gabarit
Objectif : Définir précisément les caractéristiques fréquentielles souhaitées du filtre.
Paramètres à spécifier :
• Fréquence d'échantillonnage fs (Hz)
• Type de filtre : passe-bas, passe-haut, passe-bande, coupe-bande
• Fréquences critiques :
• fp : fréquence de coupure de la bande passante (Hz)
• fa : fréquence de coupure de la bande atténuée (Hz)
• Tolérances d'amplitude :
• Ap : ondulation maximale dans la bande passante (dB)
• Aa : atténuation minimale dans la bande atténuée (dB)
Exemple pour un filtre passe-bas :
fs = 1000 Hz
fp = 100 Hz (bande passante : 0 ≤ f ≤ 100 Hz)
fa = 150 Hz (bande atténuée : f ≥ 150 Hz)
Ap = 1 dB (ondulation ≤ 1 dB en bande passante)
Aa = 40 dB (atténuation ≥ 40 dB en bande atténuée)
2. Normalisation des fréquences et établissement du gabarit normalisé
Objectif : Convertir les spécifications vers des fréquences normalisées pour faciliter la conception.
Fréquences normalisées :
• Ωp = 2πfp/fs (rad) : fréquence normalisée de la bande passante
• Ωa = 2πfa/fs (rad) : fréquence normalisée de la bande atténuée
Pré-distorsion bilinéaire : Pour compenser la distorsion introduite par la transformation bilinéaire :
ωp = (2/T) × tan(Ωp/2) [rad/s]
ωa = (2/T) × tan(Ωa/2) [rad/s]
où T = 1/fs est la période d'échantillonnage.
Gabarit normalisé : Le gabarit est maintenant exprimé en termes de fréquences analogiques normalisées
(ωp, ωa).
3. Approximation analogique
Objectif : Sélectionner et calculer une fonction de transfert analogique Ha(s) qui satisfait le gabarit
normalisé.
Types d'approximations courantes :
3.1 Approximation de Butterworth
• Caractéristiques : Réponse maximalement plate en bande passante
• Fonction de transfert :
|Ha(jω)|² = 1 / (1 + (ω/ωc)^(2n))
• Ordre requis :
n ≥ log₁₀[(10^(Aa/10) - 1)/(10^(Ap/10) - 1)] / (2×log₁₀(ωa/ωp))
3.2 Approximation de Chebyshev Type I
• Caractéristiques : Ondulation équiripple en bande passante, monotone en bande atténuée
• Fonction de transfert :
|Ha(jω)|² = 1 / (1 + ε²×Tn²(ω/ωp))
où Tn(x) est le polynôme de Chebyshev de première espèce d'ordre n.
3.3 Approximation de Chebyshev Type II
• Caractéristiques : Monotone en bande passante, ondulation équiripple en bande atténuée
3.4 Approximation elliptique (Cauer)
• Caractéristiques : Ondulation équiripple dans les deux bandes
• Avantage : Ordre minimal pour un gabarit donné
Calcul des pôles et zéros : Une fois le type et l'ordre choisis, on calcule les pôles (et éventuellement les
zéros) de Ha(s).
4. Dénormalisation de la fonction de transfert analogique
Objectif : Adapter la fonction de transfert normalisée aux spécifications fréquentielles réelles.
Transformations de fréquence :
4.1 Filtre passe-bas → passe-bas
s → s/ωc
4.2 Filtre passe-bas → passe-haut
s → ωc/s
4.3 Filtre passe-bas → passe-bande
s → (s² + ωh×ωb)/(s×B)
où B = ωh - ωb est la bande passante et ωh, ωb sont les fréquences haute et basse.
4.4 Filtre passe-bas → coupe-bande
s → s×B/(s² + ωh×ωb)
5. Transformation vers le domaine numérique
Objectif : Convertir la fonction de transfert analogique Ha(s) en fonction de transfert numérique H(z).
5.1 Transformation bilinéaire (méthode privilégiée)
s = (2/T) × (1-z⁻¹)/(1+z⁻¹)
Avantages :
• Préserve la stabilité
• Mapping biunivoque entre plans s et z
• Pas d'aliasing
Inconvénient :
• Distorsion fréquentielle (compensée par la pré-distorsion)
5.2 Méthode de l'invariance impulsionnelle
H(z) = Σ Res[Ha(s)×(1-e^(siT)z⁻¹)⁻¹, si]
où si sont les pôles de Ha(s).
Avantages :
• Préserve la forme de la réponse impulsionnelle
• Pas de distorsion fréquentielle
Inconvénient :
• Possible aliasing pour les filtres passe-haut et passe-bande
5.3 Transformation par correspondance des dérivées
Approximation de la dérivée par différences finies :
s → (1-z⁻¹)/T
Résultat final : La fonction de transfert numérique s'écrit :
H(z) = (b₀ + b₁z⁻¹ + ... + bₘz⁻ᵐ)/(1 + a₁z⁻¹ + ... + aₙz⁻ⁿ)
Considérations pratiques
Vérification de la stabilité
• Tous les pôles de H(z) doivent être à l'intérieur du cercle unité
• Vérification par calcul des racines du dénominateur
Analyse de la réponse en fréquence
• Calcul de H(e^(jω)) pour ω ∈ [0, π]
• Vérification du respect du gabarit spécifié
Optimisation des performances
• Choix de la structure de réalisation (forme directe, cascade, parallèle)
• Analyse des effets de quantification
• Considération de la complexité calculatoire
Exemple complet
Spécifications :
• Filtre passe-bas numérique
• fs = 1000 Hz, fp = 100 Hz, fa = 150 Hz
• Ap = 1 dB, Aa = 40 dB
Étape 1 : Normalisation
• Ωp = 2π×100/1000 = 0.2π rad
• Ωa = 2π×150/1000 = 0.3π rad
Étape 2 : Pré-distorsion
• ωp = 2×tan(0.1π) ≈ 0.65 rad/s
• ωa = 2×tan(0.15π) ≈ 1.0 rad/s
Étape 3 : Choix Butterworth
• n = ⌈log₁₀(999/0.259)/(2×log₁₀(1.0/0.65))⌉ = 8
Étape 4 : Calcul des coefficients Application de la transformation bilinéaire pour obtenir les coefficients
{ai, bi}.
Cette méthodologie garantit une synthèse rigoureuse et optimale des filtres RII.
