wucha.m

来自「这是一个基于峰度的ICA算法」· M 代码 · 共 9 行

function [error]=wucha(s,w,z)
[vectorSize, numSamples] = size(w);
res=w'*z;
ss=s(1,:)/norm(s(1,:),2);
for p=1:vectorSize
    sr=res(p,:)/norm(res(p,:),2);
    er(p,1)=norm(abs(ss)-abs(sr),2);
end
error=min(er);