nomax.m

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

function ys = NoMax(x)
%最大信噪比盲源分离算法
%输入：混合信号x
%输出：分离信号ys，分离矩阵W

%计算运行时间tic   t=toc
%tic;
%预处理，去均值和白化
%x=baihua(x);
[n,T] = size(x);
m=mean(x')'*ones(1,T);
x = x - m;% Remove mean
[F,D]=eig((x*x')/T);
v=F*(D^(-0.5))*F';
x=v*x;
%%%%%%%滑动平均处理%%%%%%%%%%%%
p=8;
a=ones(1,p)/p;
x=x';
S=filter(a,1,x);
%计算协方差
U=cov(S-x,1);
V=cov(x,1);
%计算特征值
[W d]=eig(V,U); 
%源信号估计
ys = (x*W)';
%t=toc;