试试这个:
t = linspace (0, 500)。^ 3;%的时间向量
s =罪(2 *π* t * 0.2);% 10 Hz正弦曲线
zx =找到(diff(签署(s)));%零交点指数“t”
为k = 1:元素个数(zx)
idx_rng = max (zx (k) 1,1): min (zx (k) + 1,元素个数(t));
tzro (k) = interp1 (s (idx_rng)、t (idx_rng), 0);%的零交叉
结束
zx_rate =梯度(tzro);%计算速度
图
yyaxis左
情节(t, s,“- b”)
持有在
情节(tzro 0(大小(tzro)),“xr”)
持有从
yyaxis正确的
情节(tzro zx_rate,“线宽”,2)
网格
传奇(“信号”,“零交叉”,讨论二阶导数过零率的)
轴(“紧”)
零交点之后,你可以做任何你想要的处理上。
。
