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

 

 
 
 

Further Reading

People are good at skipping over material they already know!

View Related Topics to







Contact Us

Name

Email *

Message *

Popular Posts

Constellation Diagrams of ASK, PSK, and FSK with MATLAB Code + Simulator

📘 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

Online Simulator for ASK, FSK, and PSK

Try our new Digital Signal Processing Simulator!   Start Simulator for binary ASK Modulation Message Bits (e.g. 1,0,1,0) Carrier Frequency (Hz) Sampling Frequency (Hz) Run Simulation Simulator for binary FSK Modulation Input Bits (e.g. 1,0,1,0) Freq for '1' (Hz) Freq for '0' (Hz) Sampling Rate (Hz) Visualize FSK Signal Simulator for BPSK Modulation ...
Read more

Fading : Slow & Fast and Large & Small Scale Fading (with MATLAB Code + Simulator)

📘 Overview 📘 LARGE SCALE FADING 📘 SMALL SCALE FADING 📘 SLOW FADING 📘 FAST FADING 🧮 MATLAB Codes 📚 Further Reading LARGE SCALE FADING The term 'Large scale fading' is used to describe variations in received signal power over a long distance, usually just considering shadowing.  Assume that a transmitter (say, a cell tower) and a receiver  (say, your smartphone) are in constant communication. Take into account the fact that you are in a moving vehicle. An obstacle, such as a tall building, comes between your cell tower and your vehicle's line of sight (LOS) path. Then you'll notice a decline in the power of your received signal on the spectrogram. Large-scale fading is the term for this type of phenomenon. SMALL SCALE FADING  Small scale fading is a term that describes rapid fluctuations in the received signal power on a small time scale. This includes multipath propagation effects as well as movement-induced Doppler fr...
Read more

Theoretical BER vs SNR for binary ASK, FSK, and PSK with MATLAB Code + Simulator

📘 Overview & Theory 🧮 MATLAB Codes 📚 Further Reading Theoretical BER vs SNR for Amplitude Shift Keying (ASK) The theoretical Bit Error Rate (BER) for binary ASK depends on how binary bits are mapped to signal amplitudes. For typical cases: If bits are mapped to 1 and -1, the BER is: BER = Q(√(2 × SNR)) If bits are mapped to 0 and 1, the BER becomes: BER = Q(√(SNR / 2)) Where: Q(x) is the Q-function: Q(x) = 0.5 × erfc(x / √2) SNR : Signal-to-Noise Ratio N₀ : Noise Power Spectral Density Understanding the Q-Function and BER for ASK Bit '0' transmits noise only Bit '1' transmits signal (1 + noise) Receiver decision threshold is 0.5 BER is given by: P b = Q(0.5 / σ) , where σ = √(N₀ / 2) Using SNR = (0.5)² / N₀, we get: BER = Q(√(SNR / 2)) Theoretical BER vs ...
Read more

What is - 3dB Frequency Response? Applications ...

📘 Overview & Theory 📘 Application of -3dB Frequency Response 🧮 MATLAB Codes 🧮 Online Digital Filter Simulator 📚 Further Reading Filters What is -3dB Frequency Response?   Remember, for most passband filters, the magnitude response typically remains close to the peak value within the passband, varying by no more than 3 dB. This is a standard characteristic in filter design. The term '-3dB frequency response' indicates that power has decreased to 50% of its maximum or that signal voltage has reduced to 0.707 of its peak value. Specifically, The -3dB comes from either 10 Log (0.5) {in the case of power} or 20 Log (0.707) {in the case of amplitude} . Viewing the signal in the frequency domain is helpful. In electronic amplifiers, the -3 dB limit is commonly used to define the passband. It shows whether the signal remains approximately flat across the passband. For example, in pulse shapi...
Read more

BER vs SNR for M-ary QAM, M-ary PSK, QPSK, BPSK, ...(MATLAB Code + Simulator)

📘 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

Pulse Shaping using Raised Cosine Filter in MATLAB

  MATLAB Code for Raised Cosine Filter Pulse Shaping clc; clear; close all ; %% ===================================================== %% PARAMETERS %% ===================================================== N = 64; % Number of OFDM subcarriers cpLen = 16; % Cyclic prefix length modOrder = 4; % QPSK oversample = 8; % Oversampling factor span = 10; % RRC filter span in symbols rolloff = 0.25; % RRC roll-off factor %% ===================================================== %% Generate Baseband OFDM Symbols %% ===================================================== data = randi([0 modOrder-1], N, 1); % Random bits txSymbols = pskmod(data, modOrder, pi/4); % QPSK modulation % IFFT to get OFDM symbol tx_ofdm = ifft(txSymbols, N); % Add cyclic prefix tx_cp = [tx_ofdm(end-cpLen+1:end); tx_ofdm]; %% ===================================================== %% Oversample the Baseband Signal %% ===============================================...
Read more

Theoretical BER vs SNR for BPSK

Theoretical Bit Error Rate (BER) vs Signal-to-Noise Ratio (SNR) for BPSK in AWGN Channel 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 N₀/2 (where N₀ 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 Rat...
Read more