echo off
SNRindB1=0:2:15;
SNRindB2=0:0.1:15;
M=16;
k=log2(M);
for i=1:length(SNRindB1),
    %仿真誤碼率
    smld_err_prb(i)=cm_sm41(SNRindB1(i));
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)));
end;
%繪圖命令
grid on;
semilogy(SNRindB1,smld_err_prb,'*');
hold
semilogy(SNRindB2,theo_err_prb);