baihua.m

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

function y=baihua(x)
%观测信号的白化处理函数
%包括均值为零处理
%去均值时，信号为列向量
[dim,sample]=size(x);
if dim<sample
    m=sample;
    x=x';
else
    m=dim;
end
xc = x - repmat(sum(x)/m,m,1);  % Remove mean
%%***************白化*********************
%白化时，信号为行向量
x=xc';
[F,D]=eig(x*x'/length(x(1,:)));
v=F*(D^-(0.5))*F';
y=v*x;