pcorrelate.m

实验发现

% 定义周期相关函数
function r=pcorrelate(a,b,t)
   n=length(a);
   h=0;
   for i=1:n-t
       h1=a(i)*conj(b(i+t));
       h=h+h1
   end
   for m=1:t
       h2=conj(b(m))*a(m+n-t);
       h=h+h2
   end
   r=abs(h);