8 Detección de Señales Binarias en Presencia de Ruido Blanco Gaussiano
8 Detección de Señales Binarias en Presencia de Ruido Blanco Gaussiano
8 Detección de Señales Binarias en Presencia de Ruido Blanco Gaussiano
s (t ) 0 t T para el 1 binario
si (t ) 1
s2(t ) 0 t T para el 0 binario
Por ejemplo, un 1 puede ser representado por una tensión eléctrica V que se mantiene
constante durante un tiempo T, y un 0 por una tensión –V que también se mantiene constante
por igual duración de tiempo. Esta señal puede ser transmitida directamente (transmisión en
banda base) o bien puede ser usada para modular una portadora. La señal recibida está
afectada por el ruido y por lo tanto existe una cierta probabilidad de que el receptor cometa un
error al decidir si se transmitió un 1 o un 0. La señal recibida por el receptor, r(t) se puede
representar por:
donde n(t) es el ruido blanco gaussiano aditivo, de media cero y varianza 2, que
interfiere sobre la señal que fue transmitida.
Trataremos de encontrar qué características debe tener nuestro receptor para que
pueda hacer una detección lo más fiel posible del bit que fue transmitido. El transmisor “sabe”
que fue transmitido un 1 o un 0, pero debido al efecto del ruido esas señales eléctricas que
representan al 1 y al 0 fueron deformadas y eso puede “confundir” al detector a la hora de
discernir qué señal se transmitió.
La señal transmitida s(t)(representada por uno de sus dos estados V), junto con el
ruido Gaussiano n(t), se la hace pasar por un integrador y luego sí es muestreada. Esto se
puede hacer con un amplificador operacional con una resistencia de entrada R y un capacitor
de realimentación C, como se muestra en la Figura 2. Por lo tanto, la salida del integrador
produce una señal que es la integral del pulso enviado, multiplicada por una constante
1/RC = . Al final de cada tiempo de bit T, una llave en paralelo con el capacitor se cierra para
descargarlo y así comenzar la siguiente integración desde cero. A la salida del integrador y
justo al instante de tomar la muestra, tenemos:
1
s(t ) n(t )dt 1 1
T T T
v0 (T )
0 0
s(t )dt
n(t )dt
0
(1)
1 T VT
s0 (T )
Vdt
0
(2)
o bien con el signo opuesto cuando la entrada del integrador es –V. La muestra de
tensión correspondiente al ruido es:
1 T
n0 (T )
n(t)dt
0
(3)
Cabe destacar que n(t) es un proceso aleatorio, mientras que no(T) es una variable
aleatoria. La varianza (potencia media) de n0(T) se expresa como:
N 0T
02 (4)
2 2
Además n0(T) tiene una fdp que es Gaussiana. El valor N0/2 es la densidad espectral del
ruido a la entrada del integrador. No demostraremos aquí la obtención de la última ecuación,
pero la misma surge de plantear que:
Figura 2. Esquema propuesto inicialmente para mejorar la relación señal a ruido, en el instante de
muestreo, de la señal recibida.
La salida del integrador, antes de la llave de muestreo, es v0(t) = s0(t) + n0(t). La señal
de salida s0(t) es una rampa para cada intervalo de tiempo T. Al final de la rampa, o sea, al
instante del muestreo T, la señal s0(t) tiene un valor que es +VT/ o –VT/, dependiendo de si
el bit transmitido fue 1 o 0. Mientras que el ruido, en el instante del muestreo, tiene un valor
aleatorio n0(T). Ambas situaciones se muestran en la Figura 3. Finalmente, el voltaje total al
momento del muestreo es:
v 0 (T ) s0 (T ) n0 (T ) (6)
Nótese que la relación señal a ruido aumenta con el incremento de la duración de bit T
y que depende de V2T que es la energía normalizada de la señal. Es decir que para mejorar la
relación señal a ruido se puede aumentar el tiempo de bit T (transmitiendo menos bits por
segundo, o sea, más lentamente), o bien aumentar la tensión V de la señal, o ambas cosas.
Aumentar la tensión implica transmitir con mayor potencia, con lo cual se necesitaría un
equipo transmisor más grande, fuentes de alimentación más grandes, quizás disipadores más
grandes, etc., aumentando el costo y el tamaño del sistema. La otra alternativa, aumentar el
tiempo de bit, implica transmitir a menor velocidad. Esto puede verse, por ejemplo, en una
comunicación vía modem (aunque en este caso no es banda base sino una señal modulada, el
ejemplo vale). Un modem que cumple con la norma de comunicación V.90 debería conectarse
normalmente a 56 Kbps. Sin embargo, suele ocurrir que la línea telefónica tenga un ruido
mayor al esperado y para mantener la relación señal a ruido en un nivel aceptable el modem
transmite más lentamente. Otro ejemplo lo muestran los satélites que se envían al espacio
(por ejemplo la NASA o la Agencia Espacial Europea) para tomar fotografías, hacer estudios,
mediciones, etc. en diversos planetas. Debido a la gran distancia que separa al satélite de la
antena terrestre, la relación señal a ruido es realmente muy baja como consecuencia de la
atenuación de la señal transmitida. Por lo tanto, para mantener una relación señal a ruido
razonable en el receptor la alternativa es aumentar T. La conclusión es que la velocidad de
transmisión es muy baja.
Hemos visto entonces, que, para el ejemplo anterior de transmisión binaria bipolar, el
filtro que nos mejora la relación señal a ruido es un integrador. Esto ha sido para este ejemplo.
En los párrafos siguientes trataremos de ver si es éste el mejor filtro que podemos poner o
existe uno mejor, y además generalizar la situación y tratarla no sólo para la señal bipolar que
se planteó sino para cualquier señal binaria.
z(T ) ai (T ) n0 (T ) i 1, 2 (8)
1 1 n 2
p(n0 ) exp 0 (9)
0 2 2 0
2
donde 0 es la varianza del ruido a la salida del filtro lineal, incluso después de la llave
de muestreo. Este ruido, si es sumado a las señales ai(T), se representa por funciones de
1 1z a
2
p( z | s1 ) exp 1
(10)
0 2 2 0
1 1z a
2
p( z | s2 ) exp 2
(11)
0 2 2 0
Una vez que la señal recibida r(t) es convertida en un número z(T), ya no importa más
la verdadera forma de la señal. Todas las formas de onda que sean transformadas al mismo
valor de z(T) son idénticas desde el punto de vista de la detección. Veremos más adelante que
un cierto tipo de filtro, llamado filtro adaptado, ubicado en el bloque 1, “mapea” las señales de
H1
z(T ) (12)
H2
donde H1 y H2 son las dos posibles hipótesis binarias. Elegir H1 es equivalente a decidir
que s1(t) se ha enviado, y elegir H2 es equivalente a decidir que s2(t) fue enviado. La
desigualdad en la relación anterior indica que si z(T) > entonces se elige H1. De lo contrario,
se elige H2. Para el caso de la igualdad la decisión es arbitraria y se ajusta al detector de
manera tal que elija una de las dos hipótesis al azar, como si tirara una moneda.
Esta comparación es similar a la que se hizo en el ejemplo al comienzo del texto. Allí,
una vez tomada la muestra, se la comparaba con cero: si el número era positivo se decidía por
un “1” binario, caso contrario se decidía por un “0” binario.
¿Qué valor debe tener ? Hasta ahora no hemos dicho nada acerca de qué valor debe
tener para que sea un nivel umbral apropiado. Un criterio para elegir se basa en la
minimización de la probabilidad de error. Es decir, se busca un valor tal de como para que el
detector se “equivoque” lo menos posible en su toma de decisiones. Este valor óptimo de que
hace mínima la probabilidad de error lo llamaremos 0.
Para hallar tal valor partimos de la siguiente relación, llamada test de relación de
verosimilitud:
H1
p(z | s1 ) P(s2 )
(13)
p(z | s2 ) P(s1 )
H2
donde P(s1) y P(s2) son las probabilidades a priori de que s1(t) y s2(t) sean transmitidos,
y H1 y H2 son las dos posibles hipótesis. Según la regla de la ecuación (13), para minimizar la
probabilidad de error en la estimación, debemos elegir la hipótesis H1 si la relación entre las
fdp condicionales es mayor que la relación entre las probabilidades a priori. Caso contrario se
elige la hipótesis H2. Esta ecuación (13) surge de la teoría de la probabilidad condicional.
H1
a1 a2
z(T ) 0 (14)
2
H2
Una vez más, volviendo a nuestro ejemplo del comienzo del texto, y para el caso
particular del circuito integrador que se ha usado, a1 sería VT/, a2 sería –VT/ y 0 sería 0.
Obviamente, si las probabilidades a priori no son iguales (ambas igual a 0,5), entonces
0 se desplaza hacia la izquierda o hacia la derecha.
Si en lugar de una transmisión binaria fuese una transmisión M-aria, debería haber M
funciones de densidad de probabilidad representando a las M señales. La decisión de máxima
verosimilitud debe ser hecha entonces según el máximo valor de probabilidad de todas las fdp.
Para el caso de una transmisión binaria, hay dos maneras por las cuales se puede
producir un error en la detección o estimación. Un error e va a ocurrir cuando se envíe una
señal s1(t) y el ruido del canal resulte tal que el valor z(T) sea menor que el umbral 0 (una vez
más la Figura 5 sirve de ayuda). En tal caso, la probabilidad de error o probabilidad de error de
bit se puede expresar como:
0
P(e | s1 ) P(H2 | s1 ) pz | s1 dz (16)
P(e | s2 ) P(H1 | s2 ) 0
p(z | s2 )dz (17)
2
PB P(e, s )
i 1
i (18)
o equivalentemente,
(Las ecuaciones (19) y (20) surgen del teorema de Bayes). Es decir, dado que fue
transmitido s1(t), ocurre un error si se elige la hipótesis H2; o dado que se transmitió s2(t),
ocurre un error si fue elegida la hipótesis H1. Para el caso en que las probabilidades a priori
sean iguales, o sea, P(s1) = P(s2) = ½, PB resulta ser:
PB 0
p(z | s2 )dz (23)
1 1z a
2
exp dz
2
PB (24)
2
0 0 2 0
2
Como ya sabemos, o es la varianza del ruido a la salida del filtro lineal. Haciendo
u = (z-a2)/0 tenemos:
du 1
dz 0
Por lo tanto,
dz 0 du
a1 a2
z 0 u a2 0
2
a1 a2
0u a2
2
a a2 a
u 1 2
2 0 0
a1 a2
u
2 0
u 1 u2 a a2
PB
u ( a1 a2 ) / 2 0 2
exp
2
du Q 1
2
0
(25)
En el ejemplo analizado al comienzo de este capítulo habíamos visto que a la salida del
filtro integrador, en el instante de muestreo, obteníamos los siguientes valores de a1 y a2:
VT
a1
VT
a2
2
Y habíamos visto que 0 era:
N 0T N 0T
02 2
0
2 2 2
a a2
PB Q 1
2 0
VT VT
PB Q
N 0T
2
2 2
V 2T 2 2
PB Q
N T 2 2
0
V 2T
PB Q
N /2
0
Hasta aquí hemos optimizado el valor del umbral haciendo mínima la probabilidad de
error PB y hallando una expresión para el cálculo de dicha probabilidad de error. Ahora
veremos cómo debe ser el filtro que, a la salida del muestreador, maximiza la relación señal a
ruido, y cómo puede usárselo para que maximice el argumento de la función Q(x) a fin de
que la probabilidad de error sea mínima (es decir, mínima no en función del umbral de
detección 0 sino en función del valor z(T) que entrega el filtro en el instante de muestreo).
Filtro adaptado
Un filtro adaptado es un filtro lineal, diseñado para que pueda dar, a su salida, la
máxima relación señal a ruido para una determinada forma de onda del símbolo transmitido.
Supongamos que una señal conocida s(t) más RBGA n(t) se presenta a la entrada de un filtro
lineal e invariante en el tiempo, seguido éste por un muestreador, como se muestra en la
Figura 4. En el tiempo t = T la salida z(T) del filtro consiste en una componente de señal a y
una componente de ruido n0. La varianza del ruido de salida (potencia media de ruido) se
2
denota por 0 , de manera que la relación entre la potencia de la señal y la potencia media de
ruido, (S/N)T en el tiempo t = T, a la salida del bloque 1 es:
S a2
2 (27)
N T 0
Lo que queremos hacer ahora es encontrar la característica que debe tener nuestro
filtro. Es decir, debemos hallar la función de transferencia óptima H0(f) que maximice la
ecuación (27). Podemos expresar la señal a(t) a la salida del filtro en términos de la función de
transferencia H(f) (antes de su optimización) y de la transformada de Fourier de la señal de
entrada:
a(t )
H(f )S(f )e j2ft df (28)
N0
2
02 H(f ) df (29)
2
2
Gno (f )df Gni (f ) H(f ) df
y además
N0
Gni (f )
2
Combinando las ecuaciones (28) y (29) con la (27), podemos expresar la relación señal
a ruido como sigue:
S
H(f )S(f )e j 2fT df
(30)
N T
2
N0 / 2 H(f ) df
Ahora tenemos que hallar un valor tal H(f) = H0(f) para el cual la relación señal a ruido
(S/N)T se hace máxima. Para ello usaremos la desigualdad de Schwarz que expresa lo
siguiente:
2
2 2
f1 (x)f2 (x)dx f1 (x) dx f2 (x) dx (31)
2
2 2
H(f )S(f )e j2fT df H(f ) df S(f ) df (32)
luego
S 2
2
S(f ) df (33)
N
T N 0
S 2E E
max (34)
N
T N 0 N 0 2
2
E S(f ) df (35)
h(t ) 1 kS (f )e j2fT (37)
ks(T t ) 0t T
h(t ) (38)
0 para cualquier otro valor
O sea, la respuesta impulsiva de este filtro que tiene a su salida la máxima relación
señal a ruido, es la imagen invertida de la señal s(t), retardada en un tiempo T que es la
duración de símbolo. En la Figura 6 se puede ver un ejemplo de composición de h(t).
t
z(t ) r(t ) h(t ) r( )h(t )d
0
(39)
t
z(t ) r( )sT t d
0
(40)
t
r( )sT t d
0
Finalmente, cuando actúa el muestreador que se encuentra a continuación del filtro z(t)
se convierte en
T
z(T ) r( )s( )d
0
(41)
La operación descripta por la ecuación (41) es la integración del producto entre la señal
recibida r(t) y una réplica de la forma de onda transmitida s(t) sobre el intervalo de un tiempo
de símbolo T. Esta operación se la conoce como correlación entre r(t) y s(t).
Anteriormente, al determinar el umbral de decisión óptimo hemos visto que para tal
caso la probabilidad de error es PB = Q[(a1-a2)/20]. El hecho de encontrar un umbral óptimo
no es suficiente para optimizar el proceso de detección. Para minimizar PB necesitamos elegir
un filtro tal que maximice el argumento de Q(x). Es decir, un filtro que maximice (a1-a2)/20, o
equivalentemente que maximice
a1 a2
2
(42)
02
S a a 2
2E d
1 22 (43)
N T 0 N0
donde N0/2 es la densidad espectral de potencia del ruido a la entrada del filtro y Ed es
la energía de la diferencia entre las señales s1(t) y s2(t) a la entrada del filtro:
T
s (t) s (t) dt
2
Ed 1 2 (44)
0
Ed
PB Q (45)
2N0
Por lo tanto, del análisis hecho en este apartado, concluimos que el filtro que se usa
previo al muestreo de la señal en un esquema de detección binaria, es un filtro adaptado a la
diferencia entre las señales binarias s1(t) y s2(t) que son transmitidas, o su equivalente
implementado con un correlador. Ambos casos se muestran en la Figura 7.
Es importante notar que, si el filtro está adaptado a s1(t) – s2(t), la máxima SNR a la
salida se produce sólo si la señal de entrada es también s1(t) – s2(t). Y en verdad, en un
sistema de comunicación, el receptor (filtro adaptado) recibe s1(t) o s2(t) y no s1(t) – s2(t), por
lo tanto la SNR no sería máxima. Sin embargo, al tratarse de señales en banda base que son
de forma rectangular, la respuesta impulsiva del filtro adaptado a s1(t) – s2(t) tiene una forma
rectangular que es de amplitud proporcional a s1(t) y s2(t), que son las señales que van a estar
pasando por él. Por lo tanto, el filtro está adaptado pero con un factor de escala que afecta por
igual tanto a la señal de entrada como a la señal de ruido y entonces se sigue cumpliendo la
máxima SNR a la salida. De todos modos, en otro capítulo veremos que el filtro adaptado
también puede usarse en la detección de señales pasabanda que no tienen forma rectangular.
Señalización unipolar.
donde A > 0 es la amplitud de la señal s1(t). Supongamos que esta señal unipolar como
la que se muestra en la Figura 8(a), más ruido blanco Gaussiano, se presenta a la entrada de
un correlador o de un filtro adaptado con tiempo de muestreo t = T. El correlador para
detectar este tipo de señal se muestra en la Figura 8(b). Este correlador multiplica e integra la
señal r(t) que llega, con la diferencia de las señales prototipo, [s1(t) – s2(t)] = A, y luego del
tiempo T compara el valor obtenido z(T) con el valor umbral 0. El valor del umbral óptimo en
este caso es 0 = (a1 + a2)/2 = (½)A2T. Si la salida del correlador es mayor que 0 entonces
estima que se recibió s1(t); caso contrario estima que se recibió s2(t).
Ed
Q E b
2
PB Q Q A T (47)
2N 0 2N N0
0
donde la energía media por bit es Eb = A2T/2 (es la energía de s1(t) más la energía de
s2(t), dividido 2).
Señalización bipolar.
z(T ) z1(T ) z2 (T )
Para señales antipodales resulta ser a1 = -a2, por lo tanto el umbral de decisión óptimo
es 0 = 0. De esta manera, si el test estadístico z(T) es positivo se decide por s1(t); caso
contrario se decide por s2(t).
2 A 2T
PB Q Q 2E b
(49)
N N0
0
Figura 10. Desempeño de tasa de error de bit para señalización unipolar y bipolar.
El filtro adaptado presenta ciertas propiedades interesantes, que pueden resultar útiles
de aplicar (quizás para comprender mejor algunos otros temas) y que describiremos a
continuación.
So (f ) Ho (f ) S(f )
S * (f ) S(f ) exp( j2fT ) (50)
2
S(f ) exp( j2fT )
so (t) Rx t T (51)
so (T ) Rx (0) E (52)
N0
2
02 H0 (f ) df
2
N0 N0
2
02 S(f ) df E (53)
2 2
E2 2E
SNRO,max (54)
N0 E / 2 N0
Aquí se ve que la SNR a la salida del filtro no depende de la forma de onda de la señal
elegida sino de su energía.
Resumen
Las señales binarias que salen de un transmisor deben ser recibidas y detectadas por el
receptor. Si el canal de comunicación no tiene ruido, esta tarea es bastante sencilla y en
general consiste en tomar una muestra de tensión cada T segundos (tiempo de bit) y
compararla contra un valor umbral para determinar si se transmitió un 1 o un 0. Si el canal
tiene ruido, va a afectar a la forma de onda de la señal transmitida y puede ocurrir que justo
en el instante de muestreo el valor instantáneo esté del “lado contrario” del umbral de
comparación y entonces el receptor va a identificar erróneamente al símbolo transmitido,
confundiendo un 1 con un 0 o viceversa. Una manera más realista de hacer la detección del
símbolo es observándolo durante todo el tiempo que dura y no sólo un instante. Esta solución
conduce al uso de un filtro, llamado filtro adaptado, seguido de una llave que se cierra cada T
segundos tomando una muestra que resulta ser proporcional al valor promedio del símbolo
durante la duración del mismo, mucho más representativo que el valor de una sola muestra.
Este valor es comparado con un umbral de decisión contra el cual el receptor decide si el