Chapter 7 - Optimal Receiver Design PDF
Chapter 7 - Optimal Receiver Design PDF
Kyung Hee
Modulation
University
z bandwidth efficient
z energy efficient
Digital Communications
Spring 2006
Modulation
z bandwidth efficient
z energy efficient
Kyung Hee
Noise Model
University
N
The noise n(t) has autocorrelation φnn ( τ) = 0 δ( τ) and
2
power spectral density Φnn ( f ) = N0 2
n( t )
Young S. Kim Wireless Communication Systems Lab. - 4 -
Kyung Hee
Problem Statement
University
Noise Model
N
The noise n(t) has autocorrelation φnn ( τ) = 0 δ( τ) and
2
power spectral density Φnn ( f ) = N0 2
n( t )
Signal Space Representation
Kyung Hee
Signal Space Representation (continued)
University
where rk = sm, k + nk
where rk = sm, k + nk
The Orthogonal Noise: n ′( t )
T T ⎛ K ⎞
∫ m
s ( t ) n ′ ( t ) dt = ∫ m
s ( t ) ⎜ n( t ) − ∑ n f ( t )
k k ⎟ dt
0 0 ⎝ k =1 ⎠
T K ⎛ K ⎞
= ∫ ∑ sm,k f k (t )⎜ n( t ) − ∑ nk f k (t )⎟ dt
0 k =1 ⎝ k =1 ⎠
K T K T
= ∑ sm, k ∫ f k (t )n(t )dt − ∑ sm,k nk ∫ f k2 (t )dt
k =1 0 k =1 0
K K
= ∑ sm,k nk − ∑ sm,k nk = 0
k =1 k =1
n
Young S. Kim Wireless Communication Systems Lab. - 8 -
Kyung Hee
The Orthogonal Noise: n ′( t )
University
T T ⎛ K ⎞
∫ m
s ( t ) n ′ ( t ) dt = ∫ m
s ( t ) ⎜ n( t ) − ∑ n f ( t )
k k ⎟ dt
0 0 ⎝ k =1 ⎠
T K ⎛ K ⎞
= ∫ ∑ sm,k f k (t )⎜ n( t ) − ∑ nk f k (t )⎟ dt
0 k =1 ⎝ k =1 ⎠
K T K T
= ∑ sm, k ∫ f k (t )n(t )dt − ∑ sm,k nk ∫ f k2 (t )dt
k =1 0 k =1 0
K K
= ∑ sm,k nk − ∑ sm,k nk = 0
k =1 k =1
n
MAP (Maximum a posteriori probability)
Decision Rule
Suppose that signal vectors {s1 , … , s M } are transmitted
with probabilities { p1,… , pm } respectively, and the signal
vector r is received
Kyung Hee
Maximum Likelihood (ML) Decision Rule
University
Kyung Hee
The Noise p.d.f
University
⎡T T ⎤
E [ni ⋅ nk ] = E ⎢ ∫ n( t ) fi ( t )dt ⋅ ∫ n( s) f k ( s)ds⎥
⎣0 0 ⎦
⎡T T ⎤ TT
= E ⎢ ∫ ∫ n(t )n( s) fi ( t ) f k ( s)dsdt ⎥ = ∫ ∫ E [n( t )n( s)] f i ( t ) f k ( s)dsdt
⎣0 0 ⎦ 00
TT TT N
= ∫ ∫ φnn (t − s) f i (t ) f k ( s)dsdt = ∫ ∫ 0 δ( t − s) f ( t ) f ( s)dsdt
i k
00 00 2
0 f ( t ) f ( t )dt = ⎧ N 0 2 , i = k
TN
=∫ i k ⎨
0 2 ⎩ 0, i≠k
⎡T T ⎤
E [ni ⋅ nk ] = E ⎢ ∫ n( t ) fi ( t )dt ⋅ ∫ n( s) f k ( s)ds⎥
⎣0 0 ⎦
⎡T T ⎤ TT
= E ⎢ ∫ ∫ n(t )n( s) fi ( t ) f k ( s)dsdt ⎥ = ∫ ∫ E [n( t )n( s)] f i ( t ) f k ( s)dsdt
⎣0 0 ⎦ 00
TT TT N
= ∫ ∫ φnn (t − s) f i (t ) f k ( s)dsdt = ∫ ∫ 0 δ( t − s) f ( t ) f ( s)dsdt
i k
00 00 2
0 f ( t ) f ( t )dt = ⎧ N 0 2 , i = k
TN
=∫ i k ⎨
0 2 ⎩ 0, i≠k
The Noise p.d.f (continued)
.
[ ]
Since E nk2 = N0 2, each noise component has a variance of N0 2
( )
K 1
=∏ exp − nk 2 N0
k =1 πN0
−K 2 ⎛ K 2 ⎞
= ( πN0) exp⎜ − ∑nk N0⎟
⎝ k =1 ⎠
Kyung Hee
Conditional pdf of Received Signal
University
⎛ K ⎞
p(r sm ) = ( πN 0 ) −K 2
(
exp⎜ − ∑ rk − sm, k
⎝ k =1
2
N0⎟
⎠
)
.
[ ]
Since E nk2 = N0 2, each noise component has a variance of N0 2
( )
K 1
=∏ exp − nk 2 N0
k =1 πN0
−K 2 ⎛ K 2 ⎞
= ( πN0) exp⎜ − ∑nk N0⎟
⎝ k =1 ⎠
⎛ K ⎞
p(r sm ) = ( πN 0 ) −K 2
(
exp⎜ − ∑ rk − sm, k
⎝ k =1
2
N0⎟
⎠
)
Structure of Optimum Receiver
⎡ ⎛ K ⎞⎤
s = arg max ln ⎢ pm ⋅ ( πN 0 )
−K 2
exp⎜ − ∑ rk − sm, k
⎝ k =1
(2
N0⎟ ⎥
⎠⎦
)
{s1,…,s M } ⎣
1 K
K
s = arg max ln[ pm ] − ln[ πN 0 ] − ∑ rk − sm, k (
2
)
{s ,…,s }
1 M
2 N 0 k =1
Kyung Hee
Structure of Optimum Receiver (continued)
University
K
s = arg max ln[ pm ] − ln[ πN 0 ]
{s1,…,s M } 2
1 ⎛ K 2 K K
2⎞
− ⎜ ∑ k r − 2 r s +
∑ k m,k ∑ m,k ⎟
s
N 0 ⎝ k =1 k =1 k =1 ⎠
⎡ ⎛ K ⎞⎤
s = arg max ln ⎢ pm ⋅ ( πN 0 )
−K 2
exp⎜ − ∑ rk − sm, k
⎝ k =1
(2
N0⎟ ⎥
⎠⎦
)
{s1,…,s M } ⎣
1 K
K
s = arg max ln[ pm ] − ln[ πN 0 ] − ∑ rk − sm, k (
2
)
{s ,…,s }
1 M
2 N 0 k =1
Kyung Hee
Structure of Optimum Receiver (continued)
University
K
s = arg max ln [pm ]− ln[ πN 0 ]
2
{s1,…,s M }
1 ⎛ K 2 K K
2⎞
− ⎜ ∑ k r − 2 r s +
∑ k m,k ∑ m,k ⎟
s
N 0 ⎝ k =1 k =1 k =1 ⎠
N0 K 1 K
s = arg max ln[ pm ] + ∑ rk sm,k − ∑ sm,k 2
{s1,…,s M } 2 k =1 2 k =1
Kyung Hee
Interpreting This Result
University
N0
ln[ pm ] weights the a priori probabilities
2
z If the noise is large, pm counts a lot
z If the noise is small, our received signal will be an accurate
estimate and pm counts less
K T
∑ rk sm, k = ∫ sm (t )r (t )dt represents the correlation with
k =1 0
the received signal
1 K 1T 2 E
∑ sm,k = ∫ sm (t )dt = m
2
represents signal energy
2 k =1 20 2
N0 K 1 K
s = arg max ln[ pm ] + ∑ rk sm,k − ∑ sm,k 2
{s1,…,s M } 2 k =1 2 k =1
K T
∑ rk sm, k = ∫ sm (t )r (t )dt represents the correlation with
k =1 0
the received signal
1 K 1T 2 E
∑ sm,k = ∫ sm (t )dt = m
2
represents signal energy
2 k =1 20 2
An Implementation of the Optimal Receiver -
Correlation Receiver
r (t ) × T
∫0 dt Σ Σ
s1(t ) − E1 2 N0
ln( p1 )
2 Choose
Largest
r (t ) × T
∫0 dt Σ Σ
sM (t ) − EM 2 N0
ln( pM )
2
Kyung Hee
Simplifications for Special Cases
University
r (t ) × T
∫0 dt Σ Σ
s1(t ) − E1 2 N0
ln( p1 )
2 Choose
Largest
r (t ) × T
∫0 dt Σ Σ
sM (t ) − EM 2 N0
ln( pM )
2
Young S. Kim Wireless Communication Systems Lab. - 19 -
r (t ) × T
∫0 dt r1
f1(t )
r = [ r1 rK ]
r (t ) × T
∫0 dt rK
f K (t )
K
∑ s1,k rk Σ Σ
k =1
r ×
− E1 2 N0
ln( p1) Choose
2
⎡ s11, sM,1 ⎤ Largest
⎢ ⎥
⎢ ⎥
⎢⎣s1, K K
sM, K ⎥⎦ ∑ Σ Σ
sM,k rk
k =1
− EM 2 N0 ln p
( M)
2
r (t ) × T
∫0 dt r1
f1(t )
r = [ r1 rK ]
r (t ) × T
∫0 dt rK
f K (t )
K
∑ s1,k rk Σ Σ
k =1
r ×
− E1 2 N0
ln( p1) Choose
2
⎡ s11, sM,1 ⎤ Largest
⎢ ⎥
⎢ ⎥
⎢⎣s1, K K
sM, K ⎥⎦ ∑ Σ Σ
sM,k rk
k =1
− EM 2 N0 ln p
( M)
2
Matched Filter Implementation
T T
Then rk = ∫ r ( t ) f k ( t )dt = ∫ r ( t ) f k (T − ( T − t ) )dt
0 0
T
= ∫ r ( t ) hk ( T − t )dt = r ( t ) ⊗ hk ( t ) t = T
0
where r(t ) ⊗ hk (t ) t =T denotes the convolution of the signals r(t)
and hk (t) evaluated at time T
t=T
r (t ) h1(t ) r1
r = [ r1 rK ]
t=T
r (t ) hK (t ) rK
T T
Then rk = ∫ r ( t ) f k ( t )dt = ∫ r ( t ) f k (T − ( T − t ) )dt
0 0
T
= ∫ r ( t ) hk ( T − t )dt = r ( t ) ⊗ hk ( t ) t = T
0
where r(t ) ⊗ hk (t ) t =T denotes the convolution of the signals r(t)
and hk (t) evaluated at time T
t=T
r (t ) h1(t ) r1
r = [ r1 rK ]
t=T
r (t ) hK (t ) rK
Example of Optimal Receiver Design
s 3 (t ) s 4 (t )
+1 +1
t t
1 2 -1 1 2
-1
f 1(t ) f 2(t )
+1 +1
t t
-1 1 2 -1 1 2
s1( t ) = 1 ⋅ f1 ( t ) + 1 ⋅ f 2 ( t ) s2 ( t ) = 1 ⋅ f1 ( t ) − 1 ⋅ f 2 ( t )
s3 ( t ) = −1 ⋅ f1 ( t ) + 1 ⋅ f 2 ( t ) s4 ( t ) = −1 ⋅ f1 ( t ) − 1 ⋅ f 2 ( t )
T=2
E1 = E2 = E3 = E4 = 2
s 3 (t ) s 4 (t )
+1 +1
t t
1 2 -1 1 2
-1
f 1(t ) f 2(t )
+1 +1
t t
-1 1 2 -1 1 2
s1( t ) = 1 ⋅ f1 ( t ) + 1 ⋅ f 2 ( t ) s2 ( t ) = 1 ⋅ f1 ( t ) − 1 ⋅ f 2 ( t )
s3 ( t ) = −1 ⋅ f1 ( t ) + 1 ⋅ f 2 ( t ) s4 ( t ) = −1 ⋅ f1 ( t ) − 1 ⋅ f 2 ( t )
T=2
E1 = E2 = E3 = E4 = 2
1st Implementation of Correlation Receiver
r (t ) × 2
∫0 dt Σ
s1(t ) N0
ln( p1 ) Choose
2
Largest
r (t ) × 2
∫0 dt Σ
s4 (t ) N0
ln( p4 )
2
r (t ) × 2
∫0 dt r1
f1(t )
r = [ r1 r2 ]
r (t ) × 2
∫0 dt r2
f2 (t )
r (t ) × 2
∫0 dt Σ
s1(t ) N0
ln( p1 ) Choose
2
Largest
r (t ) × 2
∫0 dt Σ
s4 (t ) N0
ln( p4 )
2
Young S. Kim Wireless Communication Systems Lab. - 27 -
r (t ) × 2
∫0 dt r1
f1(t )
r = [ r1 r2 ]
r (t ) × 2
∫0 dt r2
f2 (t )
Reduced Complexity Correlation Receiver -
Processing Stage
1 ⋅ r1 + 1 ⋅ r2 Σ
N0 ln( p1 ) 2
1 ⋅ r1 − 1 ⋅ r2 Σ Choose
N0 ln( p2 ) 2
Largest
−1 ⋅ r1 + 1 ⋅ r2 Σ
N0 ln( p3 ) 2
−1 ⋅ r1 − 1 ⋅ r2 Σ
N0 ln( p4 ) 2
t=2 r = [ r1 r2 ]
r (t ) h2 (t ) r2
1 ⋅ r1 + 1 ⋅ r2 Σ
N0 ln( p1 ) 2
1 ⋅ r1 − 1 ⋅ r2 Σ Choose
N0 ln( p2 ) 2
Largest
−1 ⋅ r1 + 1 ⋅ r2 Σ
N0 ln( p3 ) 2
−1 ⋅ r1 − 1 ⋅ r2 Σ
N0 ln( p4 ) 2
t=2 r = [ r1 r2 ]
r (t ) h2 (t ) r2
Summary of Optimal Receiver Design
Kyung Hee
Decision Regions
University
[ ]
N K 1 K
≥ 0 ln p j + ∑ rk s j ,k − ∑ s j ,k 2 , ∀i ≠ j
2 k =1 2 k =1
Decision Regions
[ ]
N K 1 K
≥ 0 ln p j + ∑ rk s j ,k − ∑ s j ,k 2 , ∀i ≠ j
2 k =1 2 k =1
Kyung Hee
Average Energy Per Bit: E
University b
K
Ei = ∑ si ,k 2 is the energy of the ith signal
k =1
Es
Eb = is the average energy per bit
log 2 M
Es
Eb = is the average energy per bit
log 2 M
Kyung Hee
Examples of Decision Regions - QPSK
University
0.8
0.6
0.4
0.2
-0.2
-0.4
-0.6
-0.8
-1
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
0.8
0.6
0.4
0.2
-0.2
-0.4
-0.6
-0.8
-1
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
QPSK with Unequal Signal Probabilities
0.8
0.6
0.4
0.2
-0.2
-0.4
-0.6
-0.8
-1
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
0.8
0.6
0.4
0.2
-0.2
-0.4
-0.6
-0.8
-1
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
0.8
0.6
0.4
0.2
-0.2
-0.4
-0.6
-0.8
-1
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
0.8
0.6
0.4
0.2
-0.2
-0.4
-0.6
-0.8
-1
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
Unequal Signal Powers
3.5
2.5
1.5
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Kyung Hee
Signal Constellation for 16-ary QAM
University
sigspace( [1.5 -1.5; 1.5 -0.5; 1.5 0.5; 1.5 1.5; 0.5 -1.5; 0.5 -
0.5; 0.5 0.5; 0.5 1.5;-1.5 -1.5; -1.5 -0.5; -1.5 0.5; -1.5 1.5; -
0.5 -1.5; -0.5 -0.5; -0.5 0.5; -0.5 1.5],10)
2
1.5
0.5
-0.5
-1
-1.5
-2
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2
3.5
2.5
1.5
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4
sigspace( [1.5 -1.5; 1.5 -0.5; 1.5 0.5; 1.5 1.5; 0.5 -1.5; 0.5 -
0.5; 0.5 0.5; 0.5 1.5;-1.5 -1.5; -1.5 -0.5; -1.5 0.5; -1.5 1.5; -
0.5 -1.5; -0.5 -0.5; -0.5 0.5; -0.5 1.5],10)
2
1.5
0.5
-0.5
-1
-1.5
-2
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2
Notes on Decision Regions
Signal points need not lie within their decision regions for
case of low Eb N0 and unequal probabilities.