clear all
echo off
SNRindB1=0:2:12;
SNRindB2=0:0.1:12;
for i=1:length(SNRindB1),
     smld_err_prb(i)=cm_sm34(SNRindB1(i));  %仿真誤碼率
end;
for i=1:length(SNRindB2),
     SNR=exp(SNRindB2(i)*log(10)/10);       %信噪比
     theo_err_prb(i)=exp((-2)*SNR*(sin(pi/16).^2));   %理論符號誤碼率
end;
%隨后為繪圖曲線
semilogy(SNRindB1,smld_err_prb,'*');
hold
semilogy(SNRindB2,theo_err_prb);
grid