Salta al contenuto
Note per Studenti Esercizio - campionamento e interpolazione di segnali rect, triangolari ed esponenziali (esercizi del corso, cap. 2, es. 20 e 21)

Esercizio - campionamento e interpolazione di segnali rect, triangolari ed esponenziali (esercizi del corso, cap. 2, es. 20 e 21)

In questa pagina 3

Testo (raccolta di esercizi del corso, capitolo 2, esercizi 20 e 21).

  • Es. 20. Scrivere i segnali campionati a tempo TT dei segnali: (1) x(t)=rect⁡(t3T)x(t)=\operatorname{rect}\big(\frac{t}{3T}\big); (2) x(t)=34rect⁡(t−5/23T)x(t)=\frac34\operatorname{rect}\Big(\frac{t-5/2}{3T}\Big); (3) x(t)=exp⁡(−t) 1(t)x(t)=\exp(-t)\,\mathbb 1(t); (4) x(t)=3triang⁡(t6T)+1x(t)=3\operatorname{triang}\big(\frac{t}{6T}\big)+1.
  • Es. 21. Interpolare (con T=0,5T=0{,}5):
    1. x(nT)=34rect⁡(nT−5/37T)x(nT)=\frac34\operatorname{rect}\Big(\frac{nT-5/3}{7T}\Big) con h(t)=rect⁡(tT)h(t)=\operatorname{rect}\big(\frac tT\big);
    2. x(nT)={1n=02n=1−5n=3x(nT)=\begin{cases}1&n=0\\2&n=1\\-5&n=3\end{cases} (e 00 per gli altri nn) con h(t)=triang⁡(t2T)h(t)=\operatorname{triang}\big(\frac{t}{2T}\big);
    3. lo stesso x(nT)x(nT) del punto 2 con h(t)=triang⁡(t)h(t)=\operatorname{triang}(t).

Teoria usata: Campionamento e interpolazione - mantenimento, lineare e idealeIl campionamento ricava da s(t) la sequenza s_c(mT) = s(mT); l'interpolazione ricostruisce s(t) = Σ s_c(mT) h(t − mT) con la risposta impulsiva h dell'interpolatore: mantenimento (rect su [0,T]), lineare (triangolo) o ideale h(t) = sinc(t/T). Il teorema del campionamento garantisce la ricostruzione esatta con la sinc se T è abbastanza piccolo, cioè F_s = 1/T ≥ 2B per un segnale di banda B. La coppia campionatore-interpolatore è una codifica-decodifica di sorgente senza perdite se il teorema è rispettato.Campionamento e interpolazione - mantenimento, lineare e ideale →, Segnali a tempo continuo - sinusoidi, rect e trasformazioni dell'asse dei tempiUn segnale è una funzione del tempo; è analogico se tempo e valori sono continui, digitale se entrambi sono discreti. La sinusoide A cos(2πf₀t + φ₀) ha periodo T = 1/f₀. Con rect(t) (1 per |t| ≤ 1/2) e triang(t) = 1 − |t| (|t| ≤ 1) si costruiscono impulsi; rect((t − T/2)/T) è 1 su [0,T]. Traslazione s(t − t₀) sposta a destra di t₀, scalamento s(t/A) dilata di A; in s((t − t₀)/A) si scala prima e si trasla poi. Una trasmissione è una somma Σ a_m z(t − mT). Il canale AWGN è r(t) = A·s_Tx(t − t₀) + w(t).Segnali a tempo continuo - sinusoidi, rect e trasformazioni dell'asse dei tempi →.

Campionare significa calcolare x(kT)x(kT), cioè sostituire t=kTt=kT nell'espressione: l'argomento di rect⁡\operatorname{rect} e triang⁡\operatorname{triang} diventa un'espressione in kk. Interpolare vuol dire y(t)=∑nx(nT) h(t−nT)y(t)=\sum_nx(nT)\,h(t-nT): ogni campione moltiplica una copia di hh ritardata in nTnT.

Esercizio 20: campioni

(1) x(t)=rect⁡(t3T)x(t)=\operatorname{rect}\big(\frac t{3T}\big). Con t=kTt=kT l'argomento è kT3T=k3\frac{kT}{3T}=\frac k3: x(kT)=rect⁡(k3)x(kT)=\operatorname{rect}\big(\frac k3\big). Il rect vale 11 se ∣k3∣≤12\big|\frac k3\big|\le\frac12, cioè ∣k∣≤1,5|k|\le1{,}5, per kk intero k∈{−1,0,1}k\in\{-1,0,1\}: x(kT)=1  per k=−1,0,1,x(kT)=0 altrove.x(kT)=1\ \text{ per }k=-1,0,1,\qquad x(kT)=0\ \text{altrove.} (Tre campioni, indipendenti da TT: il periodo si semplifica perché larghezza e passo sono entrambi multipli di TT.)

(2) x(t)=34rect⁡(t−5/23T)x(t)=\frac34\operatorname{rect}\big(\frac{t-5/2}{3T}\big). x(kT)=34rect⁡(kT−5/23T)x(kT)=\frac34\operatorname{rect}\Big(\frac{kT-5/2}{3T}\Big), uguale a 34\frac34 quando ∣kT−52∣≤3T2|kT-\frac52|\le\frac{3T}2. Qui TT non si semplifica e il numero di campioni dipende da TT: per T=1T=1 la condizione ∣k−2,5∣≤1,5|k-2{,}5|\le1{,}5 dà k∈{1,2,3,4}k\in\{1,2,3,4\}; per T=0,5T=0{,}5, ∣0,5k−2,5∣≤0,75|0{,}5k-2{,}5|\le0{,}75 cioè ∣k−5∣≤1,5|k-5|\le1{,}5 dà k∈{4,5,6}k\in\{4,5,6\}. In ogni caso i campioni valgono 34\frac34 dentro il rettangolo e 00 fuori.

(3) x(t)=e−t1(t)x(t)=e^{-t}\mathbb 1(t). Con t=kTt=kT: x(kT)=e−kT 1(k)x(kT)=e^{-kT}\,\mathbb 1(k), cioè e−kTe^{-kT} per k≥0k\ge0 e 00 per k<0k<0 (il gradino 1(kT)=1(k)\mathbb 1(kT)=\mathbb 1(k) perché T>0T>0). Per T=0,5T=0{,}5: 1; 0,6065; 0,3679; 0,2231; 0,1353;…1;\ 0{,}6065;\ 0{,}3679;\ 0{,}2231;\ 0{,}1353;\dots (progressione geometrica di ragione e−Te^{-T}).

(4) x(t)=3triang⁡(t6T)+1x(t)=3\operatorname{triang}\big(\frac t{6T}\big)+1. x(kT)=3triang⁡(k6)+1x(kT)=3\operatorname{triang}\big(\frac k6\big)+1. Il triangolo vale 1−∣k∣61-\frac{|k|}{6} per ∣k∣≤6|k|\le6 e 00 altrove, quindi x(kT)=3(1−∣k∣6)+1=4−∣k∣2  (∣k∣≤6),x(kT)=1  (∣k∣>6).x(kT)=3\Big(1-\frac{|k|}6\Big)+1=4-\frac{|k|}2\ \ (|k|\le6),\qquad x(kT)=1\ \ (|k|>6). Valori: 4; 3,5; 3; 2,5; 2; 1,5; 14;\ 3{,}5;\ 3;\ 2{,}5;\ 2;\ 1{,}5;\ 1 per ∣k∣=0,…,6|k|=0,\dots,6. La costante +1+1 non è nulla fuori dal supporto del triangolo: la sequenza resta pari a 11 all'infinito (segnale con energia infinita).

Esercizio 21: interpolazione

(1) h(t)=rect⁡(tT)h(t)=\operatorname{rect}\big(\frac tT\big). L'interpolatore è un mantenimento centrato (ogni campione è tenuto per TT, da nT−T2nT-\frac T2 a nT+T2nT+\frac T2; diverso da quello della lezione, con il rect su [0,T][0,T]). La formula è y(t)=∑n=−∞+∞x(nT) h(t−nT)=∑n34rect⁡(nT−5/37T)rect⁡(t−nTT).y(t)=\sum_{n=-\infty}^{+\infty}x(nT)\,h(t-nT)=\sum_n\frac34\operatorname{rect}\Big(\frac{nT-5/3}{7T}\Big)\operatorname{rect}\Big(\frac{t-nT}{T}\Big). Con T=0,5T=0{,}5 i campioni non nulli sono quelli con ∣0,5n−53∣≤1,75|0{,}5n-\frac53|\le1{,}75, cioè −0,17≤n≤6,83-0{,}17\le n\le6{,}83: n=0,1,…,6n=0,1,\dots,6, sette campioni tutti uguali a 34\frac34. I rettangoli rect⁡(t−nTT)\operatorname{rect}\big(\frac{t-nT}{T}\big) di campioni adiacenti si toccano senza sovrapporsi, quindi la somma è un unico rettangolo che va da −T2=−0,25-\frac T2=-0{,}25 a 6T+T2=3,256T+\frac T2=3{,}25: larghezza 7T=3,57T=3{,}5 e centro 3T=1,53T=1{,}5, y(t)=34rect⁡(t−3T7T)=34rect⁡(t−1,53,5).y(t)=\frac34\operatorname{rect}\Big(\frac{t-3T}{7T}\Big)=\frac34\operatorname{rect}\Big(\frac{t-1{,}5}{3{,}5}\Big).

(2) h(t)=triang⁡(t2T)h(t)=\operatorname{triang}\big(\frac t{2T}\big). Qui x(0)=1x(0)=1, x(T)=2x(T)=2, x(3T)=−5x(3T)=-5 e gli altri campioni sono 00, quindi tre termini: y(t)=1⋅triang⁡(t2T)+2triang⁡(t−T2T)−5triang⁡(t−3T2T).y(t)=1\cdot\operatorname{triang}\Big(\frac{t}{2T}\Big)+2\operatorname{triang}\Big(\frac{t-T}{2T}\Big)-5\operatorname{triang}\Big(\frac{t-3T}{2T}\Big). Con T=0,5T=0{,}5, t2T=t\frac t{2T}=t e t−T2T=t−0,5\frac{t-T}{2T}=t-0{,}5 ecc., quindi le espressioni dei punti 2 e 3 coincidono: y(t)=triang⁡(t)+2triang⁡(t−0,5)−5triang⁡(t−1,5).y(t)=\operatorname{triang}(t)+2\operatorname{triang}(t-0{,}5)-5\operatorname{triang}(t-1{,}5). (3) h(t)=triang⁡(t)h(t)=\operatorname{triang}(t): è l'espressione appena scritta (con T=0,5T=0{,}5: 2T=12T=1). Ogni termine è un triangolo di semibase 1=2T1=2T, cioè più largo del passo T=0,5T=0{,}5 dei campioni: i triangoli vicini si sovrappongono. Valori della somma:

tt −0,5-0{,}5 00 0,250{,}25 0,50{,}5 0,750{,}75 11 1,251{,}25 1,51{,}5 1,751{,}75 22 2,52{,}5
y(t)y(t) 0,50{,}5 22 2,252{,}25 2,52{,}5 0,50{,}5 −1,5-1{,}5 −3,25-3{,}25 −5-5 −3,75-3{,}75 −2,5-2{,}5 00

(per esempio y(0)=triang⁡(0)+2triang⁡(−0,5)−5triang⁡(−1,5)=1+2⋅0,5−0=2y(0)=\operatorname{triang}(0)+2\operatorname{triang}(-0{,}5)-5\operatorname{triang}(-1{,}5)=1+2\cdot0{,}5-0=2.) Il segnale è la somma di tre triangoli: salta dal valore 22 al valore −5-5 tra 0,50{,}5 e 1,51{,}5.

Grafico interattivo: Interpolazione dei campioni x(0) = 1, x(T) = 2, x(3T) = −5 con T = 0,5 e h(t) = triang(t) (semibase 2T): y(t) = triang(t) + 2 triang(t − 0,5) − 5 triang(t − 1,5)

Osservazione. In t=nTt=nT questo y(t)y(t) non riproduce i campioni: y(0)=2≠x(0)=1y(0)=2\ne x(0)=1, y(T)=2,5≠x(T)=2y(T)=2{,}5\ne x(T)=2, mentre y(3T)=−5=x(3T)y(3T)=-5=x(3T) solo perché i vicini sono nulli. Il motivo è che h(nT)≠0h(nT)\ne0 per n≠0n\ne0 (h(T)=triang⁡(12)=12h(T)=\operatorname{triang}(\frac12)=\frac12): la semibase 2T2T è troppo larga. Un vero interpolatore lineare ha h(t)=triang⁡(tT)h(t)=\operatorname{triang}\big(\frac tT\big) (semibase TT, nulla in tutti i multipli non nulli di TT) e in questo caso y(t)y(t) passerebbe per i punti (0,1)(0,1), (T,2)(T,2), (2T,0)(2T,0), (3T,−5)(3T,-5), (4T,0)(4T,0). Con T=0,5T=0{,}5 il punto 2 assume triang⁡(t2T)=triang⁡(t)\operatorname{triang}\big(\frac{t}{2T}\big)=\operatorname{triang}(t): coincidenza che rende il punto 3 uguale al punto 2.

(Verifiche con Python: numero e valori dei campioni dell'es. 20 e tabella di y(t)y(t) dell'es. 21. Le espressioni coincidono con quelle della soluzione del corso, che riporta i punti 1 e 2 dell'es. 21; il punto 3 è uguale al 2 per T=0,5T=0{,}5.)

Errori comuni

  • Scrivere x(kT)x(kT) senza sostituire t=kTt=kT anche nei centri e nelle larghezze (nel punto 2 il centro 52\frac52 resta in secondi, non diventa 52T\frac52T).
  • Dimenticare la costante +1+1 fuori dal supporto del triangolo (punto 4).
  • Scambiare il rect rect⁡(tT)\operatorname{rect}\big(\frac tT\big) (centrato in 00) con quello della lezione rect⁡(t−T/2T)\operatorname{rect}\big(\frac{t-T/2}{T}\big) (su [0,T][0,T]): cambia il ritardo di T2\frac T2.
  • Usare un triangolo di semibase 2T2T credendo che sia l'interpolatore lineare: non riproduce i campioni.

Versione ripasso

Testo. Es. 20: campioni a passo TT di rect⁡(t3T)\operatorname{rect}\big(\frac t{3T}\big), 34rect⁡(t−5/23T)\frac34\operatorname{rect}\big(\frac{t-5/2}{3T}\big), e−t1(t)e^{-t}\mathbb 1(t), 3triang⁡(t6T)+13\operatorname{triang}\big(\frac t{6T}\big)+1. Es. 21 (T=0,5T=0{,}5): interpolare (1) x(nT)=34rect⁡(nT−5/37T)x(nT)=\frac34\operatorname{rect}\big(\frac{nT-5/3}{7T}\big) con h=rect⁡(tT)h=\operatorname{rect}\big(\frac tT\big); (2) x(0)=1x(0)=1, x(T)=2x(T)=2, x(3T)=−5x(3T)=-5 con h=triang⁡(t2T)h=\operatorname{triang}\big(\frac t{2T}\big); (3) lo stesso con h=triang⁡(t)h=\operatorname{triang}(t).

(Campionamento e interpolazione - mantenimento, lineare e idealeIl campionamento ricava da s(t) la sequenza s_c(mT) = s(mT); l'interpolazione ricostruisce s(t) = Σ s_c(mT) h(t − mT) con la risposta impulsiva h dell'interpolatore: mantenimento (rect su [0,T]), lineare (triangolo) o ideale h(t) = sinc(t/T). Il teorema del campionamento garantisce la ricostruzione esatta con la sinc se T è abbastanza piccolo, cioè F_s = 1/T ≥ 2B per un segnale di banda B. La coppia campionatore-interpolatore è una codifica-decodifica di sorgente senza perdite se il teorema è rispettato.Campionamento e interpolazione - mantenimento, lineare e ideale →)

  • Campionare: t=kTt=kT. Interpolare: y(t)=∑nx(nT)h(t−nT)y(t)=\sum_nx(nT)h(t-nT).
  • 20.1 rect⁡(k3)\operatorname{rect}\big(\frac k3\big): 11 per k=−1,0,1k=-1,0,1. 20.2 34rect⁡(kT−5/23T)\frac34\operatorname{rect}\big(\frac{kT-5/2}{3T}\big): 34\frac34 se ∣kT−52∣≤3T2|kT-\frac52|\le\frac{3T}2 (T=1T=1: k=1..4k=1..4; T=0,5T=0{,}5: k=4,5,6k=4,5,6). 20.3 e−kT1(k)e^{-kT}\mathbb 1(k). 20.4 3triang⁡(k6)+1=4−∣k∣23\operatorname{triang}\big(\frac k6\big)+1=4-\frac{|k|}2 per ∣k∣≤6|k|\le6 e 11 fuori.
  • 21.1 campioni 34\frac34 per n=0..6n=0..6; rect adiacenti: y(t)=34rect⁡(t−3T7T)=34rect⁡(t−1,53,5)y(t)=\frac34\operatorname{rect}\big(\frac{t-3T}{7T}\big)=\frac34\operatorname{rect}\big(\frac{t-1{,}5}{3{,}5}\big).
  • 21.2/21.3 (2T=12T=1): y(t)=triang⁡(t)+2triang⁡(t−0,5)−5triang⁡(t−1,5)y(t)=\operatorname{triang}(t)+2\operatorname{triang}(t-0{,}5)-5\operatorname{triang}(t-1{,}5) (identici per T=0,5T=0{,}5). y(0)=2y(0)=2, y(0,5)=2,5y(0{,}5)=2{,}5, y(1,5)=−5y(1{,}5)=-5: non riproduce i campioni perché h(T)=12≠0h(T)=\frac12\ne0 (servirebbe triang⁡(tT)\operatorname{triang}\big(\frac tT\big)).
  • Valori utili: e−kTe^{-kT} con T=0,5T=0{,}5: 1; 0,6065; 0,3679; 0,2231; 0,13531;\ 0{,}6065;\ 0{,}3679;\ 0{,}2231;\ 0{,}1353. Punto 4 per ∣k∣=0,…,6|k|=0,\dots,6: 4; 3,5; 3; 2,5; 2; 1,5; 14;\ 3{,}5;\ 3;\ 2{,}5;\ 2;\ 1{,}5;\ 1. Interpolazione triangolare: y(0)=1+2⋅0,5=2y(0)=1+2\cdot0{,}5=2, y(0,25)=2,25y(0{,}25)=2{,}25, y(0,75)=0,5y(0{,}75)=0{,}5, y(1)=−1,5y(1)=-1{,}5, y(1,25)=−3,25y(1{,}25)=-3{,}25.
  • Errori: non sostituire t=kTt=kT ovunque; dimenticare il +1+1 del punto 4; rect centrato contro rect su [0,T][0,T]; semibase 2T2T scambiata per interpolatore lineare.

Teoria collegata