awgn_per.m

function[prob]=awgn_per(snr,T)
s=10^(snr/10);
ec=(2/T)*s;
n=10000;
num=0;
for i=1:n
    t=round(rand(1,1));
    if t==0;
        d=-ones(1,T);
    else
        d=ones(1,T);
    end

    redata=[d];
    
    pn=randsrc(1,T);
     
   trs=sqrt(ec)*redata.*pn;
   
   trs=trs+randn(1,T);
   de=trs.*pn;
   dc=sum(de);
   if(dc<0)
       deci=-1;
   else
       deci=1;
   end
   if(deci~=d)
       num=num+1;
   end
end
   prob=num/n;