16qam_ray_q.m

给出比较少见的16PSK的仿真算法

           
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));   
   smld_err_q_prb(i)=cm_sm53(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 on
semilogy (SNRindB1,smld_err_prb,'*');
semilogy (SNRindB1,smld_err_q_prb,'o');
grid on
xlabel('Eb/No in dB')
ylabel('error probability')
title('16QAM系统的瑞利信道仿真')
legend('理论误码率','均衡前仿真误码率','均衡后仿真误码率');