mixdata.m

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

function mixdata

下面是产生混合数据的代码
把需要的数据种类的产生代码拷贝到MATLAB命令窗口，
然后在系统窗口中点击“加载数据”按钮，把混合数据
变量输入加载数据对话框中，点击加载就可以。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
观测信号和源信号一样多的情况
%超高斯亚高斯混合信号
t=1:10000;
load chirp
y=y';
s(1,:)=y(1:10000);
s(2,:)=sin(t/300);
A=rand(2);
x=A*s;
%超高斯（语音）混合信号
s(1,:)=wavread('man.wav');
s(2,:)=wavread('women.wav');
A=rand(2);
x=A*s;
%亚高斯（正弦）信号
t=1:1000;
s(1,:)=sin(t/30);
s(2,:)=0.2*cos(t/20)*sin(t/40);
A=rand(2);
x=A*s;

观测信号大于源信号的情况
%超高斯亚高斯混合信号
t=1:10000;
load chirp
y=y';
s(1,:)=y(1:10000);
s(2,:)=sin(t/300);
A=rand(3,2);
x=A*s;
%超高斯（语音）混合信号
s(1,:)=wavread('man.wav');
s(2,:)=wavread('women.wav');
A=rand(3,2);
x=A*s;
%亚高斯（正弦）信号
t=1:1000;
s(1,:)=sin(t/30);
s(2,:)=0.2*cos(t/20)*sin(t/40);
A=rand(3,2);
x=A*s;