盲源分离。。。。。

alwaysfly

　　有用盲源分离的没

　　讨论下。。。。。。

diligency

有啊
我是啊

AaronSpark

　　利用特征值分解的盲源分离

　　function [W] = amuse(X)

　　% BSS using eigenvalue value decomposition

　　% Program written by A. Cichocki and R. Szupiluk

　　%

　　% X [m x N] matrix of observed (measured) signals,

　　% W separating matrix,

　　% y estimated separated sources

　　% p time delay used in computation of covariance matrices

　　% optimal time-delay default p= 1

　　%

　　% First stage: Standard prewhitening

　　[m,N]=size(X);

　　if nargin==1,

　　n=m; %

　　end;

　　Rxx=(X*X')/N;

　　[Ux,Dx,Vx]=svd(Rxx);

　　Dx=diag(Dx);

　　% n=xxx;

　　if n

　　%when the number of sources are known or can a priori estimated

　　Dx=Dx-real((mean(Dx(n+1:m))));

　　Q= diag(real(sqrt(1./Dx(1:n))))*Ux(:,1:n)';

　　%

　　else % under assumption of no additive noise and when the

　　% number of sources is unknown

　　n=max(find(Dx>1e-199)); %Detection the number of sources

　　Q= diag(real(sqrt(1./Dx(1:n))))*Ux(:,1:n)';

　　end;

　　%

　　% else %assumes no noise

　　% Q=inv(sqrtm(Rxx));

　　% end;

　　%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

　　% Second stage: Fast separation using sorting EVD

　　% notation the same as used in the Chapter 4

　　Xb=Q*X;

　　p=1;

　　% paramter p can take here value different than 1

　　% for example -1 or 2.

　　N=max(size(Xb));

　　Xb=Xb-kron(mean(Xb')',ones(1,N));

　　Rxbxbp=(Xb(:,1:N-1)*Xb(:,2:N)')/(N-1);

　　Rxbxbp= Rxbxbp+Rxbxbp';

　　[Vxb Dxb]=eig(Rxbxbp);

　　[D1 perm]=sort(diag(Dxb));

　　D1=flipud(D1);

　　Vxb=Vxb(:,flipud(perm));

　　W = Vxb'*Q;

　　%y = Vxb' * x1;

AaronSpark

你具体做什么问题？用什么算法实现盲源分离？

alwaysfly

　　初学
想用一种盲源分离方法
分离信号中的周期成分和
瞬态成分

找到一种
利用短时傅立叶变换的方法

可惜对信号处理一窍不通
导致进展缓慢

想请教一下
你们一般利用盲源分离处理什么问题
采用什么分离方法
是专门学信号的
还是辅修的

diligency

我也是，我是用盲源信号理论做图像处理的

kobe735

回复 5 # alwaysfly 的帖子

alwaysfly,你能否把利用短时傅立叶变换的进行盲源分离的程序发给我，先谢谢啦！邮箱caojunhong735@126.com