uniform.m

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

function [y,PI,rs]=uniform()
%
%copyright 2005.6.10
s(1,:)=wavread('man.wav');
s(2,:)=wavread('women.wav');

A=rand(2);
x=A*s;
x=baihua(x);

W=eye(2);
y=W*x;
deta=0.001;
for i=1:100
fy=atan(y);
W=W+deta*(eye(2)-fy*y')*W;
W=W/norm(W);
y=W*x;
end
PI=xnzb(W*A);
rs=xsxs(s,y);