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;