使用重新分配的谱图找到和跟踪山脊
加载一个数据文件，其中包含一只棕色大蝙蝠发出的回声定位脉冲(Eptesicus fuscus)，采样间隔为7微秒。使用信号和时间信息创建一个MATLAB®时间表。
负载batsignalt =(0:长度(batsignal) 1) * DT;sg =时间表(秒(t)”,batsignal);
开放信号分析仪并将时间表从Workspace浏览器拖到Signal表中。点击显示网格▼创建两个并排显示。选择每个显示并单击时频按钮添加声谱图视图。
将时间表拖到两个显示器上。
选择光谱图选项卡。在右边的显示器上，检查一下重新分配。为每个显示:
设置时间分辨率为280微秒，并指定相邻段重叠85%。
使用泄漏滑块增加泄漏，直到RBW约为4.5 kHz。
设置为- 45db和- 20db。
重新分配的谱图清楚地显示了三个时频脊。要跟踪脊线，请选择右边的显示。在显示选项卡上,单击生成脚本并选择光谱图脚本。脚本出现在编辑器中。
%计算光谱图%由MATLAB(R) 9.7和信号处理工具箱8.2生成。% Generated on: 26-Dec-2018 17:21:44%的参数timeLimits = seconds([0 0.002793]);%秒frequencyLimits = [0 71430 .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视图粉红色的
使用tfridge功能来跟踪山脊。
(冰箱,~,lridge) = tfridge (P, F, 0.01,“NumRidges”3,“NumFrequencyBins”10);持有在plot3(秒(T),冰箱,P (lridge),“:”,“线宽”, 3)从
感谢伊利诺伊大学贝克曼中心的柯蒂斯·康顿、肯·怀特和阿尔·冯提供的蝙蝠数据并允许在这个例子中使用。
另请参阅
应用程序
功能
你也可以从以下列表中选择一个网站:
