Skip to main content

Menu

Home 📘 Wireless Communication (Main Page) Constellation Diagrams BER vs SNR GMSK MSK Channel Impulse Response Modulation Schemes MIMO Technology OFDM MATLAB Projects Contact
Login Try App

MATLAB Code for Constellation Diagram of QAM configurations such as 4, 8, 16, 32, 64, 128, and 256-QAM


 

One of the best-performing modulation techniques is QAM [↗]. Here, we modulate the symbols by varying the carrier signal's amplitude and phase in response to the variation in the message signal (or voltage variation). So, we may say that QAM is a combination of phase and amplitude modulation. Additionally, it performs better than ASK or PSK [↗]. In fact, any constellation for any type of modulation, signal set (or, symbols) is structured in a way that prevents them from interacting further by being distinct by phase, amplitude, or frequency.


MATLAB Script (for 4-QAM)

% This code is written by SalimWirelss.Com
% This is an example of 4-QAM. Here constellation size is 4
% or total number of symbols/signals is 4
% We need 2 bits once to represent four constellation points
% QAM modulation is the combination of Amplitude modulation plus
% Phase Modulation. We map the decimal value of the input symbols, i.e.,
% 00, 01, 10, 11 to 1 + 1i, -1 + 1i, 1 - 1i, and -1 - 1i, respectively.


clc;clear all;close all;

M = 4; % Number of levels after quantization / size of signal constellation

k = log2(M); % Number of bits per symbol

rng(10) %assaining the value of seed integer

N =10000; % Number of bits to process

InputBits = randi([0 1],1,N); % Generating randon bits

InputSymbol_matrix = reshape(InputBits,length(InputBits)/k,k); % Reshape data into binary k-tuples, k = log2(M)

InputSymbols_decimal = bi2de(InputSymbol_matrix); % Convert binary to decimal

for n= 1:N/k

if InputSymbols_decimal(n)==0

QAM(n)= complex(1,1);

elseif InputSymbols_decimal(n)==1

QAM(n)= complex(-1,1);

elseif InputSymbols_decimal(n)==2

QAM(n)= complex(1,-1);

else

QAM(n)= complex(-1,-1);

end

end



%Transmission of 4QAM data over AWGN channel

snrdB = 10;

Y=awgn(QAM,snrdB); %received signal


%Threshold Detection

for n= 1:N/k

if (real(Y(n))>0 && imag(Y(n))>0)

Z(n)=complex(1,1);

elseif (real(Y(n))>0 && imag(Y(n))<0)

Z(n)=complex(1,-1);


elseif (real(Y(n))<0 && imag(Y(n))>0)

Z(n)=complex(-1,1);

else

Z(n)=complex(-1,-1);

end

end

figure(1)
scatter(real(QAM), imag(QAM))
xlim([-3, 3]);
ylim([-3, 3]);
legend('Transmitted Symbols')

figure(2)
scatter(real(Y), imag(Y))
xlim([-3, 3]);
ylim([-3, 3]);
legend('Received Symbols')
 

Output 

 
 
 
Fig 1: Constellation points of 4-QAM (Transmitted)


 
 
Fig 2: Constellation points of 4-QAM (Received)


Copy the MATLAB Code for 4-QAM


 

Another MATLAB Code (for 16-QAM)

%The code is developed by SalimWireless.Com

clc;
clear;
close all;

% Define parameters
M = 16; % Modulation order for 16-QAM
numSymbols = 10000; % Number of symbols to modulate

% Generate random data
data = randi([0 M-1], numSymbols, 1); % Ensure data is a column vector

% Modulate the data using 16-QAM
modData = qammod_custom(data, M);

snrdB = 15;
Y = awgn(modData,snrdB); %received signal

% Plot the constellation of the modulated signal
figure;
subplot(2,1,1);
scatter(real(modData), imag(modData), 'o');
grid on;
xlabel('In-phase');
ylabel('Quadrature');
title('Constellation Diagram of Modulated Signal (16-QAM)');
axis([-1.5 1.5 -1.5 1.5]); % Set axis limits for better visualization

subplot(2,1,2);
scatter(real(Y), imag(Y), 'o');
grid on;
xlabel('In-phase');
ylabel('Quadrature');
title('Constellation Diagram of Noisy Received Signal before demodulation');
axis([-1.5 1.5 -1.5 1.5]); % Set axis limits for better visualization

% Demodulate the received signal
receivedData = qamdemod_custom(modData, M);

% Ensure receivedData is a column vector for comparison
receivedData = receivedData(:);


% Custom 16-QAM Modulation Function
function modData = qammod_custom(data, M)
% QAMMOD_CUSTOM Modulate data using 16-QAM
% data - Column vector of integers (each element is between 0 and M-1)
% M - Modulation order (should be 16 for 16-QAM)

% Check if M is 16
if M ~= 16
error('This function is designed for 16-QAM modulation.');
end

% Define the 16-QAM constellation
constellation = [-3-3i, -3-1i, -1-3i, -1-1i, ...
-3+3i, -3+1i, -1+3i, -1+1i, ...
+3-3i, +3-1i, +1-3i, +1-1i, ...
+3+3i, +3+1i, +1+3i, +1+1i];

% Normalize constellation
constellation = constellation / sqrt(mean(abs(constellation).^2)); % Scale to unit average power

% Map data to constellation points
modData = constellation(data + 1);
end

% Custom 16-QAM Demodulation Function
function demodData = qamdemod_custom(modData, M)
% QAMDEMOD_CUSTOM Demodulate data using 16-QAM
% modData - Column vector of complex numbers (modulated symbols)
% M - Modulation order (should be 16 for 16-QAM)

% Check if M is 16
if M ~= 16
error('This function is designed for 16-QAM demodulation.');
end

% Define the 16-QAM constellation
constellation = [-3-3i, -3-1i, -1-3i, -1-1i, ...
-3+3i, -3+1i, -1+3i, -1+1i, ...
+3-3i, +3-1i, +1-3i, +1-1i, ...
+3+3i, +3+1i, +1+3i, +1+1i];

% Normalize constellation
constellation = constellation / sqrt(mean(abs(constellation).^2)); % Scale to unit average power

% Ensure modData is a column vector
modData = modData(:);

% Compute the distances from each modData point to all constellation points
numSymbols = length(modData);
numConstellations = length(constellation);
distances = zeros(numSymbols, numConstellations);
for k = 1:numConstellations
distances(:, k) = abs(modData - constellation(k)).^2;
end

% Find the closest constellation point for each modData point
[~, demodData] = min(distances, [], 2);

% Convert to zero-based index
demodData = demodData - 1;
end

Output  


 
 
 
 

Copy the MATLAB Code above from here (for 16-QAM)

 

MATLAB code for M-ary QAM (e.g., 4, 8, 16, 32, 64, 128, 256)

%The code is developed by SalimWireless.com
% M-ary QAM Modulation and Demodulation
clc;
clear;
close all;


% Parameters
M = 32; % Order of QAM (M-QAM)
N = 1000; % Number of symbols
SNR = 10; % Signal-to-Noise Ratio in dB


% Generate random data symbols
dataSymbols = randi([0 M-1], N, 1);


% Modulate using M-QAM
txSignal = qammod(dataSymbols, M);


% Add AWGN noise
rxSignal = awgn(txSignal, SNR, 'measured');


% Demodulate
demodulatedSymbols = qamdemod(rxSignal, M);


% Calculate symbol error rate
symbolErrors = sum(dataSymbols ~= demodulatedSymbols);
SER = symbolErrors / N;


% Display results
disp(['Symbol Error Rate (SER): ', num2str(SER)]);


% Plot constellation diagrams
figure;
subplot(2, 1, 1);
plot(real(txSignal), imag(txSignal), 'o');
grid on;
title('Transmitted Signal Constellation');
xlabel('In-Phase');
ylabel('Quadrature');


subplot(2, 1, 2);
plot(real(rxSignal), imag(rxSignal), 'o');
grid on;
title('Received Signal Constellation');
xlabel('In-Phase');
ylabel('Quadrature');

Output








Copy the MATLAB Code above from here (e.g., for QAM configurations such as 4, 8, 16, 32, 64, 128, and 256-QAM.)


MATLAB Code for BER vs SNR for 4-QAM, 16-QAM, 32-QAM, and so on

 
 
 
(Get MATLAB Code)


 Online Simulator for M-ary QAM Constellations (4-QAM, 16-QAM, 64-QAM, 256-QAM)

 

 
 
Also read about

Next>>

People are good at skipping over material they already know!

View Related Topics to







Contact Us

Name

Email *

Message *

Popular Posts

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. BER = (number of bits received in error) / (total number of tran...
Read more

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

Power Spectral Density Calculation Using FFT in MATLAB

📘 Overview 🧮 Steps to calculate the PSD of a signal 🧮 MATLAB Codes 📚 Further Reading Power spectral density (PSD) tells us how the power of a signal is distributed across different frequency components, whereas Fourier Magnitude gives you the amplitude (or strength) of each frequency component in the signal. Steps to calculate the PSD of a signal Firstly, calculate the first Fourier transform (FFT) of a signal Then, calculate the Fourier magnitude of the signal The power spectrum is the square of the Fourier magnitude To calculate power spectrum density (PSD), divide the power spectrum by the total number of samples and the frequency resolution. {Frequency resolution = (sampling frequency / total number of samples)} Sampling frequency (fs): The rate at which the continuous-time signal is sampled (in Hz). ...
Read more

OFDM Symbols and Subcarriers Explained

This article explains how OFDM (Orthogonal Frequency Division Multiplexing) symbols and subcarriers work. It covers modulation, mapping symbols to subcarriers, subcarrier frequency spacing, IFFT synthesis, cyclic prefix, and transmission. Step 1: Modulation First, modulate the input bitstream. For example, with 16-QAM , each group of 4 bits maps to one QAM symbol. Suppose we generate a sequence of QAM symbols: s0, s1, s2, s3, s4, s5, …, s63 Step 2: Mapping Symbols to Subcarriers Assume N sub = 8 subcarriers. Each OFDM symbol in the frequency domain contains 8 QAM symbols (one per subcarrier): Mapping (example) OFDM symbol 1 → s0, s1, s2, s3, s4, s5, s6, s7 OFDM symbol 2 → s8, s9, s10, s11, s12, s13, s14, s15 … OFDM sym...
Read more

Calculation of SNR from FFT bins in MATLAB

📘 Overview 🧮 MATLAB Code for Estimation of SNR from FFT bins of a Noisy Signal 🧮 MATLAB Code for Estimation of Signal-to-Noise Ratio from Power Spectral Density Using FFT and Kaiser Window Periodogram from real signal data 📚 Further Reading   Here, you can find the SNR of a received signal from periodogram / FFT bins using the Kaiser operator. The beta (β) parameter characterizes the Kaiser window, which controls the trade-off between the main lobe width and the side lobe level in the frequency domain. For that you should know the sampling rate of the signal.  The Kaiser window is a type of window function commonly used in signal processing, particularly for designing finite impulse response (FIR) filters and performing spectral analysis. It is a general-purpose window that allows for control over the trade-off between the main lobe width (frequency resolution) and side lobe levels (suppression of spectral leakage). The Kaiser window is defined...
Read more

MATLAB code for BER vs SNR for M-QAM, M-PSK, QPSk, BPSK, ...

🧮 MATLAB Code for BPSK, M-ary PSK, and M-ary QAM Together 🧮 MATLAB Code for M-ary QAM 🧮 MATLAB Code for M-ary PSK 📚 Further Reading MATLAB Script for BER vs. SNR for M-QAM, M-PSK, QPSK, BPSK % Written by Salim Wireless clc; clear; close all; num_symbols = 1e5; snr_db = -20:2:20; psk_orders = [2, 4, 8, 16, 32]; qam_orders = [4, 16, 64, 256]; ber_psk_results = zeros(length(psk_orders), length(snr_db)); ber_qam_results = zeros(length(qam_orders), length(snr_db)); for i = 1:length(psk_orders) psk_order = psk_orders(i); for j = 1:length(snr_db) data_symbols = randi([0, psk_order-1], 1, num_symbols); modulated_signal = pskmod(data_symbols, psk_order, pi/psk_order); received_signal = awgn(modulated_signal, snr_db(j), 'measured'); demodulated_symbols = pskdemod(received_signal, psk_order, pi/psk_order); ber_psk_results(i, j) = sum(data_symbols ~= demodulated_symbols) / num_symbols; end end for i...
Read more

Online Channel Impulse Response Simulator

  Fundamental Theory of Channel Impulse Response The fundamental theory behind the channel impulse response in wireless communication often involves complex exponential components such as: \( h(t) = \sum_{i=1}^{L} a_i \cdot \delta(t - \tau_i) \cdot e^{j\theta_i} \) Where: \( a_i \) is the amplitude of the \( i^{th} \) path \( \tau_i \) is the delay of the \( i^{th} \) path \( \theta_i \) is the phase shift (often due to Doppler effect, reflection, etc.) \( e^{j\theta_i} \) introduces a phase rotation (complex exponential) The convolution \( x(t) * h(t) \) gives the received signal So, instead of representing the channel with only real-valued amplitudes, each path can be more accurately modeled using a complex gain : \( h[n] = a_i \cdot e^{j\theta_i} \) 1. Simple Channel Impulse Response Simulator  (Here you can input only a unit impulse signal) Input Signal (Unit Impu...
Read more

Coherence Bandwidth and Coherence Time

🧮 Coherence Bandwidth 🧮 Coherence Time 🧮 MATLAB Code s 📚 Further Reading Coherence Bandwidth Coherence bandwidth is a concept in wireless communication and signal processing that relates to the frequency range over which a wireless channel remains approximately constant in terms of its characteristics. Coherence bandwidth is inversely related to the delay spread time (e.g., RMS delay spread). The coherence bandwidth is related to the delay spread of the channel, which is a measure of the time it takes for signals to traverse the channel due to multipath. The two are related by the following approximation: Coherence Bandwidth ≈ 1/(delay spread time) Or, Coherence Bandwidth ≈ 1/(root-mean-square delay spread time) (Coherence bandwidth in Hertz) For instance, if the root-mean-square delay spread is 500 ns (i.e., {1/(2*10^6)} seconds), the coherence bandwidth is approximately 2 MHz (1 / 500e-9) in ...
Read more