easi.m

来自「盲信号分离(BSS)是指在对彼此独立的源信号混合过程及各源信号本身均未知的情况下」· M 代码 · 共 28 行

%EASI算法
function [W ECT]=EASI(X,A,u)

N=size(X,2);
I=eye(size(A));
W=I;%W（0）

ECT=zeros(1,N);

for ind=1:N
    y=W*X(:,ind);
    
    %gy=y.^3;
    gy=y-tanh(y);
    
    W=W - u*(y*y' - I + gy*y' - y*gy')*W;
    
    C=abs(W*A);
    [height width]=size(C);
    mci=max(C,[],2);
    mcj=max(C,[],1);
    mci=repmat(mci,1,width);
    mcj=repmat(mcj,height,1);
    ECT(ind)=sum(sum(C./mci))+sum(sum(C./mcj))-height-width;
    
end

disp(strcat('EASI_ECT:',num2str(ECT(ind))));