baihua.m

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

function y=baihua(x)
%观测信号的白化处理函数
%包括均值为零处理
%去均值时，
[n,T] = size(x);

%if nargin == 1
 %   m = n;
%end

x = x - mean(x')'*ones(1,T);% Remove mean

%%***************白化*********************

[F,D]=eig((x*x')/T);
%[vals,idxs]	= sort(diag(D));
%mmvar = n-m+1:n; 
%F=F(:,idxs(mmvar))';
%D=D(idxs(mmvar),idxs(mmvar));
%v = sqrt(diag(1./diag(D)))  * F;

v=F*(D^(-0.5))*F';
y=v*x;
%修改日期：2005.8.8