bksa.m

针对线性混合信号盲分离时

%基于峭度的盲分离开关算法
%程序调用说明：该算法程序调用格式为：[y,w]=bksa(x),其中x是一个n*T的数据矩阵（n个传感器，T个采样点），
% y是估计出的源信号矩阵（n*T），w是n*n分离矩阵。

function [y,w]=bksa(x)
[n,T]=size(x);
x=x-mean(x')'*ones(1,T); %去均值
%-------白化处理-----------
[F,D]=eig((x*x')/T);
v=F*(D^(-0.5))*F';
x=v*x;
x1=x(1,:);
x2=x(2,:);
w=eye(2);

afa=0.001;
%-------主要程序-----------
for i=1:30
    y=w*x;
    y1=y(1,:);
    k1=peak(y1);
    J1=sign(k1);
    y2=y(2,:);
    k2=peak(y2);
    J2=sign(k2);
    J=[J1,J2];
    J=diag(J);  %开关矩阵
    sum=0;
 for k=1:1000
     D(:,:,k)=y(:,k)*y(:,k)'+J*tanh(y(:,k)*y(:,k)');
     af(:,:,k)=diag(diag(D(:,:,k)));
     dw=afa*(af(:,:,k)-D(:,:,k))*w;
     w=w+dw;
  end
end
w;
y=w*x;