Esercizio - cascata di filtri discreti con decimazione (prova scritta del 17 settembre 2025)
Questa pagina non ha ancora la versione ripasso: qui sotto c'è il testo completo.
In questa pagina 4
Testo (prova scritta del corso Teoria dei Segnali, UniPD, 17 settembre 2025, problema 1). Data la cascata di trasformazioni lineari
con , ms e
- Dare una rappresentazione grafica accurata, nel dominio della frequenza, dell'uscita corrispondente all'ingresso e calcolare il valore .
- Calcolare l'energia di quando l'ingresso è , con C. Che grandezza rappresenta il segnale in questo caso?
- Identificare la trasformazione equivalente alla cascata, scrivendone il nucleo, e determinare se è: a) tempo-invariante; b) causale; c) stabile.
Nota sul testo. Il tema non definisce (compare solo in ). Qui si risolve con (la lettura più naturale: diventa un semplice guadagno sul proprio reticolo ); alla fine dei punti (1) e (3) si dice cosa cambierebbe con altri valori. Questa scelta va confermata sul testo originale della prova.
Teoria usata: Segnali a tempo discretoUn segnale a tempo discreto è una funzione complessa $s(nT)$ definita sui multipli interi del quanto temporale $T$ (insieme $\mathbb Z(T)$, velocità $F_p=1/T$). Le definizioni sono quelle dei segnali continui con la somma al posto dell'integrale e il quanto $T$ al posto di $dt$: area $\sum T,s(nT)$, energia $\sum T|s(nT)|^2$, convoluzione $\sum T,x(kT)y(nT-kT)$. L'impulso ideale discreto vale $1/T$ nell'origine. Esponenziali e sinusoidi discreti sono periodici solo se $f_0/F_p$ è razionale e hanno frequenza ambigua a meno di multipli di $F_p$. I segnali periodici con periodo $NT$ sono descritti da $N$ valori e si trattano al calcolatore.Segnali a tempo discreto →, 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 →, Campionamento e ricostruzioneIl campionamento $s_c(nT)=s(nT)$ trasforma un segnale continuo in uno discreto e, in frequenza, ripete lo spettro con periodo $F_c=1/T$: $S_c(f)=\sum_kS(f-kF_c)$. Se le repliche si sovrappongono si ha aliasing e il segnale non è recuperabile. Un interpolatore $\mathbb Z(T)\to\mathbb R$ con risposta impulsiva $g$ produce $\tilde s(t)=\sum_nT g(t-nT)s(nT)$ e $\tilde S=G,S_c$. Teorema del campionamento: se $s$ ha banda $B$ e $F_c\ge2B$, l'interpolatore ideale $G=\operatorname{rect}(f/F_c)$ ricostruisce esattamente $s(t)=\sum s(nT)\operatorname{sinc}(F_c(t-nT))$. Se le ipotesi non valgono c'è un errore, in banda e fuori banda, riducibile con un prefiltro anti-aliasing.Campionamento e ricostruzione →, Trasformazioni e sistemi - proprietàUna trasformazione manda un segnale d'ingresso su I in un segnale d'uscita su U, y = φ[x]. Si classifica per invertibilità, memoria, causalità, stabilità BIBO, linearità, tempo-invarianza e realtà. Una trasformazione lineare ha un nucleo h(t,u) = risposta all'impulso applicato in u; la finestra y = w·x è lineare ma NON tempo-invariante. Le tf LTI (filtri) sono il caso con h(t,u) = g(t-u).Trasformazioni e sistemi - proprietà →, 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 →, Risposta in frequenza e filtriLa risposta in frequenza G(f) = F[g] di un filtro LTI dà Y = G·X. Gli esponenziali complessi sono autofunzioni (autovalore G(f)), quindi un ingresso sinusoidale esce sinusoidale con ampiezza moltiplicata per |G(f0)| e fase aumentata di arg G(f0). Il filtro è reale se G è hermitiana, invertibile se G non si annulla. I filtri ideali sono rect in frequenza; non distorsione secondo Heaviside: |G| costante e fase lineare.Risposta in frequenza e filtri →, ConvoluzioneLa convoluzione $xy(t)=\int x(u),y(t-u),du$ combina due segnali ribaltando e traslando il secondo, moltiplicando e integrando. È commutativa, associativa, lineare; l'area del risultato è il prodotto delle aree; l'estensione è la somma delle estensioni (estremo con estremo); l'impulso $\delta$ è l'elemento neutro; la convoluzione con il gradino integra. Per due segnali periodici di uguale periodo si usa la convoluzione ciclica (integrale su un periodo). È l'operazione del filtraggio: l'uscita di un filtro è la convoluzione dell'ingresso con la risposta impulsiva.Convoluzione →, Energia, potenza e valor medio dei segnaliSu un segnale continuo si calcolano quattro numeri riassuntivi: l'area $\int s,dt$, il valor medio (componente continua) $\lim\frac1{2T}\int_{-T}^Ts,dt$, l'energia $\int|s|^2dt$ e la potenza media $\lim\frac1{2T}\int_{-T}^T|s|^2dt$. Un segnale ad energia finita ha potenza nulla e uno a potenza finita non nulla ha energia infinita. Per un segnale periodico di periodo $T_p$ area ed energia si calcolano su un periodo; un segnale periodico è la ripetizione periodica $\sum_ku(t-kT_p)$ di un suo periodo. Valgono per traslazione l'invarianza, per scala $s(at)$ la divisione per $|a|$.Energia, potenza e valor medio dei segnali →.
Dati numerici: ms. Le frequenze di ripetizione sono Hz (reticolo ) e Hz (reticolo ).
(1) Spettro di per
Passo 1: l'ingresso. Con la convenzione del corso , la trasformata discreta è per ogni (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 →).
Passo 2: il filtro . Sul reticolo la risposta in frequenza è , cioè la ripetizione periodica, con periodo Hz, della trasformata continua di . Questa vale (dalla coppia con ), nulla per Hz. Le repliche distano Hz mentre il triangolo occupa solo Hz: non si sovrappongono, e nella banda base Hz È un filtro passa-basso triangolare, che lascia passare esattamente la banda che la decimazione per può ospitare senza aliasing. Poiché , risulta .
Passo 3: il decimatore. Il decimatore da a tiene un campione ogni quattro; in frequenza è la ripetizione periodica con il nuovo periodo Hz, cioè la somma di quattro repliche di traslate di multipli di Hz: Il triangolo largo Hz ripetuto ogni Hz dà un'onda triangolare che va da (in ) a (in ): le repliche si toccano nei punti dove valgono , senza aliasing (verifica in Hz: , che è il valore del solo triangolo centrale).
Grafico interattivo: Spettro dopo il filtro g1 (tratteggiato, un solo triangolo largo 25 Hz in banda base ±50 Hz) e dopo il decimatore (continuo): onda triangolare di periodo 25 Hz, da 1 a 0, senza sovrapposizioni
Passo 4: il filtro . Con : sul reticolo i campioni sono , e per ogni intero e per . Dunque e zero altrove, che è (l'impulso su vale in ): un semplice guadagno , per ogni . (Equivalente: la ripetizione di con periodo ricopre l'asse con valore costante .) Quindi periodica di periodo Hz, con massimo in Hz e zeri in Hz.
Grafico interattivo: Spettro dell'uscita z per x = δ: Z(f) = ½·Yc(f), onda triangolare di periodo 25 Hz con massimo 1/2 in f = 0, ±25, ±50 Hz e zeri in ±12,5, ±37,5 Hz
Passo 5: il valore . Per l'antitrasformata discreta si integra su un periodo (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 →): (L'area di un triangolo di altezza e base è .) Controllo nel tempo: (perché ), e . ✓ (confermato numericamente).
Se fosse diverso. Per risulta in banda base e : per si ottiene . Per (, guadagno ) . Questi valori sono stati verificati numericamente.
(2) Energia di e significato fisico
Passo 1: il segnale . Per il filtro dà (l'impulso ideale vale in e il peso lo compensa). Il decimatore tiene gli istanti (anche è sul reticolo , perché ):
Passo 2: l'energia. Sul reticolo l'energia è (Energia, potenza e valor medio dei segnaliSu un segnale continuo si calcolano quattro numeri riassuntivi: l'area $\int s,dt$, il valor medio (componente continua) $\lim\frac1{2T}\int_{-T}^Ts,dt$, l'energia $\int|s|^2dt$ e la potenza media $\lim\frac1{2T}\int_{-T}^T|s|^2dt$. Un segnale ad energia finita ha potenza nulla e uno a potenza finita non nulla ha energia infinita. Per un segnale periodico di periodo $T_p$ area ed energia si calcolano su un periodo; un segnale periodico è la ripetizione periodica $\sum_ku(t-kT_p)$ di un suo periodo. Valgono per traslazione l'invarianza, per scala $s(at)$ la divisione per $|a|$.Energia, potenza e valor medio dei segnali →). Ponendo : La somma si calcola così. dà ; per pari non nullo ; per dispari e . Allora usando . Perciò Controllo con Parseval. in banda base, quindi ✓ (e numericamente ).
Che grandezza è ? ha le dimensioni di e è una carica: è una corrente (C/s = A), un impulso di carica C. Il filtro ha dimensione e il peso della convoluzione è un tempo, quindi e hanno ancora le dimensioni di una corrente: è la corrente di uscita del filtro, cioè la risposta impulsiva moltiplicata per la carica dell'impulso. L'energia, in As (e C/s), vale : è l'energia normalizzata (dissipata su una resistenza da , J).
(3) Trasformazione equivalente
Passo 1: i nuclei dei blocchi. In questo corso una trasformazione lineare ha un nucleo con , integrale di Haar (somma pesata dal passo sul reticolo) (Trasformazioni e sistemi - proprietàUna trasformazione manda un segnale d'ingresso su I in un segnale d'uscita su U, y = φ[x]. Si classifica per invertibilità, memoria, causalità, stabilità BIBO, linearità, tempo-invarianza e realtà. Una trasformazione lineare ha un nucleo h(t,u) = risposta all'impulso applicato in u; la finestra y = w·x è lineare ma NON tempo-invariante. Le tf LTI (filtri) sono il caso con h(t,u) = g(t-u).Trasformazioni e sistemi - proprietà →). Qui:
- filtro : , per ;
- decimatore: , per , (con la somma pesata riproduce );
- filtro : , per .
Passo 2: la cascata. Il nucleo della cascata è (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 →). Filtro e decimatore insieme danno per : il decimatore restringe il filtro agli istanti . Allora, con e , perché, per , solo per e vale altrove. Il risultato è cioè : la trasformazione equivalente è un filtro passa-basso triangolare con guadagno seguito dal campionamento ogni . Controllo numerico del nucleo con la somma sulle : , , , uguali a quelli della formula.
a) Tempo-invarianza. La definizione del corso richiede che per ogni del dominio d'ingresso il segnale traslato produca l'uscita traslata della stessa quantità, e ha senso solo se appartiene anche al dominio d'uscita. Qui il dominio d'ingresso è e quello d'uscita : una traslazione di ms non è consentita sull'uscita. Perciò la trasformazione non è tempo-invariante: è periodicamente tempo-invariante, cioè invariante solo per (multipli di ms), come il campionatore (stessa situazione dei sistemi con decimatore).
Controesempio: per l'uscita è , con . Traslando l'ingresso di l'uscita diventa , e in vale : non è un'uscita traslata (e non appartiene nemmeno a ). Invece per , perché dipende solo da .
b) Causalità. Una trasformazione è causale se dipende solo da con , cioè se per . Qui pur essendo : l'uscita in dipende da ingressi futuri. Non causale, perché non si annulla per argomenti negativi (zeri solo in ).
c) Stabilità. Il criterio BIBO (ingresso limitato, uscita limitata) è soddisfatto se . Si ha La somma sul reticolo di passo è uguale all'integrale perché lo spettro di () è molto più stretto della frequenza di ripetizione (formula di Poisson). Dunque e la trasformazione è stabile, con guadagno massimo . (Verificato numericamente: .)
Con altri valori di . Per , ha un salto in : il nucleo decade come , non è sommabile e la trasformazione diventa instabile. Resta non causale e solo periodicamente tempo-invariante.
Errori tipici
- Dimenticare che il decimatore ripete lo spettro con il nuovo periodo (non ) e sommare solo la replica centrale.
- Scambiare il peso dell'integrale: sul reticolo si somma con il passo, per questo e .
- Dichiarare la trasformazione tempo-invariante perché dipende da : la traslazione non è nel dominio d'uscita.
Vedi anche: Esercizio - cascata interpolatore, filtro triangolare e campionatore (prova scritta del 6 febbraio 2026), Esercizio - sinc modulato, estensione spettrale e campionamento a massima efficienza (prova scritta del 3 settembre 2025).