盲源分离AMUSE算法matlab代码

function [ W,Y ,d,D] = Aamuse(X) [N,M]=size(X);

for i=1:N;

X(i,:)=X(i,:)-mean(X(i,:)); end

S=X;

r1=S*S'/M;

[v d]=eig(r1);

for i=1:N

d(i,i)=1/sqrt(d(i,i));

end

Q=d*v';

y=Q*S;

r2=y(:,1:M-10)*y(:,11:M)';

z=(r2+r2')/2;

[V D]=eig(z);

W=inv(Q)*V;

Y=inv(W)*X;

end