echo on
rho_b1=0:5:35;                                  %用于計(jì)算仿真誤碼率 
rho_b2=0:0.1:35;                                %用于計(jì)算理論誤碼率  
for i=1:length(rho_b1),
    smld_err_prb(i)=err1(rho_b1(i));
    echo off;
end;

%plot(smld_err_prb,'*');

echo on;
for i=1:length(rho_b2),
    temp=10^(rho_b2(i)/10);                     %轉(zhuǎn)換為db 
    if (temp>2)
        theo_err_rate(i)=1/(exp(1)*temp);       %p413  
    else
        theo_err_rate(i)=(1/2)*exp(-temp/2);
    end;
    echo off;
end;
echo on;



semilogy(rho_b1,smld_err_prb,'*',rho_b2,theo_err_rate,'r');