www.pudn.com > zxb.rar > ExInfomax.m

function y=ExInfomax(x) 
 
% 
%扩展信息最大化 
%编程时间：2006.1.12 
%修改：2006.3.13 
%作者：张小兵 
%*********初始化程序****************** 
%超高斯亚高斯 
%fs=8000; 
%f=14.324; 
%load chirp; 
%y=y(1:1000); 
%k=1:1000; 
%s1=sin(2*pi*f*k/fs); 
%s2=y'; 
%s=[s1;s2]; 
%加载两语音 
%load wtwo 
%s=s(:,5000:25000-1); 
%s1=s(1,:); 
%s2=s(2,:); 
%混合 
%A=[1.5 1;1 1.4]; 
%x=A*s; 
%预处理 
[n,t] = size(x); 
m=mean(x')'*ones(1,t); 
x = x - m;% Remove mean 
[F,D]=eig((x*x')/t); 
v=F*(D^(-0.5))*F'; 
x=v*x;%white 
 
x1=x(1,:); 
x2=x(2,:); 
w=eye(2); 
afa=0.001; 
B=200; 
[N,P]=size(x); 
%********迭代程序********************** 
for i=1:200 
    noblocks=fix(P/B); 
    for t=1:B:noblocks*B,   
        y=w*x(:,t:t+B-1); 
        y1=y(1,:); 
        k1=peak(y1); 
        J1=sign(k1); 
        y2=y(2,:); 
        k2=peak(y2); 
        J2=sign(k2); 
        J=[J1,J2]; 
        J=diag(J); 
        D=y*y'+J*tanh(y)*y'; 
        af=diag(diag(D)); 
        %af=eye(2); 
        dw=afa*(af-D)*w; 
        w=w+dw; 
    end    
end 
 y=w*x;