Salta al contenuto
Note per Studenti Esercizio 2 · proprietà, risposta impulsiva e risposta in frequenza di un sistema discreto (tema d'esame giugno 2023)

Esercizio 2proprietà, risposta impulsiva e risposta in frequenza di un sistema discreto (tema d'esame giugno 2023)

Esame
In questa pagina 5

Testo (tema d'esame giugno 2023, esercizio 2, 8 punti). Si consideri il sistema a tempo discreto SS con relazione ingresso-uscita y(n)=x(n)−12x(n−2)−12x(n+2),n∈Z.y(n)=x(n)-\frac12x(n-2)-\frac12x(n+2),\qquad n\in\mathbb Z. (1) Completare la tabella: statico/dinamico, causale/anticausale/non causale, lineare, tempo-invariante, stabile (BIBO). (2) Trovare la risposta impulsiva. (3) Trovare la risposta in frequenza e tracciarne il modulo. (4) Determinare l'uscita corrispondente all'ingresso x(n)=a+bsin⁡(πn+c)x(n)=a+b\sin(\pi n+c) con a,b,ca,b,c reali positivi.


Teoria usata: Sistemi e loro proprietàUn sistema trasforma un segnale di ingresso in uno di uscita, $y=\Sigma[x]$. Le proprietà da saper verificare sono: memoria (statico/dinamico), causalità, linearità (additività + omogeneità), tempo-invarianza, stabilità BIBO, invertibilità, realtà. Per negare una proprietà basta un controesempio; per affermarla serve una dimostrazione con segnali generici.Sistemi e loro proprietà →, Causalità, stabilità e connessioni dei sistemi LTIPer un sistema LTI le proprietà si leggono sulla risposta impulsiva $h$: causale se $h$ è nulla per $t<0$; reale se $h$ è reale; statico se $h=K\delta$; BIBO stabile se e solo se $h$ è assolutamente integrabile (sommabile nel discreto). Serie e parallelo corrispondono a $h_1*h_2$ e $h_1+h_2$.Causalità, stabilità e connessioni dei sistemi LTI →, 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 →, Segnali periodici, sinusoidi ed esponenziali immaginari puriUn segnale continuo è periodico di periodo T se x(t+T)=x(t); $e^{j2\pi f_0t}$ ha periodo minimo $1/|f_0|$ e una somma di periodici è periodica solo se il rapporto dei periodi è razionale. A tempo discreto $e^{j2\pi\nu n}$ è periodico solo se $\nu$ è razionale, e le pulsazioni che differiscono di $2\pi$ sono lo stesso segnale: ogni sinusoide si riporta alla forma canonica con $\omega\in[0,\pi]$.Segnali periodici, sinusoidi ed esponenziali immaginari puri →.

(1) Le proprietà

Proprietà Esito Motivo
statico falso y(n)y(n) dipende da x(n±2)x(n\pm2), non solo da x(n)x(n)
dinamico vero idem
causale falso usa x(n+2)x(n+2), cioè il futuro
anticausale falso usa anche x(n−2)x(n-2), il passato (e x(n)x(n))
non causale vero dipende sia dal passato sia dal futuro
lineare vero è una combinazione lineare di campioni dell'ingresso, con coefficienti costanti
tempo-invariante vero i coefficienti non dipendono da nn: ritardare xx di n0n_0 ritarda yy di n0n_0
stabile (BIBO) vero ∣y(n)∣≤∣x(n)∣+12∣x(n−2)∣+12∣x(n+2)∣≤2Lx\lvert y(n)\rvert\le\lvert x(n)\rvert+\frac12\lvert x(n-2)\rvert+\frac12\lvert x(n+2)\rvert\le2L_x se ∣x∣≤Lx\lvert x\rvert\le L_x

Per la tempo-invarianza (test del confronto): ritardo prima, sistema poi, dà x(n−n0)−12x(n−n0−2)−12x(n−n0+2)x(n-n_0)-\frac12x(n-n_0-2)-\frac12x(n-n_0+2); sistema prima, ritardo poi, dà y(n−n0)y(n-n_0) con nn sostituito ovunque con n−n0n-n_0: stessa espressione ✓.

(2) Risposta impulsiva

Con x=δx=\delta si ha y(n)=δ(n)−12δ(n−2)−12δ(n+2)y(n)=\delta(n)-\frac12\delta(n-2)-\frac12\delta(n+2): h(n)=−12 δ(n+2)+δ(n)−12 δ(n−2).h(n)=-\frac12\,\delta(n+2)+\delta(n)-\frac12\,\delta(n-2). Tre campioni non nulli: un filtro FIR (sempre stabile, ∑∣h∣=2\sum|h|=2).

(3) Risposta in frequenza

H(ω)=∑nh(n)e−jωn=−12ej2ω+1−12e−j2ωH(\omega)=\sum_nh(n)e^{-j\omega n}=-\frac12e^{j2\omega}+1-\frac12e^{-j2\omega}. Con la formula di Eulero ej2ω+e−j2ω2=cos⁡2ω\frac{e^{j2\omega}+e^{-j2\omega}}2=\cos2\omega: H(ω)=1−cos⁡(2ω).H(\omega)=1-\cos(2\omega). È reale, pari, periodica di periodo π\pi (e quindi anche 2π2\pi) e non negativa: ∣H(ω)∣=1−cos⁡2ω=2sin⁡2ω|H(\omega)|=1-\cos2\omega=2\sin^2\omega. Si annulla per 2ω=2πk2\omega=2\pi k, cioè ω=kπ\omega=k\pi (in particolare ω=0\omega=0 e ω=±π\omega=\pm\pi) e vale al massimo 22 per ω=±π2\omega=\pm\frac\pi2. È un filtro passa-banda che sopprime la continua e la frequenza più alta (ω=π\omega=\pi).

Grafico interattivo: |H(ω)| = 1 − cos(2ω) per il sistema y(n) = x(n) − ½x(n−2) − ½x(n+2): nullo in ω = 0 e ω = ±π, massimo 2 in ω = ±π/2

(4) Uscita per x(n)=a+bsin⁡(πn+c)x(n)=a+b\sin(\pi n+c)

Si identificano le componenti in frequenza di xx. La costante aa ha pulsazione 00. Per la sinusoide, sin⁡(πn+c)=sin⁡(πn)cos⁡c+cos⁡(πn)sin⁡c=(−1)nsin⁡c\sin(\pi n+c)=\sin(\pi n)\cos c+\cos(\pi n)\sin c=(-1)^n\sin c (perché sin⁡πn=0\sin\pi n=0 per nn intero): è la sinusoide più veloce possibile, a pulsazione π\pi. Quindi x(n)=a+bsin⁡c (−1)n,x(n)=a+b\sin c\,(-1)^n, con componenti a ω=0\omega=0 e ω=π\omega=\pi. Poiché H(0)=1−cos⁡0=0H(0)=1-\cos0=0 e H(±π)=1−cos⁡2π=0H(\pm\pi)=1-\cos2\pi=0, entrambe le componenti sono soppresse: y(n)=0∀n.y(n)=0\quad\forall n. Verifica diretta: y(n)=x(n)−12x(n−2)−12x(n+2)y(n)=x(n)-\frac12x(n-2)-\frac12x(n+2) e xx ha periodo 22 (x(n±2)=x(n)x(n\pm2)=x(n)), quindi y(n)=x(n)[1−12−12]=0y(n)=x(n)\left[1-\frac12-\frac12\right]=0 ✓. (Controllato numericamente con a=2a=2, b=3b=3, c=0,7c=0{,}7: ∣y∣<10−13|y|<10^{-13}.)

Errori comuni

  • Dire "non causale e quindi anticausale": anticausale richiede dipendenza solo da presente e futuro.
  • Dire "dinamico perché c'è una somma": è dinamico perché dipende da altri istanti.
  • Scrivere H(ω)=1−cos⁡ωH(\omega)=1-\cos\omega (dimenticando che gli sfasamenti sono di due campioni, non uno).
  • Non riconoscere che sin⁡(πn+c)\sin(\pi n+c) per nn intero è (−1)nsin⁡c(-1)^n\sin c e dire che il sistema "lascia passare una sinusoide".

Versione ripasso

Testo. y(n)=x(n)−12x(n−2)−12x(n+2)y(n)=x(n)-\frac12x(n-2)-\frac12x(n+2): proprietà, hh, H(ω)H(\omega) e ∣H∣|H|, uscita per x=a+bsin⁡(πn+c)x=a+b\sin(\pi n+c) (giugno 2023).

Teoria collegata