%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
            %   對(duì)M=8的PSK系統(tǒng)進(jìn)行蒙特卡羅仿真   %
            %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear all
echo off
SNRindB3=0:2:10;
SNRindB4=0:0.1:10;
for i=1:length (SNRindB3),
   [pb , ps]=cm_sm321 (SNRindB3 (i));       %仿真比特和符號(hào)誤碼率
   smld_bit_err_prb (i)=pb;
   smld_symbol_err_prb (i)=ps;
end;
for i=1:length (SNRindB4),
   SNR=exp(SNRindB4(i)*log(10)/10);        %信噪比
   theo_err_prb(i)=2*Qfunct(sqrt(6*SNR)*sin(pi/8))/3;    %理論比特誤碼率
end;
 %%%%%%%隨后繪圖曲線%%%%%%%
semilogy (SNRindB3,smld_bit_err_prb,'+');
hold
semilogy (SNRindB3,smld_symbol_err_prb,'^');
semilogy (SNRindB4,theo_err_prb,':');
grid on
xlabel('Eb/No in dB')
ylabel('error probability')
title('8PSK系統(tǒng)的蒙特卡羅仿真')
%gtext('（注：“--”理論誤碼率；“*”誤比特率；“o”誤符號(hào)率）')