clear
close
clc
rng(10)
%Parameters
A = 1;
rb = 10000;
Tb = 1/rb;
% Amplitude
% Bit rate
% Bit interval
fc = 20*rb;
% Carrier frequency
fs = 20*fc;
Ts = 1/fs;
% Sampling frequency
% Sampling interval
Ns = floor(Tb/Ts);
% Number of samples per bit
N = 10000;
% Number of bits to be transmitted
msg = randi([0,1],1,N);
% Message signal
[t,sig] = polar_nrz_6(msg,A,rb,fs);
% Polar non return to zero
figure(1)
plot(t,sig,’r’,’LineWidth’,2)
xlabel(‘Time’)
ylabel(‘Amplitude’)
title(‘PolarNRZ Signal’)
set(gca,’fontsize’,18)
carrier = A*cos(2*pi*fc*t);
psk = sig.*carrier;
% Carrier wave
% BPSK modulated signal
figure(2)
plot(t,psk,’r’,’LineWidth’,1.5)
xlabel(‘Time’)
ylabel(‘Amplitude’)
title(‘BPSK modulated Signal’)
set(gca,’fontsize’,18)
figure(3)
plot(real(sig),imag(sig),’ro’,’MarkerFaceColor’,’r’,’MarkerSize’,10,’M
arkerEdgeColor’,’k’)
xlim([-1.5,1.5])
ylim([-1.5,1.5])
title(‘Constellation Diagram’)
set(gca,’fontsize’,18)
pulse = ones(1,Ns); % Reference pulse
EbN0_db = 0;
% Energy per bit to noise power spectral density ratio
EbN0 = 10.^(EbN0_db/10); % Converting to linear scale
Eb = Ns*sum(abs(psk).^2)/length(psk); % Signal energy per bit
N0 = Eb/EbN0;
% Noise power spectral density
n = sqrt(N0/2)*randn(1,length(psk)); % Generating Noise
received = psk+n; % Received signal
% Correlation decoding
r_bb = received.*cos(2*pi*fc*t);
dem = filter(pulse,1,r_bb); % Demodulated signal
dem = dem(Ns:Ns:end);
% Sampling at Tb
det = dem;
% Detection
det(dem>=0)=1;
det(dem
Purchase answer to see full
attachment