bpskprob.m

来自「压缩文件中是Error Correction Coding - Mathemat」· M 代码 · 共 29 行

% Plot probability of error for BPSK

% Copyright 2004 by Todd K. Moon
% Permission is granted to use this program/data
% for educational/research only

% P0 = 0.1;  % =P(+Eb)
% P1 = 1-P0;  % P(-Eb)
startsnrdb = 0;
snrstepdb = 0.1;
endsnrdb = 10;
N0 = 2;  % sigma^2 = N0/2 = 1
sigma2 = N0/2;
sigma = sqrt(sigma2);
plist = [];
snrlist = [];
for snrdb=startsnrdb:snrstepdb:endsnrdb
  snrlist = [snrlist snrdb];
  snr = 10^(snrdb/10);   % snr not in dB
  tau = sigma2/2*log(P1/P0);
  
  Eb = N0*snr*R;
  Ebsqr = sqrt(Eb);
  
  p = qf((tau + Ebsqr)/sigma)*P1 + qf((Ebsqr - tau)/sigma)*P0;
  plist = [plist p];
end
% semilogy(snrlist,plist,'r');