16psk_ray.m

给出比较少见的16PSK的仿真算法

clear all
echo off
SNRindB1=0:2:20;
SNRindB2=0:0.1:20;

for i=1:length(SNRindB1),
   [pb,ps]=cm_sm42(SNRindB1(i));
   smld_bit_err_prb(i)=pb;
   smld_symbol_err_prb(i)=ps;
end;

for i=1:length(SNRindB2),
   SNR=exp(SNRindB2(i)*log(10)/10);
   theo_err_prb(i)=0.5*(1-1/sqrt(1+1/SNR));
 end;
semilogy(SNRindB2,theo_err_prb);
hold on
semilogy(SNRindB1,smld_bit_err_prb,'*');
semilogy(SNRindB1,smld_symbol_err_prb,'O');
axis([0,10,10^(-5),10^1]);
grid
xlabel('Eb/No in dB')
ylabel('error probability')
title('均衡前单径瑞利信道下16PSK通信系统仿真')
legend('理论误码率','仿真误比特率','仿真误符号率');