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

    clear all;
rho_b1=0:5:35;      %rho in dB代表仿真的误码率
rho_b2=0:0.1:35;    %rho in dB代表理论计算得出的误码率
for i=1:length(rho_b1)
    s_err_prb(i)=li8_16fun(rho_b1(i));     %仿真误码率
end;
for i=1:length(rho_b2)
    temp=10^(rho_b2(i)/10);
    if(temp>2)
        t_err_rate(i)=1/(exp(1)*temp);    %如果rho>2的理论误码率
    else
        t_err_rate(i)=(1/2)*exp(-temp/2); %如果rho<2的理论误码率
    end
end
semilogy(rho_b1,s_err_prb,'rp',rho_b2,t_err_rate,'-');