6 Digital Modulation Muqaibel
6 Digital Modulation Muqaibel
Modulation
and
Detection
PROF. ALI HUSSEIN MUQAIBEL
R E F E R E N C E B O O K B Y
A N D R E A G O L D S M I T H 5 . 1 U P T O T H E E N D O F 5 . 4
Geometric
Representation
of Signals
Dr. Ali Muqaibel
Ver. 3.0
Chapter 5 (Simon Haykin, Communication Systems) If you have any comments or corrections on
the slides , you may send to Dr. Ali Muqaibel,
Chapters 5 (Andrea Goldsmith) [5.1 – 5.5]
[email protected]
Outlines
• Signal-space analysis
• Euclidean distance
𝑚𝑖 𝑠𝑖 (𝑡) 𝑥(𝑡) 𝑚
ෝ
Message
Transmitter Channel Receiver
Source
𝑚:
ෝ estimate of 𝑚𝑖
• Message source emits a symbol from an alphabet of 𝑀 symbols (𝑚1 , 𝑚2 , … , 𝑚𝑀 ) every 𝑇 seconds.
𝑃𝑒 = 𝑝𝑖 𝑃 𝑚
ෝ ≠ 𝑚𝑖 |𝑚𝑖
𝑖=1
𝑃 𝑚
ෝ ≠ 𝑚𝑖 |𝑚𝑖 : conditional error probability given that 𝑚𝑖 was sent
𝑠𝑖 (𝑡) 𝑥(𝑡)
+
Transmitted Received
signal signal
𝑤(𝑡)
AWGN
• We will use geometric representation for the known set of transmitted signals {𝑠𝑖 (𝑡)}
Basis function
𝑁
0≤𝑡≤𝑇
𝑠𝑖 𝑡 = 𝑠𝑖𝑗 𝜙𝑗 (𝑡) , ቊ
𝑖 = 1,2, … , 𝑀
𝑗=1
• where the coefficients of the expansion are
𝑇
𝑖 = 1,2, … , 𝑀
𝑠𝑖𝑗 = න 𝑠𝑖 𝑡 𝜙𝑗 𝑡 𝑑𝑡 , ቊ
0
𝑗 = 1,2, … , 𝑁
• The real-valued basis functions 𝜙1 𝑡 , 𝜙2 𝑡 ,…, 𝜙𝑁 𝑡 are orthonormal
Correlator 𝑇
𝑠𝑖1 × (product-integrator) × න 𝑑𝑡 𝑠𝑖1
0
𝜙1 (𝑡) 𝜙1 (𝑡)
𝑇
𝑠𝑖2 × + 𝑠𝑖 (𝑡) 𝑠𝑖 (𝑡) × 𝑠𝑖2
න 𝑑𝑡
0
𝜙2 (𝑡) 𝜙2 (𝑡)
: :
:
𝑇
𝑠𝑖𝑁 × × 𝑠𝑖𝑁
න 𝑑𝑡
𝑁 0
0≤𝑡≤𝑇
𝑠𝑖 𝑡 = 𝑠𝑖𝑗 𝜙𝑗 (𝑡) , ቊ
𝑖 = 1,2, … , 𝑀
𝑗=1
𝜙𝑁 (𝑡) 𝜙𝑁 (𝑡)
𝑇
𝑠𝑖𝑗 = න 𝑠𝑖 𝑡 𝜙𝑗 𝑡 𝑑𝑡 , ቊ
𝑖 = 1,2, … , 𝑀
0
𝑗 = 1,2, … , 𝑁
8 Dr. Ali Hussein Muqaibel
Geometric Representation of Signals
𝑠𝑖1
𝑠𝑖2
Column vector 𝒔𝒊 = : , 𝑖 = 1,2, … , 𝑀
𝑠𝑖𝑁
• Extend the conventional notion of 2D and 3D (Euclidean space) to an 𝑁 -dimensional
Euclidean space.
→ Set of signal vectors {𝒔𝒊 |𝑖 = 1,2, . . , 𝑀} represents set of 𝑀 points in the 𝑁-dimensional
Euclidean space (signal space)
• What is the importance of representing energy signals geometrically?
• For mathematical tractable analysis and providing the basis to noise consideration in digital
communication system.
𝜙1
𝑗=1 0
𝜙1 𝑡
• The angle, 𝜃𝑖𝑘 , between two signal vectors 𝒔𝒊 and 𝒔𝒌 is given by
𝒔𝑇𝒊 𝒔𝒌
cos𝜃𝑖𝑘 =
𝒔𝒊 𝒔𝒌
• Signal vectors 𝒔𝒊 and 𝒔𝒌 are orthogonal (perpendicular) if 𝒔𝑇𝒊 𝒔𝒌
= 0 ⇒ 𝜃𝑖𝑘 = 90°
Gram-Schmidt Orthogonalization
Procedure
Dr. Ali Hussein Muqaibel
𝑠1 (𝑡) 𝑠2 (𝑡)
1 1
𝑡 𝑡
2 -1 2
𝑠3 (𝑡) 𝑠4 (𝑡)
1
𝑡 𝑡
-1 2 -1 2
Visit my Website and YouTube Channel for
more resources
𝑔𝑖 𝑡 = 𝑠𝑖 𝑡 − 𝑠𝑖𝑗 𝜙𝑗 (𝑡)
𝑗=1
• where the coefficients 𝑠𝑖𝑗 are
𝑇
𝑠𝑖𝑗 = න 𝑠𝑖 𝑡 𝜙𝑗 𝑡 𝑑𝑡 , 𝑗 = 1,2, . . , 𝑖 − 1
0
• Given 𝑔𝑖 (𝑡), then
𝑔𝑖 (𝑡)
𝜙𝑖 𝑡 = , 𝑖 = 1,2, . . , 𝑁
𝑇
0 𝑔𝑖2 𝑡 𝑑𝑡
• which forms an orthonormal set of basis functions.
• Example:
Symbol Amplitude Gray code
𝑠1 (𝑡) -3 00
𝑠2 (𝑡) -1 01
𝒔𝟏 𝒔𝟐 𝒔𝟑 𝒔𝟒 𝑠3 (𝑡) +1 11
𝜙1 (𝑡)
𝑠4 (𝑡) +3 10
-3 -1 0 +1 +3
1 1
𝑡 𝑡
2 -1 2
𝑠3 (𝑡) 𝑠4 (𝑡)
1
𝑡 𝑡
-1 2 2
-1
𝜙2 (𝑡)
𝑠4 𝑡 = 𝑠41 𝜙1 𝑡 [No new basis function]
1/ 2
𝑡
2
−1/ 2
− 2 2
𝜙1 (𝑡)
− 2
10 Dr. Ali Hussein Muqaibel
Practice: Gram Schmidt orthogonalization procedure
1. Detector 2. Decoder
This part consists of a bank of 𝑁 product- Performs maximum likelihood decision rule on
integrators (or correlators), that has input 𝑥(𝑡) and observation vector 𝒙 to obtain the estimate 𝑚 ෝ
outputs observation vector 𝒙 of transmitted symbol 𝑚𝑖 , 𝑖 = 1, . . , 𝑀 such that
the average probability of error is minimized.
𝑇
× න 𝑑𝑡 𝑥1
0
Observation vector 𝒙
𝜙1 (𝑡)
𝑇
𝑥(𝑡) × න 𝑑𝑡 𝑥2
0
𝜙2 (𝑡)
:
:
𝑇
× න 𝑑𝑡 𝑥𝑁
0
𝜙𝑁 (𝑡)
Dr. Ali Hussein Muqaibel
Equivalence of correlator and matched filter
receivers
• Detector part in correlation receiver has a set of 𝑁 correlators that
can be replaced by matched filters
• Let ℎ𝑗 (𝑡) be impulse response of a linear time-invariant (LTI) filter
𝑥(𝑡) ℎ𝑗 𝑡 = 𝜙𝑗 (𝑇 − 𝑡) 𝑦𝑗 (𝑡)
𝜙1 (𝑡)
𝑇
𝑥(𝑡) × න 𝑑𝑡 𝑥2 𝜙2 (𝑇 − 𝑡)
0
𝜙2 (𝑡)
:
:
𝑇
× න 𝑑𝑡 𝑥𝑁
0
𝜙𝑁 (𝑇 − 𝑡)
𝜙𝑁 (𝑡)
orthogonal matrix is a real square matrix whose columns/and rows are orthogonal unit vectors (orthonormal vectors)
𝐸𝑠 = 𝒔𝒊 2𝑝
𝑖
𝑖=1
where 𝑝𝑖 = 𝑃(𝑚𝑖 ).
• If the signals 𝒔𝟏 , . . , 𝒔𝑴 are translated by a vector 𝒂, then the average energy of translated
signal constellation is
𝑀
𝐸ത𝑠 = 𝒔𝒊 − 𝒂 2 𝑝𝑖
𝑖=1
• Minimum Energy Translate
• Given a signal constellation 𝒔𝒊 𝑀
𝑖=1 , the corresponding signal constellation with minimum
average energy is obtained by subtracting from each signal vector, 𝒔𝒊 , a constant 𝐸 𝒔
(B1)
0
(A1) 𝜙1
(A2) 0
𝜙1 (B2)
𝛼 2𝛼 3𝛼
Shaded region
represents 𝑃(𝐴1,2 )
𝑃(𝐴1,4 )
𝑃(𝐴1,3 )
2
Dr. Ali Hussein Muqaibel
Union Bound on Probability of Error
• There is a difference between 𝑃 𝐴𝑖,𝑘 and 𝑃(𝑚
ෝ = 𝑚𝑘 |𝑚𝑖 )
▪ 𝑃(𝑚 ෝ = 𝑚𝑘 |𝑚𝑖 ) indicates that 𝒙 lies closer to 𝒔𝒌 compared to any other vector if 𝒔𝒊 was
sent.
▪ 𝑃(𝐴𝑖,𝑘 ) depends only on two signal vectors 𝒔𝒊 and 𝒔𝒌 (pairwise error probability,
𝑃2 (𝒔𝒊 , 𝒔𝒌 ))
• Pairwise error probability: It is the probability that a receiver mistakes 𝒔𝒌
for 𝒔𝒊 assuming the system uses only pair of signals (𝒔𝒊 and 𝒔𝒌 )
𝑀
𝑃𝑒 (𝑚𝑖 ) ≤ 𝑃2 (𝒔𝒊 , 𝒔𝒌 )
𝑘=1
𝑘≠𝑖
• Recall that
2 ∞
erfc 𝑢 = න exp −𝑧 2 𝑑𝑧
𝜋 𝑣
𝜙1 (𝑡)
• Let 𝑧 = 𝑣/ 𝑁0 , then
1 𝑑𝑖𝑘
𝑃2 (𝒔𝒊 , 𝒔𝒌 ) = erfc
2 2 𝑁0
𝑃𝑒 = 𝑝𝑖 𝑃𝑒 𝑚𝑖
𝑖=1
𝑀 𝑀
1 𝑑𝑖𝑘
≤ 𝑝𝑖 erfc
2 2 𝑁0
𝑖=1 𝑘=1
𝑘≠𝑖
• There are special forms that simplifies and bound the
expression of 𝑃𝑒
𝑃𝑒 ≤ 𝑝𝑖 erfc
𝑑𝑖𝑘
2 2 𝑁0
𝑖=1 𝑘=1
𝑘≠𝑖
• Case 1: Circularly Symmetric
• If the signal constellation is circularly symmetric about the origin
→ 𝑃𝑒 𝑚𝑖 is the same for all 𝑖
𝑀 𝜙2 (𝑡)
1 𝑑𝑖𝑘
𝑃𝑒 ≤ erfc
2 2 𝑁0
𝑘=1
𝑘≠𝑖 𝜙1 (𝑡)
• Example:
• QPSK signal constellation
Probability 1
𝑃𝑒 ≤ 𝑝𝑖 erfc
2
𝑖=1 𝑘=1
𝑑𝑖𝑘
2 𝑁0
𝑘≠𝑖
• Case 2: Minimum Distance Bound
• Let the minimum distance of signal constellation, 𝑑𝑚𝑖𝑛 , be the smallest
Euclidean distance between any two signal points
𝑑𝑚𝑖𝑛 = min 𝑑𝑖𝑘
𝑖≠𝑘
• erfc(𝑢) is monotonically decreasing function of 𝑢, then
𝑑𝑖𝑘 𝑑𝑚𝑖𝑛
erfc ≤ erfc ∀𝑖, 𝑘 ∈ {1, . . , 𝑀}
2 𝑁0 2 𝑁0
• Hence,
(𝑀−1) 𝑑𝑚𝑖𝑛
𝑃𝑒 ≤ erfc
2 2 𝑁0
• Example:
• 16-QAM or signal constellations with large number of points
▪ Non-linear modulation
▪ FSK
▪ MSK
CPM : Continuous Phase Modulation
CPFSK : Continuous Phase FSK
CP-BFK : Continuous Phase-Binary FSK
MSK : Minimum Shift Keying
PSK : Phase Shift Keying
QPSK : Quadrature Phase Shift Keying
O-QPSK : Offset Quadrature Phase Shift Keying
2
𝜙𝑖 𝑡 = cos 2𝜋𝑓𝑖 𝑡 , 0 ≤ 𝑡 ≤ 𝑇𝑏
𝑇𝑏
0, elsewhere
where 𝑖 = 1,2.
▪ BFSK is a continuous phase signal
𝑇𝑏
𝑠𝑖𝑗 = න 𝑠𝑖 𝑡 𝜙𝑗 𝑡 𝑑𝑡
0 𝑇
𝒔𝟐 = 0 𝐸𝑏
𝑇𝑏
2𝐸𝑏 2
=න cos 2𝜋𝑓𝑖 𝑡 cos 2𝜋𝑓𝑗 𝑡 𝑑𝑡
0 𝑇𝑏 𝑇𝑏
𝑇
𝐸 , 𝑖=𝑗 𝒔𝟏 = 𝐸𝑏 0
=൝ 𝑏
0, 𝑖≠𝑗
1 𝐸𝑏
𝑃𝑒 = erfc
2 2𝑁0
1 𝐸𝑏
𝑃𝑒 = erfc
2 𝑁0
2𝐸𝑏
𝑠𝑖 𝑡 = cos 2𝜋𝑓𝑖 𝑡 + 𝜃(0) , 0 ≤ 𝑡 ≤ 𝑇𝑏
𝑇𝑏
0, elsewhere
• where 𝜃 0 denotes phase value at 𝑡 = 0 (summing up the past), 𝑖 = 1,2, 𝐸𝑏 is transmitted energy per bit, and
𝑇𝑏 is bit duration.
▪ Binary symbol 1 → 𝑠1 𝑡 → 𝑓1
▪ Binary symbol 0 → 𝑠2 (𝑡) → 𝑓2
2𝐸𝑏
𝑠 𝑡 = cos 2𝜋𝑓𝑐 𝑡 + 𝜃(𝑡)
𝑇𝑏
4𝜋𝑑
• where the continuous-time phase is 3𝜋𝑑
2𝜋𝑑
𝜃 𝑡 − 𝜃(0)
𝜋𝑑
𝜃 𝑡 = 𝜃 0 ± 𝑇 𝑡, 0 ≤ 𝑡 ≤ 𝑇𝑏 𝜋𝑑
𝑏 𝑡
−𝜋𝑑
• and deviation ratio −2𝜋𝑑
−3𝜋𝑑
𝑑 = 𝑇𝑏 (𝑓1 − 𝑓2 ) −4𝜋𝑑
𝜋𝑑
• Binary symbol 1 ⇒ 𝜃 𝑡 = 𝜃 0 + 𝑇 𝑡 ⇒ 𝑠 𝑡 = 𝑠1 (𝑡) 2𝑇𝑏 4𝑇𝑏 6𝑇𝑏 8𝑇𝑏
𝑏
𝜋𝑑
• Binary symbol 0 ⇒ 𝜃 𝑡 = 𝜃 0 − 𝑇 𝑡 ⇒ 𝑠 𝑡 = 𝑠2 (𝑡)
𝑏
𝜃 𝑡 − 𝜃(0)
2𝜋𝑑
𝜋𝑑 𝜋𝑑
• At 𝑡 = 𝑇𝑏 , 𝜃 𝑡 =𝜃 0 ± 𝑡, 0 ≤ 𝑡 ≤ 𝑇𝑏 𝑡
𝑇𝑏 −𝜋𝑑
−2𝜋𝑑
−3𝜋𝑑
𝜋𝑑, for symbol 1 −4𝜋𝑑
𝜃 𝑇𝑏 − 𝜃 0 = ቊ
−𝜋𝑑, for symbol 0
2𝑇𝑏 4𝑇𝑏 6𝑇𝑏 8𝑇𝑏
• Sending the symbol 1 increases the phase of CPFSK signal by 𝜋𝑑 radians; however, sending
the symbol 0 will reduce it by 𝜋𝑑
• Different variations of phase, 𝜃(𝑡), over time results in a phase tree (Phase Trellis)
1
→ CPFSK with 𝑑 = 1/2 is called minimum shift keying (MSK) , 𝑑 = 𝑇𝑏 (𝑓1 − 𝑓2 ), 𝑓1 − 𝑓2 =
2𝑇𝑏
𝜋
→ 𝜃 𝑡 = ± , 𝑡 = 𝑇𝑏
2
• Example:
• Determine the phase trellis for a binary sequence 1101000 with 𝜃 0 = 0 and 𝑑 = 1/2.
• Frequency deviation is the minimum frequency spacing between 𝑓1 and 𝑓2 to be coherently orthogonal
.
1
→ CPFSK with 𝑑 = 1/2 is called minimum shift keying (MSK) , 𝑑 = 𝑇𝑏 (𝑓1 − 𝑓2 ), 𝑓1 − 𝑓2 = 2𝑇
𝑏
𝜋
→ 𝜃 𝑡 = ±2, 𝑡 = 𝑇𝑏
𝑀 2 4 8 16 32 64
𝜌 1 1 0.75 0.5 0.312 0.187
𝑔1 𝑡 + 𝑤 𝑡 , 𝑠1 𝑡 sent
𝑥 𝑡 =ቊ
𝑔2 𝑡 + 𝑤 𝑡 , 𝑠2 𝑡 sent
where 0 ≤ 𝑡 ≤ 𝑇, 𝑤(𝑡) is additive white Gaussian noise (zero mean and PSD of 𝑁0 /2)
2𝐸𝑏
𝑠𝑖 𝑡 = cos 2𝜋𝑓𝑖 𝑡 , 0 ≤ 𝑡 ≤ 𝑇𝑏
𝑇𝑏
0, elsewhere
1 𝐸𝑏
𝑃𝑒 = exp −
2 2𝑁0 Noncoherent receiver for BFSK signals
2𝐸𝑏
cos 2𝜋𝑓𝑐 𝑡 , 0 ≤ 𝑡 ≤ 𝑇𝑏
𝑇𝑏
𝑠1 𝑡 =
2𝐸𝑏
cos 2𝜋𝑓𝑐 𝑡 , 𝑇𝑏 ≤ 𝑡 ≤ 2𝑇𝑏
𝑇𝑏
2𝐸𝑏
cos 2𝜋𝑓𝑐 𝑡 , 0 ≤ 𝑡 ≤ 𝑇𝑏
𝑇𝑏
𝑠2 𝑡 =
2𝐸𝑏
cos 2𝜋𝑓𝑐 𝑡 + 𝜋 , 𝑇𝑏 ≤ 𝑡 ≤ 2𝑇𝑏
𝑇𝑏
1 𝐸𝑏
𝑃𝑒 = exp −
2 𝑁0
{𝒃𝒌 } 1 0 0 1 0 0 1 1
{𝑑𝑘−1 } 1 1 0 1 1 0 1 1
{𝑑𝑘 } 1 1 0 1 1 0 1 1 1
𝜃 0 0 𝜋 0 0 𝜋 0 0 0
DPSK 1
exp −𝐸𝑏 /𝑁0
2
Noncoherent binary FSK 1
exp −𝐸𝑏 /2𝑁0
2
4 0.5 0.34 dB
8 0.333 3.91 dB
16 0.25 8.52 dB
32 0.2 13.52 dB
For the same Data rate
• For 𝑀 > 4
▪ M-PSK is circular
▪ M-QAM is rectangular
• Distance between points is smaller in M-PSK than M-QAM
A signal constellation diagram is shown in the figure. Let the basis functions be as follows:
2 2
1 = cos ( 2 f c t ) , 0 t T
T 2
where f c = b) 01 a) 00
2 T
2 = sin ( 2 f c t ) , 0 t T
T
1
a) What type of modulation is being used (ASK, PSK, or FSK)? Why?
c) 10 11
b) Build (sketch) the block diagram for correlators based coherent receiver.
What is the received sequence of bits? Show how you get the answer?
d) Let a=1, Assuming that all decisions where correct in the previous question. Which one of the
three symbols was exposed to the largest amount of noise (Show how you find the answer)? How
much is this noise energy?
e) Re-draw the previous constellation diagram and re-assign the bits to the constellation points in a
better way. Justify your answer?
f) A new constellation diagram is shown. Assign the bits and find the values of x and y as
function of a such that the new constellation diagram exhibit the same probability of error as the
original diagram.
2
x
1
y
-x
g) On the new constellation diagram (part f), draw the decision boundaries.
h) Calculate the average energy per bit for the original constellation diagram and for the diagram
introduced in part (f). Which one of the two constellations do you prefer? Justify your answer.
i) Show that the two basis functions are orthogonal. (show all details)
2sin cos = sin ( 2 )
________________________________________________________________________
Good Luck, Dr. Ali Muqaibel