wucha.m

这是一个基于峰度的ICA算法

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);