Matlab找到信号周期的范围

Question

我有一个信号在某些地方是周期性的而在其他地方没有，我希望能够找到它周期性的范围（及时）。我不能在这里显示我的原始信号，但我会用一个示例信号来说明我的问题：

示例信号：

vect=[randn([1,500]) sin(x) 1:500]; 
x=linspace(0, 20*Pi, 1000) 

我想要的东西会告诉我信号是x周期性的，介于500到1500之间。

我尝试使用xcorr函数（我使用xcorr(y,'unbiased')然后找到了零滞后峰值，发现其他峰值在零滞后的某个％范围内，以定义周期性区域，但我无法弄清楚如何将滞后与x范围联系起来。

编辑：我目前使用的代码

    [c, lags] = xcorr(y,'unbiased');
    lag_zero=find(lags==0)
    [peaks,locs]=findpeaks(c,'MINPEAKHEIGHT',.5*c(lag_zero)); %finding peaks 
    cindex=find(((c(lag_zero)-(0.5*c(lag_zero)))<c) & (c<(c(lag_zero)+(0.5*c(lag_zero)))));
    maxlags=lags(max(cindex));
    [c2,lags2]=xcov(y,y,maxlags,'unbiased'); %this is just to narrow the periodic part down
    plot(lags2,c2);
    period=abs(x(locs(floor((length(locs))/2)))-x(locs(floor(((length(locs))/2)-1))))

Answer 1

我认为最简单的起点是窗口数据，然后使用像this这样的技术来查找功率高于某个本底噪声的频率。这种分析很像季节性，但有不同的时期，但这不是问题。从Cross Validated中查看this link。

Answer 2

运行[~,F,T,P] = spectrogram(vect, 128, 120, 128, 1)即可获得

通过视觉检查，您可以看到元素500和1500之间有单色信号。您可以检查P以查找值非常低的矩形区域，可能更好地使用10*log10(P)，并使用F和T来提取间隔和频率。