Salta al contenuto
Note per Studenti Esercizio 32 · capacità di canali gaussiano, binario simmetrico e a cancellazione (esercizio del corso)

Esercizio 32capacità di canali gaussiano, binario simmetrico e a cancellazione (esercizio del corso)

In questa pagina 6

Testo (eserciziario di un corso UniPD equivalente, esercizio sulla capacità). Calcolare la capacità dei seguenti canali con ingresso xx e uscita yy, tempo di simbolo T=0,2T=0{,}2 s, specificando la statistica di xx che raggiunge la capacità:

  1. x∈Rx\in\mathbb R, E[x2]=1E[x^2]=1, y=3x+wy=3x+w con ww gaussiano a media nulla e varianza 44;
  2. x∈Rx\in\mathbb R, E[x2]=1E[x^2]=1, y=3x+wy=3x+w con ww gaussiano a media 44 e varianza 11;
  3. x∈{−1,1}x\in\{-1,1\}, y=D(3x+w)y=D(3x+w) con D(⋅)D(\cdot) demodulazione a minima distanza, ww gaussiano a media nulla e varianza unitaria;
  4. x∈{−2,2}x\in\{-2,2\}, P(y=x∣x)=0,9P(y=x\mid x)=0{,}9 e P(y=0∣x)=0,1P(y=0\mid x)=0{,}1.

Teoria usata: Capacità di canale - canale binario simmetrico, a cancellazione e AWGNLa capacità $C=\max_{p_x}I(x;y)$ è il massimo di informazione mutua tra ingresso e uscita del canale; per il teorema di Shannon si può comunicare con probabilità d'errore arbitrariamente piccola se e solo se il rate è minore di $C$. Per il canale binario simmetrico $C=1-H_2(p)$ bit per uso (ingresso uniforme), per il canale a cancellazione $C=1-\varepsilon$, per l'AWGN $C=\frac12\log_2(1+\text{SNR})$ per uso reale, cioè $C=B\log_2(1+\text{SNR})$ bit/s su una banda $B$. Il limite $R_b<C$ dà il minimo $\frac{E_b}{N_0}\ge\frac{2^\nu-1}\nu$ ($-1{,}59$ dB per $\nu\to0$).Capacità di canale - canale binario simmetrico, a cancellazione e AWGN →, 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 →, 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) Canale AWGN

È un canale additivo gaussiano con guadagno 33: la potenza del segnale utile all'uscita è 32⋅E[x2]=93^2\cdot E[x^2]=9 e quella del rumore 44, quindi SNR=94\text{SNR}=\frac94. La capacità si raggiunge con ingresso gaussiano a media nulla e varianza 11 (imposta dal vincolo di potenza): C=1T⋅12log⁡2(1+94)=10,4log⁡2(3,25)=2,5⋅1,700=4,25 bit/s.C=\frac1T\cdot\frac12\log_2\left(1+\frac94\right)=\frac1{0{,}4}\log_2(3{,}25)=2{,}5\cdot1{,}700=4{,}25\ \text{bit/s}. (Il fattore 12\frac12 perché xx è reale: un uso del canale ogni TT secondi, 12log⁡2(1+SNR)\frac12\log_2(1+\text{SNR}) bit per uso.)

(2) AWGN con rumore a media non nulla

Il rumore ha media 44: y=3x+4+w0y=3x+4+w_0 con w0∼N(0,1)w_0\sim\mathcal N(0,1). La costante 44 è nota e si sottrae all'uscita senza perdere informazione, quindi l'informazione mutua non cambia: stesso ingresso (gaussiano, media nulla, varianza 11) e C=10,4log⁡2(1+91)=2,5log⁡2(10)=8,30 bit/s.C=\frac1{0{,}4}\log_2\left(1+\frac91\right)=2{,}5\log_2(10)=8{,}30\ \text{bit/s}.

(3) Ingresso binario con demodulatore a soglia: canale binario simmetrico

Il demodulatore a minima distanza decide x^=sign(y)\hat x=\text{sign}(y) con y=3x+wy=3x+w: i due punti sono ±3\pm3, la soglia è in 00, il rumore ha σ=1\sigma=1. La probabilità di errore sul bit è p=P[w>3]=Q(31)=Q(3)=1,35⋅10−3.p=P[w>3]=Q\left(\frac31\right)=Q(3)=1{,}35\cdot10^{-3}. Il canale con ingresso e uscita binari e probabilità di errore pp (simmetrica per i due ingressi) è un canale binario simmetrico, con capacità C=1−H2(p)C=1-H_2(p) bit per uso, raggiunta da ingresso uniforme (P[x=1]=12P[x=1]=\frac12): H2(p)=−plog⁡2p−(1−p)log⁡2(1−p)=0,0148 bit,C=1T(1−0,0148)=5⋅0,9852=4,93 bit/s.H_2(p)=-p\log_2p-(1-p)\log_2(1-p)=0{,}0148\ \text{bit},\qquad C=\frac1T\left(1-0{,}0148\right)=5\cdot0{,}9852=4{,}93\ \text{bit/s}.

(4) Canale a cancellazione simmetrico

x∈{−2,2}x\in\{-2,2\}; l'uscita è y=xy=x con probabilità 0,90{,}9 oppure y=0y=0 ("cancellazione": il ricevitore vede che il simbolo è andato perso, non è un errore) con probabilità ε=0,1\varepsilon=0{,}1. È un canale a cancellazione binario simmetrico con C=1T(1−ε)=5⋅0,9=4,5 bit/s,C=\frac1T(1-\varepsilon)=5\cdot0{,}9=4{,}5\ \text{bit/s}, raggiunta da ingresso uniforme (P[x=2]=12P[x=2]=\frac12).

Canale Ingresso ottimo CC [bit/s]
1 AWGN, SNR=94\text{SNR}=\frac94 gaussiano N(0,1)\mathcal N(0,1) 4,254{,}25
2 AWGN, rumore con media gaussiano N(0,1)\mathcal N(0,1) 8,308{,}30
3 BSC, p=Q(3)p=Q(3) uniforme 4,934{,}93
4 cancellazione, ε=0,1\varepsilon=0{,}1 uniforme 4,504{,}50

Il canale 3 ha lo stesso rumore del canale 2 (σ2=1\sigma^2=1, guadagno 33), ma con ingresso binario e decisione a soglia non può superare 11 bit per uso, cioè 1T=5\frac1T=5 bit/s: 4,934{,}93 bit/s sono quasi il massimo del canale binario, mentre il gaussiano del canale 2 arriva a 8,308{,}30 bit/s perché ogni uso porta più di un bit.

(Verificato con Python: 4,25114{,}2511, 8,30488{,}3048, p=1,350⋅10−3p=1{,}350\cdot10^{-3}, H2(p)=0,01481H_2(p)=0{,}01481, C=4,926C=4{,}926, 4,54{,}5.)

Errori comuni

  • Dimenticare il fattore 12\frac12 nella capacità reale 12log⁡2(1+SNR)\frac12\log_2(1+\text{SNR}) (ottenendo 8,58{,}5 bit/s invece di 4,254{,}25).
  • Calcolare l'SNR con il guadagno non al quadrato (33 invece di 99).
  • Pensare che la media del rumore ridotta cambi la capacità: una costante nota si sottrae.
  • Confondere cancellazione ed errore: il canale a cancellazione ha C=1−εC=1-\varepsilon, un BSC con p=0,1p=0{,}1 ha C=1−H2(0,1)=0,531C=1-H_2(0{,}1)=0{,}531.

Versione ripasso

Testo. Capacità, T=0,2T=0{,}2 s, di: y=3x+wy=3x+w (E[x2]=1E[x^2]=1, σw2=4\sigma_w^2=4); idem con media 44 e σ2=1\sigma^2=1; x=±1x=\pm1 con demodulazione a soglia e σ2=1\sigma^2=1; x=±2x=\pm2 con cancellazione 0,10{,}1.

Teoria collegata