peak.m

针对线性混合信号盲分离时

%------计算峭度的子程序-------
function k=peak(x)
  [n,T]=size(x);
  x=x-mean(x')'*ones(1,T);
  k4=mean(x.^4);
  k2=mean(x.^2);
  k=k4/(k2^2)-3;
  k=sign(k);