主要内容

黑线图

可视化光谱峰度

描述

实例

kgram=千字图(x)快速返回黑线图kgram信号矢量x作为一个矩阵。黑线图使用归一化频率(跨越[0 π]的均匀间隔频率向量)来计算时间值。

实例

kgram=千字图(x,sampx)快速返回黑线图信号矢量x按速率或时间间隔采样sampx,作为矩阵。

kgram=千字图(xt)快速返回黑线图kgram的时间表xt作为一个矩阵。

实例

kgram=千字图(___,水平)快速返回黑线图使用指定的水平水平确定要使用的窗口分辨率的水平,因此需要计算多少谱峰情况。

实例

[kgram,F,W,足球俱乐部,厕所,体重) = kurtogram (___)快速返回黑线图连同一组参数,您可以用于后续的带通滤波器设计和光谱峰度:

  • F-频率矢量kgram

  • W窗口大小矢量kgram

  • 足球俱乐部-最大光谱峰度所在的频率

  • 厕所-峰度图上最大谱峰的窗口大小

  • 体重-最佳带通滤波器的建议带宽

您可以将此语法与前面语法中的任何输入参数一起使用。

黑线图(___)绘制kurtogram以及关键的优化参数,而不返回任何数据。您可以将此语法与以前语法中的任何输入参数一起使用。

例子

全部折叠

计算非平稳信号的峰度图。比较峰度图的不同电平设置。检查使用归一化频率的峰度图。使用峰度图提供滤波器设置,可用于对信号进行预处理,以增强瞬态检测。

产生一个带有啁啾分量和高斯白噪声的信号。

fs = 1000;t = 0:1 / fs: 10;f1 = 300;f2 = 400;xc =唧唧声(10 t f1, f2);x = xc randn(1、长度(t));

使用采样率绘制kurtogram财政司司长

kurtogram(x,fs)

图中包含一个轴。标题为K_{max}=8.5267(级别7)、最佳窗口长度=256、中心频率=388.6719 Hz、带宽=3.9062 Hz的轴包含图像类型的对象。

峰度图显示了窗口长度和频率范围内的峰度结果。高峰度水平对应于高水平的非平稳或非高斯行为。顶部文本中提供了峰值峰度,以及与其相关的窗口长度和中心频率。带宽是窗口长度的函数.

探索将最高等级降至5级的效果。

水平= 5;kurtogram (x, fs,级别)

图中包含一个坐标轴。在第5级,标题为K_{max} = 1.8848、最优窗长= 64、中心频率= 304.6875 Hz、带宽= 15.625 Hz的轴包含一个类型为image的对象。

较低的分辨率是明显的,导致较低的峰度值和偏移的中心频率。

现在在不指定采样率或时间的情况下绘制峰度。

kurtogram(x)

图中包含一个坐标轴。在第7层,标题为K_{max} = 8.5267,最优窗长= 256,中心频率= 2.4421 rad/sample,带宽= 0.024544 rad/sample的坐标轴包含一个类型为image的对象。

峰度图现在显示为归一化频率。

图顶部的参数为带通滤波器提供了建议,该滤波器可用于预滤数据并增强非平稳分量的微分。你也可以黑线图返回这些值,以便更直接地输入到滤波或光谱峰度函数中。

[kgram, f, w, fc, wc, bw) = kurtogram (x);厕所
wc = 256
足球俱乐部
fc = 2.4421
体重
bw = 0.0245

这些值与第一个绘图的最佳窗口大小、中心频率和带宽相匹配。kgram是实际的kurtogram矩阵,并且FW是伴随它的频率和窗口大小向量。

这个肺结核函数使用默认值pspectrum窗口大小(时间分辨率)。您可以指定要使用的窗口大小。在本例中,使用函数黑线图返回最佳窗口大小并将该结果用于肺结核

创建带有高斯白噪声的啁啾信号。

fs=1000;t=0:1/fs:10;f1=300;f2=400;x=chirp(t,f1,10,f2)+randn(1,长度(t));

用默认窗口大小绘制光谱峰度。

pkurtosis (x, fs)标题(“具有默认窗口大小的光谱峰度”)

图中包含一个坐标轴。标题为“光谱峰度与默认窗口大小”的轴包含2个类型为line的对象。这些对象代表光谱峰度,0.95置信区间。

现在计算最佳窗口大小使用黑线图

kurtogram(x,fs)

图中包含一个轴。标题为K_{max}=8.5267(级别7)、最佳窗口长度=256、中心频率=388.6719 Hz、带宽=3.9062 Hz的轴包含图像类型的对象。

kurtogram图还显示了300和400 Hz之间的啁啾,并显示最佳窗口大小为256.Feedw0进入肺结核

w0=256;肺结核(x,fs,w0)标题(最佳窗口尺寸为256的谱峰度)

图中包含一个坐标轴。标题为“光谱峰度与最佳窗口大小256”的轴包含2个类型为line的对象。这些对象代表光谱峰度,0.95置信区间。

主偏移具有较高的峰度值。较高的值改善了平稳分量和非平稳分量之间的区别,并增强了将非平稳分量作为特征提取的能力。

输入参数

全部折叠

时间序列信号黑线图返回指定为矢量的快速峰度图。

采样率或采样时间,指定为以下之一:

  • 正数值标量-频率单位为赫兹

  • 期间的连续样本之间的时间间隔X

  • 矢量,期间数组,或datetime的每个元素对应的时间瞬间或持续时间x

例如,请参见计算非平稳信号的峰度图

sampx表示一个时间向量,时间样本可以是非均匀的,用pspectrum中值时间间隔和平均时间间隔必须服从的约束:

1. One hundred. < 中间时间间隔 平均时间间隔 < 100

如果您指定sampx是空的,那么黑线图使用标准化频率。换句话说,它假设一个均匀分布的频率向量跨越[0π]。

信号时间表黑线图返回快速kurtogram,指定为时间表它包含单个变量和单个列。xt必须包含递增的有限行时间。如果时间表缺少或重复时间点,可以使用中的提示进行修复缺少、重复或不均匀时间的干净时间表xt可以非均匀采样,用pspectrum中值时间间隔和平均时间间隔必须服从的约束:

1. One hundred. < 中间时间间隔 平均时间间隔 < 100

最大峰度图级别,用于驱动要计算的案例数,指定为正整数。电平驱动的频率窗口大小黑线图因此,使用频率分辨率。随着频率分辨率的提高,时间分辨率降低。如果任一分辨率过低,则光谱峰度估计值将很差。当您指定水平,平衡对时间和频率分辨率的影响。

输出参数

全部折叠

快速峰度图,以矩阵形式返回,其维数由水平kgram有两个*水平行和3*2水平列。矩阵的每一行表示频率向量中每个元素的谱峰度结果,以及由行位置定义的窗口大小相对于序列:

[ 0 , 1. , 日志 2. ( 3. ) , 2. , 1. + 日志 2. ( 3. ) , 3. , 3. + 日志 2. ( 3. ) , ... , N , N + 日志 2. ( 3. ) , ... , L E v E L ] ,

一个关卡的等效窗口大小在哪里N是2N+1样品。

与之相关的频率向量kgram,作为向量返回。的长度F等于中的列数kgram

关联的窗口大小向量kgram,作为向量返回。的长度F等于中的列数kgram

最大光谱峰度值的频率kgram,作为标量返回:

  • 如果没有指定,单位是rad/秒sampx,导致黑线图使用归一化频率

  • 在赫兹,如果sampx定义

您可以使用足球俱乐部作为最优带通滤波器的中心频率,使滤波信号的包络峰度最大化。最大化包络峰度可以让你更容易地提取结果的脉冲分量作为一个特征。

最大光谱峰度值的窗口大小kgram,在示例中作为标量返回。您可以使用厕所提供的最佳窗口大小肺结核.例如,请参见使用自定义窗口大小绘制光谱峰度

最优带通滤波器的建议带宽,返回为标量:

  • 如果没有指定,单位是rad/秒sampx,导致黑线图使用归一化频率

  • 赫兹,如果你已经指定sampx

您可以使用体重创建最大化滤波信号包络峰度的滤波器。体重等于外汇/厕所,在那里外汇信号采样频率是这个吗黑线图源于sampx

更多关于

全部折叠

Kurtogram

这个黑线图函数提供执行时可以使用的关键信息谱峰态分析使用肺结核黑线图使用快速的峰度图算法计算多个窗口大小的谱峰度。随着峰度图及其相关的频率和窗口向量,黑线图返回最佳窗口大小和其他筛选器优化参数。它可以将计算结果可视化。

快速峰度图算法使用带通滤波和简化计算来近似每个窗口大小和频率的谱峰度,而不是计算短时傅里叶变换(STFT)作为高保真度肺结核所做的事。它还减少了算法需要跨越频率窗平面相对于完整的峰度图的迭代次数[1]

参考文献

[1] Antoni,J.和R.B.Randall,“用于检测瞬态故障的Kurtogram快速计算。”机械系统与信号处理.第20卷,第1期,2007,108-124页。

扩展能力

R2018a中引入