盲信号分离算法 matlab

提供一种盲分离的matlab实现，适用语音信号和一些其他的混合信号。

function y=Goge1（x）
 t=0:0.01:10; 
s1=cos（2*pi*t）; 
[NnNp]=size（s1）; 
s2=randn（1Np）; 
s=[s1;s2];
a=rand（22）;
x=a*s; 
x1=x（1:）; 
x2=x（2:）; 
figure 
subplot（211）;plot（s1）;axis tight;%title（‘混合以前的独立信号s1‘）; 
subplot（212）;plot（s2）;axis tight;%title（‘混合以前的独立信号s2‘）; 
figure 
subplot（211）;plot（x1）;axis tight;%title（‘混合以后的观测信号x1‘）; 
subplot（212）;plot（x2）;axis tight;%title（‘混合以后的观测信号x2‘）; 

[nT]=size（x）;
if T>1000
    P=1000;
end
x = x-mean（x‘）‘*ones（1T）;
[FD]=eig（（x*x‘）/T）;
v=F*（D^（-0.5））*F‘;
x=v*x;
a=ones（1P）/P;
g=filter（a1x‘）;
U=cov（x‘1）;
V=cov（g1）;
[WD]=eig（VU）;
%get independent component 
y=x‘*W;
figure 
subplot（311）;plot（y）;axis tight;
subplot（312）;plot（y（:1））;axis tight;%title（‘经过ICA后的独立成分y1‘）; 
subplot（313）;plot（y（:2））;axis tight%title（‘经过ICA后的独立成分y2‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        881  2010-05-25 16:48  Goge1.m

----------- ---------  ---------- -----  ----

                  881                    1