Ocw Uc3m TC Diapos t2 Ruido 2p
Ocw Uc3m TC Diapos t2 Ruido 2p
Ocw Uc3m TC Diapos t2 Ruido 2p
C AP ÍTULO 2
RUIDO EN LOS SISTEMAS DE COMUNICACIONES
Marcelino Lázaro
1 / 88
Índice de contenidos
Revisión de conceptos de probabilidad
I Variable aleatoria
I Procesos aleatorios
Caracterización del ruido en sistemas de comunicaciones
I Procesos blancos
I Procesos gausianos
I Suma de procesos aleatorios
I Modelo estadı́stico del ruido térmico
I Relación señal a ruido
q q
.... .... ⌦
.... ...
.... ...
q q
.. 2 4 ...
..
.. .... ....... .............................
⌦ ! IR .
... ...1...... ...
... .
...
.........................
.....
3
..........................
..... .
.. .........
.
..
.
.......
....
....
... . ......... ........ . ....
..... .
. ........ .
..... .. . ...
2 ⌦ ! X( ) 2 IR ......... .. .
.................
.. ..
.....
.
....
... ..........................................................
...... .................
. ..
.
.
...
.
...
...
..
.. ... ... ..
. ... ... ..
... ..
? ..
. ? ..
.. ? ..
..
. ?
..
.. -
X( 2 ) X( 1 ) X( 3 ) X( 4 ) R
I
Rango de X: RangoX = {x 2 IR : 9 2 ⌦, X( ) = x}
I V.a. discreta: rango formado por conjunto discreto de valores
I V.a. continua: rango continuo de valores
Descripción (probabilı́stica):
I Función de distribución: FX (x)
I Función densidad de probabilidad: fx (x)
Función de distribución
Definición
FX (x) = P(X x)
Interpretación frecuencial (probabilı́stica)
nx
FX (x) = P(X x) = lı́m
n!1 n
0.8 0.8
0.6 0.6
0.4 0.4
0.2 0.2
0 0
0 2 4 6 8 −5 0 5
Valores de la v.a.X Valores de la v.a.X
3
0.8
2
1 0.6
0
−1 0.4
−2
0.2 Teórica
−3 Estimada
−4 0
0 100 200 300 400 500 −3 −2 −1 0 1 2 3
Realizaciones de la v.a. X Valores de la v.a.X
P(a X b) = FX (b) FX (a )
P(a < X < b) = FX (b ) FX (a)
P(a X < b) = FX (b ) FX (a )
6 P(X = a) = FX (a) FX (a )
7 P(X > x) = 1 FX (x)
P(x X x + x)
fX (x) = lı́m
x !0 x
⇢
1 nx
fX (x) = lı́m lı́m
x !0 x n!1 n
Estima de la f.d.p.
4 0.5
3
0.4
2
1 0.3
0
−1 0.2
−2
0.1
−3
−4 0
0 100 200 300 400 500 −5 0 5
Realizaciones de la v.a. X Valores de la v.a.X
⇢ fX (x) 6
1
b a, a<x<b
fX (x) =
0, en otro caso 1
b a
-
a b x
Ejemplo aplicación en comunicaciones
I Fase aleatoria en una sinusoide: v.a. uniforme entre 0 y 2⇡
6f
X (x)
p1 ......
2⇡ .... ....
.
1 (x µ)2 .... .....
.
fX (x) = p ·e 2 2 .. ...
2⇡ .. ...
..... ...
...
..... ...
...
..
. . .....
......................... ........................ -
µ x
Ejemplo de aplicación en comunicaciones
I Modelado del ruido térmico
.u ..u ....u
u..u .u
..u ...uu
1
....u
. .
..... ...... ....u....
....... ....u....u.... .........
. . .
. ...... ...
....u .u
. .
...... .....u ................. .... ....
. .
.... .u ........ .... .
u
0.8
...u.....u.... u...u ......... ...... .... .... ........ .... ......... ...... .....u.... ......... ...u
. . . .
. . .
. ..
........
....... ....... ........ .... ............ ....u
0.6 ....... ........ .......... .... ... .... .. . . . .. .. . . . . .
......... ...u..u........ .... ...u.... .... .... ..u........ .... .... ...... ........ ........ ........ ..... .... ........ .......
u... .. ..
........ ... ........ ....u..... .... .... .... .u ..................... .... ................. .... .... ............ .... .... ....u .... ....u ............ .... ......... ........u......... ....u... ......... .... .....
0.4 .
.. . ...
. . . . . .u . . . ... . . . . . . . .. . . . . . . . . . . .u. . . . . . .
.... ....u...... .......... .... ..... ..... .... .... .... .......... ......... .... ..... .u.... .... ....u.... ......... ..... ..... ..... ...u .... .... ....uu.... .... .... .... .... .... ....u ...... ..... ...u.... .....
. .u . . . . . . .. u. .
0.2
.... ..u......... ....u..... ..... ........ ..... ...u.... ...u ..... .....u .... ......... ....u ..... ..... .... ......... ..... ..... ..... ..... ....u..... ..... ........ .... ........ ...u .... .....
. . . . . . .. . . . . .. . .. . .
..... ... ... ... ....... .. ... ..... ... ... u... ...... ...u.. ... ...u..u... .u...u.. ... ... ... .. ......u ..... ... ..... .... ...u...
0
...... ..... ... ..... ... ... ..... .. ... ... .. ..... ...... ... .. ...u... ...... ... .. ...u... ..u ...... ... ........ .... ..
.
....
...... ...... ... ... ... ...u ......... ... ... ..... ... ..... ...... ...u... ...u...... ... ..u... ... ... ... ...u... ....... ......
.......u .....
-0.2 .
.... ...... ... ... ...... u
...... .. .. ........... .. ..... .. .. .... ... ...... ... ...... ..u ... ... ...... ......
-0.4 .u... ... .
. . . . . . .
. .
..... ...u.. ..........
. . . . . .
.... ..... . .
. .u ... ......... ...... ..... .. ....
u....
... .. ... .. .....
... ..
... ..u .u .
.
.u .
. .
. . .
. .
. . . .
... .. ..... .... .
...u .u
.
..
. .
. .
. .
..
. ...u ...... ...u .....
...u
u
-0.6 . ..... ... .
..... .... .. ..... . ... .. ..... . . ......
. . . ... ..... ..u..
. . . . ....
......u... .u
... .. ... .. ...
......u
.....
.... ... ...... ... ..u
.u .u .u .. u ....u..u
-0.8
. .
-1 .u
0 20 40 60 80 100
Realización
c Marcelino Lázaro, 2014 Teorı́a de la Comunicación Variable Aleatoria 15 / 88
14
12
10
0
−1 −0.8 −0.6 −0.4 −0.2 0 0.2 0.4 0.6 0.8 1
200
150
100
50
0
−1 −0.8 −0.6 −0.4 −0.2 0 0.2 0.4 0.6 0.8 1
35
30
25
20
15
10
0
−3 −2 −1 0 1 2 3
600
500
400
300
200
100
0
−4 −3 −2 −1 0 1 2 3 4
µ=0 x x µ=0
c Marcelino Lázaro, 2014 Teorı́a de la Comunicación Variable Aleatoria 21 / 88
d d
µ x x µ
x µ d d
Q =Q =1 Q
2
N (µ, )
d d
x µ µ x
c Marcelino Lázaro, 2014 Teorı́a de la Comunicación Variable Aleatoria 23 / 88
2
Integrales sobre N (µ, ) en intervalos
En general se pueden escribir como sumas o diferencias de diferentes
términos involucrando integrales desde un punto a ±1, que ya hemos
visto como se obtienen utilizando la función Q(x)
Un ejemplo ilustrativo
2
N (µ, )
d1 d2
x1 µ x2
R x2 2
R1 2
R1 2
⇥ d1
⇤ ⇥ d2
⇤
x1
N (µ, )= x1
N (µ, ) x2
N (µ, )= 1 Q Q
d1 d2
x1 µ µ x2
c Marcelino Lázaro, 2014 Teorı́a de la Comunicación Variable Aleatoria 24 / 88
Funciones de una variable aleatoria
Una función Y = g(X) de una v. a. es una variable aleatoria.
Función de distribución
XNr
fX (xi )
fY (y) =
i=1
|g0 (xi )|
Momentos estadı́sticos
Valor esperado (media, o esperanza matemática)
Z 1
mX = E[X] = x · fX (x) dx
1
Momento de orden n
Z 1
mnX = xn · fX (x) dx
1
Varianza
Z 1
2
⇥ 2
⇤
X = E (x mX ) = (x mX )2 · fX (x) dx
1
h i
NOTA: X2 = E (x mX )2 = E[X 2 ] (E[X])2 = E[X 2 ] (mX )2
@2
fX,Y (x, y) = FX,Y (x, y)
@x@y
Casos particulares
I Correlación: g(X, Y) = X · Y
I Covarianza: g(X, Y) = (X mX ) · (Y mY )
Incorrelación
Coeficiente de correlación
Cov(X, Y)
⇢X,Y = , 0 |⇢X,Y | 1
X · Y
Bg,h 2
X,Y (z, w) = {(x, y) 2 IR : g(x, y) z, h(x, y) w}
1
1 (x µ)C 1 (x µ)T
fX (x1 , x2 , · · · , xn ) = p ·e 2
(2⇡)n det(C)
Señal 1 de 10
4 600
3
500
400
1
0 300
−1
200
−2
100
−3
−4 0
0 100 200 300 400 500 600 700 800 900 1000 −4 −3 −2 −1 0 1 2 3 4
Promedio de 2 señales
4 700
3 600
2
500
1
400
300
−1
200
−2
100
−3
−4 0
0 100 200 300 400 500 600 700 800 900 1000 −4 −3 −2 −1 0 1 2 3 4
Promedio de 5 señales
4 700
3 600
2
500
1
400
300
−1
200
−2
100
−3
−4 0
0 100 200 300 400 500 600 700 800 900 1000 −4 −3 −2 −1 0 1 2 3 4
Promedio de 10 señales
4 700
3 600
2
500
1
400
300
−1
200
−2
100
−3
−4 0
0 100 200 300 400 500 600 700 800 900 1000 −4 −3 −2 −1 0 1 2 3 4
Promedio de 20 senales
4 700
3 600
2
500
1
400
300
−1
200
−2
100
−3
−4 0
0 100 200 300 400 500 600 700 800 900 1000 −4 −3 −2 −1 0 1 2 3 4
Promedio de 50 senales
4 700
3 600
2
500
1
400
300
−1
200
−2
100
−3
−4 0
0 100 200 300 400 500 600 700 800 900 1000 −4 −3 −2 −1 0 1 2 3 4
3 600
2
500
1
400
300
−1
200
−2
100
−3
−4 0
0 100 200 300 400 500 600 700 800 900 1000 −4 −3 −2 −1 0 1 2 3 4
3 600
2
500
1
400
300
−1
200
−2
100
−3
−4 0
0 100 200 300 400 500 600 700 800 900 1000 −4 −3 −2 −1 0 1 2 3 4
Señal 1 de 10
400
2
350
1.8
1.6 300
1.4
250
1.2
200
1
0.8 150
0.6
100
0.4
50
0.2
0 0
0 50 100 150 200 250 300 350 400 450 500 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Promedio de 2 señales
700
2
1.8 600
1.6
500
1.4
1.2 400
1
300
0.8
0.6 200
0.4
100
0.2
0 0
0 50 100 150 200 250 300 350 400 450 500 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Promedio de 5 señales
900
2
800
1.8
700
1.6
1.4 600
1.2 500
1
400
0.8
300
0.6
200
0.4
100
0.2
0 0
0 50 100 150 200 250 300 350 400 450 500 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Promedio de 10 señales
1200
2
1.8
1000
1.6
1.4 800
1.2
600
1
0.8
400
0.6
0.4
200
0.2
0 0
0 50 100 150 200 250 300 350 400 450 500 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Promedio de 20 señales
1200
2
1.8
1000
1.6
1.4 800
1.2
600
1
0.8
400
0.6
0.4
200
0.2
0 0
0 50 100 150 200 250 300 350 400 450 500 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Promedio de 50 señales
1200
2
1.8
1000
1.6
1.4 800
1.2
600
1
0.8
400
0.6
0.4
200
0.2
0 0
0 50 100 150 200 250 300 350 400 450 500 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
1.8
1000
1.6
1.4 800
1.2
600
1
0.8
400
0.6
0.4
200
0.2
0 0
0 50 100 150 200 250 300 350 400 450 500 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
X(t) = f (t, ✓)
Descripción estadı́stica
I Completa: 8 (t1 , t2 , · · · , tn ) 2 IRn , 8 n
Estacionariedad y cicloestacionariedad
Estacionariedad en sentido estricto: 8 (t1 , t2 , · · · , tn ), 8n, 8
Cicloestacionariedad
1 mX (t + To ) = mX (t)
2 RX (t + ⌧ + To , t + To ) = RX (t + ⌧, t), para todo t y ⌧
RX ( ⌧ ) = RX (⌧ )
|RX (⌧ )| RX (0)
Ergodicidad
Promedios para un proceso X(t) y una función g(x):
1 Promedio estadı́stico
Z 1
E[g(X(t))] = g(x) · fX(t) (x) · dx
1
" Z T
#
1 2
PX = E lı́m X 2 (t) · dt
T!1 T T
2
Z T Z T
1 2 1 2
= lı́m E[X 2 (t)] · dt = lı́m RX (t, t) · dt
T!1 T T T!1 T T
2 2
En general
RX,Y (t1 , t2 ) = RY,X (t2 , t1 )
Estacionariedad conjunta: X(t) e Y(t) son conjuntamente
estacionarios si
I Ambos son individualmente estacionarios
I RX,Y (t1 , t2 ) = RX,Y (⌧ ), con ⌧ = t1 t2
F Notación alternativa: RX,Y (t + ⌧, t) = RX,Y (⌧ )
Teorema de Wiener-Khinchin
Si para cualquier valor finito ⌧ y cualquier intervalo A, de
longitud |⌧ |, la autocorrelación del proceso aleatorio cumple
Z
RX (t + ⌧, t)dt < 1
A
SX (j!) = T F[RX (⌧ )]
entonces
eX (⌧ )]
SX (j!) = T F[R
donde Z
eX (⌧ ) = 1
R RX (t + ⌧, t) · dt
To To
PX = RX (0)
I Proceso cicloestacionario
PX = R̃X (0)
X(t) - Y(t) -
h(t)
Por tanto
Z 1
mY (t) = E X(s) · h(t s) · ds
1
Z 1
= E[X(s)] · h(t s) · ds
1
Z 1
= mX · h(t s) · ds
1
Z 1
u=t s
= mX h(u) · du
1
Correlación de salida RY (⌧ )
mY = mX · H(0)
X(t) - Y(t) -
h(t)
-
SX,Y (j!)-
H ⇤ (j!)
-
SY (j!) -
|H(j!)|2
I Procesos cicloestacionarios
N 1
j! 1 X
SX (e ) = T F R̃X [k] , R̃X [k] = RX [n + k, n]
N
n=0
I Potencia
Z (
1 ⇡
RX [0], X[n] estacionario
PX = SX (ej! ) d! =
2⇡ ⇡ R̃X [0], X[n] cicloestacionario
Estadı́sticos cruzados
mX · E[Y(t)] +mX · mY
| {z }
mY
=RX,Y (⌧ ) mX · mY
Procesos incorrelados:
Por definición : Cov(X(t + ⌧ ), Y(t)) = 0, 8⌧
Consecuencia : RX,Y (⌧ ) = mX · mY
I Si al menos uno de los procesos (incorrelados) tiene media nula
RZ (⌧ ) = RX (⌧ ) + RY (⌧ )
SZ (j!) = SX (j!) + SY (j!)
F Caso habitual de suma de señal y ruido: incorrelación, ruido con media nula
c Marcelino Lázaro, 2014 Teorı́a de la Comunicación Ruido térmico 77 / 88
Proceso gausiano
Definición: X(t) es un proceso gausiano si para todo n y
todo {t1 , t2 , · · · , tn }, las variables aleatorias {X(ti )}ni=1 tienen
una distribución conjuntamente gausiana
Propiedades de los procesos gausianos
I mX (t) y RX (t1 , t2 ), proporcionan una descripción estadı́stica
completa del proceso
F Permiten calcular vector de medias y matriz de covarianzas
- Para X(ti ), ) µi = mX (ti )
- X(ti ), X(tj ), ) Ci,j =Cov(X(ti ),X(tj ))=RX (ti , tj ) mX (ti ) · mX (tj )
I Para procesos gausianos, estacionariedad en sentido
estricto y en sentido amplio son equivalentes
I Si X(t) pasa por un sistema lineal e invariante, el proceso de
salida, Y(t) es gausiano
I Para X(t) gausiano, estacionario y de media nula, una
condición suficiente para la ergodicidad de X(t) es
Z 1
|RX (⌧ )| d⌧ < 1
1
c Marcelino Lázaro, 2014 Teorı́a de la Comunicación Ruido térmico 78 / 88
Procesos aleatorios conjuntamente gausianos
Definición: Los procesos X(t) e Y(t) son conjuntamente
gausianos, si para todo n, m, y todo {t1 , t2 , · · · , tn } y
{⌧1 , ⌧2 , · · · , ⌧m }, las variables aleatorias
Proceso blanco
Un proceso es blanco si su densidad espectral de potencia
es constante para todas las frecuencias
SX (j!) = C
SX (j!) 6
C
-
!
I Consecuencias
F Función de autocorrelación de un proceso blanco
estacionario
1
RX (⌧ ) = T F {C} = C · (⌧ )
F La potencia de un proceso blanco es infinita
Z 1 Z 1
1 1
PX = SX (j!) d! = C d! = 1
2⇡ 1 2⇡ 1
c Marcelino Lázaro, 2014 Teorı́a de la Comunicación Ruido térmico 80 / 88
Filtrado de un proceso blanco
Proceso X(t) blanco con SX (j!) = C se filtra (h(t) / H(j!))
Densidad espectral de potencia a la salida del filtro (Y(t))
SY (j!) = SX (j!) · |H(j!)|2 = C · |H(j!)|2
Ruido térmico
Densidad espectral de potencia del ruido térmico (mecánica cuántica)
h!
Sn (j!) = h!
1) 4⇡(e 2⇡kT
8
>
> h: Constante de Planck (6.6⇥10 34 Julios ⇥ segundo)
<
k: Constante de Boltzmann (1.38⇥10 23 Julios/o Kelvin)
donde
>
> T: Temperatura en grados Kelvin
:
!: Pulsación (2⇡ veces la frecuencia) en rad/s
6Sn (j!)
............................................2........6
Sn (j!)
..2............................................. .................................................
kT kT
..
...
...
...
...
...
...
...
...
... ............
............ 0,8 kT 0,8 kT
2 2
0,6 kT
2
0,6 kT
2
0,4 kT
2
0,4 kT
2
0,2 kT
2
0,2 kT
2
-2000 -1500 -1000 -500 ! 0 500 1000 1500 2000 -200 -150 -100 -50 ! 0 50 100 150 200
2⇡
(GHz) 2⇡
(GHz)
Estadı́stica gausiana
c Marcelino Lázaro, 2014 Teorı́a de la Comunicación Ruido térmico 82 / 88
Modelo de ruido térmico
Proceso aleatorio n(t) blanco, gausiano, estacionario,
ergódico
I Media nula (mn = 0)
I Función de autocorrelación
N0
Rn (⌧ ) =
· (⌧ )
2
I Densidad espectral de potencia
N0
Sn (j!) =
2
Sn (j!) 6
N0 /2
-
!
Valor de la constante N0
N0 = k · T a Watt/Hz
c Marcelino Lázaro, 2014 Teorı́a de la Comunicación Ruido térmico 83 / 88
|H(j!)| 6
.
. ..........Hmax
. ....
.. .. ...
...
... ....
... ....
.... ......
.................... ................ -
2⇡ · Beq !
c Marcelino Lázaro, 2014 Teorı́a de la Comunicación Ruido térmico 87 / 88
PZ = N0 · B PZ = N0 · B · G
F Filtro con ancho de banda equivalente Beq y ganancia Geq
PZ = N0 · Beq · Geq