Try this:
t = linspace(0, 5, 500).^3;
s = sin(2*pi*t*0.2);
zx = find(diff(sign(s)));
for k = 1:numel(zx)
idx_rng = max(zx(k)-1, 1):min(zx(k)+1,numel(t));
tzro(k) = interp1(s(idx_rng), t(idx_rng), 0);
end
zx_rate = gradient(tzro);
figure
yyaxis left
plot(t, s, '-b')
hold on
plot(tzro, zeros(size(tzro)), 'xr')
hold off
yyaxis right
plot(tzro, zx_rate, 'LineWidth',2)
grid
legend('Signal','Zero-Crossings','Zero-Crossing Rate')
axis('tight')
Once you have the zero-crossings, you can do whatever processing on them that you want.
.
Best Answer