16qam_ray.m

来自「给出比较少见的16PSK的仿真算法」· M 代码 · 共 24 行

           
clear all
echo off
SNRindB1=0:2:10;
SNRindB2=0:0.1:10;
M=16;
k=log2(M);
for i=1:length (SNRindB1),
   smld_err_prb(i)=cm_sm52(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 (SNRindB2,theo_err_prb);
hold
semilogy (SNRindB1,smld_err_prb,'*');
grid on
xlabel('Eb/No in dB')
ylabel('error probability')
title('16QAM系统的瑞利信道仿真')
legend('理论误符号率','仿真误符号率');