16psk_ray_q.m

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

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

for i=1:length(SNRindB1),
   [pb,ps]=cm_sm42(SNRindB1(i));
   [pb1,ps]=cm_sm43(SNRindB1(i));
   smld_bit_err_prb(i)=pb;
   smld_bit_err_q_prb(i)=pb1;
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_bit_err_q_prb,'O');
axis([0,10,10^(-5),10^1]);
grid
xlabel('Eb/No in dB')
ylabel('error probability')
title('单径瑞利信道下16PSK通信系统仿真')
legend('理论误码率','均衡前仿真误码率','均衡后仿真误码率');
echo off;