function [ys,W]=SNR_Max(x)
[n,T]=size(x);
x=x-mean(x')'*ones(1,T);%去均值
%------------------白化处理-------------------------------
[F,D]=eig((x*x')/T);
v=F*(D^(-0.5))*F';
x=v*x;
p=80;
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)';%原信号估计
mplot(x);
figure(2);
mplot(ys);
评论0