kurt.m

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

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