% 定義非周期相關(guān)函數(shù)，t值可以為負(fù)值
function r=correlate(a,b,t)
   n=length(a);
   h=0;
   if (t>=0)
      for i=1:n-t
           h1=a(i)*conj(b(i+t));
           h=h+h1;
      end
   else
       for i=1:n+t
           h1=conj(b(i))*a(i-t);
           h=h+h1;
       end
   end    
   r=h;