Salta al contenuto
Note per Studenti Esercizio - M-PAM a 4 Mbit/s con impulso a due tratti, energia necessaria e attenuazione della fibra (primo compitino del 29 novembre 2024, es. 183)

Esercizio - M-PAM a 4 Mbitᐟs con impulso a due tratti, energia necessaria e attenuazione della fibra (primo compitino del 29 novembre 2024, es. 183)

Esame
In questa pagina 6

Testo (primo compitino del 29 novembre 2024, corso di Fondamenti di Telecomunicazioni, Ing. Informatica UniPD). Si vuole trasmettere un flusso di bit con velocità (bit rate) di 44 Mbit/s utilizzando una modulazione MM-PAM con periodo di simbolo T=1 μT=1\ \mus, su canale AWGN con potenza di rumore σI2=10−2\sigma_I^2=10^{-2}.

  1. Calcolare il valore di MM (cardinalità della costellazione) e specificare la mappatura da bit a simbolo che minimizza la probabilità di errore.
  2. Si consideri come forma d'onda prototipo (o impulso fondamentale) hTx(t)={A tt∈[0,T/2]A t/2t∈(T/2,T]0altrove.h_{Tx}(t)=\begin{cases}A\,t&t\in[0,T/2]\\ A\,t/2&t\in(T/2,T]\\ 0&\text{altrove.}\end{cases} Calcolare un valore di AA perché l'energia di hTx(t)h_{Tx}(t) valga 10−1910^{-19}.
  3. Usando ora una forma d'onda prototipo hTx(t)h_{Tx}(t) rettangolare, calcolare l'energia media in ricezione per avere una probabilità di errore sul simbolo di 10−310^{-3}.
  4. Con hTx(t)h_{Tx}(t) del punto precedente, disegnare il ricevitore ottimo con il minor numero di filtri (o correlatori), specificando il valore di tutti i parametri.
  5. Il segnale modulato viene ora trasmesso con potenza 1010 W su una fibra ottica. Sapendo che all'uscita di un tratto di fibra lungo 150150 km la potenza ricevuta è 1,71{,}7 mW, calcolare l'attenuazione della fibra in dB/km.

Teoria usata: Modulazione PAM - livelli, energia media, codifica di Gray e probabilità d'erroreNella PAM i segnali sono un impulso $h_{Tx}(t)$ moltiplicato per un'ampiezza $\alpha_m=2m-1-M$: $s_m(t)=\alpha_mh_{Tx}(t)$. La base è $\phi_1=\frac{h_{Tx}}{\sqrt{E_h}}$, i punti sono $\alpha_m\sqrt{E_h}$ su una retta, $d_{\min}=2\sqrt{E_h}$, $E_s=\frac{M^2-1}{3}E_h$. I due simboli esterni sbagliano da una parte sola, gli $M-2$ interni da due: $P[E]=2\left(1-\frac1M\right)Q\left(\sqrt{\frac{E_s}{\sigma_I^2}\frac{3}{M^2-1}}\right)$. Con la codifica di Gray (simboli adiacenti differiscono per un bit) $P_{bit}\simeq\frac{P[E]}{\log_2M}$.Modulazione PAM - livelli, energia media, codifica di Gray e probabilità d'errore →, Ricevitore ottimo - banco di correlatori e filtri adattatiIl ricevitore ottimo calcola le coordinate $r_i[k]=\int r(t),\phi_i(t-kT),dt$ con un banco di $I$ correlatori (moltiplicatore, integratore sul simbolo, campionatore) e le passa a un blocco di decisione MAP, ML o MD. Ogni correlatore è equivalente a un filtro adattato con risposta impulsiva $\phi_i(t_0-t)$, causale se $t_0$ è la fine del supporto, campionato in $t_0+kT$. Alternativa senza base: $M$ filtri adattati ai segnali $s_m(T-t)$ e decisione $\arg\max_m\left[\langle r,s_m\rangle-\frac{E_m}2\right]$ (MD), che conviene se $M=I$. Per $M=2$ basta un solo ramo, lungo $s_1-s_2$, con una soglia.Ricevitore ottimo - banco di correlatori e filtri adattati →, Canali di trasmissione - cavo, fibra ottica e radio in spazio liberoIl canale moltiplica il segnale per un fattore $A$ e aggiunge il rumore, che nasce al ricevitore: $r(t)=A,s_{Tx}(t)+w(t)$, con $\frac{E_r}{E_{Tx}}=A^2$. Se $A<1$ c'è attenuazione $a=\frac1A$, in dB $(a^2){dB}=10\log{10}a^2$. Cavo e fibra: attenuazione proporzionale alla lunghezza, $(a^2){dB}=\tilde a\cdot l{km}$ (cavo in rame: $\tilde a$ cresce con $\sqrt f$; fibra circa $0{,}2$ dB/km). Radio in spazio libero: $(a){dB}=32{,}4+20\log{10}d_{km}+20\log_{10}f_{MHz}-g_{Tx,dB}-g_{Rx,dB}$. L'energia al ricevitore $E_s=A^2E_{Tx}$ è quella da usare nelle formule di $P[E]$; poiché $P=\frac{E}{T}$, vale $\frac{E_r}{E_{Tx}}=\frac{P_r}{P_{Tx}}$.Canali di trasmissione - cavo, fibra ottica e radio in spazio libero →, Variabili aleatorie continue - uniforme, gaussiana, esponenziale e funzione QUna v.a. continua ha funzione di distribuzione F_x(a) = P(x ≤ a) e densità f_x(a) = dF_x/da ≥ 0 con integrale 1; P(A ≤ x ≤ B) = F_x(B) − F_x(A) = ∫_A^B f_x e P(x = A) = 0. Uniforme U([a,b]): densità 1/(b−a), media (a+b)/2, varianza (b−a)²/12. Esponenziale unilatera di media m: (1/m)e^(−a/m), a ≥ 0. Gaussiana N(m, σ²): P(x ≥ A) = Q((A − m)/σ), con Q(A) coda della gaussiana standard. Somme di gaussiane indipendenti sono gaussiane (medie e varianze si sommano); teorema del limite centrale.Variabili aleatorie continue - uniforme, gaussiana, esponenziale e funzione Q →.

(1) Cardinalità e mappatura

Una MM-PAM porta log⁡2M\log_2M bit per simbolo e si trasmette un simbolo ogni TT: il bit rate è Rb=log⁡2MTR_b=\frac{\log_2M}{T}. Quindi log⁡2M=RbT=4⋅106⋅10−6=4 ⇒ M=24=16.\log_2M=R_bT=4\cdot10^6\cdot10^{-6}=4\ \Rightarrow\ M=2^4=16. Mappatura. Il canale è AWGN e il ricevitore decide a minima distanza (simboli equiprobabili): l'errore più probabile è finire nella regione adiacente. La mappatura che minimizza la probabilità d'errore sul bit è la codifica di Gray, nella quale simboli adiacenti differiscono per un solo bit, così un errore di simbolo comporta quasi sempre un solo bit sbagliato. Per M=16M=16 si costruisce per riflessione, dal livello più basso a quello più alto: 0000, 0001, 0011, 0010, 0110, 0111, 0101, 0100, 1100, 1101, 1111, 1110, 1010, 1011, 1001, 10000000,\,0001,\,0011,\,0010,\,0110,\,0111,\,0101,\,0100,\,1100,\,1101,\,1111,\,1110,\,1010,\,1011,\,1001,\,1000 (ogni parola differisce dalla precedente per un bit; formula: g(n)=n XOR ⌊n/2⌋g(n)=n\ \mathrm{XOR}\ \lfloor n/2\rfloor). Con Gray Pbit≃P[E]4P_{bit}\simeq\frac{P[E]}{4}.

(2) Valore di AA per Eh=10−19E_h=10^{-19}

L'energia di hTxh_{Tx} è l'integrale del suo quadrato, tratto per tratto: Eh=∫0T/2A2t2 dt+∫T/2T(A t2)2dt=A2(T/2)33+A24⋅T3−(T/2)33.E_h=\int_0^{T/2}A^2t^2\,dt+\int_{T/2}^{T}\left(\frac{A\,t}{2}\right)^2dt=A^2\frac{(T/2)^3}{3}+\frac{A^2}{4}\cdot\frac{T^3-(T/2)^3}{3}.

  • Primo tratto: A2T324A^2\frac{T^3}{24}.
  • Secondo tratto: A212(T3−T38)=A212⋅7T38=7A2T396\frac{A^2}{12}\left(T^3-\frac{T^3}{8}\right)=\frac{A^2}{12}\cdot\frac{7T^3}{8}=\frac{7A^2T^3}{96}.
  • Somma: Eh=A2T3(496+796)=1196A2T3E_h=A^2T^3\left(\frac{4}{96}+\frac{7}{96}\right)=\frac{11}{96}A^2T^3.

Con T=10−6T=10^{-6} s, T3=10−18T^3=10^{-18}: Eh=1,1458⋅10−19 A2E_h=1{,}1458\cdot10^{-19}\,A^2. Imponendo Eh=10−19E_h=10^{-19}: A2=10−191,1458⋅10−19=0,8727 ⇒ A=0,934A^2=\frac{10^{-19}}{1{,}1458\cdot10^{-19}}=0{,}8727\ \Rightarrow\ A=0{,}934 (in V/s: h(t)=Ath(t)=At ha le unità di una tensione se AA è una tensione divisa per un tempo; il corso usa il numero senza unità). Verificato anche con integrazione numerica (2⋅1062\cdot10^6 punti: Eh=1,1458⋅10−19A2E_h=1{,}1458\cdot10^{-19}A^2).

(3) Energia media in ricezione per P[E]=10−3P[E]=10^{-3}

Con M=16M=16 la probabilità d'errore della PAM è P[E]=2(1−1M)Q(EsσI2⋅3M2−1)=158 Q(3Es255 σI2),P[E]=2\left(1-\frac1M\right)Q\left(\sqrt{\frac{E_s}{\sigma_I^2}\cdot\frac{3}{M^2-1}}\right)=\frac{15}{8}\,Q\left(\sqrt{\frac{3E_s}{255\,\sigma_I^2}}\right), con EsE_s energia media al ricevitore (Modulazione PAM - livelli, energia media, codifica di Gray e probabilità d'erroreNella PAM i segnali sono un impulso $h_{Tx}(t)$ moltiplicato per un'ampiezza $\alpha_m=2m-1-M$: $s_m(t)=\alpha_mh_{Tx}(t)$. La base è $\phi_1=\frac{h_{Tx}}{\sqrt{E_h}}$, i punti sono $\alpha_m\sqrt{E_h}$ su una retta, $d_{\min}=2\sqrt{E_h}$, $E_s=\frac{M^2-1}{3}E_h$. I due simboli esterni sbagliano da una parte sola, gli $M-2$ interni da due: $P[E]=2\left(1-\frac1M\right)Q\left(\sqrt{\frac{E_s}{\sigma_I^2}\frac{3}{M^2-1}}\right)$. Con la codifica di Gray (simboli adiacenti differiscono per un bit) $P_{bit}\simeq\frac{P[E]}{\log_2M}$.Modulazione PAM - livelli, energia media, codifica di Gray e probabilità d'errore →). Si impone P[E]=10−3P[E]=10^{-3}: Q(x)=10−315/8=5,333⋅10−4 ⇒ x=Q−1(5,333⋅10−4)=3,2723,Q(x)=\frac{10^{-3}}{15/8}=5{,}333\cdot10^{-4}\ \Rightarrow\ x=Q^{-1}(5{,}333\cdot10^{-4})=3{,}2723, e x2=3Es255σI2x^2=\frac{3E_s}{255\sigma_I^2} dà Es=x2⋅255 σI23=10,708⋅255⋅10−23=9,10 V2s.E_s=\frac{x^2\cdot255\,\sigma_I^2}{3}=\frac{10{,}708\cdot255\cdot10^{-2}}{3}=9{,}10\ \text{V}^2\text{s}. (Il rapporto è EsσI2=910\frac{E_s}{\sigma_I^2}=910, cioè 29,629{,}6 dB: per 1616 livelli serve molto SNR. La soluzione del corso riporta Es=9,01E_s=9{,}01; con questo valore si otterrebbe P[E]=1,06⋅10−3P[E]=1{,}06\cdot10^{-3}, quindi 9,109{,}10 è quello che dà esattamente 10−310^{-3}.)

(4) Ricevitore ottimo con il minimo numero di filtri

La PAM ha una sola dimensione (I=1I=1): basta un filtro (o un correlatore). Con hTxh_{Tx} rettangolare di durata TT, la base è ϕ1(t)=hTx(t)Eh=1Trect⁡(t−T/2T)\phi_1(t)=\frac{h_{Tx}(t)}{\sqrt{E_h}}=\frac1{\sqrt T}\operatorname{rect}\left(\frac{t-T/2}{T}\right).

(Nota: EsEh=M2−13=85\frac{E_s}{E_h}=\frac{M^2-1}{3}=85.)

(5) Attenuazione della fibra

Il segnale parte con potenza 1010 W e arriva con 1,71{,}7 mW dopo 150150 km. Il rapporto di potenze è PTxPRx=101,7⋅10−3=5882\frac{P_{Tx}}{P_{Rx}}=\frac{10}{1{,}7\cdot10^{-3}}=5882 e in dB (a)dB=10log⁡105882=37,70 dB.(a)_{dB}=10\log_{10}5882=37{,}70\ \text{dB}. Poiché l'attenuazione in dB è proporzionale alla lunghezza (Canali di trasmissione - cavo, fibra ottica e radio in spazio liberoIl canale moltiplica il segnale per un fattore $A$ e aggiunge il rumore, che nasce al ricevitore: $r(t)=A,s_{Tx}(t)+w(t)$, con $\frac{E_r}{E_{Tx}}=A^2$. Se $A<1$ c'è attenuazione $a=\frac1A$, in dB $(a^2){dB}=10\log{10}a^2$. Cavo e fibra: attenuazione proporzionale alla lunghezza, $(a^2){dB}=\tilde a\cdot l{km}$ (cavo in rame: $\tilde a$ cresce con $\sqrt f$; fibra circa $0{,}2$ dB/km). Radio in spazio libero: $(a){dB}=32{,}4+20\log{10}d_{km}+20\log_{10}f_{MHz}-g_{Tx,dB}-g_{Rx,dB}$. L'energia al ricevitore $E_s=A^2E_{Tx}$ è quella da usare nelle formule di $P[E]$; poiché $P=\frac{E}{T}$, vale $\frac{E_r}{E_{Tx}}=\frac{P_r}{P_{Tx}}$.Canali di trasmissione - cavo, fibra ottica e radio in spazio libero →), a~=37,70 dB150 km=0,251 dB/km.\tilde a=\frac{37{,}70\ \text{dB}}{150\ \text{km}}=0{,}251\ \text{dB/km}. (Ordine di grandezza tipico di una fibra: 0,20{,}2 dB/km.)

Errori comuni

  • Scrivere M=4M=4 (44 Mbit/s sembra "quattro"): i bit per simbolo sono RbT=4R_bT=4, quindi M=24=16M=2^4=16.
  • Usare la mappatura naturale (binario) invece di Gray: a parità di P[E]P[E] raddoppia circa PbitP_{bit}.
  • Dimenticare il fattore 2(1−1M)=1582\left(1-\frac1M\right)=\frac{15}{8} nell'inversione della QQ.
  • Confondere EsE_s con EhE_h: Es=M2−13Eh=85EhE_s=\frac{M^2-1}{3}E_h=85E_h.
  • Usare 20log⁡1020\log_{10} per un rapporto di potenze (è 10log⁡1010\log_{10}).

Versione ripasso

Testo. 44 Mbit/s con MM-PAM, T=1 μT=1\ \mus, σI2=10−2\sigma_I^2=10^{-2}. (1) MM e mappatura ottima; (2) AA perché Eh=10−19E_h=10^{-19} con h=Ath=At su [0,T2][0,\frac T2], At2\frac{At}2 su (T2,T](\frac T2,T]; (3) EsE_s in ricezione per P[E]=10−3P[E]=10^{-3} con impulso rettangolare; (4) ricevitore con il minimo numero di filtri; (5) fibra: 1010 W →\to 1,71{,}7 mW dopo 150150 km, dB/km.

(1) Modulazione PAM - livelli, energia media, codifica di Gray e probabilità d'erroreNella PAM i segnali sono un impulso $h_{Tx}(t)$ moltiplicato per un'ampiezza $\alpha_m=2m-1-M$: $s_m(t)=\alpha_mh_{Tx}(t)$. La base è $\phi_1=\frac{h_{Tx}}{\sqrt{E_h}}$, i punti sono $\alpha_m\sqrt{E_h}$ su una retta, $d_{\min}=2\sqrt{E_h}$, $E_s=\frac{M^2-1}{3}E_h$. I due simboli esterni sbagliano da una parte sola, gli $M-2$ interni da due: $P[E]=2\left(1-\frac1M\right)Q\left(\sqrt{\frac{E_s}{\sigma_I^2}\frac{3}{M^2-1}}\right)$. Con la codifica di Gray (simboli adiacenti differiscono per un bit) $P_{bit}\simeq\frac{P[E]}{\log_2M}$.Modulazione PAM - livelli, energia media, codifica di Gray e probabilità d'errore →

  • log⁡2M=RbT=4\log_2M=R_bT=4, M=16M=16; mappatura di Gray (adiacenti differiscono per un bit), Pbit≃P[E]4P_{bit}\simeq\frac{P[E]}4.

(2) Eh=A2T324+7A2T396=1196A2T3=1,1458⋅10−19A2⇒A=0,934E_h=\frac{A^2T^3}{24}+\frac{7A^2T^3}{96}=\frac{11}{96}A^2T^3=1{,}1458\cdot10^{-19}A^2\Rightarrow A=0{,}934.

(3) P[E]=158Q(3Es255σI2)=10−3P[E]=\frac{15}8Q\left(\sqrt{\frac{3E_s}{255\sigma_I^2}}\right)=10^{-3}: Q(x)=5,333⋅10−4Q(x)=5{,}333\cdot10^{-4}, x=3,272x=3{,}272; Es=x2⋅255σI23=9,10E_s=\frac{x^2\cdot255\sigma_I^2}3=9{,}10 V2^2s (EsσI2=910\frac{E_s}{\sigma_I^2}=910, 29,629{,}6 dB). (Il corso: 9,019{,}01, che darebbe 1,06⋅10−31{,}06\cdot10^{-3}.)

(4) Ricevitore ottimo - banco di correlatori e filtri adattatiIl ricevitore ottimo calcola le coordinate $r_i[k]=\int r(t),\phi_i(t-kT),dt$ con un banco di $I$ correlatori (moltiplicatore, integratore sul simbolo, campionatore) e le passa a un blocco di decisione MAP, ML o MD. Ogni correlatore è equivalente a un filtro adattato con risposta impulsiva $\phi_i(t_0-t)$, causale se $t_0$ è la fine del supporto, campionato in $t_0+kT$. Alternativa senza base: $M$ filtri adattati ai segnali $s_m(T-t)$ e decisione $\arg\max_m\left[\langle r,s_m\rangle-\frac{E_m}2\right]$ (MD), che conviene se $M=I$. Per $M=2$ basta un solo ramo, lungo $s_1-s_2$, con una soglia.Ricevitore ottimo - banco di correlatori e filtri adattati →

  • Un solo filtro adattato hRc=ϕ1(t0−t)=1Trect⁡h_{Rc}=\phi_1(t_0-t)=\frac1{\sqrt T}\operatorname{rect}, t0=Tt_0=T; Eh=Es85=0,1071E_h=\frac{E_s}{85}=0{,}1071, Eh=0,327\sqrt{E_h}=0{,}327; soglie 2jEh2j\sqrt{E_h}, j=−7,…,7j=-7,\dots,7; poi bit di Gray.

(5) Canali di trasmissione - cavo, fibra ottica e radio in spazio liberoIl canale moltiplica il segnale per un fattore $A$ e aggiunge il rumore, che nasce al ricevitore: $r(t)=A,s_{Tx}(t)+w(t)$, con $\frac{E_r}{E_{Tx}}=A^2$. Se $A<1$ c'è attenuazione $a=\frac1A$, in dB $(a^2){dB}=10\log{10}a^2$. Cavo e fibra: attenuazione proporzionale alla lunghezza, $(a^2){dB}=\tilde a\cdot l{km}$ (cavo in rame: $\tilde a$ cresce con $\sqrt f$; fibra circa $0{,}2$ dB/km). Radio in spazio libero: $(a){dB}=32{,}4+20\log{10}d_{km}+20\log_{10}f_{MHz}-g_{Tx,dB}-g_{Rx,dB}$. L'energia al ricevitore $E_s=A^2E_{Tx}$ è quella da usare nelle formule di $P[E]$; poiché $P=\frac{E}{T}$, vale $\frac{E_r}{E_{Tx}}=\frac{P_r}{P_{Tx}}$.Canali di trasmissione - cavo, fibra ottica e radio in spazio libero →

  • 101,7⋅10−3=5882=37,70\frac{10}{1{,}7\cdot10^{-3}}=5882=37{,}70 dB; 37,70150=0,251\frac{37{,}70}{150}=0{,}251 dB/km.

Errori tipici

  • M=4M=4 invece di 1616 (log⁡2M=RbT\log_2M=R_bT).
  • Mappatura naturale al posto di Gray.
  • Dimenticare 2(1−1M)2\left(1-\frac1M\right) nell'inversione.
  • EsE_s confuso con EhE_h (Es=85EhE_s=85E_h); 20log⁡1020\log_{10} per potenze.

Teoria collegata