sync_dl=[1 0 1 1 0 0 1 1 1 0 1 0 0 1 1 1 1 1 0 0 1 1 0 0 0 0 0 0 0 1 0 1 1 0 1 0 1 0 0 1 1 0 0 0 0 1 1 0 1 0 0 0 1 0 0 0 1 1 1 0 0 1 0 0];
sync_ul=randsrc(1,128);
gp1=zeros(1,32);
gp2=zeros(1,96);
m=1;
%%%%%%%%%%%%%%%%%%%
for i=1:64;
    sync_dl(i)= sync_dl(i)+(sync_dl(i)==0)*(-1);
end
%%%%%%%%%%%%%%%%%%%%%%%%
d=[gp1,sync_dl,gp2,sync_ul,gp1];
q=randsrc(1,32000);
k=100;
%%%%%%%%%%%%%
for i=0:4
    a=k+i*6400;
for j=1:352;
    q(a+j)=d(j);
end
end
%%%%%%%%%%%%%%
for snr=-0.8:0.01:0.8
%%%%%%%%
 b=0;
for n=1:500;
 q2=awgn(q,snr);
 q3=reshape(q2,6400,5);
 y=zeros(1,6400);
%%%%%%
 for i=1:5
     y=y+q3(:,i)';
 end
%%%%%%
 for i=1:6273;
    p1=0;
    p2=0;
    p3=0;
  for j=0:31;%
     p1=p1+y(i+j)^2;%
     p2=p2+y(i+96+j)^2;%
  end%
    for t=32:95;
      p3=p3+y(i+t)^2;
    end
 z(i)=(p1+p2)/p3; 
 end
%%%%%%%%%%%%%
[c,d]=min(z);
  if (d==k+1)
    b=b+1;
  else
    b=b;
  end
end
p(m)=b/500;
m=m+1;
end
snr=-0.8:0.01:0.8;
plot(snr,p);