主要内容

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

加载包含由大棕色蝙蝠发出的回声机脉冲的数据文件(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)持有离开

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

感谢Curtis Condon,Ken White和Beckman中心的Al Feng在伊利诺伊大学的BAT数据和许可在此示例中使用它。

也可以看看

应用

职能

相关例子

更多关于