comp_exam8_1.m

来自「the source code of Ziemer s book Digital」· M 代码 · 共 29 行

% Bit error probabilities -  M-ary PSK
%
clf
z_dB = 5:.1:40;
z = 10.^(z_dB/10);
for j = 1:8
   M=2^j;
   A=2/log2(M);
   k=(sin(pi/M))^2*log2(M);
   if j == 1 | j == 2
      Pb = qfn(sqrt(2*z));
   elseif j >= 3
      Pb = A*qfn(sqrt(k*2*z));
   end
   semilogy(z_dB,Pb),axis([5 40 10^(-6) 1]),xlabel('E_b/N_0, dB'),ylabel('P_E_, _b_i_t'),...
   if j == 1 | j == 2
      text(z_dB(50)+.2, Pb(50), 'M = 2 & 4')
   elseif j == 3 | j == 4
      text(z_dB(50)+.2, Pb(50), ['M = ', num2str(M)])
   elseif j > 4
      text(z_dB(180)+.2, Pb(180), ['M = ', num2str(M)])
   end
   if j==1
      hold on
      grid
   end
end
title('Bit error probability for MPSK')