main.m

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

function PI=main(s)
%
%2005.10.31
%load wtwo;
tic;

[dim,sample]=size(s);
B=8000;
N=sample/B;
%A=rand(dim);
%for p=2:80
for i=1:N
    u=s(:,(i-1)*B+1:i*B);
    x=rand(dim)*u;
    y=NoMax(x);
    r=xsxs(u,y);
    PI(p,i)=xnzb(r);
end
%end
t=toc
%plot(PI);