selfcorr.m

来自「本压缩包括最小均方滤波器的源码」· M 代码 · 共 12 行

function h=selfcorr(y,n)
%计算自相关矩阵
h=zeros(n,n);
m=n-length(y);
x=[y zeros(1,m)];
for i=1:n
   for j=1:n
       for k=1:n-abs(i-j)
           h(i,j)=h(i,j)+x(k)*x(k+abs(i-j));
       end
   end
end