findzc.m

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

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) = [];