bksa.m

来自「针对线性混合信号盲分离时」· M 代码 · 共 33 行

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);
w=[-0.1698 0.5757;0.3878 -0.1883]; %[-0.1690 0.5753 0.3875 -0.1871]
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;