getsignconsec.m

来自「一种新的时频分析方法的matlab源程序。」· M 代码 · 共 22 行

function [y,yidx]=getSignConsec(x)
% y = getSignConsec(x) - returns vector of consective sign counts
%  by Kenneth Arnold <[e-mail]>

y = [];
yidx = [];

x=sign(x);

curSign = x(1);
signStart = 1;
for i=2:max(size(x))
    if x(i) ~= curSign
        if signStart ~= 1 & curSign ~= 0
            y = [y i-signStart];
            yidx = [yidx (signStart+i)/2];
        end
        curSign = x(i);
        signStart = i;
    end
end