qammotocalo.m

16QAM通信系统的MonteCalo仿真程序

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);  %转化信噪比为数值表示
    theo_err_prb(i)=4*Qfunct(sqrt(3*k*SNR/(M-1)));
end;
semilogy(SNRindB1,smld_bit_err_prb,'*');%用对数坐标作出实际信噪比——误比特率点
hold                                    %将上一曲线保持住
semilogy(SNRindB2,theo_err_prb);        %作出对数坐标理论的信噪比——误比特率曲线