盲源分离AMUSE算法matlab代码
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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