2 - Detection of Digital Signals
2 - Detection of Digital Signals
y(Tb) is the detector output at t=Tb( bit duration for both "0" and "1"
signals ). A comparator with threshold voltage Vth is used to decide if the
y(Tb)
x(t)+n(t) Detector
(binary) + data
Vth
The problem is to find the best choice of Vth such that the overall net error
probability(average error prob) of "1" & "0" is minimum. This optimum
Vth depends on type of the binary signals and the PDF of n(t).
Unipolar binary signals: First we discuss the simplest case if x(t)=0 volt
for logic"0" and x(t)=+A volt for logic"1" over the bit duration Tb.
2
+A
x(t) 0 0
1 1
0
Tb
p(1T)
p(0T)
Note that for noiseless case (y(t)=x(t)) then the PDF of y(t) is discrete
similar to the PDF of x(t).
When n(t)≠0, the PDF of y(t) will be the same as the PDF of n(t) (f(n)) but
d.c. shifted by 0 volt or +A volts corresponding to x(t). Let f1(y) be the
PDF of y when "1" is transmitted and fo(y) be the PDF of y when "0" is
transmitted.
0T 0R
p(0R/1T)
p(1R/0T)
1T
1R
y2
1
pe fo( y)dy
Vth A/ 2 2
e 2 2
dy . Note that this definite integral can
x
Vth 0.912
p(1R / 0T )
Vth
fo( y )dy Q(
) Q(
0.4
) Q(2.28) 0.01125
and
A Vth 2 0.912
Vth
p(0 R / 1T )
f 1( y)dy Q(
) Q(
0.4
) Q(2.72) 3.28 *103
Homework: Show that for bipolar signals A over AWGN channel with
variance 2 and q=p(0T)/p(1T), then:
2 ln q
Vth | optimum and if q=1, then Vth=0, and pe=Q(A/)
2A
1 0.5
p(1R / 0T ) 0.75(1 y )dy 0.156 p(0 R / 1T ) 0.75[1 ( y 1) 2 ]dy
2
0.5 0
Ternary signals: Suppose that the ternary voltage levels –A, 0, +A volts
are used to represent logic states 0,1,2 respectively. The decision rules will
be:
p (0 ) f 1(Vth 2)
Vth2 can be found from p (1 ) fo(Vth 2)
T
T
9
1-pe
X1 Y1
X2 pe
e
pe
1-2pe
X2 Y2
X3
pe
pe
X3 Y3
1-pe
Ex: Equiprobable ternary signals -2A, 0, +3A volts are affected by zero
mean AWGN with variance 2=0.16A2. Find optimum threshold settings
and the net error prob.
Solution: Since equiprobable, then: Vth1=[0-2A]/2=-A volts,
Vth2=[0+3A]/2=1.5A. Also =0.4A, then:
pe1=p(0R/1T)=p(1R/0T)=Q[d1/(2)]=Q[2A/(2*0.4*A)]=Q(2.5)=6.2*10-3
pe2=p(2R/1T)=p(1R/2T)=Q[d2/(2)]=Q[3A/(2*0.4*A)]=Q(3.75)=8.9*10-5
and the net error prob=(1/3)*pe1+(1/3)*pe1+(1/3)*pe2+(1/3)*pe2
=(2/3)*pe1+(2/3)*pe24.19*10-3
11
errors
5*4
[1]p(a block with 2 errors)= C25 pe2 (1 pe) 3 * (0.02) 2 (0.98) 3 0.0037
2
11
Matched filter
Once we studied how to choose the best threshold to minimize the net
error prob, next, we consider the detection problem, i.e. the detection of
the signal x(t) embedded in AWGN noise n(t) (only AWGN case is
considered here). The output of this detector is y(Tb) which was the signal
used at the decision block. The device that carries out the function of
detection is called a Matched Filter . This matched filter is simply a linear
system with impulse response h(t) (or frequency response H()) that acts
| y (Tb ) | 2
as a filter. The job of this filter is to maximize the ratio n 2 (t ) . This
o
ratio is in fact the ratio of the signal power to noise power at the output of
this filter at the instants of sampling at Tb(bit duration). The input
Matched filter
x(t)+ni(t) H() y(t)+no(t)
h(t)
2
1
And no2(t)=o/p noise power= 2 Gn ( ) | H ( ) | d where Gn() is the
2
input noise power spectral density which is flat for AWGN and is equal to
o/2 (two-sided spectrum), then:
12
1
no (t)= 2 (o / 2) | H ( ) | d
2
2
The ratio to be maximized will be:
1
| y (Tb ) | 2
|
2
H ( ) X ( )e jwTbd |2
n o2 (t ) 1 then using Shwartz inequality
2
( o / 2) | H ( ) | 2
d
f ( x) f ( x)dx | ( | f1 ( x) | dx ) ( | f 2 ( x) | 2 dx ) where
2 2
that states that: | 1 2
*
equal sign holds if f1(x)=k f2 (x) (*=conjugate), k is any constant that
may represent the gain, say k=1.
Now if f1(x)=|H()|, f2(x)=X()ejwTb, then applying this inequality on the
numerator of above ratio get:
1
( | H ( ) | d ) ( | X ( )e jwTb |2 d )
2
| y (Tb ) |2 4
2
n o2 (t ) (o / 2)
2
| H ( ) | 2
d
| y (Tb ) |2 1
d
2
| X ( ) |
o
jwTb 2
And since |e | =1, then: n 2 (t ) and since:
o
| y (Tb ) |2 E
| X () | d E =energy of x(t), then: n 2o (t ) o or:
2
| y (Tb ) |2 E
[ 2 ] max
n o (t ) o note that both E and o have the same units
Watts/(rad/sec). But if expressed in Watt/Hz, then:
| y (Tb ) |2 E
[ 2 ] max
n o (t ) (o / 2) =Max SNR at matched filter output.
o=one-sided AWGN spectral density(Watt/Hz).
To find the condition for maximization, recall that we used equal sign in
Schwartz inequality if f1(x)=f*2(x) , or H()=(X()ejwTb)*=X*()e-jwTb, or
in time domain:
h(t)=x(Tb-t).
13
Hence the impulse response of a matched filter matched to the signal x(t)
is the negative time of x(t) shifted by Tb.
Ex: Find the impulse response and the max SNR at the o/p of a matched
filter used to detect the signal so(t) shown embedded in AWGN with one-
sided spectral density of o.
h(t) so(t)
V V
Tb t Tb t
t
0 0
2V 2Tb
(SNR)max= 3 .
o
Practical Implementation of the Matched filter:
We know that h(t)=x(Tb-t), and since y(t)=h(t)x(t) (=convolution)
t Tb
y (t ) x( )h(t )d and at t=Tb y (Tb ) x( )h(Tb )d
0 0
Tb Tb
(t )dt =Energy of x(t) over Tb
2
Or y (Tb ) x ( ) x (Tb Tb ) d x
0 0
Above equation indicates that if the i/p x(t) is noise free then y(Tb) is E,
but in general if the i/p is x(t)+n(t), then y(Tb) will be a Gaussian random
variable (n(t) is Gaussian) with mean E and variance that is linearly
proportional to o. The practical implementation can be deduced as the
product of the input x(t)+n(t) and x(t) then integrate the result over t=0 to
t=Tb. The figure below shows the practical implementation of the matched
14
filter.
Tb t=Tb
x(t)+n(t)
( )dt
0 Y(Tb)
x(t)
Matched
filter yo(Tb)
ho(t)
s(t)+n(t)
data
Matched
filter y1(Tb)
h1(t)
Tb
Recall that for equiprobable case p(0T)=p(1T)=0.5, and for Gaussian case,
then:
Tb
[s o (t ) s1 (t )]2 dt
Eo E1 2 E01
pe Q( 0
) = Q( )
2 o 2 o
Ex: Matched filter is used to detect the two equiprobable signals shown in
AWGN. Find the error prob in terms of average signal power, bit rate, and
noise spectral density
o.
s1(t)
so(t) Tb
V
t
-V
Tb t
(V )(V ) 2
Tb
V2
Eo1 2
t dt Tb
0
T b 3
2-BPSK signals(Binary Phase Shift Keying): Here the two signals are
so(t)=Acost, s1(t)=-Acost over bit duration Tb. Again, using the same
general equation, then:
Eo=E1=[A2Tb]/2 since both so(t) and s1(t) are two sinusoids, also:
Tb
0
Also S=[0.5Eo+0.5E1]/Tb=A2/2=average signal power, then:
f1 fc f2
fd=f
A2 sin d Tb
Tb
A2
2 0
E01 [cos t cos 2 t ] dt
d
d c
2
A 2Tb sin d Tb
Or E o1 2 d Tb which is a sinc function as shown:
sin d Tb sin d Tb
[ A2Tb (1 ( )] [ S (1 ( )]
d Tb d Tb
pe Q( ) Q( )
2 o o Rb
S
Where S=[0.5Eo+0.5E1]/Tb=0.5A , and if Rb , then:
2
o
sin d Tb
pe Q( (1 ) which is a general formula used for matched
d Tb
filter detection of FSK signals.
Noncoherent detection:
Previous equations for error prob of digital carrier systems using matched
filter under AWGN can not be used for noncoherent detection(not
matched filter). In fact, exact derivation for noncoherent case is
mathematically complicated. Approximate results are given without
derivations:
1-noncoherent ASK(OOK):
pe=0.5[e-/2+Q()] 0.5 e-/2 for >>1, since here e-/2 >>Q()
2-noncoherent FSK:
pe 0.5 e-/2 for >>1(the same as OOK)
3-noncoherent PSK=DPSK:
pe 0.5 e- for >>1 (better than OOK or FSK)
19
Which is worst than BPSK for the same .( remember that both BPSK and
FSK have the same average power of A2 /2, but the performance of BPSK
is better at the expense of the requirement of carrier recovery at the
receiver).