covar.asv

这是盲信号的代码 都已经通过编译了 做这方面的同仁可以参考一下 我觉得蛮惯用的

function C=covar(x,y)
%
%2005.8.27
%C(i)=E(x(t)*y(t+i))-E(x(t))*E(y(t+i));
lx=length(x);
ly=length(y);
if lx ~= ly
    fprintf('length of vector must be equal\n');
    return ;
end
Y=[y,y];
for i=1:lx
    y=Y(i:lx+i-1);