kurt.m

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

function k=kurt(x)
%如果x是一个列向量，则返回该列向量的峰度，
%如果x是一个矩阵，则返回其列向量的的峰度。
[m,n]=size(x);
k=[];
for v=1:n
    a=x(:,v);
    k(v)=mean(a.^4)-3*(mean(a.^2))^2;
end