Introduction aux filtres de réponse à l'impulsion à l'infini (RII)

Les filtres Infinite Impulse Response (IIR) sont une classe fondamentale de filtres numériques caractérisés par un mécanisme de rétroaction qui permet à la sortie de dépendre des entrées actuelles et précédentes ainsi que des sorties précédentes. Cette structure récursive permet aux filtres IIR d'effectuer des transitions de fréquence pointues avec des coefficients beaucoup moins élevés que leurs homologues de la réponse Finite Impulse (FIR).

La conception des filtres IIR commence généralement par un filtre prototype analogique, souvent le type Butterworth, Chebyshev ou elliptique (Cauer), qui est ensuite transformé en domaine numérique par des techniques comme la transformation bilinéaire ou l'invariance des impulsions. Chaque prototype offre des compromis spécifiques entre la planéité des bandes passantes, l'atténuation des bandes stoppées, la bande passante et la linéarité de phase.

Concepts de base dans la conception de filtres IIR

Fonction de réponse et de transfert de fréquence

Un filtre IIR est décrit par sa fonction de transfert dans le domaine z:

Les coefficients (feedback) et (feedforward) déterminent la magnitude du filtre et la réponse de phase. L'ordre du filtre N (le nombre de pôles) influence fortement la raideur du roulis de la bande passante à la bande stop.

Spécifications de conception clés

  • Frequence de bord de bande passante (-p[):[la fréquence supérieure (ou inférieure) où la réponse est dans une tolérance spécifiée.
  • Frequence de bord de bande ( .s[[[[FLT:]][[FLT:]][[FLT:]][[FLT:]][FLT:]][[FLT:]][FLT:]][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:]][FLT:][FLT:]][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:]][FLT:][F][FLT:][F][F
  • Onduleur de bande passante (Rp[):[ l'écart maximal admissible par rapport au gain idéal dans la bande passante (en dB).
  • Atténuation de bande d'arrêt (As):[Atténuation minimale requise dans le bouchon (en dB).
  • Ordre des filtres (N): détermine la complexité du calcul et la raideur de la transition.

Prototype analogique et transformation numérique

La plupart des méthodes de conception de l'IIR commencent par un filtre prototype analogique à passe-bas qui répond aux spécifications de magnitude. La fonction de transfert analogique H(s) est ensuite convertie en un H(z) en utilisant une transformation de magnitude. La méthode la plus courante pour préserver la réponse de magnitude est la transformation bilinéaire, qui mastique l'ensemble s-plan au z-plan par ]s = (2/T) * (1 - z^-1)/(1 + z^-1).

Filtres Butterworth: réponse maximale à plat

Caractéristiques

Les filtres Butterworth sont conçus pour fournir une réponse de fréquence aussi plate que possible dans la bande passante, sans ondulations. La réponse au carré de magnitude est donnée par:

où - c est la fréquence de coupure (généralement le point -3 dB) et -N est l'ordre du filtre.

Avantages

  • Place de sortie:[ Idéal pour les applications qui doivent préserver l'amplitude du signal sur une large gamme de fréquences, telles que les systèmes audio haute fidélité et la surveillance médicale.
  • Réponse de phase prévisible:[ Le déplacement de phase varie linéairement avec la fréquence à l'intérieur de la bande passante, minimisant la distorsion de la forme d'onde pour les signaux transitoires.
  • Équations de conception simples:[ Les emplacements des pôles sont facilement calculés sur un cercle de rayon --c dans le plan --, permettant des conceptions simples à haut ordre.

Limitations

  • Modérer le roulage:[ Pour obtenir une transition nette, il faut un ordre de filtre élevé, augmentant le coût de calcul et le retard du groupe.
  • Peu de pente d'atténuation du bandeau d'arrêt:[ La pente asymptotique est -20N dB/décennie, qui est moins agressive que les filtres Chebyshev ou elliptiques pour le même ordre.

Étapes de conception

  1. Préciser le bord de la bande passante, le bord de la bande arrêt, l'ondulation de la bande passante (généralement 0,5–3 dB) et l'atténuation minimale de la bande arrêt.
  2. Calculer l'ordre de filtre requis N en utilisant .
  3. Déterminer les emplacements des pôles analogiques : les pôles sont uniformément espacés sur la moitié gauche du cercle de rayon - - c.
  4. Convertir en filtre numérique en utilisant une transformation bilinéaire avec pré-warping.
  5. Mettre en place comme cascade de sections de deuxième ordre (SOS) pour la stabilité numérique.

Des tableaux détaillés et des formules de calcul en forme fermée pour les coefficients de Butterworth sont largement disponibles; voir par exemple la référence autorisée Dispositifs en analogie -Guide de conception du filtre IIR.

Filtres Chebyshev: Décollage de Steeper avec des ripelles

Type I et type II

Les filtres Chebyshev échangent la planéité de la bande passante pour une transition plus raide. Il y a deux variantes:

  • Chebyshev Type I: Les ipples sont autorisés dans la bande passante, tandis que le stopband est monotonique. La réponse de magnitude est , où T N est le polynôme de l'ordre de Chebyshev N et ε est un facteur d'ondulation dérivé de l'ondulation de la bande passante autorisée R p.
  • Chebyshev Type II (inverse Chebyshev):[ La bande de passe est plate, mais des ondulations apparaissent dans le bouchon. Ce dessin est moins fréquent parce que les nuls de la bande de stop peuvent être indésirables dans certaines applications.

Avantages

  • Faster roll-off que Butterworth: Les filtres Chebyshev obtiennent la même pente de transition avec un ordre inférieur, réduisant la charge de calcul.
  • Bien adapté aux spécifications où une petite ondulation de bande passante est acceptable:[ Les valeurs typiques d'ondulation de bande passante varient de 0,1 à 3 dB. De nombreux systèmes de communication tolèrent ces ondulations en échange d'un filtrage plus net.

Limitations

  • Réponse de phase non linéaire dans la bande passante:[ Les ondulations introduisent une distorsion de phase, qui peut dégrader la qualité des signaux qui dépendent de la cohérence de phase (p. ex., formes d'onde modulées).
  • La magnitude du radeau doit être soigneusement choisie :[ Une ondulation plus élevée donne un roll-off plus raide, mais augmente la variation de retard du groupe et l'anneauage potentiel dans le domaine temporel.

Approche de conception

  1. Définir le bord du bandeau, le bord du bandeau, l'onduleur du bandeau (p. ex. 1 dB) et l'atténuation du bandeau.
  2. Calculer l'ordre du filtre : .
  3. Trouvez les pôles du prototype analogique : les pôles sont situés sur une ellipse dans le plan s, avec l'axe principal le long de l'axe réel pour le type I.
  4. Carte numérique en utilisant la transformation bilinéaire avec pré-warping de fréquence.
  5. Se décomposer en sections SOS pour éviter la sensibilité aux coefficients.

Pour un exemple pratique utilisant la bibliothèque Python , consultez SciPy , documentation de conception de filtre IIR.

Filtres elliptiques: Bande passante de transition optimale

Caractéristiques

Les filtres elliptiques (ou Cauer) permettent le roll-off le plus abrupt possible entre le bandeau et le bouchon pour un ordre de filtre donné en permettant des ondulations dans les deux bandes. La réponse carrée de magnitude implique des fonctions elliptiques jacobins, ce qui rend le design mathématiquement plus impliqué.

L'avantage clé est que pour un ordre fixe N, le filtre elliptique offre la bande de transition la plus courte. Cette propriété en fait le filtre de choix dans les applications nécessitant des coupures extrêmement pointues, comme les filtres anti-aliasing dans les récepteurs de communication ou la sélection de canaux à bande étroite.

Avantages

  • Roulements les plus agressifs :[ Atteint une bande passante donnée avec l'ordre le plus bas possible par rapport à Butterworth ou Chebyshev.
  • Compensation flexible:[ Permet de spécifier indépendamment le taux d'ondulation de la bande passante, l'atténuation de la bande arrêt et le rapport de transition.

Limitations

  • Les ipples dans la bande passante et la bande stop: peuvent introduire une distorsion inacceptable du signal dans les applications de mesure audio ou de précision à haute fidélité.
  • Conception complexe: Nécessite une évaluation de la fonction elliptique et des algorithmes spécialisés. La plupart des concepteurs comptent sur un logiciel de conception de filtre (MATLAB, Octave, ou des outils dédiés).
  • Plus grande sensibilité aux coefficients:[ La transition étroite peut entraîner des effets de quantification lors de la mise en œuvre sur les processeurs à points fixes.

Procédure de conception

  1. Spécifiez les quatre paramètres : - p, - s, R p, A s.
  2. Calculer le facteur de sélectivité k = γ p/γ s et le facteur de discrimination k1 = ε/ √(10^(As/10)-1).
  3. Trouver la commande N[ en utilisant les fonctions intégrales elliptiques: , où K[ est l'intégrale elliptique complète du premier genre.
  4. Déterminez les emplacements des pôles et zéro des fonctions elliptiques jacobins; les zéros sont situés sur l'axe imaginaire dans le bouchon pour les filtres elliptiques de type I.
  5. Appliquer la transformation bilinéaire pour obtenir la fonction de transfert numérique.
  6. Mettre en œuvre comme sections SOS, en veillant à ce que les zéros (qui sont réels pour la version numérique) soient correctement appariés pour éviter de prendre un pic.

Un traitement complet de la conception du filtre elliptique peut être trouvé dans MATLAB=s documentation de fonction ellip, qui fournit des routines intégrées pour générer des coefficients.

Étapes pratiques de conception communes à tous les types

1. Collecte des spécifications

Commencez par recueillir les exigences d'application : fréquence d'échantillonnage, bande passante et fréquences de bande arrêt (normalisées à la moitié du taux d'échantillonnage), ondule de bande passante admissible et atténuation requise du bande arrêt. Par exemple, un filtre anti-aliasage dans un système audio de 48 kHz peut nécessiter une ondule de bande passante inférieure à 0,1 dB et une atténuation du bande arrêt de 80 dB à 24 kHz.

2. Génération de prototypes analogiques

En utilisant des outils ou des calculs manuels, générer le prototype de passe-bas analogique qui correspond aux spécifications de magnitude aux fréquences normalisées. Pré-chauffer les fréquences critiques parce que la transformation bilinéaire compresse l'axe de fréquence.

3. Conversion numérique

Pour les filtres à haut ordre (N > 2), convertissez la fonction de transfert en cascade de sections de second ordre (SOS). Cette décomposition améliore grandement la stabilité numérique dans les implémentations à points fixes et flottants.

4. Mise en œuvre et vérification

Implémentez le filtre dans votre environnement cible (C, Python, HDL, etc.). Vérifiez la réponse de fréquence au moyen d'un outil d'analyse spectrale et assurez-vous que l'onduleur de bande passante, l'atténuation de bande stop et la bande passante de la transition répondent aux spécifications originales.

Comparaison des filtres Butterworth, Chebyshev et Elliptic

Property Butterworth Chebyshev (Type I) Elliptic
Passband Response Monotonic, maximally flat Ripples present Ripples present
Stopband Response Monotonic Monotonic Ripples present
Roll-off Sharpness Lowest (N-dependent) Medium Highest
Filter Order for Given Specs Highest required Lower than Butterworth Lowest required
Phase Linearity Best (near linear) Moderate Poor (significant nonlinearity)
Complexity of Design Low Moderate High

Outils et logiciels pour la conception de filtres IIR

Les flux de travail modernes du traitement des signaux numériques reposent sur des outils de haut niveau qui automatisent la génération de coefficients de filtrage.

  • Python avec SciPy:[ Fonctions telles que , , et retournent directement les matrices SOS à fréquences normalisées. Exemple : .
  • MATLAB / Octave: Commandes comme , , et calculateur et dénominateur polynômes. La fonction se convertit en SOS.
  • Des packs de conception de filtres spécialisés : des outils gratuits et commerciaux tels que FilterPro (Texas Instruments), Nuhertz Filter Solutions ou AADE Filter Design offrent des spécifications graphiques interactives et une sortie de coefficient réelle.

Pour une comparaison exhaustive de ces environnements de conception, veuillez consulter une enquête de 2021 sur les outils de conception de filtres numériques (ResearchGate).

Conseils d'applications pour le monde réel

Le choix des trois types de filtres dépend des contraintes spécifiques du système:

  • Processus audio haute fidélité:[ Utilisez des filtres Butterworth pour préserver la forme du signal d'origine avec une distorsion de phase minimale. Le léger compromis dans le roulage est compensé par une qualité sonore transparente.
  • Choix des canaux de télécommunications: Les filtres elliptiques sont préférés lorsque le rejet strict des canaux adjacents doit être réalisé avec des ressources matérielles limitées. La distorsion de phase est tolérable en signaux modulés en amplitude.
  • Conditionnement du signal biomédical (p. ex., ECG):[ Les filtres Chebyshev de type II offrent une bande passante plate pour éviter de déformer le complexe QRS, tout en rejetant efficacement les interférences de 50/60 Hz.
  • Systèmes de contrôle: Les filtres Butterworth sont souvent choisis pour leur réponse de phase prévisible et leur manque de dépassement, ce qui est critique pour les marges de stabilité du système.

Considérations avancées

Effets finaux de la longueur de mots

Lors de la mise en œuvre des filtres IIR sur les processeurs à points fixes (par exemple 16 bits ou 32 bits), la quantification des coefficients peut déplacer les emplacements des pôles, ce qui peut entraîner une instabilité ou des performances dégradantes.

Conceptions de passe haut, passe-bande et arrêt-bande

Les trois types de prototypes peuvent être transformés de lowpass en d'autres réponses sélectives par des transformations spectrales (replacer la z-variable avec une fonction rationnelle appropriée). Par exemple, une transformation lowpass-to-bandpass double l'ordre du filtre et est facilement implémentée dans les outils logiciels.

Comparaison avec les filtres FIR

Les filtres RIP offrent une phase linéaire et une stabilité garantie, mais nécessitent des commandes beaucoup plus élevées pour correspondre aux performances de roulement de l'IIR. Dans les applications où la phase n'est pas critique, le coût de calcul plus faible des filtres RII – en particulier elliptiques – en fait le choix préféré.

Conclusion

Les filtres Butterworth offrent la conception la plus simple et la réponse la plus prévisible en phase, ce qui en fait un choix sûr pour les applications audio et basse fréquence. Les filtres Chebyshev offrent un roll-off amélioré avec un ondulateur contrôlé, particulièrement utile dans les communications. Les filtres elliptiques offrent la sélectivité de fréquence la plus agressive pour l'ordre le plus bas, mais avec une complexité de conception accrue et un ondulateur de bande.

L'ingénieur moderne peut utiliser des bibliothèques logicielles sophistiquées pour prototyper rapidement ces filtres, laissant le temps de se concentrer sur les compromis et la vérification au niveau du système. Comprendre les bases mathématiques garantit que le filtre sélectionné répond à la fois aux spécifications de la fréquence-domaine et aux contraintes du monde réel de votre système de traitement de signaux numériques.

Pour plus de détails sur la théorie du design et la dérivation de ces filtres, voir Julius O. Smith="s online book on digital filters et les textes classiques d'Oppenheim et Schafer. De plus, le rapport d'application Texas Instruments ="Design of IIR Digital Filters=" fournit un guide pratique orienté vers la mise en œuvre.