cm_sm32.m

16psk通信系统的信噪比与误码率计算并绘图

echo on 
SNRindB1=0:2:10;
SNRindB2=0:0.1:10 ;
for i=1:length(SNRindB1);
   [pb,ps]=cm_sm(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)=2*Qfunct(sqrt(8*SNR)*sin(pi/16))/4; 
 end; 
semilogy(SNRindB1,smld_bit_err_prb,'*'); 
hold on 
semilogy(SNRindB1,smld_symbol_err_prb,'O'); 
semilogy(SNRindB2,theo_err_prb); 
axis([0,10,10^(-5),10^0]); 
grid 
xlabel('Eb/No in dB') 
ylabel('error probability') 
title('16PSK通信系统的蒙特卡洛仿真') 
gtext('（注：“—”理论误码率；“*”误比特率；“o”误符号率 ）')