主要内容

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

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

负载batsignalt=(0:长度(信号)-1)*DT;sg=时刻表(秒(t)”,信号灯;

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

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

选择光谱图标签。在右侧显示屏上,选中重新分配. 对于每个显示器:

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

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

  • 将功率限制设置为–45 dB和–20 dB。

重新分配的频谱图清楚地显示了三个时频脊。要跟踪脊线,请选择右侧的显示。上陈列选项卡,单击生成脚本选择谱图脚本. 脚本将显示在编辑器中。

%计算谱图%由MATLAB(R)9.7和信号处理工具箱8.2生成。%生成日期:2018年12月26日17:21:44%的参数时间限制=秒([0 0.002793]);%秒频率限制=[0 71428.57];%赫兹泄漏= 0.9;timeResolution = 0.00028;%秒百分比=85;reassignFlag=true;%%%索引到感兴趣的信号时间区域sg_batsignal_ROI=sg(:,“蝙蝠信号”); sg_batsignal_ROI=sg_batsignal_ROI(时间范围(时间限制)(1),时间限制(2),“关闭”), 1);%计算谱估计%在不带输出参数的情况下运行下面的函数调用以绘制结果[P F T] = pspectrum (sg_batsignal_ROI,...“光谱图”,...“频率限制”,频率限制,...“漏”泄漏量...“TimeResolution”,时间分辨率,...“重叠百分比”overlapPercent,...“再分配”,重新分配旗帜);

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

网格(秒(T),F,P)xlabel(“时间”) ylabel (“频率”)轴心视图(2)彩色地图粉红色

图中包含一个轴对象。“轴”对象包含“曲面”类型的对象。

使用tfridge函数跟踪脊线。

[fridge,~,lridge]=t桥(P,F,0.01,“NumRidges”,3,“NumFrequencyBins”,10); 持有图3(秒(T)、冰箱、P(lridge),“:”,“线宽”(3)保持

图中包含一个轴对象。轴对象包含面、线四个对象。

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

另见

应用程序

功能

相关实例

更多关于