Salta al contenuto
Note per Studenti Esercizio 6 · segnale esponenziale, entropie estreme e quantizzatore a 3 bit (temi d'esame giugno 2025 e giugno 2026)

Esercizio 6segnale esponenziale, entropie estreme e quantizzatore a 3 bit (temi d'esame giugno 2025 e giugno 2026)

Esame
In questa pagina 5

Testo (temi d'esame giugno 2025 e giugno 2026, esercizio 1). Si consideri un segnale con distribuzione esponenziale ps(x)=λe−λxp_s(x)=\lambda e^{-\lambda x} per x≥0x\ge0 (nulla altrove), con λ=2,30 V−1\lambda=2{,}30\ \text{V}^{-1}.

  1. Se il segnale fosse quantizzato con un quantizzatore (non necessariamente uniforme) con L=16L=16 livelli, quali sono l'entropia massima e minima in uscita? (Suggerimento: esiste un quantizzatore a L=16L=16 livelli con entropia in uscita di 1010 bit?)
  2. Si consideri un segnale con NN componenti {s1(t),…,sN(t)}\{s_1(t),\dots,s_N(t)\} da quantizzare con un banco di NN quantizzatori a 1616 livelli. Quali sono l'entropia massima e minima in uscita? (giugno 2025)
  3. Si progetti un quantizzatore uniforme a 33 bit tale che la probabilità di saturazione sia al più 10−510^{-5} (giugno 2025) oppure inferiore a 10−310^{-3} (giugno 2026); si rappresenti la funzione caratteristica e si determini l'SNR.
  4. Quanto vale l'entropia in uscita dal quantizzatore? È ottima? Si suggerisca un modo per migliorarla. (giugno 2025)

Teoria usata: Informazione ed entropiaL'informazione di un evento di probabilità $p$ è $\log_2\frac1p$ bit; l'entropia $H(x)=\sum p\log_2\frac1p$ è l'informazione media e misura l'incertezza della sorgente: $0\le H\le\log_2M$, con il massimo quando i simboli sono equiprobabili. Per più simboli: $H(x,y)\le H(x)+H(y)$ (uguaglianza se indipendenti), $H(x|y)=H(x,y)-H(y)$. Per una sorgente con $F_s$ simboli al secondo il rate di informazione è $F_sH_s$, il rate nominale $F_s\log_2M$ e l'efficienza $\eta=\frac{H_s}{\log_2M}$.Informazione ed entropia →, Quantizzatore uniforme - livelli, mid-riser ed erroriUn quantizzatore mappa i campioni reali su $L=2^b$ livelli. Quello uniforme (PCM) sceglie un range dinamico $[-V_{sat},V_{sat}]$ e un passo $\Delta=\frac{2V_{sat}}L$; nel tipo mid-riser i livelli sono $\pm\frac\Delta2,\pm\frac{3\Delta}2,\dots$ e non c'è lo zero. L'errore $e_q=a_q-a$ ha una parte granulare (in $[-\frac\Delta2,\frac\Delta2]$, circa uniforme, potenza $\frac{\Delta^2}{12}$) e una di saturazione (fuori range). Per renderlo piccolo servono $P_{sat}$ piccola e $L$ grande.Quantizzatore uniforme - livelli, mid-riser ed errori →, SNR di quantizzazione e progetto del quantizzatoreL'SNR di quantizzazione è $\Lambda_q=\frac{M_a}{M_e}$. Con errore granulare uniforme e saturazione trascurabile vale $\Lambda_q=\frac{\sigma^2}{\Delta^2/12}=3\frac{\sigma^2}{V_{sat}^2},2^{2b}$, cioè $[\Lambda_q]{dB}=6{,}02,b+4{,}77+20\log{10}\frac\sigma{V_{sat}}$: ogni bit in più dà $+6$ dB. Per progettare: $V_{sat}$ dalla probabilità di saturazione ($V_{sat}=\sigma,Q^{-1}\left(\frac{P_{sat}}2\right)$ per un gaussiano), poi $b$ dall'SNR richiesto, arrotondando per eccesso.SNR di quantizzazione e progetto del quantizzatore →, Codici di Shannon-Fano e di HuffmanIn un codice ottimo le parole più probabili non sono più lunghe di quelle meno probabili e le due parole più lunghe differiscono solo per l'ultimo simbolo. Shannon-Fano costruisce l'albero dall'alto dividendo ripetutamente i simboli in due gruppi di probabilità quasi uguali; Huffman lo costruisce dal basso unendo ogni volta i due simboli meno probabili ed è sempre ottimo tra i codici a prefisso. La lunghezza media $L_y$ è la somma delle probabilità dei nodi uniti, l'efficienza è $\eta=\frac{H}{L_y}$.Codici di Shannon-Fano e di Huffman →, Variabili aleatorie e vettori aleatori per le comunicazioniUna variabile aleatoria è descritta da una PMD (discreta) o da una PDF (continua) e dalla funzione di ripartizione; per la gaussiana $P[x>a]=Q\left(\frac{a-m}{\sigma}\right)$ con la funzione coda $Q$. Media $m_x$, varianza $\sigma_x^2$ e potenza statistica $M_x=\sigma_x^2+m_x^2$. Per un vettore aleatorio contano l'indipendenza, le probabilità condizionate (totali e di Bayes) e la correlazione; combinazioni lineari di gaussiane sono gaussiane, e gaussiane scorrelate sono indipendenti.Variabili aleatorie e vettori aleatori per le comunicazioni →.

(1) Entropia di un quantizzatore a L=16L=16 livelli

L'uscita del quantizzatore è una variabile discreta con alfabeto di 1616 valori. Per la proposizione sull'entropia (Informazione ed entropiaL'informazione di un evento di probabilità $p$ è $\log_2\frac1p$ bit; l'entropia $H(x)=\sum p\log_2\frac1p$ è l'informazione media e misura l'incertezza della sorgente: $0\le H\le\log_2M$, con il massimo quando i simboli sono equiprobabili. Per più simboli: $H(x,y)\le H(x)+H(y)$ (uguaglianza se indipendenti), $H(x|y)=H(x,y)-H(y)$. Per una sorgente con $F_s$ simboli al secondo il rate di informazione è $F_sH_s$, il rate nominale $F_s\log_2M$ e l'efficienza $\eta=\frac{H_s}{\log_2M}$.Informazione ed entropia →) 0≤H≤log⁡2L0\le H\le\log_2L:

  • massima Hmax=log⁡216=4H_{max}=\log_216=4 bit, quando i 1616 livelli sono equiprobabili. Con un quantizzatore non uniforme è possibile ottenerla anche per un segnale esponenziale (scegliendo soglie xkx_k tali che P[xk<x<xk+1]=116P[x_{k}<x<x_{k+1}]=\frac1{16}, cioè xk=−1λln⁡(1−k16)x_k=-\frac1\lambda\ln\left(1-\frac k{16}\right));
  • minima Hmin=0H_{min}=0: se un livello ha probabilità quasi 11 (per esempio una sola soglia molto in alto: quasi tutti i campioni cadono nella prima regione). Con un quantizzatore uniforme non si raggiunge esattamente, ma con un non uniforme sì, al limite.

Suggerimento del testo: un'entropia di 1010 bit richiederebbe L≥210=1024L\ge2^{10}=1024 livelli; con L=16L=16 è impossibile (10>410>4).

(2) NN quantizzatori in parallelo

Il vettore di uscita ha NN componenti (alfabeto 16N16^N). L'entropia congiunta soddisfa max⁡iH(qi)≤H(q)≤∑iH(qi)≤Nlog⁡216\max_iH(q_i)\le H(\mathbf q)\le\sum_iH(q_i)\le N\log_216 (Informazione ed entropiaL'informazione di un evento di probabilità $p$ è $\log_2\frac1p$ bit; l'entropia $H(x)=\sum p\log_2\frac1p$ è l'informazione media e misura l'incertezza della sorgente: $0\le H\le\log_2M$, con il massimo quando i simboli sono equiprobabili. Per più simboli: $H(x,y)\le H(x)+H(y)$ (uguaglianza se indipendenti), $H(x|y)=H(x,y)-H(y)$. Per una sorgente con $F_s$ simboli al secondo il rate di informazione è $F_sH_s$, il rate nominale $F_s\log_2M$ e l'efficienza $\eta=\frac{H_s}{\log_2M}$.Informazione ed entropia →):

  • massima Hmax=Nlog⁡216=4NH_{max}=N\log_216=4N bit: componenti indipendenti e ciascuna a livelli equiprobabili;
  • minima Hmin=0H_{min}=0: tutte le uscite quasi deterministiche (e, a parità di quantizzatori, la dipendenza tra le componenti — per esempio componenti identiche — riduce la congiunta fino a H(q1)H(q_1)).

(3) Quantizzatore uniforme a 3 bit

Range. Il segnale è solo positivo: un range simmetrico [−Vsat,Vsat][-V_{sat},V_{sat}] sprecherebbe metà dei livelli sul semiasse negativo, dove non cade mai un campione. Si usa [0,Vsat][0,V_{sat}] con L=23=8L=2^3=8 livelli, Δ=VsatL\Delta=\frac{V_{sat}}L, e la probabilità di saturazione è la coda dell'esponenziale: Psat=P[s>Vsat]=e−λVsatP_{sat}=P[s>V_{sat}]=e^{-\lambda V_{sat}}, da cui Vsat=1λln⁡1Psat.V_{sat}=\frac1\lambda\ln\frac1{P_{sat}}.

Giugno 2025 (Psat≤10−5P_{sat}\le10^{-5}) Giugno 2026 (Psat<10−3P_{sat}<10^{-3})
VsatV_{sat} ln⁡1052,3=5,006\frac{\ln10^5}{2{,}3}=5{,}006 V ln⁡1032,3=3,003\frac{\ln10^3}{2{,}3}=3{,}003 V
Δ=Vsat8\Delta=\frac{V_{sat}}8 0,62570{,}6257 V 0,37540{,}3754 V
livelli (i+12)Δ(i+\frac12)\Delta 0,313; 0,939; 1,564; 2,190; 2,816; 3,441; 4,067; 4,6930{,}313;\ 0{,}939;\ 1{,}564;\ 2{,}190;\ 2{,}816;\ 3{,}441;\ 4{,}067;\ 4{,}693 V 0,188; 0,563; 0,939; 1,314; 1,689; 2,065; 2,440; 2,8160{,}188;\ 0{,}563;\ 0{,}939;\ 1{,}314;\ 1{,}689;\ 2{,}065;\ 2{,}440;\ 2{,}816 V
soglie iΔi\Delta 0,626, 1,251,…, 4,3800{,}626,\ 1{,}251,\dots,\ 4{,}380 V 0,375, 0,751,…, 2,6280{,}375,\ 0{,}751,\dots,\ 2{,}628 V

Funzione caratteristica (giugno 2025): 88 gradini larghi Δ=0,626\Delta=0{,}626 V sull'asse positivo, il livello Qi=(i+12)ΔQ_i=\left(i+\frac12\right)\Delta per iΔ≤s<(i+1)Δi\Delta\le s<(i+1)\Delta, e l'ultimo gradino (4,6934{,}693 V) si estende a +∞+\infty (saturazione); per s<0s<0 non serve (probabilità nulla).

Grafico interattivo: Funzione caratteristica del quantizzatore a 3 bit per segnale positivo, range [0, 5,006] V (Δ = 0,626 V): 8 gradini di livelli 0,313; 0,939; …; 4,693 V; oltre 5,006 V satura

SNR. La potenza statistica del segnale (non la varianza: la media non è nulla) è Ms=E[s2]=2λ2=0,378M_s=E[s^2]=\frac2{\lambda^2}=0{,}378 V². Con Me=Δ212M_e=\frac{\Delta^2}{12} (errore granulare): Λq=MsΔ2/12=0,3780,62572/12=11,59 ⟹ 10,64 dB(giugno 2026: 0,3780,37542/12=32,2→15,08 dB).\Lambda_q=\frac{M_s}{\Delta^2/12}=\frac{0{,}378}{0{,}6257^2/12}=11{,}59\ \Longrightarrow\ 10{,}64\ \text{dB}\quad(\text{giugno 2026: }\frac{0{,}378}{0{,}3754^2/12}=32{,}2\to15{,}08\ \text{dB}). Includendo anche l'errore di saturazione (integrazione numerica) l'SNR vale 10,3610{,}36 dB (2025) e 14,7714{,}77 dB (2026): 0,30{,}3 dB in meno. (Con un range simmetrico [−5,006,5,006][-5{,}006,5{,}006] V e 88 livelli si avrebbe Δ=1,25\Delta=1{,}25 V e solo 4,64{,}6 dB: metà dei livelli sarebbe sprecata.)

(4) Entropia in uscita

Le probabilità dei livelli sono quelle dei tratti [iΔ,(i+1)Δ)[i\Delta,(i+1)\Delta), con l'ultimo che comprende la coda: Pi=e−λiΔ−e−λ(i+1)Δ=e−λiΔ(1−e−λΔ)P_i=e^{-\lambda i\Delta}-e^{-\lambda(i+1)\Delta}=e^{-\lambda i\Delta}\left(1-e^{-\lambda\Delta}\right), con e−λΔ=e−1,439=0,2372e^{-\lambda\Delta}=e^{-1{,}439}=0{,}2372 (giugno 2025). Sono in progressione geometrica di ragione 0,23720{,}2372:

livello 00 11 22 33 44 55 66 77 (satur.)
PiP_i 0,76290{,}7629 0,18090{,}1809 0,04290{,}0429 0,01020{,}0102 0,00240{,}0024 0,000570{,}00057 0,000140{,}00014 0,000040{,}00004

H=∑Pilog⁡21Pi=1,036 bit(giugno 2026: 1,694 bit).H=\sum P_i\log_2\frac1{P_i}=1{,}036\ \text{bit}\quad(\text{giugno 2026: }1{,}694\ \text{bit}). Non è ottima: il massimo per 88 livelli sarebbe 33 bit e l'efficienza della codifica a lunghezza fissa 33 bit è η=1,0363=0,35\eta=\frac{1{,}036}3=0{,}35 (2026: 0,560{,}56): il quantizzatore uniforme su un segnale esponenziale usa male i livelli (quelli alti sono rarissimi). Come migliorarla: (a) codificare la sorgente a lunghezza variabile (Huffman, Codici di Shannon-Fano e di HuffmanIn un codice ottimo le parole più probabili non sono più lunghe di quelle meno probabili e le due parole più lunghe differiscono solo per l'ultimo simbolo. Shannon-Fano costruisce l'albero dall'alto dividendo ripetutamente i simboli in due gruppi di probabilità quasi uguali; Huffman lo costruisce dal basso unendo ogni volta i due simboli meno probabili ed è sempre ottimo tra i codici a prefisso. La lunghezza media $L_y$ è la somma delle probabilità dei nodi uniti, l'efficienza è $\eta=\frac{H}{L_y}$.Codici di Shannon-Fano e di Huffman →): lunghezze 1,2,3,4,5,6,7,71,2,3,4,5,6,7,7 e Ly=1,311L_y=1{,}311 bit, η=0,79\eta=0{,}79 (2026: Ly=1,725L_y=1{,}725, η=0,98\eta=0{,}98); (b) usare un quantizzatore non uniforme con soglie scelte per rendere più uniformi le probabilità.

(Verificato con Python: Vsat=5,0056V_{sat}=5{,}0056 V, H=1,0359H=1{,}0359 bit, Huffman Ly=1,3108L_y=1{,}3108; 2026: Vsat=3,0034V_{sat}=3{,}0034 V, H=1,6944H=1{,}6944, Ly=1,7251L_y=1{,}7251.)

Errori comuni

  • Dire che l'entropia massima di un quantizzatore a 1616 livelli è 1616 bit: è log⁡216=4\log_216=4 bit.
  • Usare un range simmetrico per un segnale che assume solo valori positivi.
  • Calcolare PsatP_{sat} con 2Q(⋅)2Q(\cdot) (gaussiano): per l'esponenziale è e−λVsate^{-\lambda V_{sat}} (una sola coda).
  • Usare la varianza 1λ2\frac1{\lambda^2} al posto della potenza 2λ2\frac2{\lambda^2} nell'SNR.

Versione ripasso

Testo. Esponenziale λ=2,3 V−1\lambda=2{,}3\ \text{V}^{-1}: entropia max/min per L=16L=16 (e 1010 bit?); banco di NN quantizzatori; quantizzatore uniforme a 33 bit con Psat≤10−5P_{sat}\le10^{-5} (<10−3<10^{-3} nel 2026), funzione caratteristica, SNR; entropia in uscita, ottima? (giugno 2025 e 2026).

Teoria collegata