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 = 1:length(qam_orders)
    qam_order = qam_orders(i);
    ber_qam_results(i, :) = berawgn(snr_db, 'qam', qam_order);
end
figure;
semilogy(snr_db, ber_psk_results(1, :), 'o-', 'DisplayName', 'BPSK'); hold on;
for i = 2:length(psk_orders)
    semilogy(snr_db, ber_psk_results(i, :), 'o-', 'DisplayName', sprintf('%d-PSK', psk_orders(i)));
end
for i = 1:length(qam_orders)
    semilogy(snr_db, ber_qam_results(i, :), 'x-', 'DisplayName', sprintf('%d-QAM', qam_orders(i)));
end
title('BER vs SNR for Various PSK and QAM Schemes');
xlabel('SNR (dB)');
ylabel('Bit Error Rate (BER)');
legend('Location', 'best'); grid on; ylim([1e-5 1]); hold off;
web('https://www.salimwireless.com/search?q=m-ary%20PSK%20and%20QAM', '-browser');

Output

BER vs SNR for PSK & QAM

Fig: BER vs SNR graph for Various PSK and QAM

MATLAB Code for BER vs SNR for m-ary QAM

MATLAB Code for BER vs SNR for m-ary PSK

clc; clear all; close all; rng(0);
num_symbols = 1e5;
snr_db = 0:2:20;
psk_orders = [2, 4, 8, 16, 32];
ber_results = zeros(length(psk_orders), length(snr_db));
for i = 1:length(psk_orders)
    M = psk_orders(i);
    k = log2(M);
    for j = 1:length(snr_db)
        data_symbols = randi([0, M-1], 1, num_symbols);
        tx_bits = de2bi(data_symbols, k, 'left-msb');
        modulated_signal = pskmod(data_symbols, M, pi/M, 'gray');
        received_signal = awgn(modulated_signal, snr_db(j), 'measured');
        demod_symbols = pskdemod(received_signal, M, pi/M, 'gray');
        rx_bits = de2bi(demod_symbols, k, 'left-msb');
        ber_results(i, j) = sum(tx_bits(:) ~= rx_bits(:)) / (num_symbols * k);
    end
end
figure; semilogy(snr_db, ber_results(1, :), 'o-', 'DisplayName', 'BPSK'); hold on;
for i = 2:length(psk_orders)
    semilogy(snr_db, ber_results(i, :), 'o-', 'DisplayName', sprintf('%d-PSK', psk_orders(i)));
end
title('Bit Error Rate vs. SNR for Various PSK Schemes');
xlabel('SNR (dB)'); ylabel('Bit Error Rate (BER)');
legend('Location', 'best'); grid on; hold off;
web('https://www.salimwireless.com/search?q=m-ary%20PSK%20and%20QAM', '-browser');

Search This Blog