decorrection.m

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

function y=decorrection(x)
%2006.3.13
%copyright:xiaobing zhang

%预处理
[n,t] = size(x);
m=mean(x')'*ones(1,t);
x = x - m;% Remove mean
[F,D]=eig((x*x')/t);
v=F*(D^(-0.5))*F';
x=v*x;%white

i=22;
%main code 
Q=x(:,1:t-i)*x(:,i+1:t)'+x(:,i+1:t)*x(:,1:t-i)';
[W,D]=eig(x*x',Q);
%separating
y=W'*x;