mixed.m

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

function [s,x]=mixed()
%
%
%
s1=wavread('why.wav');
s2=wavread('man.wav');
s3=wavread('women.wav');
s(:,1)=s1(1:5000);
s(:,2)=s2(1:5000);
s(:,3)=s3(1:5000);
ss=s';
for i=1:3
	s=ss(i,:);
	s=s-mean(s);
	s=s/std(s,1);
	ss(i,:)=s;
end; s=ss';
[dim,sample]=size(s);
A=randn(sample);
x=s*A;