selfcorrelate.m

滤波器的设计以及零极点对应频率响应分析

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