www.gusucode.com > matlab通信工程仿真源码（张德丰等编著）程序书籍 > matlab_code/matlab通信工程仿真源码（张德丰等编著）/第8章/li8_14.m

     clear all;
SNRindB1=0:2:15; SNRindB2=0:1:15;
M=16;k=log2(M);
for i=1:length(SNRindB1)
    s_err_prb(i)=Qxiufun(SNRindB1(i));
end
for i=1:length(SNRindB2)
    SNR=exp(SNRindB2(i)*log(10)/10);
    t_err_prb(i)=Qxiufun(sqrt(3*k*SNR/(M-1)));
end
semilogy(SNRindB1,s_err_prb,'rp'); %用对数坐标作出实际信噪比-误比特率的点
hold on;
semilogy(SNRindB2,t_err_prb); %用对数坐标作出理论信噪比-误比特率曲线