nomax.asv

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

function [ys,W] = NoMax(x)
%最大信噪比盲源分离算法
%输入：混合信号x
%输出：分离信号ys，分离矩阵W
%preprocess
%tic;
x=baihua(x);
%%%%%%%moving average coefficient%%%%%%%%%%%%
p=80;
a=ones(1,p)/p;
x=x';
S=filter(a,1,x);
%computer covariance
U=cov(S-x,1);
V=cov(x,1);
%computer eigenvector
[W d]=eig(V,U); 
%recover signal
ys = (x*W)';
%t=toc