在MATLAB中,给定数据向量,在阈值内找到过零点

时间:2016-05-31 11:27:03

标签: matlab time signal-processing feature-detection

我在MATLAB中编写了一个函数来计算给定信号数据矢量的过零次数。如果我发现零交叉,我还要检查所涉及的两个矢量索引之间的绝对差值是否大于阈值 - 这是为了减少信号噪声的影响。

    zc = [];
    thresh = 2;
    for i = 1:length(v)-1
        if ( (v(i)>0 && v(i+1)<0) || (v(i)<0 && v(i+1)>0) ) && abs(v(i)-v(i+1)) >= thresh
            zc = [zc; i+1];
        end
    end
    zcCount = length(zc);

我使用了零交叉函数中的向量来测试它:http://hips.seas.harvard.edu/content/count-zero-crossings-matlab

    A = [-0.49840598306643,
         1.04975509964655,
        -1.67055867973620,
        -2.01437026154355,
         0.98661592496732,
        -0.06048256273708,
         1.19294080740269,
         2.68558025885591,
         0.85373360483580,
         1.00554850567375];

它似乎工作正常,但是有更有效的方法来实现相同的结果吗?例如。在上面的网页上,他们只需使用以下行来计算过零点:

    z=find(diff(v>0)~=0)+1;

有没有办法将阈值检查纳入同样有效的方法?

1 个答案:

答案 0 :(得分:2)

怎么样

zeroCrossIndex=diff(v>0)~=0
threshholdIndex = diff(v) >= thresh;
zcCount = sum(zeroCrossIndex & threshholdIndex)