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

打开实时脚本

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

负载蝙蝠信号t =(0:长度(batsignal) 1) * DT;sg =时间表(秒(t)”,batsignal);

开放信号分析仪然后将时间表从工作区浏览器拖动到信号表。单击显示网格创建两个并排显示。选择每个显示并单击时频按钮，以添加光谱图视图。

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

选择光谱图选项卡。在右边的显示器上，检查重新分配．为每个显示:

将时间分辨率设置为280微秒，并指定相邻线段之间85%的重叠。
使用泄漏滑块增加泄漏，直到RBW约为4.5 kHz。
功率限制设置为- 45db和- 20db。

重新分配的谱图清楚地显示了三个时频脊。要跟踪山脊，请选择右边的显示器。在显示选项卡上,单击生成脚本并选择光谱图脚本．脚本将出现在编辑器中。

%计算光谱图%由MATLAB(R) 9.7和信号处理工具箱8.2生成。%生成时间:2018年12月26日17:21:44%参数时间限制=秒（[0.002793]）；%秒frequencyLimits = [0 71428.57];%赫兹泄漏=0.9；时间分辨率=0.00028；%秒overlapPercent = 85;reassignFlag = true;％％索引为感兴趣的信号时间区域sg_batsignal_ROI = sg (:,“batsignal”）;sg_batsignal_ROI = sg_batsignal_ROI (timerange(时限(1),期限(2),“关闭”), 1);计算谱估计%运行下面不带输出参数的函数调用来绘制结果[P，F，T]=pspectrum（sg_-batsignal_-ROI，．..的谱图，．..“FrequencyLimits”frequencyLimits,．..“漏”泄漏,．..“时间分辨率”timeResolution,．..“OverlapPercent”百分之，．..“重新分配”, reassignFlag);

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

网格(秒(T), F P)包含(“时间”)伊拉贝尔(“频率”)轴牢固的(2) colormap视图粉红色的

图中包含一个轴对象。axis对象包含一个类型为surface的对象。

使用特弗里奇追踪山脊的功能。

(冰箱,~,lridge) = tfridge (P, F, 0.01,“NumRidges”3,“NumFrequencyBins”10);持有在…上plot3(秒(T),冰箱,P (lridge),':'，“线宽”, 3)从

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

感谢Curtis Condon，Ken White和Al Feng的伊利诺伊大学贝克曼中心的蝙蝠数据，并允许在这个例子中使用它。

另请参阅

应用程序

信号分析仪

功能

pspectrum|特弗里奇|时间表

信号处理工具箱文档

金宝app

基于MATLAB的信号处理深度学习

下载白皮书