signal.m

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

function s = signal()
%
%
%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);