kurtscum.m

来自「混沌生成的很多形式」· M 代码 · 共 22 行

function k = kurtscum(x)
[row, col] = size(x);

if row == 1 & col == 1
    k = 0;
else
   mu=mean2(x);
   m1=mean2((x-mu).^1);
   m2=mean2((x-mu).^2);
   m3=mean2((x-mu).^3);
   m4=mean2((x-mu).^4);
   %m1 = moment(x,1);
   %m2 = moment(x,2);
   %m3 = moment(x,3);
   %m4 = moment(x,4);
   sm2 = sqrt(m2);
   k1 = m4-3*(m2).^2;
   k2 = 12*m1.^2.*m2;
   k3 = 4*m1.*m3+6*m1.^4;
   k =(k1+k2-k3)./(sm2+eps).^4;
end