correlate.m

来自「实验发现」· M 代码 · 共 19 行

% 定义非周期相关函数，t值可以为负值
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;