peak.m

这是盲信号的代码 都已经通过编译了 做这方面的同仁可以参考一下 我觉得蛮惯用的

function k=peak(x)    
%author lucky zhang
%copyright 2005.6.21
%
[n,T] = size(x);
x = x - mean(x')'*ones(1,T);% Remove mean


    k4=mean(x.^4);
    k2=mean(x.^2);
    k=k4/(k2^2)-3;
   k=sign(k);