% MATLAB script for Illustrated Problem 5, Chapter 5.
clear all
echo on
SNRindB1=0:1:10;
SNRindB2=0:10;
for i=1:length(SNRindB1),
  % simulated error rate
  smld_err_prb(i)=smldPe55(SNRindB1(i));   
end;
for i=1:length(SNRindB2),
  SNR=exp(SNRindB2(i)*log(10)/10);   
  % theoretical error rate      
  theo_err_prb(i)=Qfunct(sqrt(2*SNR));     
end;
% Plotting commands follow
%hold on
semilogy(SNRindB1,smld_err_prb,'*');
hold on
semilogy(SNRindB2,theo_err_prb);
xlabel('Eb/N0 (dB)');
ylabel('Probability of Bit Error');
legend('Simulation','Analytical');
grid