加载一个数据文件,其中包含一个大棕蝙蝠发出的回声定位脉冲(褐斑伊蚊)并以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提供的蝙蝠数据,以及在本例中使用蝙蝠的许可。