%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;    enddisp(strcat('EASI_ECT:',num2str(ECT(ind))));