Scribd
Upload
33 views1 page

% Calculation of Error Probabilit Using Monte Carlo Simulation

This document contains MATLAB code that simulates communication over a channel to calculate bit error rate (BER) using Monte Carlo simulation. It generates random transmitted symbols with noise, calculates the received signal, detects the transmitted bit, and compares with actual bit to determine error. It runs this simulation over a range of signal-to-noise ratio (SNR) values to generate a BER curve. It also calculates the theoretical BER using a formula to plot both simulation results and theoretical curve on the same graph.

Uploaded by

Nida Waheed
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
33 views1 page

% Calculation of Error Probabilit Using Monte Carlo Simulation

This document contains MATLAB code that simulates communication over a channel to calculate bit error rate (BER) using Monte Carlo simulation. It generates random transmitted symbols with noise, calculates the received signal, detects the transmitted bit, and compares with actual bit to determine error. It runs this simulation over a range of signal-to-noise ratio (SNR) values to generate a BER curve. It also calculates the theoretical BER using a formula to plot both simulation results and theoretical curve on the same graph.

Uploaded by

Nida Waheed
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 1

clc

clear
D=2;
sigma = 1;
Eb = 1/sqrt(2);
EbNo_r_per_ch_dB = 5:5:30;
EbNo_r_per_ch = 10.^(EbNo_r_per_ch_dB/10);
No = Eb*2*sigma^2*10.^(-EbNo_r_per_ch_dB/10);
BER = zeros(1,length(No));
SNR_rx_per_b_per_ch = zeros(1,length(No));
% Calculation of error probabilit using Monte Carlo simulation:
for i = 1:length(No)
no_bits = 0;
no_errors = 0;
P_r_t = 0;
P_n_t = 0;
r = zeros(2,2);
R = zeros(1,2);
% Total rxd power
% Total noise power
% Assumption: m = 1 (All one codeword is transmitted):
while no_errors <= 2
no_bits = no_bits + 1;
u = rand(1,2);
alpha = sigma*sqrt(-2*log(u));
phi = 2*pi*rand(1,2);
noise = sqrt(No(i)/2)*(randn(2,2) + 1i*randn(2,2));
r(1,1) = alpha(1)*sqrt(Eb)*exp(1i*phi(1))+noise(1,1);
r(1,2) = noise(1,2);
r(2,1) = alpha(2)*sqrt(Eb)*exp(1i*phi(2))+noise(2,1);
r(2,2) = noise(2,2);
R(1) = abs(r(1,1))^2 + abs(r(2,1))^2;
R(2) = abs(r(1,2))^2 + abs(r(2,2))^2;
if R(1) <= R(2)
m_h =0;
else
m_h = 1;
end
P_n_t = P_n_t + No(i);
P_r_t = P_r_t + 0.5*(abs(r(1))^2 + abs(r(2))^2);
no_errors = no_errors + (1-m_h);
end
SNR_rx_per_b_per_ch(i) = (P_r_t-P_n_t)/P_n_t;
BER(i) = no_errors/no_bits;
end
% Calculation of error probabilit using the theoretical formula:
rho = EbNo_r_per_ch;
rho_dB = 10*log10(rho);
rho_b = D*rho;
rho_b_dB = 10*log10(rho_b);
K_D = factorial((2*D-1))/factorial(D)/factorial((D-1));
P_2 = K_D./rho.^D;
% Plot the results:
semilogy(rho_b_dB,BER,'-*',rho_b_dB,P_2,'-o')
xlabel('Average SNR/bit (dB) '); ylabel('BER')
legend('Monte Carlo simulation','Theoretical value')

You might also like