selfcorrelate.m

来自「滤波器的设计以及零极点对应频率响应分析」· M 代码 · 共 8 行

function r=selfcorrelate(s)
l=length(s);
tmp=zeros(size(s));
rs=[s;s];
for i=0:floor(l/2)
  r(i+1)=s(1:l)'*rs((1+i):(l+i));
end
end