example1.asv

contemporarycommunicatonsystem－using MATLAB经典书目源程序！

SNRindB1=0:1:20;	  		
SNRindB2=0:0.1:20;		
d=1:0.1:2.0;

echo on;
for i=1:length(d),
s1(i)=(d(i)-1)^2/(sqrt(d(i)^2+1)*(2-d(i)));
s2(i)=sqrt(8)*(d(i)-1)^2/(sqrt(d(i)^2+1)*(4*d(i)-2));

for j=1:length(SNRindB2),
  % signal to noise ratio
  SNR_per_bit=exp(SNRindB2(j)*log(10)/10);  
  % theoretical error rate
  theo_err_prb(j)=0.5*Qfunct(s1*sqrt(SNR_per_bit))+0.5*Qfunct(s2*sqrt(SNR_per_bit));  
  echo off;
end;
end;

echo on;
% Plotting commands follow
semilogy(SNRindB2,theo_err_prb(1),'r');
hold;