echo on
SNRindB1=0:2:15;   %定義信噪比序列，共有八個值
SNRindB2=0:0.1:15; %掃描用的信噪比序列，采用dB作單位
M=16;  k=log2(M);
for i=1:length(SNRindB1)              %計算誤比特率
   
    smld_bit_err_prb(i)=qammoto(SNRindB1(i));
   
end;
for i=1:length(SNRindB2)
    SNR=exp(SNRindB2(i)*log(10)/10);  %轉(zhuǎn)化信噪比為數(shù)值表示
    theo_err_prb(i)=4*Qfunct(sqrt(3*k*SNR/(M-1)));
end;
semilogy(SNRindB1,smld_bit_err_prb,'*');%用對數(shù)坐標(biāo)作出實(shí)際信噪比——誤比特率點(diǎn)
hold                                    %將上一曲線保持住
semilogy(SNRindB2,theo_err_prb);        %作出對數(shù)坐標(biāo)理論的信噪比——誤比特率曲線