maxinfo.m

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

function [y,w]=maxinfo(x)
%2006年1月11号
%修改版
[N,P]=size(x);
mixes=x;
mx=mean(mixes'); 
c=cov(mixes');
x=x-mx'*ones(1,P);         % subtract means from mixes
wz=2*inv(sqrtm(c));        % get decorrelating matrix
x=wz*x;                    % decorrelate mixes so cov(x')=4*eye(N);
L=0.0001;
B=30;
M=N;
w=rand(M,N);
noblocks=fix(P/B);
Id=eye(N);
%BI=B*Id;
BI=Id;
%for i=1:150
%    for t=1:B:noblocks*B,
%      u=w*x(:,t:t+B-1); 
%      w=w+L*(BI+(1-2*(1./(1+exp(-u))))*u')*w;
%    end
%end

    for t=1:P
      u=w*x(:,t); 
      w=w+L*(BI+(1-2*(1./(1+exp(-u))))*u')*w;
    end

%mixes=a*sources;     
%plot(mixes(1,:))     
y=w*wz*mixes; 
w=w*wz;
%plot(y(1,:))   
%xlabel('y1');