使用重新分配的光谱图寻找和跟踪脊
载入一个包含由一只大棕蝠发出的回声定位脉冲的数据文件(Eptesicus fuscus),并以7微秒的采样率进行测量。使用信号和时间信息创建MATLAB®时间表。
负载batsignalt =(0:长度(蝙蝠信号)-1)*DT;Sg =时间表(秒(t)',蝙蝠信号);
开放信号分析仪并将时间表从工作空间浏览器信号表。点击显示网格创建两个并排的显示。选择每个显示,并在显示选项卡上,单击时频添加光谱图视图。
将时间表拖到两个显示器上。
选择光谱图选项卡。对于每个显示:
将功率限制设置为- 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;timerresolution = 0.00028;%秒overlapPercent = 85;reassignFlag = true;% %索引到感兴趣的信号时间区域sg_batsignal_ROI = sg(:,“batsignal”);sg_batsignal_ROI = sg_batsignal_ROI(timerange(timeLimits(1),timeLimits(2),“关闭”), 1);计算光谱估计运行下面不带输出参数的函数调用来绘制结果[P,F,T] = pspectrum(sg_batsignal_ROI;...的谱图,...“FrequencyLimits”frequencyLimits,...“漏”泄漏,...“TimeResolution”timeResolution,...“OverlapPercent”overlapPercent,...“再分配”, reassignFlag);
运行脚本。绘制重新分配的频谱图。
网格(秒(T), F P)包含(“时间”) ylabel (“频率”)轴紧(2) colormap视图粉红色的
使用tfridge
函数跟踪脊线。
[冰箱,~,lridge] = tfridge(P,F,0.01, num脊=3,NumFrequencyBins=10);持有在plot3(秒(T),冰箱,P (lridge),":"线宽= 3)从
感谢伊利诺伊大学贝克曼中心的Curtis Condon、Ken White和Al Feng提供的蝙蝠数据,并允许在这个例子中使用它。
另请参阅
应用程序
功能
相关的例子
- 找出相关信号之间的延迟
- 通过改变窗户泄漏来解决音调问题
- 使用持久谱发现干扰
- 利用复包络的调制与解调
- 从音乐信号中提取声音
- 非均匀采样信号的重采样与滤波
- 使用自己的函数衰减饱和信号
- 计算振动信号的包络谱
- 从鲸鱼的歌声中提取感兴趣的区域