加载包含由大棕色蝙蝠发出的回声机脉冲的数据文件(Eptesicus fuscus.)并以7微秒的采样间隔测量。使用信号和时间信息创建MATLAB®时间表。
加载击球路t =(0:长度(击球路)-1)* dt;sg =时间表(秒(t)',batsignal);
打开信号分析仪并将时间表从工作区浏览器拖动到信号表。点击显示网格▼创建两个并排显示。选择每个显示,然后单击时频按钮添加频谱图视图。
将时间表拖到两个显示器上。
选择谱图标签。在右侧显示屏上,检查重新分配。对于每个显示:
将时间分辨率设置为280微秒,并在邻接段之间指定85%重叠。
使用泄漏滑块增加泄漏直到RBW约为4.5 kHz。
将电源限制设置为-45 db和-20 dB。
重新分配的频谱图清楚地显示了三个时频脊。要跟踪脊,请在右侧选择显示。在这一点展示选项卡,单击生成脚本并选择谱图脚本
。脚本出现在编辑器中。
%计算频谱图MATLAB(R)9.7生成的%和信号处理工具箱8.2。%生成:26-of-2018 17:21:44% 参数timelimits =秒([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,......'重叠的',重叠,......'重新分配',重新译名;
运行脚本。绘制重新分配的频谱图。
网格(秒(t),f,p)xlabel('时间')ylabel('频率')轴紧的查看(2)ColorMap粉色的
使用TFRIDGE.
跟踪脊的功能。
[冰箱,〜,leridge] = tfridge(p,f,0.01,'numridges',3,'numfrequencybins',10);抓住在Plot3(秒(T),冰箱,P(液晶),':'那'行宽',3)持有离开
感谢Curtis Condon,Ken White和Beckman中心的Al Feng在伊利诺伊大学的BAT数据和许可在此示例中使用它。