amari.m

来自「互信息盲源分离」· M 代码 · 共 19 行

function out=amari(C,A)

[b,a]=size(C);

dummy=pinv(A)*C;
dummy=sum(ntu(abs(dummy)))-1;

dummy2=pinv(C)*A;
dummy2=sum(ntu(abs(dummy2)))-1;

out=(sum(dummy)+sum(dummy2))/(2*a*(a-1));


function CN=ntu(C)
[m n]=size(C);
for t=1:n,
 CN(:,t)=C(:,t)./max(abs(C(:,t)));
end