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;