加载包含由大棕蝠发出的回声定位脉冲数据文件(大棕鲇)中并用7微秒的采样间隔进行测量。创建使用信号和时间信息的MATLAB®时间表。
加载batsignalT =(0:长度(batsignal)-1)* DT;SG =时间表(秒(t)的”,batsignal);
打开信号分析仪并将时间表从工作区浏览器拖动到信号表。点击显示网格▼创建两个并排侧显示器。选择每个展示和点击时频按钮添加一个频谱图。
拖动时间表两个显示器。
选择光谱图标签。在右边的显示,请检查重新分配。对于每个显示:
将时间分辨率设置为280微秒,并指定相邻段之间85%的重叠。
使用泄漏滑块增加泄漏,直到RBW约为4.5 kHz。
将电源限制-45 dB和-20分贝。
重新分配的频谱图清楚地显示出三个时间频脊。为了跟踪脊,选择右边的显示。在显示选项卡,单击生成脚本并选择谱图脚本
。该脚本编辑器中显示。
%计算频谱%的MATLAB(R)9.7和信号处理工具箱8.2生成。26日 - 12月2018 17时21分44秒:%上生成%的参数时间限制=秒([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(timeLimits(1),期限(2),'关闭'),1);%计算谱估计%运行下面的函数调用,而不输出参数绘制结果[P F T] = pspectrum (sg_batsignal_ROI,...“频谱”,...'FrequencyLimits',frequencyLimits,...'泄漏',泄漏,...“TimeResolution”,timeResolution,...'OverlapPercent'overlapPercent,...“再分配”,reassignFlag);
运行该脚本。剧情重新分配的频谱。
目(秒(T),F,P)xlabel('时间')ylabel (“频率”)轴紧视图(2)的颜色表粉
使用tfridge
功能跟踪脊。
[冰箱,〜,lridge] = tfridge(P,F,0.01,'NumRidges',3,'NumFrequencyBins',10);保持在plot3(秒(T),冰箱,P(lridge),“:”,“线宽”3)保持离
感谢伊利诺斯大学贝克曼中心的柯蒂斯·康顿、肯·怀特和艾尔·冯提供蝙蝠数据,并允许在本例中使用。