signintermittency.m

一种新的时频分析方法的matlab源程序。

function y=signIntermittency(x,maxConsec)
% y = signIntermittency(x,maxConsec) - remove consecutive signs greater than max
%  by Kenneth Arnold <[e-mail]>

y=x;
x=sign(x);

curSign = x(1);
signStart = 1;
for i=2:max(size(x))
    if x(i) ~= curSign
        if signStart ~= 1 & curSign ~= 0 & i-signStart >= maxConsec
            y(signStart:i) = 0;
        end
        curSign = x(i);
        signStart = i;
    end
end