signal.m

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

function signal(A)
s1=wavread('man.wav');
s2=wavread('women.wav');
s=[s1';s2'];
[n,m]=size(s);
A=rand(n);
x=A*s;
[y,PI]=based_neg(x);

figure(1);
plot(PI);
figure(2);
for i=1:n
    subplot(n,1,i);
    plot(y(i,:));
end