findzc.m

这是一个关于hht变换很有用的工具箱

function zcx = findzc(in)

npts = length(in);
zcx = zeros(1,npts/2);
wasPositive = in(1) > 0;
idx = 1;

for i=2:npts
    isPositive = in(i) > 0;
    
    if wasPositive ~= isPositive
        zcx(idx) = i + in(i)/(in(i-1)-in(i));
        idx = idx+1;
        wasPositive = isPositive;
    end
end

% clear unused points
zcx(idx:end) = [];