谢谢你引用我的代码!
我不确定你打算通过“中位数”零交叉。简单的计算更精确的零交叉使用简单的线性插值。(您可以使用
interp1
为此,但是我的代码可能会更有效的计算对于这个简单的问题。)
的代码
- - - - - -
t = linspace(0, 20 *π,5000);
y =罪(2 *π* t);
zci = @ (v)找到(v (:)。* circshift (v (:), [1 0]) < = 0);%返回零交点参数向量的指标
zx = zci (y);
为k1 = 1:元素个数(zx) 1%插入计算精确的零交叉
电视= t (zx (k1): zx (k1) + 1);%独立变量向量
青年志愿= y (zx (k1): zx (k1) + 1);%因变量向量
b =[[1; 1]、电视()):\青年志愿(:);%线性插值
ti (k1) = b - b (1) / (2);%实际零交点(独立变量)的值
结束
图(1)
情节(t, y)
持有在
情节(t (zx), y (zx),“pg”)
情节(ti, 0(大小(ti)),“+ r”)
持有从
图(2)%放大图(可选)
情节(t, y)
持有在
情节(t (zx), y (zx),“pg”)
情节(ti, 0(大小(ti)),“+ r”)
持有从
轴([60 63 ylim])
