two_signal.m

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

function PI = two_signal()

s1=wavread('man.wav');
s2=wavread('women.wav');
s=[s1';s2'];
%A=[1.5 1;1 1.4];
[dim,sample]=size(s);
n=1:sample;
t=n/16000;
A=rand(2);
x=A*s;
[y,B,W]=fastica(x);
fprintf('dim is %d\n',dim);
fprintf('sample is %d',sample);
G=W*A;
PI=xnzb(G);
subplot(221);
plot(t,s1');
subplot(223);
plot(t,s2');
xlabel('s');
subplot(222);
plot(t,y(1,:));
subplot(224);
plot(t,y(2,:));
xlabel('y');