查找并跟踪脊使用重分配谱图

加载包含由大棕蝠发出的回声定位脉冲数据文件(大棕鲇)中并用7微秒的采样间隔进行测量。创建使用信号和时间信息的MATLAB®时间表。

加载batsignalT =(0:长度(batsignal)-1)* DT;SG =时间表(秒(t)的”,batsignal);

打开信号分析仪并将时间表从工作区浏览器拖动到信号表。点击显示网格▼创建两个并排侧显示器。选择每个展示和点击时频按钮添加一个频谱图。

拖动时间表两个显示器。

选择光谱图标签。在右边的显示,请检查重新分配。对于每个显示:

  • 将时间分辨率设置为280微秒,并指定相邻段之间85%的重叠。

  • 使用泄漏滑块增加泄漏,直到RBW约为4.5 kHz。

  • 将电源限制-45 dB和-20分贝。

重新分配的频谱图清楚地显示出三个时间频脊。为了跟踪脊,选择右边的显示。在显示选项卡,单击生成脚本并选择谱图脚本。该脚本编辑器中显示。

%计算频谱%的MATLAB(R)9.7和信号处理工具箱8.2生成。26日 -  12月2018 17时21分44秒:%上生成%的参数时间限制=秒([0 0.002793]);%秒frequencyLimits = [0 71428.57]。%赫兹泄漏= 0.9;timeResolution = 0.00028;%秒overlapPercent = 85;reassignFlag = TRUE;%%%索引关心的信号时区sg_batsignal_ROI = SG(:,'batsignal');sg_batsignal_ROI = sg_batsignal_ROI(TIMERANGE(timeLimits(1),期限(2),'关闭'),1);%计算谱估计%运行下面的函数调用,而不输出参数绘制结果[P F T] = pspectrum (sg_batsignal_ROI,...“频谱”...'FrequencyLimits',frequencyLimits,...'泄漏',泄漏,...“TimeResolution”,timeResolution,...'OverlapPercent'overlapPercent,...“再分配”,reassignFlag);

运行该脚本。剧情重新分配的频谱。

目(秒(T),F,P)xlabel('时间')ylabel (“频率”)轴视图(2)的颜色表

使用tfridge功能跟踪脊。

[冰箱,〜,lridge] = tfridge(P,F,0.01,'NumRidges',3,'NumFrequencyBins',10);保持plot3(秒(T),冰箱,P(lridge),“:”“线宽”3)保持

感谢伊利诺斯大学贝克曼中心的柯蒂斯·康顿、肯·怀特和艾尔·冯提供蝙蝠数据,并允许在本例中使用。

也可以看看

应用

职能

相关的例子

更多关于