bdt_f533.m

现代通信系统书第五章程序（电子工业出版社）

% MATLAB script that generates the probability of error versus the signal to noise ratio.
initial_snr=0;
final_snr=12;
snr_step=0.75;			
tolerance=eps;			   	% Tolerance used for the integration
plus_inf=20;			   	% This is practically infinity
snr_in_dB=initial_snr:snr_step:final_snr;
for i=1:length(snr_in_dB),
   snr=10^(snr_in_dB(i)/10);
   Pe_2(i)=1-quad8('bdt_int2',0,plus_inf,tolerance,[],snr,2);
   Pe_4(i)=1-quad8('bdt_int2',0,plus_inf,tolerance,[],snr,4);
   Pe_8(i)=1-quad8('bdt_int2',0,plus_inf,tolerance,[],snr,8);
   Pe_16(i)=1-quad8('bdt_int2',0,plus_inf,tolerance,[],snr,16);
   Pe_32(i)=1-quad8('bdt_int2',0,plus_inf,tolerance,[],snr,32);
   echo off;
end;
% Plotting commands follow