用重新分配的光谱图找到并追踪脊线

打开直播脚本

加载一个包含大棕蝠发出的回声定位脉冲的数据文件(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生成。% Generated on: 26-Dec-2018 17:21:44％ 参数timelimits =秒（[0 0.002793]）;%秒frequencyLimits =[0 71428.57]; / /频率限制%赫兹泄漏= 0.9;timeresolution = 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'timeResolution,…“OverlapPercent”，重叠，…'重新分配', reassignFlag);

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

网格(秒(T), F P)包含(“时间”）ylabel（'频率')轴紧的(2) colormap视图粉红色的

使用TFRIDGE.功能，以跟踪脊。

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

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

另请参阅

应用程序

信号分析仪

功能

pspectrum|TFRIDGE.|时间表

信号处理工具箱文档

金宝app

用MATLAB进行深度学习信号处理

下载白皮书