x=rand(1,100000);

pds=[50,200,400];
for p=1:length(pds)
  pd=pds(i);
  x(pd:pd:end)=10;
end


[ac,lag]=xcorr(x,'coeff');

plot(ac,lag)
