%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
            %   對(duì)M=64的QAM系統(tǒng)進(jìn)行蒙特卡羅仿真   %
            %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear all
echo off
SNRindB1=0:2:20;
SNRindB2=0:0.1:20;
M=64;
k=log2(M);
for i=1:length (SNRindB1),
   smld_err_prb(i)=cm_sm41(SNRindB1(i));       %仿真誤符號(hào)率
end;
for i=1:length (SNRindB2),
   SNR=exp(SNRindB2(i)*log(10)/10);        %信噪比
   theo_err_prb(i)=4*Qfunct(sqrt(3*k*SNR/(M-1)));    %理論誤符號(hào)率
end;
 %%%%%%%隨后繪圖曲線(xiàn)%%%%%%%
semilogy (SNRindB1,smld_err_prb,'*');
hold
semilogy (SNRindB2,theo_err_prb);
grid on
xlabel('Eb/No in dB')
ylabel('error probability')
title('64QAM系統(tǒng)的蒙特卡羅仿真')
%gtext('（注：“—”理論誤符號(hào)率；“*”仿真誤符號(hào)率）')