La Struttura della Retroazione
Un filtro FIR (finite impulse response - risposta all'impulso finita) calcola ogni campione di output come una somma ponderata degli input attuali & passati solamente. Nessuna retroazione. La risposta all'impulso ha una durata finita.
Un filtro IIR (infinite impulse response - risposta all'impulso infinita) reimmette gli output precedenti nel calcolo:
y_n = Σ b_k · x_{n−k} − Σ a_k · y_{n−k}
Il termine ricorsivo −Σ a_k · y_{n−k} crea la retroazione. Un singolo impulso all'input ecoeggia attorno al ciclo di retroazione indefinitamente (decadendo geometricamente se stabile).
Perché Usare la Retroazione?
Un filtro IIR può ottenere una selettività in frequenza accentuata con molti meno coefficienti rispetto a un filtro FIR. Un IIR a 2 poli può approssimare quello che richiede un FIR a 50 coefficienti per la stessa attenuazione della banda di arresto.
Il prezzo: potenziale instabilità. I poli di H(z) determinano la stabilità. Tutti i poli devono trovarsi rigorosamente all'interno del cerchio unitario.
La Storia della Doccia con Retroazione di Hamming
Hamming ha usato una vivida storia personale per illustrare l'instabilità della retroazione.
Stava nella stessa stanza d'hotel ripetutamente perché la familiarità lo aiutava a orientarsi quando era stanco. L'idraulico aveva installato tubi dell'acqua calda di grande diametro nella doccia. Questi creavano un significativo ritardo tra l'aggiustamento del rubinetto e la sensazione del cambio della temperatura dell'acqua.
Ogni mattina, Hamming seguiva lo stesso schema: acqua troppo fredda → aumenta l'acqua calda → ancora fredda → aumenta ancora di più → improvvisamente bollente → salta fuori → abbassa → ripeti.
Il ritardo nel percorso della retroazione significava che le sue correzioni oltrepassavano sempre il bersaglio. Non poteva adattarsi al ritardo, neanche dopo molte ripetizioni.
La lezione di ingegneria: l'instabilità nasce da un guadagno eccessivo nel percorso della retroazione O da un ritardo eccessivo nel percorso della retroazione. Entrambi si manifestano come lo stesso comportamento di ricerca. In termini di filtro: poli al livello o al di fuori del cerchio unitario producono esattamente questa risposta oscillatoria o divergente.
Caratterizzazione dell'Instabilità
Hamming ha osservato che la stessa instabilità della doccia poteva essere analizzata in due modi:
1. La sua risposta era troppo forte (guadagno eccessivo nell'azione di correzione).
2. Il suo rilevamento era troppo ritardato (troppo frettoloso nell'entrare nella vasca prima che il sistema si stabilizzasse).
Entrambe le descrizioni producono lo stesso risultato matematico: il polo del ciclo di retroazione si è spostato al di fuori del cerchio unitario.
Le Quattro Famiglie Classiche
La teoria dei filtri analogici si è sviluppata attorno a quattro famiglie di design classiche, ognuna rappresentante un diverso compromesso. Queste famiglie si trasformano al tempo discreto tramite la trasformazione bilineare o l'invarianza dell'impulso.
Butterworth (Massimamente Piatto)
Risposta della banda passante: |H(jω)|² = 1 / (1 + (ω/ω_c)^{2N}). Decrescente monotonicamente. Nessuna ondulazione nella banda passante o nella banda di arresto. I poli giacciono su un cerchio di raggio ω_c nel piano s (o cerchio trasformato nel piano z). La banda passante più piatta possibile per un dato ordine N.
Chebyshev Tipo I
Ondulazione uguale nella banda passante, monotono nella banda di arresto. Per un dato ordine N e livello di ondulazione, raggiunge un taglio più accentuato rispetto a Butterworth. I poli giacciono su un'ellisse (nel piano s).
Chebyshev Tipo II
Ondulazione uguale nella banda di arresto, monotono nella banda passante. Immagine speculare del Tipo I nel dominio della frequenza.
Ellittico (Cauer)
Ondulazione uguale in ENTRAMBI banda passante e banda di arresto. Per un dato ordine N e livelli di ondulazione, raggiunge la transizione più accentuata possibile da banda passante a banda di arresto. Utilizza funzioni ellittiche per posizionare i poli & gli zeri in modo ottimale. Hamming: il nome deriva dal fatto che le funzioni ellittiche sono utilizzate nella derivazione.
Il Compromesso Fondamentale
Tutte e quattro le famiglie raggiungono lo stesso compromesso di base diversamente: un ordine N più alto dà una transizione più accentuata. Consentendo ondulazione (Chebyshev, ellittico) raggiunge una transizione più accentuata per lo stesso N. L'ellittico raggiunge la transizione assolutamente più accentuata per qualsiasi N e specifiche di ondulazione dati.
Scelta Tra le Famiglie di Filtri
La scelta tra le famiglie dipende da ciò che l'applicazione tollera.
Messa in Discussione della Rivendicazione dell'Esperto
Hamming ricordava che certi esperti avevano affermato che tutti i filtri IIR (ricorsivi) possedevano una particolare proprietà. Si è chiesto se questo era davvero vero — e ha trovato un contro-esempio.
Il suo punto: gli esperti spesso portano avanti affermazioni che hanno assorbito a scuola senza mai riesaminarle nel contesto dei problemi attuali. Se ti chiedi se quello che ti viene detto è davvero vero, è sorprendente quanto puoi trovare che è, o sfiora di essere, falso, anche in un campo ben sviluppato.
Il contro-esempio non era il tipo di filtro che normalmente progetteresti, ma ha provato che l'affermazione era superficiale. Un singolo contro-esempio è sufficiente per confutare un'affermazione universale.
Design IIR in Pratica
Hamming ha notato che aveva sviluppato indipendentemente gran parte della teoria dei filtri IIR mentre risolveva un problema diverso: derivare formule di correzione stabili per equazioni differenziali ordinarie numeriche.
The corrector formula form: y_n = Σ a_k · y_{n−k} + Σ b_k · f(y_{n−k})
La retroazione appare sia nei termini y (retroazione lineare) che nei termini f(y) (retroazione nonlineare attraverso l'equazione differenziale). La stabilità per i filtri IIR è un caso speciale del problema più generale della stabilità per i risolutori ODE numerici.
Connessione della Retroazione Tra i Domini
La stessa struttura matematica — retroazione, poli, confine di stabilità — appare nei filtri digitali, nei risolutori ODE numerici, nei sistemi di controllo, nei ritmi biologici, & nei modelli economici.
In ogni dominio: un ciclo di retroazione calcola un nuovo stato dai stati precedenti. La stabilità richiede che la retroazione non amplifichi indefinitamente le perturbazioni.
Il confine di stabilità del cerchio unitario nel piano Z corrisponde a: l'asse immaginario nel piano s di Laplace (tempo continuo), la condizione del raggio spettrale ρ(A) < 1 per le iterazioni lineari, & la condizione dell'esponente di Lyapunov λ < 0 per i sistemi nonlineari.