Vai al contenuto
Teoria dei Segnali

Processi aleatori, densità spettrale e rumore termico

Stazionarietà ed ergodicità, densità spettrale di potenza, passaggio in un filtro con dimostrazione, rumore bianco, rumore termico, rumore kT/C con dimostrazione, shot e 1/f.

Una variabile aleatoria descrive un numero incerto; un processo aleatorio descrive un intero segnale incerto, un valore aleatorio per ogni istante. Il rumore di un circuito, il parlato, i dati trasmessi sono processi aleatori. Non si può prevedere la loro forma d’onda, ma si possono calcolare la potenza, come è distribuita in frequenza e come la modifica un filtro. In questa lezione vediamo questi strumenti e li applichiamo al rumore termico.

Processi aleatori

Un processo X(t) si può pensare come l’insieme di tutte le forme d’onda che un esperimento potrebbe produrre (le realizzazioni). Fissato un istante, X(t_1) è una variabile aleatoria. Le grandezze principali sono

    \[ \mu_X(t) = E[X(t)],\qquad R_X(t_1, t_2) = E[X(t_1)\,X(t_2)] \]

la media e l’autocorrelazione, che dice quanto sono legati i valori in due istanti diversi.

Stazionarietà. Un processo è stazionario in senso lato se la media è costante e l’autocorrelazione dipende solo dalla distanza tra gli istanti: R_X(t + \tau, t) = R_X(\tau). Le sue proprietà statistiche non cambiano nel tempo. R_X(0) = E[X^2] è la potenza media.

Ergodicità. Un processo è ergodico se le medie temporali su una sola realizzazione lunga coincidono con le medie d’insieme. È l’ipotesi che si fa sempre in pratica: si misura il rumore di un amplificatore per un tempo lungo, non di mille amplificatori nello stesso istante.

La densità spettrale di potenza

Definizione (Wiener-Khinchin). Per un processo stazionario, la densità spettrale di potenza è la trasformata di Fourier dell’autocorrelazione:

    \[ S_X(f) = \int_{-\infty}^{+\infty} R_X(\tau)\,e^{-j2\pi f\tau}\,d\tau \]

Antitrasformando in \tau = 0: E[X^2] = R_X(0) = \int S_X(f)\,df. La potenza totale è l’area sotto S_X, e la potenza in una banda è l’area su quella banda: S_X si misura in W/Hz (o V²/Hz).

Il passaggio in un filtro

Teorema. Se un processo stazionario X attraversa un sistema LTI stabile con risposta in frequenza H(f), l’uscita Y è stazionaria e

    \[ S_Y(f) = |H(f)|^2\,S_X(f) \]

Dimostrazione. L’uscita è Y(t) = \int h(\alpha)X(t - \alpha)\,d\alpha. Allora

    \[ R_Y(\tau) = E[Y(t + \tau)Y(t)] = \iint h(\alpha)h(\beta)\,E[X(t + \tau - \alpha)X(t - \beta)]\,d\alpha\,d\beta = \iint h(\alpha)h(\beta)\,R_X(\tau - \alpha + \beta)\,d\alpha\,d\beta \]

che dipende solo da \tau. L’espressione è la convoluzione di R_X con h(\tau) e con h(-\tau). Trasformando, la convoluzione diventa prodotto e la trasformata di h(-\tau) è H^*(f) (risposta reale): S_Y = H\,H^*\,S_X = |H|^2 S_X. ∎

La fase del filtro non conta: per la potenza conta solo il modulo.

Il rumore bianco

Un processo con densità spettrale costante,

    \[ S_N(f) = \frac{N_0}{2} \]

si chiama rumore bianco, per analogia con la luce bianca che contiene tutti i colori. La sua autocorrelazione è \frac{N_0}{2}\delta(\tau): due valori in istanti diversi, anche vicinissimi, sono scorrelati. È un’idealizzazione (avrebbe potenza infinita), ma è un ottimo modello finché la banda del rumore è molto più larga di quella del sistema che lo riceve. Se è anche gaussiano si parla di rumore bianco gaussiano additivo (AWGN), il modello standard dei canali di comunicazione.

Rumore bianco filtrato: densità spettrale d'ingresso, risposta del filtro e densità spettrale d'uscita
A sinistra la densità spettrale costante del rumore bianco. Al centro il modulo al quadrato della risposta di un filtro passa basso. A destra il loro prodotto, la densità spettrale del rumore all’uscita: la potenza ora è finita e concentrata nella banda del filtro.

All’uscita di un filtro la potenza del rumore è \frac{N_0}{2}\int|H(f)|^2\,df. Si definisce banda equivalente di rumore B_N la banda del filtro ideale con lo stesso guadagno massimo che lascia passare la stessa potenza: per un passa basso RC, B_N = \frac{\pi}{2}f_{-3dB}.

Il rumore termico

Gli elettroni di un resistore si agitano per effetto della temperatura, e ai suoi capi compare una tensione di rumore. Il risultato di Johnson e Nyquist (1928) è che, fino a frequenze altissime, è un rumore bianco gaussiano con

    \[ \overline{v_n^2} = 4kTRB \]

in una banda B, dove k = 1{,}38\cdot 10^{-23} J/K è la costante di Boltzmann e T la temperatura assoluta.

Esempio. Un resistore da 1 kΩ a 300 K, nella banda audio di 20 kHz: v_n = \sqrt{4\cdot 1{,}38\cdot 10^{-23}\cdot 300\cdot 1000\cdot 20\,000} \approx 0{,}58 µV. Sembra poco, ma è il limite al di sotto del quale nessun microfono o sensore collegato a quella resistenza può essere amplificato senza rumore.

Teorema (rumore kT/C). La tensione di rumore ai capi del condensatore di un circuito RC ha potenza

    \[ \overline{v_C^2} = \frac{kT}{C} \]

indipendente da R.

Dimostrazione. Il rumore del resistore, con densità 4kTR sulle sole frequenze positive, passa attraverso il partitore RC: |H(f)|^2 = 1/(1 + (2\pi fRC)^2). Integrando:

    \[ \overline{v_C^2} = \int_0^\infty \frac{4kTR}{1 + (2\pi fRC)^2}\,df = 4kTR\cdot\frac{1}{2\pi RC}\cdot\frac{\pi}{2} = \frac{kT}{C} \]

usando \int_0^\infty \frac{du}{1 + u^2} = \frac{\pi}{2}. ∎

Una resistenza più grande fa più rumore, ma restringe anche la banda nella stessa proporzione. Con C = 1 pF, \sqrt{kT/C} \approx 64 µV: è il rumore fondamentale dei circuiti a capacità commutate e dei pixel dei sensori d’immagine, e il motivo per cui i condensatori di campionamento dei convertitori di precisione non possono essere troppo piccoli.

Altri rumori

  • Rumore shot: la corrente in una giunzione è fatta di cariche discrete che la attraversano a caso; è bianco con \overline{i_n^2} = 2qIB.
  • Rumore 1/f (flicker): densità spettrale che cresce alle basse frequenze, dovuta a difetti nei materiali; domina nei MOSFET sotto qualche kilohertz.
  • Processi di Markov: in molti sistemi (code di pacchetti, canali con errori a raffica) il futuro dipende solo dallo stato presente; le catene di Markov sono il modello di base per descriverli.
Subscribe
Notify of
guest
0 Comments
Oldest
Newest Most Voted
0
Would love your thoughts, please comment.x
()
x