inde.m

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

function y=Inde(x)
%copyright 2005.7.5
%A New Approach to Blind Source Separation with Global Optimal Property
%matlab code author: 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

[N,P]=size(x);
if P>10000
    P=10000;
end

a=ones(1,P)/P;
g=filter(a,1,x');
U=cov(x',1);
V=cov(g,1);
[W,D]=eig(V,U);
y=x'*W;
y=y';