www.pudn.com > fhss02.rar > main.m

function main; 
echo on 
rho_b1=0:5:35; 
rho_b2=0:0.1:35; 
for i=1:length(rho_b1) 
    smld_err_prb(i)=ss_Pe96(rho_b1(i)); 
end 
for i=1:length(rho_b2) 
    temp=10^(rho_b2(i)/10); 
    if temp>2 
        theo_err_rate(i)=1/(exp(1)*temp); 
    else 
        theo_err_rate(i)=1/2*exp(-temp/2); 
    end 
end 
semilogy(rho_b1,smld_err_prb,'-b*',rho_b2,theo_err_rate,'-r');