switchica.m

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

function y=switchICA(x)
%改进开关算法
%author： xiaobing zhang
%copyright 2006.3.13


%**** prewhite the data
%预处理
[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

P=t;
w=eye(2);
afa=0.001;
%********迭代程序**********************
for i=1:300
y=w*x;
y1=y(1,:);
J1=peak(y1);
y2=y(2,:);
J2=peak(y2);
J=[J1,J2];
J=diag(J);
sum=0;
B=fix(P/1000);
for k=1:1000
    D=y(:,(k-1)*B+1:k*B)*y(:,(k-1)*B+1:k*B)'+J*tanh(y(:,(k-1)*B+1:k*B)*y(:,(k-1)*B+1:k*B)');
    D=D/B;
    af=diag(diag(D));
    dw=afa*(af-D)*w;
    w=w+dw;
end
W(:,:,i)=w;
end
y=w*x;
y=y';
temp=repmat(mean(y),[P 1]'); 	y=y-temp;
temp=repmat(std(y,1),[P 1]');	y=y./temp;
[dim,sample]=size(y);
if dim>sample
    y=y';
end