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

载入一个包含由一只大棕色蝙蝠(Eptesicus fuscus)，并以7微秒的采样率进行测量。使用信号和时间信息创建MATLAB®时间表。

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

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

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

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

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

对于右边的显示，选中重新分配．

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

%计算光谱图%由MATLAB(R) 9.13和信号处理工具箱9.1生成。%生成于:2022年6月15日12:02:38%的参数timeLimits = seconds([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”overlapPercent,.．.“再分配”, reassignFlag);

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

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

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

使用tfridge函数跟踪山脊。

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

图中包含一个axes对象。axis对象包含4个类型为surface、line的对象。

感谢伊利诺伊大学贝克曼中心的Curtis Condon、Ken White和Al Feng提供的蝙蝠数据，并允许我们在本例中使用这些数据。

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

另请参阅

应用程序

功能

相关的例子

更多关于