Skip to main content

MATLAB Code for QPSK Modulation and Demodulation


 

Quadrature Phase Shift Keying (QPSK) is a digital modulation scheme that conveys two bits per symbol by changing the phase of the carrier signal. Each pair of bits is mapped to one of four possible phase shifts: 0°, 90°, 180°, or 270°

00  ===> 0 degree phase shift of carrier signal

01  ===> 90 degree

11  ===> 180 degree

10  ===> 270 degree

 

MATLAB Script


 

Result

data
1
0
2
2
0
2
1
.
.
.

modData
-1.00000000000000 + 1.22464679914735e-16i
-1.83697019872103e-16 - 1.00000000000000i
1.00000000000000 + 0.00000000000000i
6.12323399573677e-17 + 1.00000000000000i
6.12323399573677e-17 + 1.00000000000000i
-1.83697019872103e-16 - 1.00000000000000i
1.00000000000000 + 0.00000000000000i
-1.00000000000000 + 1.22464679914735e-16i
6.12323399573677e-17 + 1.00000000000000i
-1.83697019872103e-16 - 1.00000000000000i
.
.
.


Fig 1: Constellation Diagram of Transmitted QPSK (or 4 PSK) Signal




Fig 2: Constellation Diagram of Received QPSK (or 4 PSK) Signal through Noisy Channel


Further Reading

  1. MATLAB Code for Constellation Diagram of QPSK
  2. Constellation Diagram of BPSK

QPSK (Passband) Signal Generation 

Choose QPSK as the modulation type in the simulation to visualize the QPSK signal.

Try the Online Simulator for QPSK Signal Generator (Click here)


Spectral Efficiency in QPSK

M = 4 Modulation

In Quadrature Phase Shift Keying (QPSK), the system transmits two bits per symbol. This allows for a significantly higher Spectral Efficiency (η) compared to binary modulation, doubling the data rate within the same frequency allocation.

Fundamental Measurement

η = Rb / B

For QPSK (M = 4), each symbol carries 2 bits. Therefore:

Rb = 2 × Rs

Where Rb is Bit Rate and Rs is Symbol Rate

Efficiency Benchmarks (η)

Nyquist Limit (B = Rs) 2.0 bps/Hz
Null-to-Null (B = 2Rs) 1.0 bps/Hz

The QPSK Calculation Logic

The measurement of efficiency in a QPSK system is driven by three core technical factors:

  1. 1
    Bit-to-Symbol Mapping: Because QPSK utilizes four distinct phases, it maps bit pairs (00, 01, 10, 11) to a single symbol. This means the Symbol Rate (Rs) is exactly half of the Bit Rate (Rb).
  2. 2
    Bandwidth Efficiency: The Required Bandwidth (B) of a signal is physically determined by the Symbol Rate. Since QPSK has a lower symbol rate than BPSK for the same bit rate, it consumes 50% less bandwidth to move the same amount of data.
  3. 3
The Sweet Spot Result: Under identical pulse-shaping conditions, QPSK provides an Efficiency (η) that is exactly double that of BPSK, reaching a theoretical max of 2 bits/s/Hz.

Why Engineers Choose QPSK

While 16-QAM or 64-QAM offer even higher efficiency, they are extremely fragile in noisy environments. QPSK is the industry's robust middle-ground, used extensively in Satellite Links (DVB-S) and LTE Control Channels where both speed and reliability are critical.

Read more: about Spectral Efficiency for various Modulation Schemes

QPSK Performance Analysis: Bit and Symbol Error Rates

Quadrature Phase Shift Keying (QPSK) is widely used due to its ability to carry twice the data rate of BPSK within the same bandwidth while maintaining the same bit error performance.

1. QPSK Bit Error Rate (BER)

For Gray-coded QPSK, the probability of bit error is identical to BPSK, as the quadrature components are independent:

\[ P_{b, QPSK} = Q\left(\sqrt{\frac{2E_b}{N_0}}\right) \]

2. QPSK Symbol Error Rate (SER)

A QPSK symbol consists of two bits. An error occurs if either (or both) bits are detected incorrectly. The approximation for SER is:

\[ P_{s, QPSK} \approx 2Q\left(\sqrt{\frac{2E_b}{N_0}}\right) \]

Note: The exact expression is \( P_s = 1 - [1 - P_b]^2 \).


Comparative BER Performance

How QPSK compares to other modulation formats in terms of \(P_b\):

BPSK (Same as QPSK)

\[ P_b = Q\left(\sqrt{\frac{2E_b}{N_0}}\right) \]

8-PSK (Higher Order)

\[ P_b \approx \frac{2}{3} Q\left( \sqrt{\frac{2E_s}{N_0}} \sin\left(\frac{\pi}{8}\right) \right) \]

16-QAM (Gray Coded)

\[ P_b \approx \frac{3}{4} Q\left( \sqrt{\frac{4}{5}\frac{E_b}{N_0}} \right) \]

64-QAM / 256-QAM

\[ P_{b, 64} \approx \frac{7}{12} Q\left( \sqrt{\frac{2}{7}\frac{E_b}{N_0}} \right), \quad P_{b, 256} \approx \frac{15}{32} Q\left( \sqrt{\frac{8}{85}\frac{E_b}{N_0}} \right) \]

Comparative Symbol Error Rate (SER)

Performance benchmarks for symbol detection:

8-PSK Approximation

\[ P_s \approx 2Q\left( \sqrt{\frac{2E_s}{N_0}} \sin\left(\frac{\pi}{8}\right) \right) \]

16-QAM Approximation

\[ P_s \approx 3Q\left( \sqrt{\frac{4}{5}\frac{E_b}{N_0}} \right) \]

Higher Order QAM (64 & 256)

\[ P_{s, 64} \approx \frac{7}{2} Q\left( \sqrt{\frac{2}{7}\frac{E_b}{N_0}} \right), \quad P_{s, 256} \approx \frac{15}{4} Q\left( \sqrt{\frac{8}{85}\frac{E_b}{N_0}} \right) \]
Technical Deep Dive: Q-function & SNR Calculation


Contact Us

Name

Email *

Message *

Popular Posts

UGC NET Electronic Science Previous Year Question Papers with Solutions

Home / Engineering & Other Exams / UGC NET 2026 PYQ ⬇️ Download Papers and Solutions 📋 Exam Pattern 💡 Preparation Tips ❓ FAQs 📊 Exam Highlights: Electronic Science (88) Feature Details Junior Research Fellowship (JRF) ₹37,000 + HRA per month Eligibility M.Sc/M.Tech in Electronics (55%) Validity of Certificate JRF (3 Years) | Lectureship (Lifetime) 📥 Download UGC NET Electronics PDFs Complete collection of previous year question papers, answer keys and explanations for Subject Code 88. Start Downloading 📂 View All Question Papers June 2025 - Question Paper Download PDF June 2025 - Solved Paper + Explanation ...

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...

Design of CMOS XOR/XNOR Gates

Design of CMOS XOR/XNOR Gates The semiconductor industry has experienced rapid integration of multimedia applications into mobile electronics, leading to very high integration density in CMOS VLSI. As operating frequencies increase, power consumption, speed, silicon area, and reliability become critical considerations. The XOR-XNOR circuits are fundamental building blocks in arithmetic circuits (Full Adders, Multipliers), compressors, comparators, parity checkers, code converters, error-detecting/correcting codes, and phase detectors. Their performance directly impacts the complex circuits they are used in. Design goals include full output voltage swing, low power consumption, reduced transistor count, minimal delay, and simultaneous non-skewed outputs. Static Logic (Static CMOS) Stat...

Online Simulator for ASK, FSK, and PSK Signal Generation

Interactive Digital Signal Processing (DSP) Tutorial and Simulator for ASK, FSK, and BPSK modulation techniques. Try our new Digital Signal Processing Simulator!   •   Interactive ASK, FSK, and BPSK tools updated for 2025. Start Now Digital Modulation Visualizer: ASK, FSK, & BPSK Simulator Learn and visualize binary modulation techniques (ASK, FSK, BPSK) in real-time with adjustable carrier and sampling parameters. Perfect for DSP students and engineers. 📡 ASK Simulator 📶 FSK Simulator 🎚️ BPSK Simulator 📚 More Topics ASK Modulator FSK Modulator BPSK Modulator More Topics 1. ASK (Amplitude Shift Keying) Simulat...

MATLAB Code for MUSIC

  MATLAB Code clc; clear; close all ; %% Step 1: Define Parameters M = 8; % Number of array sensors d = 0.5; % Sensor spacing (lambda/2) K = 2; % Number of signals N = 200; % Number of snapshots theta = [-20 30]; % True signal angles (degrees) SNR = 10; % Signal-to-noise ratio (dB) fprintf( 'Step 1: Parameters Initialized\n' ); %% Step 2: Generate Signal Sources t = 1:N; s1 = exp(1j*2*pi*0.05*t); s2 = exp(1j*2*pi*0.1*t); S = [s1; s2]; figure; plot(real(S(1,:))) title( 'Signal 1 (Real Part)' ) xlabel( 'Samples' ) ylabel( 'Amplitude' ) figure; plot(real(S(2,:))) title( 'Signal 2 (Real Part)' ) xlabel( 'Samples' ) ylabel( 'Amplitude' ) fprintf( 'Step 2: Source Signals Generated\n' ); %% Step 3: Construct Steering Matrix A = zeros(M,K); for k = 1:K A(:,k) = exp(-1j*2*pi*d*(0:M-1)'*sin(theta(k)*pi/180)); end fprintf( 'Step 3: Steering Matr...

Interactive Eye Diagram & ISI Simulator: Visualizing Signal Integrity

Interactive Eye Diagram & Noise Margin Filter Roll-off (α) α = 0.4 Noise Level Bitstream (TX) Stage 1: Components (Click Legend to Hide/Show) Stage 2: Composite Signal Stage 3: The Eye (Annotated Analysis) Simulation Workflow 1. Pulse Mapping: Each bit a[k] is mapped to a pulse shape p(t) . We use the Raised Cosine filter, which is the standard for bandwidth-limited communication. Clicking the legend above removes a specific bit's contribution to show how it affects the neighbors. 2. Linear Superposition: The total signal x(t) is the sum of all individual pulses shifted by the symbol period T . This "Combined Waveform" shows how pulses "bleed" into each other, cre...

Direction of Arrival (DoA) Online Simulator (using MUSIC)

Interactive DOA Simulator X-axis XY angle (deg): 45 XZ angle (deg): 30 Noise: 0.05 Y-axis XY angle (deg): 60 YZ angle (deg): 45 Noise: 0.05 Z-axis XZ angle (deg): 60 YZ angle (deg): 30 Noise: 0.05 Estimated DOA (deg): 0 Simulation Workflow and Mathematical Background This simulator demonstrates Direction of Arrival (DOA) estimation using three-axis sensor signals (X, Y, Z), Maximal Ratio Combining (MRC) , and the MUSIC algorithm . It allows interactive control of signal angles and noise for teaching purposes. 1. Signal Generation A pure sinewave signal of frequency f is projected onto three axes using user-defined angles in different planes: X-axis: θ XY , θ XZ Y-axis: θ XY , θ YZ Z-axis: θ XZ , θ YZ Mathematically, for each time sample t : x(t) = s(t) * cos(θ_xy_x) * cos(θ_xz_x) + n_x(t) y(t) = s(t) * sin(θ_xy_y) * cos(θ_yz_y) + n_y(t) z(t) = s(t) * sin(θ_xz_z) * sin(θ_yz_z) + n_z(t) wh...