OFDM and Cyclic Prefix – a handson demonstration using Matlab

(9 votes, average: 4.44 out of 5)

This tutorial is a part of the website http://www.gaussianwaves.com and created under Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License. You must credit the author in your work if you remix, tweak, and build upon the work below. More such tricks can be found in the ebook : Simulation of Digital communication systems using Matlab by Mathuranathan Viswanathan}

Introduction

A OFDM transceiver architecture is typically implemented using IFFT and FFT blocks. In a OFDM transmitter, the modulated information symbols are assigned to invidual subcarriers, forming a single OFDM symbol in frequency domain. The OFDM symbol in frequency domain is converted to time domain using IFFT block. The time domain samples are then transmitted across a channel that is modeled using Channel Impulse Response (CIR). On the otherhand, the receiver converts the time domain samples of the received OFDM symbols and converts them to frequency domain for further demodulation of the information in the individual carriers.

In reality, a multipath channel or any channel in nature acts as a $$linear$$ filter on the transmitted OFDM symbols. Mathematically, the transmitted OFDM symbol (denoted as $$s[n]$$) gets $$linearly$$ convolved with the CIR $$h[n]$$ and gets corrupted with additive white gaussian noise – designated as $$w[n]$$. Denoting $$linear\;convolution$$ as ‘$$\ast$$’ and using ‘$$N$$’-point DFT, received signal in discrete-time can be represented as

$$r[n] = h[n] \ast s[n] + w[n] \;\;\; (\ast \rightarrow linear\;convolution ) \;\;\;\;\;\; (1)$$

The idea behind using OFDM is to combat frequency selective fading , where different frequency components of a transmitted signal can undergo different levels of fading. The OFDM divides the available spectrum in to small chunks called sub-channels. Within the subchannels, the fading experienced by the individual modulated symbol can be considered flat. This opens-up the possibility of using a simple frequency domain equalizer to neutralize the channel effects in the individual subchannels.

Circular convolution and designing a simple frequency domain equalizer

From one of the DFT properties, we know that the $$circular\;convolution$$ of two sequences in time domain is equivalent to the multiplication of their individual responses in frequency domain.

Let $$s[n]$$ and $$h[n]$$ are two sequences of length $$N$$ with their DFTs denoted as $$S[k]$$ and $$H[k]$$ respectively. Denoting $$circular\;convolution$$ as $$\circledast$$,

$$h[n] \circledast s[n] \;\;\;\; \equiv H[k] S[k] \;\;\;\;\;\; (2)$$

If we ignore channel noise in the OFDM transmission, the received signal is written as
$$r[n] = h[n] \ast s[n] \;\;\;\;\;\;\;\;\; (\ast \rightarrow linear\;convolution ) \;\;\;\;\;\; (3)$$

Note that the channel performs a linear convolution operation on the transmitted signal. Instead, if the channel performs circular convolution, then the equation would have been
$$r[n] = h[n] \circledast s[n] \;\;\;\;\;\;\;\;\; (\circledast \rightarrow circular\;convolution ) \;\;\;\;\;\; (4)$$

Then, by applying the DFT property given in equation $$2$$,
$$r[n] = h[n] \circledast s[n] \;\;\;\; \equiv R[k] = H[k] S[k] \;\;\;\;\;\; (5)$$

As a consequence, the channel effects can be neutralized at the receiver using a simple frequency domain equalizer (actually this is a zero-forcing equalizer when viewed in time domain) that just inverts the $$estimated$$ channel response and multiplies it with the frequency response of the received signal to obtain the estimates of the OFDM symbols in the subcarriers as

$$\hat{S}[k] = \frac{R[k]}{H[k]} \;\;\;\;\;\; (6)$$

Demonstrating the role of Cyclic Prefix

The simple frequency domain equalizer shown in equation$$6$$ is possible only if the channel performs circular convolution. But in nature, all channels perform linear convolution. The linear convolution can be converted into circular convolution by adding Cyclic Prefix (CP) in the OFDM architecture. Let’s understand this by demonstration.

To simply stuffs, we will create two randomized vectors for $$s[n]$$ and $$h[n]$$. $$s[n]$$ is of length $$8$$, the channel impulse response $$h[n]$$ is of length $$3$$ and we intend to use $$N=8$$-point DFT/IDFT wherever applicable.

Now, convolve the vectors $$\mathbf{s}$$ and $$\mathbf{h}$$ linearly and circularly. The outputs are plotted in Figure $$1$$. We note that the linear convolution and the circular convolution does not yield identical results.

Let’s append a cyclic prefix to the signal $$\mathbf{s}$$ by copying last $$N_{cp}$$ symbols from $$\mathbf{s}$$ and pasting it to its front. Since the channel delay is 3 symbols (CIR of length 3), we need to add \emph{atleast} 2 CP symbols.

Lets assume that we send the cyclic-prefixed OFDM symbol $$\mathbf{s_{cp}}$$ through the channel which perform linear filtering

Compare the outputs due to cir_s_h $$=\mathbf{s} \circledast \mathbf{h}$$ and lin_scp_h $$=\mathbf{s}_{cp} \ast \mathbf{h}$$. We can immediately recognize that that the middle section of lin_scp_h is exactly same as the cir_s_h vector. This is shown in the figure $$2$$.

So far we have faked the channel to perform circular convolution by adding a cyclic extension to the OFDM symbol. We are not done yet. At the receiver, we are only interested in the the middle section of the lin_scp_h which is the channel output. The first block in the OFDM receiver removes the additional symbols from the front and back of the channel output. This becomes the received symbol $$\mathbf{r}$$ after the removal of cyclic prefix in the receiver.

Verifying DFT property

The DFT property given in equation $$5$$ can be re-written as
\begin{align} r[n] & = IDFT \left\{ R[k] \right\} \\ & = IDFT \left\{H[k] S[k] \right\} \\ & = IDFT \left\{DFT \left( h[n] \right) DFT \left( s[n] \right) \right\} \\ & = h[n] \circledast s[n] \end{align}\;\;\;\;\;\; (7)

To verify this in Matlab, take N-point DFTs of the received signal and CIR. Then, we can see that the IDFT of product of DFTs of $$\mathbf{s}$$ and $$\mathbf{h}$$ will be equal to the N-point circular convolution of $$\mathbf{s}$$ and $$\mathbf{h}$$

Recommended Signal Processing Books

• SmartGeek

It is great article. Thank you very much.

• Thanks

Sir, actually I have a query. In your matlab book, the Rayleigh channel generation has been done taking the Jakes’s model into consideration. But the Rician channel generation has been done in a static scenario, unlike the Rayleigh case. It would have been really helpful for me if you could do the same (i.e., scenarios where Rician channel is generated taking into account the doppler effect and the inter-sampling duration) and upload a post regarding the same or share the Matlab code.

• Karan Arya

hello sir

I want your help. i want matlab code of ber vs snr plot with bpsk,qpsk and 16qam modulation with mrc, egc and sc in mmse sic receiver

• I don’t have the code for this specific scenario. May be I will reserve it for the future.

• Karan Arya

Sir i have code of sic receiver but i dont know how to add mrc,sc,egc and 16qam,gmsk,qpsk

• Karan Arya

i have the code but i dont know how to add gmsk,qpsk, 16 qam and selection combing, equal gain combining