Modulation

Symbol Error Rate for 16PSK

Krishna Sankar17 years ago12 years ago314 mins

In this post, let us try to derive the symbol error rate for 16-PSK (16-Phase Shift Keying) modulation.

Consider a general M-PSK modulation, where the alphabets,

$alpha_{16PSK} =\sqrt{E_s}\left\{1,\ e^{\frac{j2\pi}{M}},\ e^{\frac{j4\pi}{M}},\ \ldots,\ e^{\frac{j2\pi(M-1)}{M}} \right}$ are used.

(Refer example 5-38 in [DIG-COMM-BARRY-LEE-MESSERSCHMITT])

Figure: 16-PSK constellation plot

Deriving the symbol error rate

Let us the consider the symbol on the real axis, i.e

$s_0=\sqrt{E_s}$ .

The received symbol $y=\sqrt{E_s} + n$ .

Where the additive noise $n$ follows the Gaussian probability distribution function,

$p(x) = \frac{1}{\sqrt{2\pi\sigma^2}}e^{\frac{-(x-\mu)^2}{2\sigma^2}$ with $\mu=0$ and $\sigma^2 = \frac{N_0}{2}$ .

The conditional probability distribution function (PDF) of received symbol $y$ given $s_0$ was transmitted is:

$p(y|s_0)=\frac{1}{\sqrt{\pi N_0}}e^{-\frac{(x-\sqrt{E_s})^2}{N_0}}$ .

As can be seen from the figure above, due to the addition of noise, the transmitted symbol gets spreaded. However, if the received symbol is present with in the boundary defined by the magenta lines, then the symbol will be demodulated correctly.

To derive the symbol error rate, the objective is to find the probability that the phase of the received symbol lies within this boundary defined by the magenta lines i.e. from $-\frac{\pi}{M}$ to $+\frac{\pi}{M}$ .

For simplifying the derivation, let us make the following assumptions:

(a) The signal to noise ratio, $\frac{Es}{N_0}$ is reasonably high.

For a reasonably high value of $\frac{Es}{N_0}$ , then the real part of the received symbol is not afected by noise i.e.,

$\Re{y}\approx\sqrt{E_s}$ and

the imaginary part of the received symbol is equal to noise, i.e.
$\Im{y}=n$ .

(b) The value of M is reasonably high (typically M >4 suffice)

For a reasonably high value of M, the constellation points are closely spaced. Given so, the distance of the constellation point $s_0$ to the magenta line can be approximated as $\sqrt{E_s}sin(\frac{\pi}{M})$ .

Figure: Distance between constellation points

Given the above two assumptions, it can be observed that the symbol $s_0$ will be decoded incorrectly, if the imaginary component of received symbol $y$ isgreater than $\sqrt{E_s}sin(\frac{\pi}{M})$ . The probability of $y$ being greater than $\sqrt{E_s}sin(\frac{\pi}{M})$ is,

$p\left(\Im(y)>\sqrt{E_s}sin(\frac{\pi}{M})\right)=\frac{1}{\sqrt{\pi N_0}}\int_{\sqrt{E_s}sin(\frac{\pi}{M})}^{\infty}e^{-\frac{{\Im{y}}^2}{N_0}}dy$ .

Changing the variable to $u=\frac{\Im{y}}{\sqrt{N_0}}$ ,

$p\left(\Im(y)>\sqrt{E_s}sin(\frac{\pi}{M})\right)=\frac{1}{\sqrt{\pi}}\int_{\sqrt{\frac{E_s}{N_0}}sin(\frac{\pi}{M})}^{\infty}e^{-u^2}du=\frac{1}{2}erfc\left[\sqrt{\frac{E_s}{N_0}}sin(\frac{\pi}{M})\right]$ .

Note: The complementary error function, $erfc(x) = \frac{2}{\sqrt{\pi}}\int_x^\infty e^{-x^2}dx$ .

Similarly, the symbol $s_0$ will be decoded incorrectly, if the imaginary component of received symbol $y$ is less than $-\sqrt{E_s}sin(\frac{\pi}{M})$ . The probability of $y$ being less than $-\sqrt{E_s}sin(\frac{\pi}{M})$ is,

$p(\Im(y) <-\sqrt{E_s}sin(\frac{\pi}{M})=\frac{1}{\sqrt{\pi N_0}}\int_{-\infty}^{-\sqrt{E_s}sin(\frac{\pi}{M})}e^{\frac{{-\Im y}^2}{N_0}}dy=\frac{1}{2}erfc\left[\sqrt{\frac{E_s}{N_0}}sin(\frac{\pi}{M})\right]$ .

The total probability of error given $s_0$ was transmittd is,

$p(e|s_0) =erfc\left[\sqrt{\frac{E_s}{N_0}}sin(\frac{\pi}{M})\right]$ .

Total symbol error rate

The symbol will be in error, if atleast one of the symbol gets decoded incorrectly. Hence the total symbol error rate from M-PSK modulation is,

$P_{e,MPSK}=erfc\left[\sqrt{\frac{E_s}{N_0}}sin(\frac{\pi}{M})\right]$ .

Simulation model

Simple Matlab/Octave script for simulating transmission and recepetion of an M-PSK modulation is attached. It can be observed that the simulated symbol error rate compares well with the theoretical symbol error rate.

Click here to download: Matlab/Octave script for simulating symbol error rate curve for 16 PSK modulation

Figure: Symbol Error rate curve for 16PSK modulation

Hope this helps.

Krishna

References

[DIG-COMM-BARRY-LEE-MESSERSCHMITT] Digital Communication: Third Edition, by John R. Barry, Edward A. Lee, David G. Messerschmitt

31 thoughts on “Symbol Error Rate for 16PSK”

Pingback: Download free e-book on error probability in AWGN
abdelfatah elkhalifa abdelati says:

April 26, 2012 at 12:25 pm

please help me in MATLAB simulink for rain attenuation and graphiical user interface . at least the block diagram . am MSC student

Reply
1. Krishna Sankar says:
  
  April 27, 2012 at 5:44 am
  
  @abdelfatah : sorry, i do not have simulink and have not also not studied about rain attenuation
  
  Reply
imran says:

April 25, 2012 at 9:41 pm

please help me by giving me the Symbol Error Rate (SER)/ Bit error rate formula/Equation for 16-APSK modulation. because i need the theoretical equation for plotting and compare simulation result.
waiting someones feedback

Reply
1. Krishna Sankar says:
  
  April 26, 2012 at 5:53 am
  
  @imran: sorry, i have not discussed 16APSK schemes in the posts
  
  Reply
Ravi says:

March 25, 2012 at 9:38 pm

Dear Mr. Krishna
I want to develop a MATLAB code for MIMO using MMSE equalizer for 4X4 antenna configuration and using BPSK, QPSK and 16-QAM modulation. Suggest me any MATLAB code for reference.

I hope you will reply me.

Reply
1. Krishna Sankar says:
  
  March 26, 2012 at 5:36 am
  
  @Ravi: You can check out posts on MIMO at
  https://dsplog.com/tag/mimo
  Most of them discuss 2Tx-2Rx IID Rayleigh channel system with BPSK case
  
  Reply
NH says:

February 4, 2012 at 6:07 am

Hi,
Why did you use N=10^5 a big number?

Reply
1. Krishna Sankar says:
  
  February 4, 2012 at 11:15 am
  
  @NH: To get a good average in symbol error rate….
  
  Reply
abhijeet says:

October 29, 2011 at 9:44 pm

sir,
i m doing b.tech in ece branch from bit mesra. i need a code for finding symbol error rate of bfsk signal in presence of additive gaussian noise.

please help me out in this.

Reply
1. Krishna Sankar says:
  
  October 30, 2011 at 7:23 pm
  
  @abhijeet: Does this post on frequency shift keying with coherent demodulation help?
  https://dsplog.com/2007/08/30/bit-error-rate-for-frequency-shift-keying-with-coherent-demodulation/
  
  Reply
lyla says:

February 22, 2010 at 9:38 am

Hello, I have a problem that I need to solve about a
modulation algorithm. I’m working on a assignment and need a code in Matlab to the
in-phase and quadrature components of a 16 QAM. Afterwards, I will run the
program in a Vector Signal Analyzer and I will have to see the I and Q waveforms
and the constellation.

I have been trying out by creating an array with 16 values (0 to 15) and
converting them to binary (0000 to 1111), afterwards, I have created the I and Q
signals, as a sine or cosine of the array values mentioned respectively using a
for loop, thus creating the different values for the phase, but I don’t know how
to express the amplitude levels (-3, -1, 1, 3).

I would appreciate if someone could lend me a hand in this problem.

Thanks in advance

Reply
1. Krishna Sankar says:
  
  March 31, 2010 at 5:32 am
  
  @lyla: You can generate random 16QAM constellation points, pass them through a pulse shaping filter and then feed to the VSA
  
  Reply
Moongeun says:

August 19, 2009 at 1:12 pm

hi~~Krishna

i have a question~~~

in the matbla source code~

ipPhasehat=2*pi/M*round(opPhase/(2*pi/M))

i don`t know why you use this code~~~

anyway thank you for your efforts~

Reply
1. Krishna Sankar says:
  
  August 22, 2009 at 4:56 am
  
  @Moongeun: To round the received phase to the nearest constellation point.
  
  Reply
Tukuben says:

April 19, 2009 at 6:31 am

Hi Krishna,
i`m an MSc communication systems engineering student in the university of portsmouth, i would love to be a member of this forum and there is this problem i want to solve:
Analysis & Simulation of 16PSK Using
MATLAB.

OBJECTIVES:
· Demonstrate the ability to learn MATLAB for the purpose of this
coursework
· Write a technical report describing the principle, theory and applications
of 16PSK. Your system will be tested over the AWGN Channel.
· Implement a MATLAB simulation of 16PSK modulator and
demodulator system.
· Measurement of the symbol error rate and bit error rate as a function of
the ratios Eb / No and Es / No achieved by the system in the presence of
additive white Gaussian Noise. Your measurements should be compared
to the theoretical estimates that you should have derived from first
principles.
· Show that you are able to identify and use technical references.

But for now am learning MATLAB and the problem is giving me though time, i wish you can help me by putting me through with the rest apart from the MATLAB. Please sir i need your help.
Regards.
Tukuben

Reply
1. Krishna Pillai says:
  
  April 19, 2009 at 6:54 am
  
  @Tukuben: The following posts maybe helpful to you
  https://dsplog.com/2008/03/18/symbol-error-rate-for-16psk/
  https://dsplog.com/2008/05/18/bit-error-rate-for-16psk-modulation-using-gray-mapping/
  
  Reply
abdulmaleek says:

January 29, 2009 at 3:40 pm

I am an MSc mobile, personal and satellite communication student in university of westminster, london. i like to be member of this forum

Reply
1. Krishna Pillai says:
  
  January 30, 2009 at 5:44 am
  
  @abdulmaleek: There are two options:
  
  1. You can subscribe to the email newsletter by providing you email address in the form on the top right of the page. Email subscription will entitle you to receive the free e-book on error rates in AWGN.
  
  2. Further, you can join via Google FriendConnect option provided in the bottom of the page. Google FriendConnect enables you to share the contents of this blog to your friends/colleagues in an easier fasion.
  
  Hope this helps.
  
  Reply
gautam says:

January 19, 2009 at 1:26 pm

sir
i am doing ME in digital communication. please refer me any project in MATLAB coding related to my subject.also suggest any reference material for MATLAB study

Reply
1. Krishna Pillai says:
  
  January 22, 2009 at 3:24 am
  
  @gautam: Maybe you can try modeling the transmitter, receiver and channel per the 3GPP LTE specifications? It includes MIMO, advanced coding scheme, OFDMA etc. Quite a lot to learn. 🙂
  
  Reply
Krishna says:

April 10, 2008 at 4:33 am

@Rajesh:
1. If the device is non-linear, one would expect that the signal at the output will have more frequencies than the signal at the input (generation of harmonics dye to x^2, x^3 terms etc).

The out-of-band frequencies maybe filtered out, however nothing much can be done about the harmonics in the in-band. It may be that the measured power includes both out-band and in-band harmonics (along with the desired signal power), hence the C/N does not degrade. However, since the in-band harmonics are not deisrable for demodulation, there is degradation in the MER.

2. When you say ‘devices not matched’, I would think that you are referring to impedence mismatch between the devices (resulting in reflection…). May I suggest the generation of in-band and out-band frequencies due to impedence mismatch as a probable reason for MER degradation.

Does this perspective makes sense? Thanks for these nice questions… 🙂

Reply
Rajesh Kher says:

April 8, 2008 at 4:58 pm

Dear Krishna

Thanks for the reply. While you are right, however we will have to add the noise level. Now the noise level at room temp will be about 3 dBuV while that for QPSK with MER of 15 dB at Signal strength of 86dBuV will be approx 71 dBuV. So that the additive White noise is negligible. We also find that C/N does not get deterioted as we pass thru varios devices.

However we have noticed that MER doe reduce by 0.7-1.0dB.

Question is

1. If thye device contains diodes then if at the point of operation diode is non linear it can add to the reduction of MER but not of C/N. How MER is related to non-linear CSO and CRB.
2. If the device is poorly matched then agian MER measurement will get affected. How does the mismatch between the device and the instruments affects the MER?

Reagards

Rajesh Kher

Reply
Krishna says:

April 3, 2008 at 8:30 pm

@ Rajesh:
1. By definition modulation error ratio (and not rate) is,
MER, dB = 10*log10(average transmit symbol power/average error power) where
error is the distance between the transmitted and received constellation.

If we consider that the received symbol Y is
Y = S + N, where S is the transmitted symbol, N the noise.
then error = Y-S = N.

From the above equation, it is reasonably intuitive that modulation error ratio (MER) is comparable to signal to noise ratio (SNR) in noise only scenario.

In the presence of other impairments, MER will reflect the distortions in the constellation due to phase noise, IQ imbalance, frequency offset etc.

Slides 20-22 from http://chapters.scte.org/cascade/SCTE%20CNR%20vs%20SNR.pdf will be useful

2. Is your question the following: If we pass a received symbol with a noise of xdB
through a splitter having resistive loss of 1dB, how much will be the resultant noise
power. Did I interpret correctly?
I would say the resultant noise will be 10*log10(10^(-x/10) + 10^(-1/10))
Did the linear addition of two noise components, do you agree?

Regards,
Krishna

Reply
Rajesh Kher says:

April 3, 2008 at 12:01 pm

Krishna

1. The symbol error rate should also be related to Modulation Error Rate (MER). Since MER is what most Meters measure a correlation will be most helpful.

2. Most of the signal are distributed over a network. What one also needs is the MER impairments due to network elements.

For example if the signal is split into two parts. That splitter will have an insertion loss which will have a spatial component. For example a 2 way lossless splitter will reduce the signal and noise by 3 dB and hence will not impair the signal. However there is an associated resistive loss aboput 1 dB. How much will the MER deteoriate?

Regards

Rajesh Kher

Reply
Krishna says:

March 28, 2008 at 10:20 am

@mahesh:
Yeah… i agree. it takes a bit of time to get used to 20*log10() or 10*log10() … 🙂
My rule of thumb is as follows:
(a) for voltage signals use 20*log10()
(b) for power signals use 10*log10().

In the code, as you observed,
(a) for scaling noise, which is a voltage signal, I used the 20*log10().
(b) for finding the theoretical symbol error rate, the Es/No in dB (which is signal power by noise power), the conversion used is 10*log10().

Hope this helps.

Regards,
Krishna

Reply
mahesh says:

March 28, 2008 at 8:25 am

In your simulation, you have added noise voltage/level i.e taken 20log10(ns) whereas for calculating theoretical BER , you have used noise power. (10log10(ns^2)). I always tend to get confused between these. Do you know any study material clarifying these and some examples?
thanks
mahesh

Reply
Krishna says:

March 28, 2008 at 4:30 am

@mahesh: true. Infact in the transmit and receive chain, there are various sources which can introduce DC (LO leakage etc). Given so, it makes it difficult to recover information present on DC subcarrier.

Reply
mahesh says:

March 27, 2008 at 9:15 am

@ Rajesh, How do you get sync information from DC?

Reply
Krishna says:

March 21, 2008 at 11:51 am

@ rajesh:
well, typically dc subcarrier is not used (as far as I have seen). Can you please point to the section in the spec (and the spec version).

Thanks,
krishna

Reply
RAJESH.N says:

March 20, 2008 at 11:06 am

hi i am doing M.E wireless tech. now i am doing my final sem project is implimentation of synchronization algorithem for wimax ieee 802.16e standard.

in wimax ieee802.16e standard
for synchronization purpose they are using preamble .in that preamble they using dc subcarrier .why they are using dc subcarrier… please if any one knows clarify my dought…its urgent

with regards
rajesh neelakandan

Reply