four_sig.m

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

function [rs,PI] = four_sig()
%
%
%two supergauss signal
%two subgauss signal
s1=wavread('man.wav');
s1=s1';
s2=wavread('women.wav');
s2=s2';
t=0:pi/20000:8*pi;
t=t(2:160001);
s3=sin(t);
s4=randn(1,160000);
s4=sort(s4);
s=[s1;s2;s3;s4];
s=baihua(s);
[dim,sample]=size(s);
if dim>sample
    s=s';
end
%A=randn(4);
A =[0.6947   -0.5441    0.5320    0.2603;...
   -0.3038   -0.4080   -0.4754   -1.0665;...
    0.6142    0.0096    0.9295    1.2096;...
    0.5773   -1.1813   -1.0454   -0.1655];
x=A*s;
%*****************************************
[y, A, W] = fastica(x);
%*******************************
rs=xsxs(s,y);
PI = xnzb(rs);