%基于峭度的盲分離開關算法
%程序調用說明：該算法程序調用格式為：[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;