Risposta in frequenza e filtri
In questa pagina 8
Qui si usano i filtri LTI della nota Sistemi lineari tempo-invarianti e risposta impulsivaUna tf lineare e tempo-invariante (LTI, filtro) ha nucleo h(t,u) = g(t-u): l'uscita è la convoluzione y = gx con la risposta impulsiva g (uscita all'impulso ideale nell'origine). Causale se e solo se g è causale; stabile BIBO se e solo se g è assolutamente integrabile (sommabile); reale se e solo se g è reale. Cascata: g = g2g1; parallelo: g1+g2; retroazione: Ge = G/(1+HG) in frequenza.Sistemi lineari tempo-invarianti e risposta impulsiva → e la trasformata di Fourier (Trasformata di FourierLa trasformata di Fourier $S(f)=\int s(t)e^{-i2\pi ft}dt$ associa a un segnale continuo (anche aperiodico) la sua rappresentazione in frequenza; l'antitrasformata $s(t)=\int S(f)e^{i2\pi ft}df$ lo ricostruisce, perché gli esponenziali $e^{i2\pi ft}$ sono ortogonali su tutto $\mathbb R$ ($\int e^{i2\pi ft}dt=\delta(f)$). Per un segnale reale $S(-f)=S^(f)$. Si calcola per i segnali notevoli (rect $\leftrightarrow$ sinc, $e^{-\alpha t}\mathbf 1(t)\leftrightarrow\frac1{\alpha+i2\pi f}$, gaussiana, $\delta\leftrightarrow1$, $1\leftrightarrow\delta$, gradino) e per i segnali periodici, la cui trasformata è un treno di impulsi di area $S_n$ in $nF$.Trasformata di Fourier →, Proprietà della trasformata di FourierLe regole della trasformata di Fourier trasformano operazioni sui segnali in operazioni sulle trasformate: linearità, ribaltamento, coniugio, traslazione nel tempo ($\times e^{-i2\pi ft_0}$) e in frequenza, convoluzione $\leftrightarrow$ prodotto, cambio di scala $s(at)\to\frac1{|a|}S(f/a)$, derivazione ($\times i2\pi f$), integrazione, regola di simmetria ($S(t)\to s(-f)$). Area $S(0)=\int s$, teorema di Parseval $\int|s|^2=\int|S|^2$. Durata e banda sono inversamente legati e un segnale non può avere durata e banda entrambe limitate; la banda del prodotto è la somma delle bande. Con queste regole si ricavano quasi tutte le trasformate senza integrare.Proprietà della trasformata di Fourier →, Trasformata di Fourier a tempo discretoLa trasformata di Fourier di un segnale discreto è $S(f)=\sum_nT,s(nT),e^{-i2\pi fnT}$: una funzione continua e periodica in $f$ di periodo $F_p=1/T$. L'antitrasformata è l'integrale su un periodo, $s(nT)=\int_0^{F_p}S(f)e^{i2\pi fnT}df$. Le regole sono quelle del caso continuo (traslazione, convoluzione $\leftrightarrow$ prodotto, Parseval $\sum T|s|^2=\int_0^{F_p}|S|^2$), con incremento e somma corrente al posto di derivata e integrale. Per segnali reali $S(f)=S^(-f)$, quindi basta $[0,F_p/2]$. Esempi: $\delta\to1$, $1\to\delta_{F_p}$, rect $\to$ sinc periodico, $a^n\mathbf 1_0\to\frac T{1-ae^{-i2\pi fT}}$.Trasformata di Fourier a tempo discreto →). Le tf LTI si chiamano anche filtri perché, come si vede, discriminano le frequenze.
Risposta in frequenza
L'uscita di un filtro è . La convoluzione diventa un prodotto in frequenza:
Formula (risposta in frequenza). è la risposta in frequenza; il suo modulo è la caratteristica d'ampiezza, il suo argomento la caratteristica di fase. Per un filtro continuo è definita su ; per un filtro discreto su è periodica di periodo .
Esempio (FIR discreto). Il filtro ha, trasformando i due membri e ricordando la regola del ritardo ( per ogni di ritardo), , quindi : periodica di periodo , nulla in .
Le trasformate di Laplace e zeta danno relazioni analoghe e con le funzioni di trasferimento (Trasformata zeta e sistemi a tempo discretoUn sistema LTI discreto è descritto da un'equazione alle differenze Σ a_i y(n-i) = Σ b_i x(n-i) con condizioni iniziali. L'uscita è evoluzione libera (dalle condizioni iniziali) più risposta forzata gx. Si calcola in tre modi: soluzione dell'equazione, risposta impulsiva con segnale fittizio, trasformata zeta. FIR: memoria finita, sempre stabile; IIR: stabile se tutti i poli hanno modulo minore di 1. In frequenza G(f) = Σ b_i e^{-i2πfi} / Σ a_i e^{-i2πfi}.Trasformata zeta e sistemi a tempo discreto →, Sistemi continui descritti da equazioni differenzialiUn'equazione differenziale lineare a coefficienti costanti Σ a_i y^(i) = Σ b_i x^(i) descrive un sistema LTI continuo con stato (condizioni iniziali y(0), y'(0), ...). L'uscita è evoluzione libera + risposta forzata gx. Si calcola risolvendo l'equazione (omogenea + particolare) oppure con Laplace: G(p) = Σ b_i p^i / Σ a_i p^i. Il sistema è stabile se tutte le radici del polinomio caratteristico hanno parte reale negativa.Sistemi continui descritti da equazioni differenziali →).
Autofunzioni
Teorema (autofunzioni dei filtri). Se all'ingresso di un filtro LTI si applica l'esponenziale complesso (su o su ), l'uscita è lo stesso esponenziale, moltiplicato per il numero complesso (l'autovalore). Più in generale dà e, su , dà (per e nella regione di convergenza).
Perché: : il segnale "esce" dall'integrale come fattore costante. Un filtro non cambia la dipendenza dal tempo dell'esponenziale, ne cambia solo ampiezza e fase: questo è il motivo per cui Fourier è lo strumento giusto per i sistemi LTI.
Esempio. Il filtro RC ha (trasformata di ). L'esponenziale di frequenza esce moltiplicato per : al crescere di il modulo diminuisce, le alte frequenze sono attenuate.
Regime sinusoidale
Sia con . Per Eulero (Numeri complessi, formula di Eulero ed esponenziali complessiUn numero complesso si scrive in forma cartesiana a+jb o polare |x|e^{jφ}; il prodotto moltiplica i moduli e somma le fasi. La formula di Eulero e^{jα}=cos α + j sin α lega esponenziali e sinusoidi e permette di trattare tutti i segnali del corso come somme di esponenziali complessi e^{(σ+jω)t}.Numeri complessi, formula di Eulero ed esponenziali complessi →) è la somma di due esponenziali di ampiezze complesse e frequenze . Per linearità l'uscita è Se il filtro è reale, e i due addendi sono uno il complesso coniugato dell'altro:
Formula (ingresso sinusoidale, filtro reale). L'uscita è una sinusoide della stessa frequenza, con ampiezza moltiplicata per la caratteristica d'ampiezza e fase aumentata della caratteristica di fase, entrambe valutate in .
Esempio numerico. Filtro RC con , ingresso . Si calcola , quindi e . L'uscita è (verificato con una convoluzione numerica a regime). La frequenza di taglio a è : a il modulo vale e la fase ; per la fase tende a .
Grafico interattivo: Filtro RC con T = 1 ms: caratteristica d'ampiezza |G(f)| (f in Hz); a f_c = 159 Hz vale 1/√2
Grafico interattivo: Filtro RC con T = 1 ms: caratteristica di fase β(f) = -arctan(2πfT) in radianti; a f_c vale -π/4, per f → ∞ tende a -π/2
Grafico interattivo: Diagramma di Bode d'ampiezza del filtro RC con T = 1 ms: |G| in dB con asse delle frequenze logaritmico (-3 dB a f_c, poi -20 dB per decade)
Proprietà lette sulla risposta in frequenza
Filtro reale. Il filtro è reale se e solo se è reale, cioè se e solo se ha simmetria hermitiana, : il modulo è pari e la fase dispari (Proprietà della trasformata di FourierLe regole della trasformata di Fourier trasformano operazioni sui segnali in operazioni sulle trasformate: linearità, ribaltamento, coniugio, traslazione nel tempo ($\times e^{-i2\pi ft_0}$) e in frequenza, convoluzione $\leftrightarrow$ prodotto, cambio di scala $s(at)\to\frac1{|a|}S(f/a)$, derivazione ($\times i2\pi f$), integrazione, regola di simmetria ($S(t)\to s(-f)$). Area $S(0)=\int s$, teorema di Parseval $\int|s|^2=\int|S|^2$. Durata e banda sono inversamente legati e un segnale non può avere durata e banda entrambe limitate; la banda del prodotto è la somma delle bande. Con queste regole si ricavano quasi tutte le trasformate senza integrare.Proprietà della trasformata di Fourier →).
Esempio. Il filtro discreto , , ha . È reale solo per reale. Per reale e si ha ; con in vale e la fase è (verificato).
Invertibilità. Un filtro LTI è invertibile se e solo se per ogni . Allora da si ha e l'inverso è il filtro LTI con Esempio. Il filtro RC ( sempre) ha inverso , cioè (la moltiplicazione per è una derivata). Il filtro discreto ha inverso , cioè (la trasformata di è ), un FIR di due termini (per e la convoluzione è proprio l'impulso: , verificato). Un filtro che si annulla in una frequenza non è invertibile: ad esempio il FIR si annulla in .
Inversa di una finestra. La finestra è invertibile se e solo se per ogni e l'inversa è la finestra .
Filtri ideali
Definizione (filtro ideale). Un filtro ideale ha nella banda passante e nella banda oscura . La risposta impulsiva è . È reale se è simmetrica rispetto a ; non è mai causale.
| filtro (continuo) | ||
|---|---|---|
| passa-basso di banda | ||
| passa-alto (taglio ) | ||
| passa-banda (centro , banda ) | ||
| passa-banda complesso | (complessa: non reale) |
Il passa-basso ideale discreto su ha risposta in frequenza periodica con , (la banda è una frazione del periodo ), e risposta impulsiva Per si ottiene : il filtro passa-tutto (l'identità). Per i valori sono verificati numericamente con l'integrale.
Grafico interattivo: Passa-basso ideale di banda B = 1: risposta in frequenza rect(f/2B) e risposta impulsiva g(t) = 2B sinc(2Bt) (non causale: oscilla anche per t < 0)
Passa-alto e passa-banda ideali discreti. Si procede come per il passa-basso: la risposta in frequenza è periodica di periodo , quindi si assegna in e la si ripete; la risposta impulsiva è l'antitrasformata su quel solo periodo.
- Passa-alto di taglio (): , cioè "tutto tranne il passa-basso". Poiché l'identità ha risposta impulsiva ,
- Passa-banda reale di centro e larghezza (con e ): e sono i campioni della risposta impulsiva del filtro continuo con la stessa in , perché l'integrale è lo stesso e la banda sta tutta in un periodo (per il passa-basso ritorna con ).
Esempio (, quindi ). Passa-basso con (banda ): vale in , in , in , in , in . Passa-alto con taglio (la banda complementare): vale in , in , in , in , in : è proprio volte quella del passa-basso, perché moltiplicare per trasla lo spettro di e il passa-basso diventa il passa-alto (verificato integrando su un periodo). Passa-banda con , : vale in , in , in , in , in (coseno che si annulla sui dispari, inviluppo sinc).
Con una sinusoide la frequenza si confronta con la banda modulo : con fuori da si ripiega prima ( equivale a ). Con il passa-alto di taglio azzera e (entrambe sono a , in banda oscura) e lascia passare e (cioè e , in banda passante); il passa-basso fa l'opposto. Verificato filtrando in frequenza campioni: la potenza d'uscita è o (la potenza d'ingresso) secondo i casi.
Grafico interattivo: Risposta impulsiva del passa-alto ideale discreto con taglio F_p/4 (T = 1): g(n) = (−1)^n · ½ sinc(n/2); vale 0,5 in n = 0 e si annulla per n = ±2, ±4, …
Cosa fanno. Applicato a , il passa-basso ideale lo lascia inalterato se e lo azzera se . In generale vale in banda passante e fuori: dove è l'estensione spettrale di (l'insieme dove ). Ne seguono due casi:
- Irrilevanza del filtro. Se , allora : il filtro è trasparente (irrilevante) per quel segnale.
- Se , .
In tutti gli altri casi il segnale è distorto. Con ( su tutto ) si ottiene : l'impulso diventa un sinc.
Esempio. ha . Applicato al passa-basso di banda : se si ha ; se si ha , cioè .
Per i filtri discreti vale lo stesso ragionamento sull'intersezione dell'estensione di con .
Larghezza di banda e banda. La larghezza di banda di un filtro è la misura dell'estensione di ; per un filtro reale si conta soltanto la parte con e si parla di banda . Così il passa-basso ideale ha banda (larghezza ), il passa-banda ideale reale ha banda e il passa-banda complesso larghezza di banda .
Interpretazione della trasformata di Fourier con un filtro passa-banda
Si applichi a un passa-banda complesso stretto, . L'uscita (tanto meglio quanto più è piccola) è la componente di nell'intorno di . In energia, per Parseval, : quindi è la densità di energia per unità di frequenza.
Condizioni di equalizzazione e non distorsione
In un sistema di telecomunicazioni si vorrebbe che il segnale finale fosse uguale al segnale originario , condizione di perfetta equalizzazione . Questo richiede l'identità, cioè il filtro passa-tutto , . Se la condizione si chiede solo per i segnali a banda limitata con basta qualunque filtro con per ed arbitrario fuori, e il più semplice (e in pratica il migliore, perché elimina i disturbi fuori banda) è il passa-basso ideale .
Si accetta spesso che il segnale arrivi amplificato/attenuato e ritardato:
Condizione (non distorsione secondo Heaviside). per ogni . Equivale a un filtro con cioè modulo costante e fase lineare (ritardo uguale per tutte le frequenze).
Una fase non lineare ritarda le varie componenti in modo diverso e distorce la forma d'onda anche se il modulo è costante. In pratica le condizioni si impongono solo nella banda utile (un amplificatore audio ad alta fedeltà: -). Per i segnali a banda limitata si richiedono e fase lineare solo in . Le stesse idee sono riprese in Interpolazione LTI e teorema del campionamentoIl campionatore R→Z(T) ripete lo spettro con periodo Fc = 1/T (Sc(f) = Σ S(f - kFc), senza fattore 1/T). Il filtro interpolatore Z(T)→R ha y(t) = Σ x(nT) g0(t-nT) con g0 = T g e in frequenza Y = G·X. Se S è nulla fuori da (-B,B) e Fc ≥ 2B, con g0(t) = sinc(Fc t) si ricostruisce esattamente s(t) dai campioni. Altrimenti c'è un errore (in banda per l'aliasing, fuori banda per la parte tagliata), ridotto da un prefiltro anti-aliasing.Interpolazione LTI e teorema del campionamento →.
Errori comuni
- Applicare la formula dell'ingresso sinusoidale a un filtro non reale: serve .
- Confondere ampiezza della sinusoide in uscita (moltiplicata per ) con la sua potenza ().
- Dimenticare che del filtro discreto è periodica: la banda si limita a un periodo.
- Dire che un filtro ideale è realizzabile: non è causale.
Vedi anche: Risposta in frequenza e risposta a ingresso sinusoidaleGli esponenziali complessi $e^{j\omega t}$ sono autofunzioni dei sistemi LTI stabili: escono moltiplicati per la risposta in frequenza $H(\omega)=\int h(\tau)e^{-j\omega\tau}d\tau$. Per un sistema reale, $A\cos(\omega_0t+\varphi)$ esce come $A|H(\omega_0)|\cos(\omega_0t+\varphi+\arg H(\omega_0))$: la frequenza non cambia, cambiano ampiezza e fase.Risposta in frequenza e risposta a ingresso sinusoidale → e Filtri, modulazione e trasmissione AMIn frequenza un filtro LTI moltiplica: $Y(\omega)=H(\omega)X(\omega)$ (serie, parallelo, retroazione diventano $H_1H_2$, $H_1+H_2$, $\frac{H_1}{1+H_1H_2}$). I filtri ideali (passa-basso, passa-alto, passa-banda) hanno $h$ non causale e instabile; il filtro RC è il passa-basso reale. La modulazione $x(t)\cos\omega_0t$ trasla lo spettro in $\pm\omega_0$ e si demodula con la stessa portante e un passa-basso di guadagno 2.Filtri, modulazione e trasmissione AM → nella materia gemella; Circuiti del primo ordine (RC e RL)Una rete con un solo condensatore o un solo induttore è del primo ordine: la variabile di stato x(t) soddisfa dx/dt + x/T = b(t). Con generatori costanti la soluzione è x(t) = x∞ + [x(0) - x∞] e^(-t/T): valore finale (regime stazionario per t > 0) più un transitorio esponenziale. La costante di tempo è T = R C per il condensatore e T = L / R per l'induttore, dove R è la resistenza di Thévenin vista dal componente dinamico. Dopo 4-5 costanti di tempo il transitorio è esaurito. Carica RC: vC = E(1 - e^(-t/RC)); scarica: vC = V0 e^(-t/RC); duali per RL.Circuiti del primo ordine (RC e RL) → per il filtro RC come circuito.
Versione ripasso
- Formula: con risposta in frequenza (la convoluzione diventa prodotto). è la caratteristica d'ampiezza, la caratteristica di fase. Nel discreto è periodica di periodo .
- Esempio FIR: dà , nulla in .
- Autofunzioni: l'ingresso dà l'uscita : lo stesso esponenziale, con autovalore . Su e per Laplace e zeta valgono relazioni analoghe (Trasformata zeta e sistemi a tempo discretoUn sistema LTI discreto è descritto da un'equazione alle differenze Σ a_i y(n-i) = Σ b_i x(n-i) con condizioni iniziali. L'uscita è evoluzione libera (dalle condizioni iniziali) più risposta forzata gx. Si calcola in tre modi: soluzione dell'equazione, risposta impulsiva con segnale fittizio, trasformata zeta. FIR: memoria finita, sempre stabile; IIR: stabile se tutti i poli hanno modulo minore di 1. In frequenza G(f) = Σ b_i e^{-i2πfi} / Σ a_i e^{-i2πfi}.Trasformata zeta e sistemi a tempo discreto →, Sistemi continui descritti da equazioni differenzialiUn'equazione differenziale lineare a coefficienti costanti Σ a_i y^(i) = Σ b_i x^(i) descrive un sistema LTI continuo con stato (condizioni iniziali y(0), y'(0), ...). L'uscita è evoluzione libera + risposta forzata gx. Si calcola risolvendo l'equazione (omogenea + particolare) oppure con Laplace: G(p) = Σ b_i p^i / Σ a_i p^i. Il sistema è stabile se tutte le radici del polinomio caratteristico hanno parte reale negativa.Sistemi continui descritti da equazioni differenziali →).
- Sinusoide, filtro reale: . Il filtro è reale se .
- Esempio RC con ms, ingresso : , , uscita . Il filtro RC ha e frequenza di taglio Hz, dove il modulo vale e la fase (Circuiti del primo ordine (RC e RL)Una rete con un solo condensatore o un solo induttore è del primo ordine: la variabile di stato x(t) soddisfa dx/dt + x/T = b(t). Con generatori costanti la soluzione è x(t) = x∞ + [x(0) - x∞] e^(-t/T): valore finale (regime stazionario per t > 0) più un transitorio esponenziale. La costante di tempo è T = R C per il condensatore e T = L / R per l'induttore, dove R è la resistenza di Thévenin vista dal componente dinamico. Dopo 4-5 costanti di tempo il transitorio è esaurito. Carica RC: vC = E(1 - e^(-t/RC)); scarica: vC = V0 e^(-t/RC); duali per RL.Circuiti del primo ordine (RC e RL) →).
- Filtro reale: reale , cioè modulo pari e fase dispari (Proprietà della trasformata di FourierLe regole della trasformata di Fourier trasformano operazioni sui segnali in operazioni sulle trasformate: linearità, ribaltamento, coniugio, traslazione nel tempo ($\times e^{-i2\pi ft_0}$) e in frequenza, convoluzione $\leftrightarrow$ prodotto, cambio di scala $s(at)\to\frac1{|a|}S(f/a)$, derivazione ($\times i2\pi f$), integrazione, regola di simmetria ($S(t)\to s(-f)$). Area $S(0)=\int s$, teorema di Parseval $\int|s|^2=\int|S|^2$. Durata e banda sono inversamente legati e un segnale non può avere durata e banda entrambe limitate; la banda del prodotto è la somma delle bande. Con queste regole si ricavano quasi tutte le trasformate senza integrare.Proprietà della trasformata di Fourier →). Esempio: ha , reale solo per reale.
- Invertibilità: un filtro LTI è invertibile per ogni , e l'inverso ha . Esempio: il filtro RC ha inverso , cioè . Un filtro che si annulla in una frequenza non è invertibile.
- Filtri ideali: nella banda passante , altrove; . Non sono mai causali. Passa-basso di banda : , . Passa-alto di taglio : .
- Passa-basso discreto: con e , ; per è l'identità.
- Passa-alto e passa-banda discreti: si assegna in e si ripete. Passa-alto di taglio : ; passa-banda reale: . Le sinusoidi si confrontano con la banda modulo . Esempio (, taglio ): , con e .
- Azione del passa-basso ideale: lascia inalterata una sinusoide con e la annulla con . Se il filtro è irrilevante per (); se , . Esempio: con diventa .
- Interpretazione con passa-banda stretto: , quindi è la densità di energia per unità di frequenza.
- Non distorsione (Heaviside): equivale a : modulo costante e fase lineare . Per segnali con basta in , e il passa-basso ideale è la scelta migliore.
Errori tipici:
- Applicare la formula della sinusoide in uscita a un filtro non reale: serve .
- Confondere l'ampiezza in uscita con la potenza .
- Dimenticare che discreta è periodica: la banda si conta su un solo periodo.
- Ritenere realizzabile un filtro ideale: non è causale.
- Usare costante come unica condizione di non distorsione: serve anche la fase lineare.
Esercizi su questo argomento
- Esercizio - cascata di filtri discreti con decimazione (prova scritta del 17 settembre 2025)
- Esercizio - cascata interpolatore, filtro triangolare e campionatore (prova scritta del 6 febbraio 2026)
- Esercizio - connessioni in cascata e in parallelo di trasformazioni LTI
- Esercizio - corpo che cade nell'atmosfera - equazione differenziale del primo ordine
- Esercizio - filtraggio di un impulso rettangolare con convoluzione e trasformate in Python (laboratorio 5)
- Esercizio - filtri FIR e IIR e risposta in frequenza in Python (laboratorio 5)
- Esercizio - filtro passa-alto ideale e segnale sincn attraverso un passa-basso ideale
- Esercizio - interpolatori in cascata e interpolazione parabolica a tratti
- Esercizio - inversa di un filtro e verifica di g convoluta con g_i uguale a delta
- Esercizio - massa e molla - equazione differenziale del secondo ordine
- Esercizio - processi gaussiani con correlazione modulata e filtro passa-basso (prova scritta del 17 settembre 2025)
- Esercizio - proprietà della convoluzione
- Esercizio - rumore bianco gaussiano e filtro passa-basso in Python (laboratorio 6)
- Esercizio - sinc modulato, estensione spettrale e campionamento a massima efficienza (prova scritta del 3 settembre 2025)
Teoria collegata
- Risposta in frequenza dei sistemi FIR
- Campionamento e ricostruzione
- Convoluzione
- Processi aleatori attraverso sistemi LTI
- Proprietà della trasformata di Fourier
- Segnali a tempo continuo - definizioni e trasformazioni
- Serie di Fourier
- Sistemi continui descritti da equazioni differenziali
- Sistemi lineari tempo-invarianti e risposta impulsiva
- Trasformazioni e sistemi - proprietà