加载包含由大棕色蝙蝠发出的回声机脉冲的数据文件(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)持有离开
感谢伊利诺伊大学贝克曼中心的Curtis Condon、Ken White和Al Feng提供的蝙蝠数据,以及在本例中使用蝙蝠的许可。