bksa.m

独立分量分析的开关算法代码

function y=bksa(x)
%
%基于峰度的盲源开关算法
%系统仿真学报 2005.1
%copyrigth：2005.4.13
%author：lucky zhang

%*********init******************
%fs=8000;
% f=14.324;
% load chirp;
% y=y(1:1000);
% k=1:1000;
%s1=sin(2*pi*f*k/fs);
%s2=y';
%s=[s1;s2];
%A=[1.5 1;1 1.4];
%A=randn(2);
%x=A*s;


x=baihua(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
%%%%%%%%%%%性能指标%%%%%%%%%%%%%%%%%%%%%%


%%P=w*A;
%n=length(P);
%c=0;
%for j=1:2
%    a=0;
%    b=0;
%    for k=1:2
%    a=a+abs(P(j,k))/max(abs(P(j,:)));
%    b=b+abs(P(k,j))/max(abs(P(:,j)));
%    end
%    a=a-1;
%    b=b-1;
%   c=c+a+b;
%end
%PI(i)=c/(n*(n-1));
%xn(i)=xnzb(P);
end
y=w*x;

%*********快速ICA算法*******************


%********画图程序**********************

%subplot(421);
%plot(s1);
%ylabel('s1');
%subplot(423);
%plot(s2);
%ylabel('s2');
%subplot(425);
%plot(x1);
%ylabel('x1');
%subplot(427);
%plot(x2);
%ylabel('x2');
%subplot(422);
%plot(y1);
%ylabel('y1');
%text(400,1, '本算法')
%subplot(424);
%plot(y2);
%ylabel('y2');