主要内容

使用重新分配频谱图查找和跟踪脊

加载包含由大棕色蝙蝠发出的回声机脉冲的数据文件(Eptesicus fuscus.)并以7微秒的采样间隔测量。使用信号和时间信息创建MATLAB®时间表。

加载batsignalT =(0:长度(Batsignal)-1)* DT;SG =时间表(秒(t)',击球率);

打开信号分析仪并将时间表从Workspace浏览器拖到Signal表。点击显示网格▼创建两个并排显示。选择每个显示,然后单击时频按钮添加频谱图视图。

将时间表拖到两个显示器上。

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

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

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

  • 将电源限制设置为-45 db和-20 dB。

重新分配的谱图清楚地显示了三个时频脊。要跟踪脊,请在右侧选择显示。在这方面展示选项卡,单击生成脚本并选择谱图脚本。脚本出现在编辑器中。

%计算频谱图MATLAB(R)9.7生成的%和信号处理工具箱8.2。%生成:26-DEC-2018 17:21:44%的参数时间限制=秒([0 0.002793]);%秒surformlimits = [0 71428.57];%Hz.泄漏= 0.9;timeResolution = 0.00028;%秒重叠的= 85;REARSIGNFLAG = TRUE;%%%指数进入信号时间区域sg_batsignal_roi = sg(:,'击球路');sg_batsignal_roi = sg_batsignal_roi(timeRange(timelimits(1),timelimits(2),'关闭'),1);%计算谱估计%运行下面的函数调用,没有输出参数绘制结果[P F T] = pspectrum (sg_batsignal_ROI,......'谱图'......'surformlimits',频率,......'泄漏',泄漏,......“TimeResolution”,timeresolution,......'重叠的'overlapPercent,......“再分配”,重新开始);

运行脚本。绘制重新分配的频谱图。

网格(秒(t),f,p)xlabel('时间') ylabel (“频率”)轴查看(2)ColorMap粉色的

图包含轴对象。轴对象包含类型表面的对象。

使用tfridge追踪脊的功能。

[冰箱,〜,leridge] = tfridge(p,f,0.01,'numridges',3,'numfrequencybins'10);抓住Plot3(秒(t),冰箱,p(leridge),“:”“线宽”,3)持有离开

图包含轴对象。轴对象包含4个类型表面的物体,线。

感谢伊利诺伊大学贝克曼中心的Curtis Condon、Ken White和Al Feng提供的蝙蝠数据,以及在本例中使用蝙蝠的许可。

也可以看看

应用

职能

相关例子

更多关于