Introduzione ai filtri di risposta dell'impulso infinito (IIR)

I filtri Infinite Impulse Response (IIR) sono una classe fondamentale di filtri digitali caratterizzati da un meccanismo di feedback che permette all'output di dipendere sia dagli input attuali che da quelli precedenti, sia dalle uscite precedenti. Questa struttura ricorsiva consente ai filtri IIR di raggiungere transizioni di frequenza taglienti con meno coefficienti rispetto ai loro omologhi Finite Impulse Response (FIR) e quindi i filtri IIR sono elaborati in tempo reale come strumenti di elaborazione audio efficiente e ampiamente implementati in tempo reale.

Il design dei filtri IIR inizia tipicamente con un filtro analogico del prototipo – spesso il tipo Butterworth, Chebyshev o ellittico (Cauer) – che viene poi trasformato nel dominio digitale attraverso tecniche come la trasformazione bilineare o l'invarianza degli impulsi.

Concetti core in IIR Filter Design

Funzione di risposta e trasferimento di frequenza

Un filtro IIR è descritto dalla sua funzione di trasferimento nel z]-dominio:

I coefficienti (feedback) e [ (feedforward) determinano la magnitudine del filtro e la risposta di fase. L'ordine del filtro N[]] (il numero di poli) influenza pesantemente la ripidazza del roll-off da passband a stopband.

Specifiche di progettazione chiave

  • Frequenza del bordo della banda ([[]]ω]p]][]] la frequenza superiore (o inferiore) in cui la risposta è all'interno di una tolleranza specificata.
  • Frequenza del bordo della fascia ([[]]][]]]]][[]]]]]][[]]]]] la frequenza in cui l'attenuazione incontra un livello minimo richiesto.
  • Spaccamento del passabanda ([[]R]][]]]]][] la deviazione massima consentita dal guadagno ideale nella banda passante (in dB).
  • A]][]][[[]]]]]]][] l'attenuazione minima necessaria nella banda di arresto (in dB).
  • Ordine di deposito (N):] determina la complessità computazionale e la ripidità di transizione.

Prototipo analogico e trasformazione digitale

La maggior parte dei metodi di progettazione IIR inizia con un filtro analogico a basso passaggio che soddisfa le specifiche di magnitudo. La funzione di trasferimento analogico H(s)] viene poi convertita in un digitale H(z)] usando una trasformazione a frequenza mappante.

Filtri Butterworth: Risposta massimamente piatta

Caratteristiche

I filtri Butterworth sono progettati per fornire una risposta di frequenza che è il più piatta possibile nella banda passante, senza rigorosi. La risposta quadrata di magnitudo è data da:

dove ω c[] è la frequenza di taglio (tipicamente il punto -3 dB) e ]N[] è l'ordine del filtro. Come ]]N[]] aumenta, la transizione diventa più ripida, ma la banda passante rimane monotonica.

Vantaggi

  • Flat passband:[] Ideale per applicazioni che devono preservare l'ampiezza del segnale attraverso un ampio range di frequenza, come sistemi audio ad alta fedeltà e monitoraggio medico.
  • Risposta di fase prevedibile:[ Il turno di fase varia linearmente con la frequenza all'interno della banda passante, minimizzando la distorsione delle forme d'onda per i segnali transitori.
  • Equazioni di design semplici:[] Le posizioni dei polacchi sono facilmente calcolate su un cerchio di raggio ω c nel []s]-plano, consentendo disegni di ordine elevato semplici.

Limitazioni

  • Moderate roll-off:[] Con una transizione affilata richiede un alto ordine di filtro, aumentando i costi computazionali e il ritardo di gruppo.
  • La pendenza di attenuazione della banda di arresto della pozza: Il pendio asintotico è -20N dB/decensione, che è meno aggressivo di Chebyshev o filtri ellittici per lo stesso ordine.

Passi di progettazione

  1. Specificare bordo passband, bordo stopband, banda passante (tipicamente 0.5–3 dB), e l'attenuazione minima della banda di stopband.
  2. Compiti l'ordine del filtro richiesto ]N] usando .
  3. Determinare le posizioni dei pali analogici: i poli sono uniformemente distanziati sulla metà sinistra del cerchio del raggio ω c.
  4. Convertirsi in filtro digitale utilizzando la trasformazione bilineare con pre-riscaldamento.
  5. Implementazione come cascata di sezioni di secondo ordine (SOS) per stabilità numerica.

Tavoli di progettazione dettagliate e formule chiuse per i coefficienti Butterworth sono ampiamente disponibili; vedi ad esempio il riferimento autorevole [ Guida di progettazione filtro IIR di Analog Devices[.

Filtri Chebyshev: Steeper Roll-Off con Ripples

Tipo I e Tipo II

I filtri Chebyshev commerciano la flatnesse passband per una transizione più ripida. Ci sono due varianti:

  • Chebyshev Tipo I:] I ricci sono ammessi nella banda passante, mentre la banda ferma è monotonica. La risposta di magnitudo è , dove T N è il polinomio di Chebyshev N[Fpleband:
  • Chebyshev Tipo II (inverso Chebyshev): La banda passante è piatta, ma le ricrescenze appaiono nella banda di arresto. Questo disegno è meno comune perché la banda di stopband non è indesiderabile in alcune applicazioni.

Vantaggi

  • La maggior parte dei filtri di Butterworth:[ Chebyshev raggiunge la stessa pendenza di transizione con un ordine inferiore, riducendo il carico computazionale.
  • Ben adattato per specifiche in cui una piccola fascia passante è accettabile:[] I valori tipici della banda passante coprono variano da 0,1 a 3 dB. Molti sistemi di comunicazione tollerano queste increspature in cambio di filtraggio più affilato.

Limitazioni

  • Risposta di fase non lineare nel passband:[ Le ondulazioni introducono distorsioni di fase, che possono degradare la qualità dei segnali che dipendono dalla coerenza di fase (ad esempio, forme d'onda modulate).
  • La magnitudine del rifornimento deve essere accuratamente scelta:[ L'increspatura più alta produce un roll-off più ripido ma aumenta la variazione del ritardo del gruppo e il potenziale ringing nel dominio del tempo.

Approccio di progettazione

  1. Definire il bordo del passband, il bordo del stopband, il passband ripple (ad esempio, 1 dB), e l'attenuazione della banda di stopband.
  2. Calcola l'ordine del filtro: .
  3. Trova i poli del prototipo analogico: i poli si trovano su un ellisse nel ]]-plano, con l'asse maggiore lungo l'asse reale per il tipo I.
  4. Mappa al digitale utilizzando la trasformazione bilineare con la frequenza pre-riscaldamento.
  5. Decomporre in sezioni SOS per evitare la sensibilità del coefficiente.

Per un esempio pratico utilizzando la libreria SciPy di Python, consultare []SciPy’s IIR filter design documentazione[.

Filtri ellittici: Larghezza di banda di trasmissione ottimale

Caratteristiche

I filtri ellittici (o Cauer) raggiungono il più ripido possibile roll-off tra il passband e il stopband per un dato ordine di filtro, permettendo le increspature in entrambe le bande. La risposta quadrata di magnitudine coinvolge le funzioni ellittiche giaconiche, rendendo il design matematicamente più coinvolto.

Il vantaggio chiave è che per un ordine fisso N], il filtro ellittico offre la banda di transizione più breve. Questa proprietà lo rende il filtro di scelta in applicazioni che richiedono tagli estremamente taglienti, come filtri anti-aliasing nei ricevitori di comunicazione o la selezione di canali a banda stretta.

Vantaggi

  • La maggior parte dei roll-off aggressivi:] Consegui una data larghezza di banda di transizione con il più basso ordine possibile rispetto a Butterworth o Chebyshev.
  • Commercio flessibile:[ Permette specifiche indipendenti di ondulazione del passband, attenuazione della stopband e rapporto di transizione.

Limitazioni

  • Ripples in passband e stopband:[[] Può introdurre una distorsione del segnale inaccettabile in applicazioni di misurazione audio ad alta fedeltà o di precisione.
  • Design complesso:[[] Richiede una valutazione delle funzioni ellittiche e algoritmi specializzati. La maggior parte dei progettisti si affidano al software di progettazione dei filtri (MATLAB, Octave, o strumenti dedicati).
  • Sensibilità del coefficiente di gravità:[ La transizione stretta può portare ad effetti di quantizzazione quando si implementano su processori a punto fisso.

Procedura di progettazione

  1. Specificare tutti e quattro i parametri: ω p, ω s, R p, A s.
  2. Compiti il fattore di selettività [k = ω p/ω s[ e il fattore di discriminazione k1 = ε/√(10^(As/10)-1)]].
  3. Trovare l'ordine N[]] utilizzando funzioni integrali ellittiche: [], dove [K]] è l'integrale ellittico completo del primo tipo.
  4. Determinare il polo e zero posizioni dalle funzioni ellittiche giacobiche; gli zeri si trovano sull'asse immaginario nella banda di arresto per i filtri ellittici di tipo I.
  5. Applicare la trasformazione bilineare per ottenere la funzione di trasferimento digitale.
  6. Implementa le sezioni SOS, assicurando che gli zeri (che sono reali per la versione digitale) siano adeguatamente abbinati per evitare di ottenere il picco.

Un trattamento completo del design dei filtri ellittici si può trovare nella documentazione della funzione ellip di MATLAB[], che fornisce routine integrate per la generazione dei coefficienti.

Pratico Design Passi comuni a tutti i tipi

1. Specificazione della collezione

Iniziare raccogliendo i requisiti applicativi: frequenza di campionamento, frequenza di banda passante e frequenza di stopband (normalizzata a metà della frequenza di campionamento), ripple passband ammissibile e attenuazione della banda di stopband richiesta. Ad esempio, un filtro anti-aliasing in un sistema audio 48 kHz potrebbe richiedere una banda passante al di sotto di 0,1 dB e un'attenuazione della banda di 80 dB a 24 kHz.

2. Generazione di prototipi analogici

Utilizzando strumenti o calcolo manuale, generare il prototipo analogico di passaggio basso che corrisponde alle specifiche di magnitudo alle frequenze normalizzate. Pre-warp le frequenze critiche perché il bilineare trasforma compressa l'asse di frequenza.

3. Conversione digitale

Per i filtri di alto ordine (N > 2), converte la funzione di trasferimento in una cascata di sezioni di secondo ordine (SOS). Questa decomposizione migliora notevolmente la stabilità numerica nelle implementazioni a punto fisso e galleggiante.

4. Attuazione e verifica

Verificare la risposta di frequenza utilizzando uno strumento di analisi spettrale, e garantire che la banda passante si increspa, attenuazione della banda ferma e larghezza di banda di transizione soddisfano le specifiche originali.

Confronto di Butterworth, Chebyshev e Filtri Ellittici

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

Strumenti e software per la progettazione di filtri IIR

I flussi di lavoro di elaborazione digitale del segnale moderni si basano su strumenti di alto livello che automatizzano la generazione di coefficienti di filtro.

  • Python with SciPy:[] Funzioni come [], , [, e restituire direttamente le matrici SOS date frequenze normalizzate. Esempio: .
  • MATLAB / Octave:[] Comandamenti come , , , e calcolo numeratore e denominatore polinomi. La funzione converte in SOS.
  • Pacchetti di progettazione del filtro specifici:[[] strumenti gratuiti e commerciali come FilterPro (Texas Instruments), Nuhertz Filter Solutions, o AADE Filter Design offrono specifiche grafiche interattive e output coefficiente effettivo.

Per un confronto esaustivo di questi ambienti di progettazione, fare riferimento a a 2021 indagine sugli strumenti di progettazione digitale dei filtri[] (ResearchGate).

Guida all'applicazione reale

La scelta tra i tre tipi di filtro dipende dai vincoli specifici del sistema:

  • Elaborazione audio ad alta fedeltà:[] Utilizzare filtri Butterworth per preservare la forma del segnale originale con una distorsione di fase minima. Il leggero compromesso in roll-off è compensato dalla qualità del suono trasparente.
  • Selezione dei canali di telecomunicazioni:[[] I filtri ellittici sono preferiti quando il rigetto adiacente-canale rigido deve essere raggiunto con risorse hardware limitate.
  • I filtri di Chebyshev Type II offrono una banda passante piana per evitare distorcere il complesso QRS, rifiutando comunque efficacemente le interferenze di 50/60 Hz. Molti front-end ECG utilizzano un notch di tipo II di 4 ° ordine.
  • Sistemi di controllo:[] I filtri Butterworth sono spesso scelti per la loro risposta di fase prevedibile e la mancanza di overshoot, che è fondamentale per i margini di stabilità del sistema.

Considerazioni avanzate

Effetti finiti di Word-Length

Quando si implementano i filtri IIR su processori a punto fisso (ad esempio, 16-bit o 32-bit), la quantizzazione del coefficiente può spostare le posizioni dei poli, potenzialmente causando instabilità o prestazioni degradanti.

Design ad alta velocità, Bandpass e Bandstop

Tutti e tre i tipi di prototipi possono essere trasformati da un passaggio basso ad altre risposte selettive di frequenza tramite trasformazioni spettrali (rimodellando [z]]-variabile con una funzione razionale appropriata). Ad esempio, una trasformazione di passaggio-basso raddoppia l'ordine del filtro e viene implementata prontamente in strumenti software.

Confronto con filtri FIR

I filtri FIR offrono fase lineare e stabilità garantita ma richiedono ordini molto più elevati per soddisfare le prestazioni di roll-off IIR. Nelle applicazioni in cui la fase non è critica, il costo computazionale inferiore dei filtri IIR - soprattutto ellittiche - li fa la scelta preferita.

Conclusioni

I filtri IIR con Butterworth, Chebyshev e caratteristiche ellittiche comporta il bilanciamento degli obiettivi concorrenti di flatness tra passband, la ripidità di transizione e la complessità di implementazione. I filtri Butterworth forniscono il design piÃ1 semplice e la risposta di fase piÃ1 prevedibile, rendendoli una scelta sicura per applicazioni audio generiche e di bassa frequenza.

L'ingegnere moderno può sfruttare sofisticate librerie software per prototipi rapidamente questi filtri, lasciando il tempo per concentrarsi sui trade-off di livello di sistema e la verifica.

Per ulteriori informazioni sulla teoria del design e sulla derivazione di questi filtri, vedere []Il libro on line di Julius O. Smith sui filtri digitali e i testi classici di Oppenheim e Schafer. Inoltre, il rapporto sull'applicazione Texas Instruments “Design of IIR Digital Filters”] fornisce una guida pratica e orientata all'implementazione.