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

function [p]=ss_Pe96(rho_in_dB); 
 
rho=10^(rho_in_dB/10); 
Eb=rho; 
if rho>2 
    alpha=2/rho; 
else 
    alpha=1; 
end 
sgma=sqrt(1/2/alpha); 
N=10000; 
for i=1:N 
    temp=rand; 
    if temp<0.5 
        data(i)=1; 
    else 
        data(i)=0; 
    end 
end 
 
for i=1:N 
    if data(i)==0 
        r1c(i)=sqrt(Eb); 
        r1s(i)=0; 
        r2c(i)=0; 
        r2s(i)=0; 
    else 
        r1c(i)=0; 
        r1s(i)=0; 
        r2c(i)=sqrt(Eb); 
        r2s(i)=0; 
    end 
    if randr2 
        decis=0; 
    else 
        decis=1; 
    end 
    if decis~=data(i) 
        num_of_err=num_of_err+1; 
    end 
end 
p=num_of_err/N;