test.m

来自「这是盲信号的代码 都已经通过编译了 做这方面的同仁可以参考一下 我觉得蛮惯用的」· M 代码 · 共 26 行

function [y,PI]=test()
%preprocess

load wtwo;
A=randn(2);
x=A*s;

x=baihua(x);
%%%%%%%moving average coefficient%%%%%%%%%%%%
x=x';
for p=2:500
tic;
a=ones(1,p)/p;
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)';
r=xsxs(s,ys);
PI(p-1)=xnzb(r);
t=toc;
end
y=ys;