Efecto del Ruido en Sistemas de Comunicación
Efecto del Ruido en Sistemas de Comunicación
Como se dijo previamente, el ruido constituye uno de los principales problemas a resolver
en el diseño de un sistema de comunicaciones. Se analizará el caso para señales binarias
y M-arias, o multinivel. En esta primera sección se analiza el efecto del ruido como si la
interferencia intersimbolica no existiera. Luego se considerará la ISI anulando el efecto del
ruido.
y(tk)
x(t)
Filtro LPF Muestra xc(t)
H(f) y retención + Umbral
- V
Sincronización: D V
Figura 2.1.
La señal mas el ruido ingresan al receptor. La primera etapa consta de un filtro pasabajos
que elimina parte del ruido entrante sin producir interferencia intersimbolica. Estudiaremos
el efecto del ruido suponiendo que la ISI provocada es nula. El valor detectado por el
receptor luego de la operación de muestreo y retención es:
14
H1 : Es la hipótesis de que ha sido enviado un uno ‘1’: ak = A; Y = A + n.
y(t) 1 0 1 0 0
A
tk Tb
y(tk) A
xc(t) A
0
Tb/2 tk Tb+Tb/2
Figura 2.2.
15
pY(Y/H0) pY(Y/H1)
0 V A
Pe1 Pe0
Figura 2.3.
En la figura 2.3 puede verse que las regiones sombreadas corresponden a los valores
considerados dentro del error para cada hipótesis. Así, se supone que sucedió la
hipótesis H 0 , cuando Y < V , y que sucedió la hipótesis H1 cuando Y > V . Las
probabilidades de error están asociadas con las áreas bajo las funciones densidad de
probabilidad están representadas en la figura 2.3:
∞
Pe0 = P(Y>V / H 0 ) = ∫ V
p y (y/H 0 )dy (2.4)
V
Pe1 = P(Y<V / H1 ) = ∫
−∞
p y (y/H1 )dy (2.5)
donde:
P0 = P(H 0 ) ; P1 = P(H1 )
16
Si se supone que los unos 1’s y los ceros 0’s son igualmente probables:
1
P0 = P1 = (2.9)
2
Luego:
1
Pe = (Pe 0 +Pe1 ) (2.10)
2
Esto significa que si los símbolos son igualmente probables entonces, el umbral debe
estar ubicado en el punto medio que separa los símbolos, para el caso de la figura 2.3,
Vopt = A/2.
α4
V Vopt
α1 α2 α3
Figura 2.4.
Si V=Vopt Si V≠Vopt
----------- -----------
Como se analiza en la figura 2.4, la probabilidad de error es mínima para V = Vopt = A/2.
17
La función densidad de probabilidad del ruido corresponde al caso del ruido Gaussiano,
que se caracteriza por poseer valor medio cero, y dispersión σ 2 . La función densidad de
probabilidad (fdp) esta expresada en la ecuación:
y2
1 -
2σ 2
pN (y) = .e (2.12)
2.π .σ 2
Para el caso general en que la fdp se encuentra desplazada a un valor medio m , distinto
de cero:
(y-m)2
1 -
2σ2
pN (y) = .e (2.13)
2
2.π .σ
pN(y)
Q(k)
m m+kσ
Figura 2.5.
(y-m)2
1 ∞ −
P(Y>m+kσ) =
2.π .σ 2 ∫
m + kσσ
e 2 .σ 2
.dx (2.14)
( λ )2
1 ∞ −
Q(k) =
2.π
∫ k
e 2.
.dλ (2.15)
Siendo:
y-m
λ= (2.16)
σ
18
y2
∞ 1 ∞ −
Pe0 = ∫V
pN (y)dy =
2πσ 2
∫
V
e 2 .σ 2
dy = Q(V/σ) (2.17)
e igualmente:
(y-A)2
V 1 V −
Pe1 = ∫-∞
pN (y-A)dy =
2.π .σ 2 ∫ −∞
e 2 .σ 2
dy = Q((A-V)/σ) (2.18)
1 |t| < τ / 2
p(t) = (2.20)
0 |t| > τ / 2
Siendo τ ≤ D .
1 1 1 T/ 2 k=N/ 2
( ∑ ak .p(t-kD))2 dt = ∑ ak2 p 2 (t-kD)dt = ∑a
T/ 2 T/ 2
∫ ∫ ∫
2
SR = k p 2 (t-kD)dt
T −T/ 2
k T −T/ 2
k T −T/ 2
k=-N/ 2
1 D/ 2 2 2 N 0 D/ 2 2 2 N1 D/ 2 2 2
SR = ∑ ∫ ∫ N.D ∫−D/ 2
ak p (t)dt = a0 p (t)dt+ a1 p (t)dt
k ND N.D
− D/ 2 − D/ 2
1 τ/ 2 2 2 1 τ/ 2
SR = P0 ∫
D −τ/ 2
a0 p (t)dt+ P1 ∫ a12 p 2 (t)dt
D −τ/ 2
(2.21)
19
SR = A2/4 Polar NRZ
2S R Unipolar
A= (2.22)
4S R pola r
2
A A2 (1 / 2 )(S/N)R Unipolar
= = (2.23)
2σ 4N R (S/N)R Polar
N R = η B ≥ η rb / 2 (2.24)
El filtro pasabajos que se encuentra a la entrada del receptor tiene por objeto limitar el
ruido entrante al sistema. Como se ha explicado, el filtrado excesivo podría provocar ISI,
por lo que existirá un compromiso entre el filtrado de ruido y la interferencia que este
provoca. Si los pulsos son limitados en el tiempo, y el ruido presente es blanco y
Gaussiano, el filtro de características optimas en términos de ruido e ISI es el filtro
adaptado.
El filtro adaptado maximiza la relación señal-ruido, siendo esta en realidad, la relación que
existe entre el valor muestreado en el instante optimo y la potencia de ruido que ingresa a
través del filtro.
En el caso analizado, el pulso entrante tiene forma conocida p( t ) .
este pulso es una señal de energía que tiene una transformada de Fourier:
con energía:
∞
E R = AR2 ∫ | P(f) | 2 df (2.27)
−∞
20
El filtro óptimo debe reducir el ruido y maximizar el valor de la amplitud a ser detectada en
el instante óptimo t = t 0 + t d .
Gn(f)
xR(t) yD(t)
H(f)
AR+nD
t0+td
Figura 2.6.
La amplitud del pulso a la salida del filtro puede ser obtenida a partir de la función
transferencia del mismo y el espectro del pulso recibido de la siguiente forma:
∞ ∞
A= ∫ -∞
H(f).AR P(f).e -jωt0 .e jω( t0+td ) df = AR ∫ H(f).P(f).e jωt d df
-∞
(2.29)
∞
AR2 ∫ H(f).P(f).e jωt d df 2
2
A -∞
= ∞
(2.31)
σ
∫
−∞
| H(f) |2 Gn (f)df
21
∞ 2
∫−∞
V(f).W*(f)df ∞
≤ ∫ W(f) df
2
∞
(2.32)
∫
2 −∞
V(f) df
−∞
AR .P(f).e jωtd
W*(f) = (2.34)
Gn (f)
KW(f)
V(f) = (2.35)
AR
entonces, se obtiene el máximo valor
2 2
A ∞ P(f)
= AR .∫
2
df (2.36)
σ max − ∞ G (f)
n
KP*(f).e -jωtd
H opt (f) = (2.37)
Gn (f)
de la expresión se observa que el filtro enfatiza las frecuencias para las que la señal tiene
un valor alto, y de-enfatiza las frecuencias para las que el ruido es alto.
2
A 2.AR2 ∞ 2E R
∫
2
= . P(f) df = (2.38)
σ max η −∞ η
2K 2K ∞
η ∫−∞
hopt (t)= F -1 [ .P*(f).e -jωt d ] = . P*(f).e − jω(td −t) df (2.39)
η
hopt (t)=
2K ∞
η ∫−∞
[
. P(f).e jω( td −t) df ] = 2ηK .p(t
*
d − t)* =
2K
η
.p(t d − t) (2.40)
22
2.2.2 Filtro Adaptado para un pulso de forma definida
p(t) , x(t) ak
kD-τ/2 kD kD+τ/2
Figura 2.7.
∞ 2K 2K +∞
ak = ak ∫ .P*(f)P(f)e jωtd .e -jωtd df ak = ak ∫ p 2 (t)dt (2.44)
−∞ η η −∞
η
K= + ∞ 2 (2.45)
∫ p (t)dt
2
−∞
1
h(t) = .p(t d -t) (2.46)
τ eq
donde:
∞
τ eq = ∫−∞
p 2 (t)dt (2.47)
23
ak
τeq
Figura 2.8.
El mínimo valor de t d es τ / 2 .
Las formas de onda típicas en la recepción de un pulso cuadrado se ven en la figura 2.9.
x(t)
t0 t0+τeq/2 t0+τeq
h(t)
1/τeq
τeq
ak
tk-τeq tk tk+τeq
Figura 2.9.
2
∞ η ∞ η ∞ p(t) η
∫ ∫−∞ h(t) 2 dt = ∫
2 2
NR = σ = 2
H(f) df = dt (2.49)
−∞ 2 −∞ τ eq2 2
24
η
N R =σ 2 = (2.50)
2τ eq
Este es el valor mínimo de potencia de ruido que puede ingresar por el filtro. El valor
mínimo siempre se verifica para el filtro adaptado. En general se puede decir que la
potencia de ruido entrante a un filtro dado será:
η
N R =σ 2 ≥ (2.51)
2τ eq
SR
Eb = Energía promedio por bit (2.52)
rb
SR E
γb = = b Relación de la energía por bit a la potencia de ruido (2.53)
ηηr b η
La energía promedio por bit para los valores discretos de la secuencia expresada en
ecuación (1.1) puede ser calculada como:
∞ ∞
E b = E [a k2 ∫−∞
p 2 (t-kD)dt ] = E [a k2 ∫
−∞
p 2 (t)dt ] = ak2 .τ eq (2.54)
A2
2
E[a k ] = ∑ ak2 .P(ak ) = (1 / 2 )A 2 + (1 / 2 )0 2 =
k 2
Unipolar
a k2 = (2.55)
E[a 2 ] = A2
k ∑k k k
a 2
.P(a ) = (1 / 2 )(A/ 2 ) 2
+ (1 / 2 )( − A/ 2 ) 2
=
4
Polar
2 a k2 2 E b 2τ eq E
2
= = b = γ b Unipolar
4σ 4 ητ eq η η
2
A
= 2 (2.56)
2σ 4 a k = E b 2τ eq = 2 E b = 2γ Polar
4σ 2 ητ eq η η
b
25
Pero como se ha evaluado en la sección anterior, la probabilidad de error depende del
termino ( A / 2 σ ) 2 , con lo que para los formatos Unipolar y Polar en señalización binaria
se tiene:
Q( γ b ) Unipolar
Pe = (2.57)
Q( 2γ b ) Polar
p(t)
(a)
Tb
h(t) = p(Tb-t)
(b)
Tb
Figura 2.10.
La forma del pulso es revertida y desplazada para conseguir la transferencia del filtro
adaptado h( t ) , como se ve en la figura 2.10(b). Cuando un pulso de la forma p( t )
presentado en la figura 2.10 es recibido por un filtro de transferencia h( t ) se produce a la
salida una forma de onda que es la integración del producto de p( t ) por h( t ) , revertido
en el tiempo. La situación para τ = 0 y τ > 0 se ve en la figura 2.11:
26
h(-t) p(t)
Tb
h(t1-t) p(t)
Tb
figura 2.11.
Para el caso en que el ruido es blanco, el filtro adaptado puede configurarse empleando la
correlación de la entrada con una forma de onda replica de la señal de entrada.
p(t) y(tk)
Integrador
tk
∫t k −τ eq
p(t)s(t)dt
s(t)
Figura 2.12.
Para verificar que esta expresión es equivalente a la de la operación del filtro adaptado
puede verse que :
tk
y(t k ) = p(t k )*h(t k ) = ∫ p(λ )h(t k − λ )dλ (2.59)
−∞
pero:
1
p(t k -t) 0 ≤ t ≤ τ eq
h(t) = τ eq (2.60)
0 otro valor de t
27
luego:
tk tk
y(t k ) = ∫
t k −τ eq
p(λ)s(t k -(t k -λ ) ).dλ = ∫
t k −τ eq
p(λ )s(λ)dλ (2.61)
Una implementación práctica del filtro adaptado para pulsos cuadrados lo constituye el
filtro de integración y descarga.
La señalización binaria provee la mayor aislación contra el ruido debido a que todo el
rango dinámico de tensiones es utilizado solo por dos niveles. El empleo de la
señalización M-aria o multinivel, tiene por objeto aumentar la velocidad de señalización
manteniendo el ancho de banda constante, pero dicho objetivo es logrado a expensas de
la performance de error, o bien para conservar la misma probabilidad de error, a expensas
de mayor potencia.
Se calculará la probabilidad de error en presencia de ruido blanco y Gaussiano, de valor
medio cero y potencia de ruido σ 2 .
1
Pe = .(Pe0 +Pe1+...+PeM-1 ) (2.63)
M
El esquema de funciones densidad de probabilidad para una señal M-aria de por ejemplo
4 niveles se presenta en la figura 2.13:
Figura 2.13.
Los umbrales de decisión se determinan suponiendo que todos los valores son
igualmente probables de forma que se encuentran ubicados en el punto medio de los
valores medios de cada función densidad de probabilidad. Los umbrales se encuentran en
−A, 0 y +A.
Los símbolos no tienen la misma probabilidad de error. Para los símbolos de los extremos
izquierdo y derecho:
28
A
Pe0 = Pe 3 = Q (2.64)
2σ
A
Pe1 = Pe 2 = 2Q (2.65)
2σ
1 A A 3 A
Pe=
4 2Q 2 σ + 2 x 2Q 2 σ = 2 Q 2 σ (2.66)
1 A A 1 A
Pe=
M 2Q 2 σ + 2(M − 2 )Q 2 σ = 2 1- M Q 2 σ (2.68)
[ ∞
−∞
]
E M = E ∫ a k2 p 2 (t)dt = a k2 ∫
−∞
∞
p 2 (t)dt = a k2 τ eq (2.69)
2
1 M/ 2
A M 2 −1 2
a = ∑ a P(a k ) = 2
2
k
2
k ∑ 2
( 2 i-1 ) = A (2.70)
k M i =1 2 12
12 E M
A2 = (2.71)
τ eq (M 2 − 1 )
Por analogía con la expresión:
SR =E b r b (2.72)
Se tendrá:
S R =E M r (2.73)
Luego:
29
2
A 3E M 3 SR
= 2 2
= 2
(2.74)
2σ (M -1 )τ eq σ (M -1 )τ eq r N R
2
A 6 SR
≤ 2
(2.75)
2σ (M - 1) ηr
η
NR =
2τ eq
Pe
Pbe ≅ (2.76)
log 2 M
2(M - 1) A
Pbe ≅ Q (2.77)
M log 2 M 2 σ
Suponga que se tiene un canal con velocidad de señalización fija e igual a 3000 baud, y
una relación señal ruido fija en 400 veces, es decir aproximadamente 20 dB. Considérese
filtro adaptado y uso de pulsos rectangulares NRZ, con lo que r .τ eq = 1 .
2
A 3 S 6log M
= 2 = 2 2 γb
2σ M −1 N R M −1
La siguiente tabla muestra los valores calculados para r = 3000 baud y (S / N )R = 400 ;
30
M rb Kbps A / 2σ Pbe η B ,bits/Hz
-89
2 3 20.0 3.10 1
4 6 8.9 1.10-19 2
8 9 4.4 4.10-6 3
16 12 2.2 7.10-3 4
32 15 1.1 6.10-2 5
M r (kbaud) γ b
2 9.00 10
4 6.00 24
8 3.00 67
16 2.25 200
32 1.80 620
Donde puede verse que para aumentar el numero de niveles con la probabilidad error
binaria y velocidad binaria constante, se precisa aumentar la potencia 60 veces si M
varia de M = 2 a M = 32 .
31