Salta al contenuto
Note per Studenti Esercizio 5 · finestra rettangolare discreta e derivata in frequenza della TFtd (tema d'esame giugno 2024)

Esercizio 5finestra rettangolare discreta e derivata in frequenza della TFtd (tema d'esame giugno 2024)

Esame
In questa pagina 5

Testo (tema d'esame giugno 2024, esercizio 2, 7 punti). Sia u(n)u(n) la sequenza gradino discreta (u(n)=0u(n)=0 per n<0n<0, u(n)=1u(n)=1 per n≥0n\ge0). Siano definiti i segnali x1(n)=u(25−n2),x2(n)=−n⋅x1(n),y(n)=x1∗x2(n).x_1(n)=u(25-n^2),\qquad x_2(n)=-n\cdot x_1(n),\qquad y(n)=x_1*x_2(n). (1) Per ognuno dei tre segnali determinare il supporto e l'appartenenza agli spazi ℓ1(Z)\ell^1(\mathbb Z), ℓ2(Z)\ell^2(\mathbb Z), ℓ∞(Z)\ell^\infty(\mathbb Z). (2) Tracciare x1x_1 e x2x_2. (3) Mostrare che la TFtd della finestra rettangolare wN(n)=rect⁡N(n)w_N(n)=\operatorname{rect}_N(n) (non nulla per ∣n∣≤N|n|\le N) è WN(ω)=sin⁡(2N+12ω)sin⁡ω2W_N(\omega)=\frac{\sin\left(\frac{2N+1}2\omega\right)}{\sin\frac\omega2}. (4) Calcolare X1(ω)X_1(\omega), X2(ω)X_2(\omega), Y(ω)Y(\omega) in termini di WN(ω)W_N(\omega) e WN′(ω)W_N'(\omega), con NN opportuno da determinare.


Teoria usata: Segnali - supporto, area, valor medio, energia e potenzaUn segnale è una funzione del tempo (continuo $t$ o discreto $n$). Si descrive con pochi numeri: estensione, area, valor medio, energia $\int|x|^2$ e potenza (energia media). Energia finita implica potenza nulla; potenza finita non nulla implica energia infinita; per i segnali periodici tutto si calcola su un periodo.Segnali - supporto, area, valor medio, energia e potenza →, Spazi di segnali L1, L2 e L-infinitoGli spazi $L^1$ (area del modulo finita), $L^2$ (energia finita) e $L^\infty$ (segnale limitato) classificano i segnali e dicono quali teoremi si possono usare: $L^1$ per l'esistenza della trasformata di Fourier e la stabilità BIBO, $L^2$ per Parseval. Nel discreto valgono le inclusioni $\ell^1 \subset \ell^2 \subset \ell^\infty$, nel continuo no.Spazi di segnali L1, L2 e L-infinito →, Trasformata di Fourier a tempo discreto (TFtd)La TFtd di una sequenza è $X(\omega)=\sum_nx(n)e^{-j\omega n}$, funzione continua e periodica di periodo $2\pi$; si inverte con $x(n)=\frac1{2\pi}\int_{-\pi}^{\pi}X(\omega)e^{j\omega n}d\omega$. Ha le stesse proprietà della TF continua (la convoluzione diventa prodotto, $n,x(n)\leftrightarrow jX'$). È la risposta in frequenza dei sistemi discreti; con la TFD e lo zero-padding se ne ottengono campioni arbitrariamente fitti.Trasformata di Fourier a tempo discreto (TFtd) →, Calcolo della convoluzione e sue proprietàIl supporto della convoluzione è la somma dei supporti, $\operatorname{rect}*\operatorname{rect}=\Lambda$, e due esponenziali causali danno $(e^{-bt}-e^{-at})/(a-b)$. Si calcola con il metodo grafico a casi (ribaltare, traslare, individuare gli intervalli di sovrapposizione). Proprietà: lineare, commutativa, associativa, $\delta$ è l'elemento neutro, la traslazione si somma, l'area è il prodotto delle aree.Calcolo della convoluzione e sue proprietà →.

(1) Supporti e spazi

x1(n)=u(25−n2)x_1(n)=u(25-n^2) vale 11 se 25−n2≥025-n^2\ge0, cioè n2≤25n^2\le25, cioè ∣n∣≤5|n|\le5, e 00 altrimenti: x1=rect⁡5x_1=\operatorname{rect}_5, 11 campioni uguali a 11 (n=−5,…,5n=-5,\ldots,5), supporto [−5,5][-5,5].

x2(n)=−n x1(n)x_2(n)=-n\,x_1(n) ha lo stesso supporto [−5,5][-5,5] (il fattore −n-n annulla solo n=0n=0). Valori: x2(n)=−nx_2(n)=-n per ∣n∣≤5|n|\le5, cioè 5,4,3,2,1,0,−1,−2,−3,−4,−55,4,3,2,1,0,-1,-2,-3,-4,-5 per n=−5,…,5n=-5,\ldots,5.

y=x1∗x2y=x_1*x_2 ha supporto la somma dei supporti (Calcolo della convoluzione e sue proprietàIl supporto della convoluzione è la somma dei supporti, $\operatorname{rect}*\operatorname{rect}=\Lambda$, e due esponenziali causali danno $(e^{-bt}-e^{-at})/(a-b)$. Si calcola con il metodo grafico a casi (ribaltare, traslare, individuare gli intervalli di sovrapposizione). Proprietà: lineare, commutativa, associativa, $\delta$ è l'elemento neutro, la traslazione si somma, l'area è il prodotto delle aree.Calcolo della convoluzione e sue proprietà →): [−5−5, 5+5]=[−10,10][-5-5,\,5+5]=[-10,10], quindi 21 campioni.

Tutti e tre i segnali hanno supporto finito e valori finiti: quindi appartengono a ℓ1\ell^1, ℓ2\ell^2 e ℓ∞\ell^\infty (somma di un numero finito di termini finiti).

(2) Grafici

x1x_1: 1111 campioni di valore 11 tra n=−5n=-5 e n=5n=5. x2x_2: la rampa discreta −n-n nello stesso intervallo, dispari, da +5+5 a −5-5.

Grafico interattivo: x₂(n) = −n·x₁(n): rampa discreta decrescente da +5 (in n = −5) a −5 (in n = 5), nulla fuori da |n| ≤ 5

(3) La TFtd di wNw_N

WN(ω)=∑n=−NNe−jωnW_N(\omega)=\sum_{n=-N}^Ne^{-j\omega n}: somma geometrica di 2N+12N+1 termini (vedi 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 →). Si porta l'indice a partire da 00 con m=n+Nm=n+N: WN(ω)=ejωN∑m=02Ne−jωm=ejωN 1−e−jω(2N+1)1−e−jω.W_N(\omega)=e^{j\omega N}\sum_{m=0}^{2N}e^{-j\omega m}=e^{j\omega N}\,\frac{1-e^{-j\omega(2N+1)}}{1-e^{-j\omega}}. Si estrae e−jω(2N+1)/2e^{-j\omega(2N+1)/2} dal numeratore e e−jω/2e^{-j\omega/2} dal denominatore: =ejωN e−jω2N+12(ejω2N+12−e−jω2N+12)e−jω/2(ejω/2−e−jω/2)=ejω[N−2N+12+12]2jsin⁡(2N+12ω)2jsin⁡ω2.=e^{j\omega N}\,\frac{e^{-j\omega\frac{2N+1}2}\left(e^{j\omega\frac{2N+1}2}-e^{-j\omega\frac{2N+1}2}\right)}{e^{-j\omega/2}\left(e^{j\omega/2}-e^{-j\omega/2}\right)}=e^{j\omega\left[N-\frac{2N+1}2+\frac12\right]}\frac{2j\sin\left(\frac{2N+1}2\omega\right)}{2j\sin\frac\omega2}. L'esponente N−N−12+12=0N-N-\frac12+\frac12=0, quindi i fattori di fase spariscono: WN(ω)=sin⁡(2N+12ω)sin⁡ω2.W_N(\omega)=\frac{\sin\left(\frac{2N+1}2\omega\right)}{\sin\frac\omega2}. (Il passaggio vale per ω≠2πk\omega\ne2\pi k; in ω=2πk\omega=2\pi k la somma vale 2N+12N+1, e l'espressione si prolunga per continuità.)

(4) Le trasformate

Controllo numerico in ω=0,8\omega=0{,}8: la somma diretta ∑nx2(n)e−jωn\sum_nx_2(n)e^{-j\omega n} dà 1,45077j1{,}45077j e −jW5′(0,8)=1,45077j-jW_5'(0{,}8)=1{,}45077j; la somma diretta di yy dà −3,54517j-3{,}54517j e −jW5W5′-jW_5W_5' dà lo stesso valore.

Forma chiusa di yy (come verifica). Si riconosce W5W5′=12ddω(W52)W_5W_5'=\frac12\frac{d}{d\omega}\left(W_5^2\right), e W52W_5^2 è la TFtd di x1∗x1x_1*x_1, il triangolo 11−∣n∣11-|n| per ∣n∣≤10|n|\le10. Dalla regola n z(n)↔jZ′n\,z(n)\leftrightarrow jZ' con Z=W52Z=W_5^2: F[n (x1∗x1)]=j⋅2W5W5′\mathcal F[n\,(x_1*x_1)]=j\cdot2W_5W_5', quindi W5W5′=12jF[n(x1∗x1)]W_5W_5'=\frac{1}{2j}\mathcal F[n(x_1*x_1)] e Y=−jW5W5′=−12F[n(x1∗x1)]Y=-jW_5W_5'=-\frac12\mathcal F[n(x_1*x_1)]: y(n)=−n2 (11−∣n∣)per ∣n∣≤10.y(n)=-\frac n2\,(11-|n|)\quad\text{per }|n|\le10. Per esempio y(1)=−5y(1)=-5: direttamente y(1)=∑k=−45x2(1−k)⋅1=−∑k=−45(1−k)=−(10−5)=−5y(1)=\sum_{k=-4}^{5}x_2(1-k)\cdot1=-\sum_{k=-4}^5(1-k)=-(10-5)=-5 ✓; i valori vanno da y(−10)=5y(-10)=5 a y(−5)=15y(-5)=15 (massimo) e y(5)=−15y(5)=-15 (minimo): dispari. (Verificato con np.convolve.)

Grafico interattivo: y(n) = x₁*x₂(n) = −(n/2)(11−|n|) per |n| ≤ 10: dispari, con massimo 15 in n = −5 e minimo −15 in n = 5

Errori comuni

  • Dimenticare che u(25−n2)u(25-n^2) è una finestra rect⁡5\operatorname{rect}_5 (11 campioni) e non una semplice "metà retta".
  • Sbagliare il numero di campioni della convoluzione: 11+11−1=2111+11-1=21, supporto [−10,10][-10,10].
  • Dimenticare il segno meno: x2=−n x1x_2=-n\,x_1 dà X2=−jW′X_2=-jW' e non +jW′+jW'.
  • Dire che yy appartiene a ℓp\ell^p solo per alcuni pp: ha supporto finito, quindi a tutti.

Versione ripasso

Testo. x1=u(25−n2)x_1=u(25-n^2), x2=−n x1x_2=-n\,x_1, y=x1∗x2y=x_1*x_2: supporti e ℓp\ell^p; grafici; TFtd di wNw_N; X1X_1, X2X_2, YY (giugno 2024).

Teoria collegata