盲源分离。。。。。
本帖最后由 VibInfo 于 2016-10-11 15:02 编辑有用盲源分离的没
讨论下。。。。。。
本帖最后由 VibInfo 于 2016-10-11 15:02 编辑
有啊
我是啊
回复：（alwaysfly）盲源分离
本帖最后由 VibInfo 于 2016-10-11 15:02 编辑利用特征值分解的盲源分离
function = amuse(X)
% BSS using eigenvalue value decomposition
% Program written by A. Cichocki and R. Szupiluk
%
% X 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
=size(X);
if nargin==1,
n=m; %
end;
Rxx=(X*X')/N;
=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';
=eig(Rxbxbp);
=sort(diag(Dxb));
D1=flipud(D1);
Vxb=Vxb(:,flipud(perm));
W = Vxb'*Q;
%y = Vxb' * x1;
回复：（alwaysfly）盲源分离
你具体做什么问题？用什么算法实现盲源分离？ 本帖最后由 VibInfo 于 2016-10-11 15:03 编辑初学
想用一种盲源分离方法
分离信号中的周期成分和
瞬态成分
找到一种
利用短时傅立叶变换的方法
可惜对信号处理一窍不通
导致进展缓慢
想请教一下
你们一般利用盲源分离处理什么问题
采用什么分离方法
是专门学信号的
还是辅修的
本帖最后由 VibInfo 于 2016-10-11 15:03 编辑
我也是，我是用盲源信号理论做图像处理的 回复 5 # alwaysfly 的帖子
alwaysfly,你能否把利用短时傅立叶变换的进行盲源分离的程序发给我，先谢谢啦！邮箱caojunhong735@126.com
页:
[1]