Skip to main content
Login Home Wireless Communication Modulation MATLAB Beamforming Project Ideas MIMO Computer Networks

OFDM in MATLAB


 

MATLAB Script

% The code is written by SalimWireless.Com

1. Initialization

clc;
clear all;
close all;


2. Generate Random Bits

% Generate random bits
numBits = 100;
bits = randi([0, 1], 1, numBits);


3. Define Parameters

% Define parameters
numSubcarriers = 4; % Number of subcarriers
numPilotSymbols = 3; % Number of pilot symbols
cpLength = ceil(numBits / 4); % Length of cyclic prefix (one-fourth of the data length)


4. Add Cyclic Prefix

% Add cyclic prefix
dataWithCP = [bits(end - cpLength + 1:end), bits];


5. Insert Pilot Symbols

% Insert pilot symbols
pilotSymbols = ones(1, numPilotSymbols); % Example pilot symbols (could be any pattern)
dataWithPilots = [pilotSymbols, dataWithCP];

 

6. Perform OFDM Modulation (IFFT)

% Perform OFDM modulation (IFFT)
dataMatrix = reshape(dataWithPilots, numSubcarriers, []);
ofdmSignal = ifft(dataMatrix, numSubcarriers);
ofdmSignal = reshape(ofdmSignal, 1, []);


7. Display the Generated Data

% Display the generated data
disp("Original Bits:");
disp(bits);
disp("Data with Cyclic Prefix and Pilots:");
disp(dataWithPilots);
disp("OFDM Signal:");
disp(ofdmSignal);

%%%%%%%%%%%%%%%%%%%%%%%%%%% Demodulation %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


8. Demodulation

% Perform FFT on the received signal
%ofdmSignal = awgn(ofdmSignal, 1000);
ofdmSignal = reshape(ofdmSignal, numSubcarriers, []);
rxSignal = fft(ofdmSignal, numSubcarriers);
%rxSignal = [rxSignal(1,:) rxSignal(2,:) rxSignal(3,:) rxSignal(4,:)];


9. Remove Cyclic Prefix

% Remove cyclic prefix
rxSignalNoCP = rxSignal(cpLength + 1:end);


10. Extract Data Symbols and Discard Pilot Symbols

% Extract data symbols and discard pilot symbols
dataSymbols = rxSignalNoCP(numPilotSymbols + 1:end);


11. Demodulate the Symbols Using Thresholding

% Demodulate the symbols using thresholding
threshold = 0;
demodulatedBits = (real(dataSymbols) > threshold);


12. Plot the Original and Received Bits

figure(1)
stem(bits);
legend("Original Information Bits")

figure(2)
stem(demodulatedBits);
legend("Received Bits")

Output

 

 
 
Fig 1: Original Information Bits
 
 
 
 
 
Fig 2: OFDM Signal
 
 
 
 
 
Fig 3: Received Demodulated Bits

 

Copy the MATLAB Code above from here

 

 

MATLAB Code for OFDM using QPSK

% The code is written by SalimWireless.Com
clc;
clear all;
close all;

% Generate random bits (must be even for QPSK)
numBits = 20;
if mod(numBits, 2) ~= 0
numBits = numBits + 1; % Make even
end
bits = randi([0, 1], 1, numBits);

% QPSK Modulation (2 bits per symbol)
bitPairs = reshape(bits, 2, []).';
qpskSymbols = (1/sqrt(2)) * ((2*bitPairs(:,1)-1) + 1j*(2*bitPairs(:,2)-1)); % Gray-coded

% Parameters
numSubcarriers = 4; % Number of OFDM subcarriers
numPilotSymbols = 3; % Number of pilot symbols
cpLength = ceil(length(qpskSymbols) / 4); % Cyclic prefix length

% Insert pilot symbols
pilotSymbols = ones(1, numPilotSymbols); % Example pilot symbols (BPSK pilots)
dataWithPilots = [pilotSymbols, qpskSymbols.'];

% Add cyclic prefix
dataWithCP = [dataWithPilots(end - cpLength + 1:end), dataWithPilots];

% Reshape and perform IFFT (OFDM modulation)
dataMatrix = reshape(dataWithCP, numSubcarriers, []);
ofdmSignal = ifft(dataMatrix, numSubcarriers);
ofdmSignal1 = reshape(ofdmSignal, 1, []);

% Display
disp("Original Bits:");
disp(bits);
disp("QPSK Symbols:");
disp(qpskSymbols.');
disp("Data with CP and Pilots:");
disp(dataWithCP);
disp("OFDM Signal:");
disp(ofdmSignal1);

%%%%%%%%%%%%%%%%%%%%%%%%%%% Demodulation %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Reshape back to subcarrier-wise blocks and FFT
ofdmRxMatrix = reshape(ofdmSignal1, numSubcarriers, []);
rxSignal = fft(ofdmRxMatrix, numSubcarriers);

% Remove cyclic prefix
rxSignal1D = reshape(rxSignal, 1, []);
rxSignalNoCP = rxSignal1D(cpLength + 1:end);

% Remove pilots
rxDataSymbols = rxSignalNoCP(numPilotSymbols + 1:end);

% QPSK Demodulation
demodBits = zeros(1, 2*length(rxDataSymbols));
demodBits(1:2:end) = real(rxDataSymbols) > 0;
demodBits(2:2:end) = imag(rxDataSymbols) > 0;

%%%%%%%%%%%%%%%%%%%%%%%%%%% Plotting %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

figure(1)
stem(bits);
title("Original Bits");
xlabel("Bit Index"); ylabel("Bit Value");
legend("Original Bits");

figure(2)
hReal = stem(real(ofdmSignal1), 'r', 'DisplayName', 'Real Part');
hold on;
hImag = stem(imag(ofdmSignal1), 'b', 'DisplayName', 'Imaginary Part');
set(hReal, 'Marker', 'o', 'LineWidth', 1.5);
set(hImag, 'Marker', 'x', 'LineWidth', 1.5);
grid on;
title('OFDM Signal (Time Domain)');
xlabel('Sample Index');
ylabel('Amplitude');
legend;
hold off;

figure(3)
stem(demodBits);
title("Demodulated Bits");
xlabel("Bit Index"); ylabel("Bit Value");
legend("Demodulated Bits");

% Optional: Calculate BER
numErrors = sum(bits ~= demodBits);
ber = numErrors / numBits;
fprintf("Bit Error Rate (BER): %.4f\n", ber);



Output


 

 
 
 
 

 
 
  
 

 
 
 

 

MATLAB Code for OFDM Subcarriers (using 16-QAM)

clc;
clear;
close all;

% OFDM System with 16-QAM and Cooley-Tukey FFT/IFFT

% Parameters
N = 64; % Number of OFDM subcarriers
M = 16; % Modulation order (16-QAM -> M = 16)
nSymbols = 100;% Number of OFDM symbols
Ncp = 16; % Length of cyclic prefix

% Generate random data for transmission (0 to M-1 for 16-QAM)
data = randi([0 M-1], nSymbols, N);

% 16-QAM modulation of the data using custom function
modData = zeros(nSymbols, N);
for i = 1:nSymbols
modData(i, :) = qammod(data(i, :), M);
end

% Perform IFFT using Cooley-Tukey to generate the time domain OFDM signal
ofdmTimeSignal = zeros(size(modData));
for i = 1:nSymbols
ofdmTimeSignal(i, :) = ifft(modData(i, :));
end

% Add cyclic prefix
cyclicPrefix = ofdmTimeSignal(:, end-Ncp+1:end); % Extract cyclic prefix
ofdmWithCP = [cyclicPrefix ofdmTimeSignal]; % Add cyclic prefix to the signal

%% Plot Subcarriers in Frequency Domain (before IFFT)
figure;
stem(0:N-1, abs(modData(100, :))); % Plot absolute value of the subcarriers for the first symbol
title('Subcarriers in Frequency Domain for 1st OFDM Symbol (Before IFFT)');
xlabel('Subcarrier Index');
ylabel('Magnitude');

%% Plot Time Domain OFDM Signal (after IFFT)
figure;
plot(real(ofdmTimeSignal(1, :))); % Plot real part of the OFDM time signal for the first symbol
title('OFDM Signal in Time Domain for 1st OFDM Symbol (Without CP)');
xlabel('Time Sample Index');
ylabel('Amplitude');

%% Plot Time Domain OFDM Signal with Cyclic Prefix
figure;
plot(real(ofdmWithCP(1, :))); % Plot real part of the OFDM time signal with CP for the first symbol
title('OFDM Signal in Time Domain for 1st OFDM Symbol (With Cyclic Prefix)');
xlabel('Time Sample Index');
ylabel('Amplitude');

%% Receiver Side - Remove Cyclic Prefix and Demodulate
% Remove cyclic prefix
receivedSignal = ofdmWithCP(:, Ncp+1:end); % Remove cyclic prefix

% Apply FFT using Cooley-Tukey to recover the received subcarriers (back to frequency domain)
receivedSubcarriers = zeros(size(receivedSignal));
for i = 1:nSymbols
receivedSubcarriers(i, :) = fft(receivedSignal(i, :));
end

% 16-QAM Demodulation of the received subcarriers using custom function
receivedData = zeros(nSymbols, N);
for i = 1:nSymbols
receivedData(i, :) = qamdemod(receivedSubcarriers(i, :), M);
end

% Calculate symbol errors
numErrors = sum(data(:) ~= receivedData(:));
fprintf('Number of symbol errors: %d\n', numErrors);

%% Plot Received Subcarriers in Frequency Domain (after FFT at the receiver)
figure;
stem(0:N-1, abs(receivedSubcarriers(100, :))); % Plot absolute value of received subcarriers for the first symbol
title('Received Subcarriers in Frequency Domain for 1st OFDM Symbol (After FFT)');
xlabel('Subcarrier Index');
ylabel('Magnitude');

%% Plot Transmitted Data Constellation (Before IFFT)
figure;
scatterplot(modData(1, :)); % Plot for the first OFDM symbol
title('Transmitted 16-QAM Symbols for 1st OFDM Symbol');
xlabel('In-phase');
ylabel('Quadrature');

%% Plot Received Data Constellation (After Demodulation)
receivedModData = qammod(receivedData(1, :), M); % Map back for plotting
figure;
scatterplot(receivedModData);
title('Received 16-QAM Symbols for 1st OFDM Symbol');
xlabel('In-phase');
ylabel('Quadrature');

 Output

 
















Copy the MATLAB code above from here

 

Read more about

[1] OFDM in details

[2] Structure of an OFDM packet

People are good at skipping over material they already know!

View Related Topics to







Admin & Author: Salim

s

  Website: www.salimwireless.com
  Interests: Signal Processing, Telecommunication, 5G Technology, Present & Future Wireless Technologies, Digital Signal Processing, Computer Networks, Millimeter Wave Band Channel, Web Development
  Seeking an opportunity in the Teaching or Electronics & Telecommunication domains.
  Possess M.Tech in Electronic Communication Systems.


Contact Us

Name

Email *

Message *

Popular Posts

Constellation Diagrams of ASK, PSK, and FSK

📘 Overview of Energy per Bit (Eb / N0) 🧮 Online Simulator for constellation diagrams of ASK, FSK, and PSK 🧮 Theory behind Constellation Diagrams of ASK, FSK, and PSK 🧮 MATLAB Codes for Constellation Diagrams of ASK, FSK, and PSK 📚 Further Reading 📂 Other Topics on Constellation Diagrams of ASK, PSK, and FSK ... 🧮 Simulator for constellation diagrams of m-ary PSK 🧮 Simulator for constellation diagrams of m-ary QAM BASK (Binary ASK) Modulation: Transmits one of two signals: 0 or -√Eb, where Eb​ is the energy per bit. These signals represent binary 0 and 1.    BFSK (Binary FSK) Modulation: Transmits one of two signals: +√Eb​ ( On the y-axis, the phase shift of 90 degrees with respect to the x-axis, which is also termed phase offset ) or √Eb (on x-axis), where Eb​ is the energy per bit. These signals represent binary 0 and 1.  BPSK (Binary PSK) Modulation: Transmits one of two signals...
Read more

Comparisons among ASK, PSK, and FSK | And the definitions of each

📘 Comparisons among ASK, FSK, and PSK 🧮 Online Simulator for calculating Bandwidth of ASK, FSK, and PSK 🧮 MATLAB Code for BER vs. SNR Analysis of ASK, FSK, and PSK 📚 Further Reading 📂 View Other Topics on Comparisons among ASK, PSK, and FSK ... 🧮 Comparisons of Noise Sensitivity, Bandwidth, Complexity, etc. 🧮 MATLAB Code for Constellation Diagrams of ASK, FSK, and PSK 🧮 Online Simulator for ASK, FSK, and PSK Generation 🧮 Online Simulator for ASK, FSK, and PSK Constellation 🧮 Some Questions and Answers Modulation ASK, FSK & PSK Constellation MATLAB Simulink MATLAB Code Comparisons among ASK, PSK, and FSK    Comparisons among ASK, PSK, and FSK Comparison among ASK, FSK, and PSK Parameters ASK FSK PSK Variable Characteristics Amplitude Frequency ...
Read more

BER vs SNR for M-ary QAM, M-ary PSK, QPSK, BPSK, ...

📘 Overview of BER and SNR 🧮 Online Simulator for BER calculation of m-ary QAM and m-ary PSK 🧮 MATLAB Code for BER calculation of M-ary QAM, M-ary PSK, QPSK, BPSK, ... 📚 Further Reading 📂 View Other Topics on M-ary QAM, M-ary PSK, QPSK ... 🧮 Online Simulator for Constellation Diagram of m-ary QAM 🧮 Online Simulator for Constellation Diagram of m-ary PSK 🧮 MATLAB Code for BER calculation of ASK, FSK, and PSK 🧮 MATLAB Code for BER calculation of Alamouti Scheme 🧮 Different approaches to calculate BER vs SNR What is Bit Error Rate (BER)? The abbreviation BER stands for bit error rate, which indicates how many corrupted bits are received (after the demodulation process) compared to the total number of bits sent in a communication process. It is defined as,  In mathematics, BER = (number of bits received in error / total number of transmitted bits)  On the other hand, SNR ...
Read more

Constellation Diagrams of M-ary QAM | M-ary Modulation

📘 Overview of QAM 🧮 MATLAB Code for m-ary QAM (4-QAM, 16-QAM, 32-QAM, ...) 🧮 Online Simulator for M-ary QAM Constellations 📚 Further Reading 📂 Other Topics on Constellation Diagrams of QAM configurations ... 🧮 MATLAB Code for 4-QAM 🧮 MATLAB Code for 16-QAM 🧮 MATLAB Code for m-ary QAM (4-QAM, 16-QAM, 32-QAM, ...) 🧮 Simulator for constellation diagrams of m-ary PSK 🧮 Simulator for constellation diagrams of m-ary QAM 🧮 Overview of Energy per Bit (Eb / N0) 🧮 Online Simulator for constellation diagrams of ASK, FSK, and PSK 🧮 Theory behind Constellation Diagrams of ASK, FSK, and PSK 🧮 MATLAB Codes for Constellation Diagrams of ASK, FSK, and PSK QAM Unlike M-ary PSK, where the signal is modulated with diffe...
Read more

MATLAB Code for Zero-Forcing (ZF) Beamforming in 4×4 MIMO Systems

MATLAB Code for Zero-Forcing (ZF) Beamforming in 4×4 MIMO Systems clc; clear; close all; %% Parameters Nt = 4; % Transmit antennas Nr = 4; % Receive antennas (must be >= Nt for ZFBF) numBits = 1e4; % Number of bits per stream SNRdB = 0; % SNR in dB numRuns = 100; % Number of independent runs for averaging %% Precompute noise standard deviation noiseSigma = 10^(-SNRdB / 20); %% Accumulator for total errors totalErrors = 0; for run = 1:numRuns % Generate random bits: [4 x 10000] bits = randi([0 1], Nt, numBits); % BPSK modulation: 0 → +1, 1 → -1 txSymbols = 1 - 2 * bits; % Rayleigh channel matrix: [4 x 4] H = (randn(Nr, Nt) + 1j * randn(Nr, Nt)) / sqrt(2); %% === Zero Forcing Beamforming at Transmitter === W_zf = pinv(H); % Precoding matrix: [Nt x Nr] txPrecoded = W_zf * txSymbols; % Apply ZF precoding % Normalize transmit power (optional but useful) txPrecoded = txPrecoded / sqrt(mean(abs(txPrecoded(:)).^2)); %% Channel transmission with AWGN noise = noiseSigma * (randn(...
Read more

Theoretical BER vs SNR for BPSK

Let's simplify the explanation for the theoretical Bit Error Rate (BER) versus Signal-to-Noise Ratio (SNR) for Binary Phase Shift Keying (BPSK) in an Additive White Gaussian Noise (AWGN) channel.  Key Points Fig 1: Constellation Diagrams of BASK, BFSK, and BPSK [↗] BPSK Modulation: Transmits one of two signals: +√Eb ​ or -√Eb , where Eb​ is the energy per bit. These signals represent binary 0 and 1 . AWGN Channel: The channel adds Gaussian noise with zero mean and variance N0/2 (where N0 ​ is the noise power spectral density). Receiver Decision: The receiver decides if the received signal is closer to +√Eb​ (for bit 0) or -√Eb​ (for bit 1) . Bit Error Rate (BER) The probability of error (BER) for BPSK is given by a function called the Q-function. The Q-function Q(x) measures the tail probability of the normal distribution, i.e., the probability that a Gaussian random variable exceeds a certain value x.  Understanding the Q...
Read more

Constellation Diagram of FSK in Detail

📘 Overview 🧮 Simulator for constellation diagram of FSK 🧮 Theory 🧮 MATLAB Code 📚 Further Reading   Binary bits '0' and '1' can be mapped to 'j' and '1' to '1', respectively, for Baseband Binary Frequency Shift Keying (BFSK) . Signals are in phase here. These bits can be mapped into baseband representation for a number of uses, including power spectral density (PSD) calculations. For passband BFSK transmission, we can modulate signal 'j' with a lower carrier frequency and signal '1' with a higher carrier frequency while transmitting over a wireless channel. Let's assume we are transmitting carrier signal fc1 for the transmission of binary bit '1' and carrier signal fc2 for the transmission of binary bit '0'. Simulator for 2-FSK Constellation Diagram Simulator for 2-FSK Constellation Diagram SNR (dB): ...
Read more

MATLAB Code for ASK, FSK, and PSK

📘 Overview & Theory 🧮 MATLAB Code for ASK 🧮 MATLAB Code for FSK 🧮 MATLAB Code for PSK 🧮 Simulator for binary ASK, FSK, and PSK Modulations 📚 Further Reading ASK, FSK & PSK HomePage MATLAB Code MATLAB Code for Amplitude Shift Keying (ASK) % The code is written by SalimWireless.Com % Clear previous data and plots clc; clear all; close all; % Parameters Tb = 1; % Bit duration fc = 10; % Carrier frequency N = 10; % Number of bits % Generate carrier signal t = 0:Tb/100:1; carrier_signal = sqrt(2/Tb) * sin(2*pi*fc*t); % Generate message signal rng(10); % Set random seed for reproducibility binary_data = rand(1, N); % Generate random binary data t_start = 0; t_end = Tb; for i = 1:N t = [t_start:0.01:t_end]; % Generate message signal if binary_data(i) > 0.5 binary_data(i) = 1; message_signal = ones(1, length(t)); ...
Read more