主要内容

使用重新分配谱图发现和跟踪脊

打开生活的脚本

加载数据文件包含一个蝙蝠回声定位脉冲发出的(Eptesicus fuscus)和测量的采样率7微秒。创建一个MATLAB®时间表使用信号和时间信息。

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

开放信号分析仪并拖动的时间表工作空间浏览器表的信号。点击显示网格创建两个并排显示。选择每个显示,显示选项卡上,单击时频添加一个谱图视图。

把时间表显示。

选择光谱图选项卡。为每个显示:

  • 将权限设置为-45分贝和-20分贝。

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

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

显示在右边,检查重新分配。

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

%计算光谱图%由MATLAB 9.13 (R)和信号处理工具箱9.1。%上生成:15 - 2022年6月——12:02:38%的参数期限=秒(0.002793 [0]);%秒frequencyLimits = 71428.57 [0];%赫兹泄漏= 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”overlapPercent,…“再分配”,reassignFlag);

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

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

图包含一个坐标轴对象。坐标轴对象包含时间,ylabel频率包含一个类型的对象的表面。

使用tfridge函数跟踪脊。

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

图包含一个坐标轴对象。坐标轴对象包含时间,ylabel频率包含4对象类型的表面。

由于柯蒂斯Condon肯白、阿尔贝克曼中心的冯伊利诺伊大学的蝙蝠数据和权限在这个例子中使用它。

另请参阅

应用程序

功能

相关的例子

更多关于