Skip to main content

MATLAB Code for Pulse Width Modulation (PWM) and Demodulation


 

MATLAB Code for Analog Pulse Width Modulation (PWM)

clc;
clear all;
close all;
fs=30; %frequency of the sawtooth signal
fm=3; %frequency of the message signal
sampling_frequency = 10e3;
a=0.5; % amplitide

t=0:(1/sampling_frequency):1; %sampling rate of 10kHz


sawtooth=2*a.*sawtooth(2*pi*fs*t); %generating a sawtooth wave


subplot(4,1,1);
plot(t,sawtooth); % plotting the sawtooth wave
title('Comparator Wave');

msg=a.*sin(2*pi*fm*t); %generating message wave

subplot(4,1,2);
plot(t,msg); %plotting the sine message wave
title('Message Signal');


for i=1:length(sawtooth)
if (msg(i)>=sawtooth(i))
pwm(i)=1; %is message signal amplitude at i th sample is greater than
%sawtooth wave amplitude at i th sample
else
pwm(i)=0;
end
end

subplot(4,1,3);
plot(t,pwm,'r');
title('PWM');
axis([0 1 0 1.1]); %to keep the pwm visible during plotting.

%% Demodulation
% Demodulation: Measure the pulse width to reconstruct the signal
demodulated_signal = zeros(size(msg));

for i = 1:length(pwm)-1
if pwm(i) == 1
% Measure the time until the next falling edge
j = i + 1;
while pwm(j) == 1 && j < length(pwm)
j = j + 1;
end
% Reconstruct the analog value based on pulse width
demodulated_signal(i) = mean(msg(i:j-1));
end
end

% Low-Pass Filter Design
Fs = 1 / (t(2) - t(1)); % Sampling frequency
Fc = 5; % Cutoff frequency (adjust based on your signal)
[b, a] = butter(4, Fc / (Fs / 2), 'low'); % 4th-order Butterworth filter

% Apply the Low-Pass Filter
filtered_signal = filtfilt(b, a, demodulated_signal);

% Plot the demodulated and filtered signal for comparison
subplot(4,1,4);
plot(t, filtered_signal, 'r', 'LineWidth', 1.5); % Filtered signal in red
title('Demodulated Signal (Filtered)');
xlabel('Time');
ylabel('Amplitude');
grid on;
 

Output 



 MATLAB Code for Digital Pulse Width Modulation (PWM)


% This code is developed by SalimWireless.Com
clc; clear; close all;
% Digital SPWM Generator using Square Wave Carrier in MATLAB

% === PARAMETERS ===
fs_carrier = 20;       % Carrier frequency in Hz
f_signal = 5;           % Message signal frequency in Hz
sampleRate = 50000;      % Samples per second
duration = 1;            % Duration in seconds

% === TIME VECTOR ===
t = linspace(0, duration, sampleRate * duration);

% === MESSAGE SIGNAL (SINE WAVE) ===
signal = sin(2 * pi * f_signal * t);

% === NORMALIZE SIGNAL TO 0–1 FOR DUTY CYCLE ===
normalizedSignal = (signal + 1) / 2;  % Scale from [-1, 1] to [0, 1]

% === PWM GENERATION BASED ON SQUARE CARRIER PERIODS ===
samplesPerCarrierPeriod = floor(sampleRate / fs_carrier);
pwm = zeros(1, length(t));

% Generate PWM: For each carrier cycle, set ON time based on message amplitude at start
for i = 1:samplesPerCarrierPeriod:length(t)
    startIndex = i;
    if startIndex > length(t)
        break;
    end
    
    % Duty cycle at start of period
    duty = normalizedSignal(startIndex);
    onSamples = floor(samplesPerCarrierPeriod * duty);
    
    % Set PWM high for onSamples
    endIndex = min(startIndex + samplesPerCarrierPeriod - 1, length(t));
    onEndIndex = min(startIndex + onSamples - 1, endIndex);
    
    pwm(startIndex:onEndIndex) = 1;
end

% === GENERATE SQUARE CARRIER FOR REFERENCE PLOTTING ===
carrierSquare = double(mod(t * fs_carrier, 1) < 0.5);

% === TRIM TO FIRST 3 CYCLES OF MESSAGE SIGNAL FOR VISUALIZATION ===
samplesToPlot = floor(3 * (sampleRate / f_signal));
t_plot = t(1:samplesToPlot);
signal_plot = signal(1:samplesToPlot);
carrier_plot = carrierSquare(1:samplesToPlot);
pwm_plot = pwm(1:samplesToPlot);

% === PLOTTING ===
figure('Name', 'PWM with Square Wave Carrier', 'Color', 'w');
hold on;
plot(t_plot, signal_plot, 'b', 'LineWidth', 1.2);
plot(t_plot, carrier_plot, 'g--', 'LineWidth', 1);
stairs(t_plot, pwm_plot, 'r', 'LineWidth', 1.2);
hold off;

xlabel('Time (s)');
ylabel('Amplitude');
title('PWM Output with Square Wave Carrier');
legend('Message Signal (Sine)', 'Square Carrier', 'PWM Output', 'Location', 'southoutside', 'Orientation', 'horizontal');
grid on;
web('https://www.salimwireless.com/search?q=pwm%20pulse%20modulation', '-browser');

Output

 
Parameter PAM PWM PPM DM PCM
Parameter varied Signal amplitude Pulse duration Pulse timing Sample difference (delta) Digital code
Pulse duration Fixed Adjustable Fixed Fixed Fixed
Resistance to noise Poor Average Good Average Good
Bandwidth requirement Low Moderate High Low High
Implementation complexity Low Medium High Low High
MATLAB implementation PAM Script PWM Script PPM Script DM Script PCM Script
Further reading PAM PWM PPM DM PCM

PWM Signal Generation

 

 
 
 
Return to Home →

Further Reading

📊 MATLAB & ONLINE SIMULATIONS (Web-Based)
🎓 GATE & UGC - NET SOLUTIONS
📡 WIRELESS & SIGNAL PROCESSING


Contact Us

Name

Email *

Message *

Popular Posts

Hybrid Beamforming | Page 1

Beamforming Techniques Hybrid Beamforming... Page 1 | Page 2 | Hybrid Beamforming: Hybrid beam formation was developed to address some of the limitations of digital pre-coding approaches. Every antenna element is connected to an RF chain in digital pre-coding (beam forming) method. We also know that each RF chain is in charge of providing a separate data stream between the transmitter and the receiver. We know that a larger number of independent data streams leads to higher data rates. It has a spatial multiplexing feature for MIMO. As a result, we may assume that switching from MIMO to massive MIMO will benefit us more in terms of spatial multiplexing in massive MIMO, where each antenna is coupled to a single RF chain. We'll proceed with a definition of hybrid beam forming. Overview of hybrid beam forming with example: Unlike digital beam forming, more than one antenna element is connected to a single RF chain in hybr...
Read more

MATLAB Code for 8-PSK, 16-PSK, ...

📘 Overview & Theory 🧮 MATLAB Code for BPSK, QPSK, 8-PSK, 16-PSK, 32-PSK 🧮 Simulator for m-ary PSK 📚 Further Reading   MATLAB Code for BPSK, QPSK, 8-PSK, 16-PSK, 32-PSK clc; clear all; close all; rng(10) M = 8; % M = 2, 4, 8, 16, 32, etc. N_Bits = 2520; Phase = 0; data_info_bit = randi([0,1],N_Bits,1); data_temp = bi2de(reshape(data_info_bit,N_Bits/log2(M),log2(M))); modData = pskmod(data_temp,M,Phase); figure(1); scatterplot(modData); channelAWGN = 15; rxData2 = awgn(modData, channelAWGN); figure(2); scatterplot(rxData2); demodData = pskdemod(rxData2,M,Phase);   for BPSK, Constellation Size, M = 2 for QPSK, M = 4 for 8-PSK, M = 8, and so on    Output Figure: 8-PSK Modulation Figure: 8-PSK Demodulation after adding AWGN Noise Using the above MATLAB code you'll able be to modulate and demodulate 2-PSK, 4-PSK, 8-PSK, 16-PSK, 32-PSK and so on.  16-PSK   Fig: 16-PSK In this above code ' M ' is the number of the conste...
Read more

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

How Windowing Affects Your Periodogram

The windowed periodogram is a widely used technique for estimating the Power Spectral Density (PSD) of a signal. It enhances the classical periodogram by mitigating spectral leakage through the application of a windowing function. This technique is essential in signal processing for accurate frequency-domain analysis.   Power Spectral Density (PSD) The PSD characterizes how the power of a signal is distributed across different frequency components. For a discrete-time signal, the PSD is defined as the Fourier Transform of the signal’s autocorrelation function: S x (f) = FT{R x (Ï„)} Here, R x (Ï„)}is the autocorrelation function. FT : Fourier Transform   Classical Periodogram The periodogram is a non-parametric PSD estimation method based on the Discrete Fourier Transform (DFT): P x (f) = \(\frac{1}{N}\) X(f) 2 Here: X(f): DFT of the signal x(n) N: Signal length However, the classical periodogram suffers from spectral leakage due to abrupt truncation of the ...
Read more

Advanced M-ary Modulation Simulator: Constellation, min dist, Efficiency, SER, EVM (RMS)

Advanced M-ary Communication Lab Analytical & Statistical Performance of Digital Modulation Theoretical Probability of Error (\(P_s\)) \[ P_s = Q\left(\sqrt{\frac{2 E_b}{N_0}}\right) \] Modulation (M-ary) BPSK (M=2) QPSK (M=4) 8-PSK (M=8) 16-QAM (M=16) 64-QAM (M=64) 256-QAM (M=256) SNR (\(E_b/N_0\)): 12 dB Efficiency 2 bps/Hz Min Dist (\(d_{min}\)) 1.41 Symbol Error 1.2e-5 EVM (RMS) 0.0% Constellation Diagram Noise PDF & Decision Tail 1. Geometric Mapping ...
Read more

Frequency Shift Keying (FSK) Modulation & Demodulation (with Simulation)

Frequency Shift Keying (FSK) Theoretical Foundations: Frequency Shift Keying (FSK) is a discrete frequency modulation scheme wherein the digital information is encoded via instantaneous shifts in the carrier signal's frequency. The fundamental implementation is Binary FSK (BFSK), which maps binary data onto two distinct, discrete spectral states. A binary '1' (the "mark" state) is represented by a carrier frequency \( f_1 \), while a binary '0' (the "space" state) corresponds to frequency \( f_2 \). Each symbol is sustained for a bit interval denoted by \( T_b \). FSK Transmitter Characterization: The mathematical model for the modulated BFSK output \( s(t) \) is defined as: \[ s(t) = \begin{cases} A_c \cos(2\pi f_1 t), & \text{for } m = 1 \\ A_c \cos(2\pi f_2 t), & \text{for } m = 0 \end{cases} \] ...
Read more

Galois Fields: GF(2) and GF(2m) and Primitive Polynomial

Galois Fields: GF(2) and GF(2 m ) 1. What is a Galois Field (GF)? A Galois Field (GF) is a finite set of elements in which the four basic arithmetic operations—addition, subtraction, multiplication, and division (except by zero)— are all well defined and closed. GF(q) ⇒ a field with exactly q elements 2. The Simplest Field: GF(2) GF(2) is the smallest possible finite field and forms the foundation of all digital systems. GF(2) = {0, 1} Addition in GF(2) Addition is performed modulo 2 (XOR operation): + 0 1 0 0 1 1 1 0 Multiplication in GF(2) × 0 1 0 0 0 1 0 1 GF(2) is used in binary logic, XOR operations, and simple error-control codes. 3. Meaning of GF(2 m ) GF(2 m ) is a finite field containing exactly 2 m elements . Each element represents an m-bit symbol . Field Number of Elements GF(2) 2 GF(2²) 4 GF(2³) 8 GF(2⁸) 256 Important: GF(2 m ) is not integer arithmetic modulo 2 m . It is polynomial-based arithmetic. 4....
Read more