kurtogram

可视化的光谱峰度

折叠所有页面

句法

的kgram = kurtogram（x）的
的kgram = kurtogram（X，sampx）
的kgram = kurtogram（XT）
的kgram = kurtogram（___，水平）
[的kgram，F，W，FC，WC，BW] = kurtogram（___）
kurtogram（___）

描述

例

的kgram= kurtogram（X）返回快kurtogram的kgram信号矢量的X作为基质。kurtogram使用归一化频率（跨越均匀间隔的频率矢量[0π]）来计算时间值。

例

的kgram= kurtogram（X，sampx）返回快kurtogram信号矢量的X在速率或时间间隔采样的sampx作为基质。

的kgram= kurtogram（XT）返回快kurtogram的kgram时间表XT作为基质。

例

的kgram= kurtogram（___，水平）返回快kurtogram使用指定的水平。水平确定窗口的分辨率，以利用水平，因此多少光谱峭度的情况下进行计算。

例

[的kgram，F，w ^，FC，厕所，BW] = kurtogram（___）返回快kurtogram伴随着一组可用于后续的带通滤波器的设计和光谱峰度参数：

  • F- 频率矢量的kgram

  • w ^- 为窗口大小矢量的kgram

  • FC- 频率其中的最大光谱峰度位于

  • 厕所- 窗口大小，其中最大频谱上kurtogram峰度位于

  • BW- 为对最佳带通滤波器带宽的建议

您可以使用此语法与任何以前的语法输入参数。

kurtogram（___）地块kurtogram，与主要的关键优化参数一起，不返回任何数据。您可以使用此语法与任何以前的语法输入参数。

例子

全部收缩

开立真实脚本

计算非平稳信号的kurtogram。对于kurtogram比较不同的级别设置。检查使用标准化频率的kurtogram。使用kurtogram以提供可用于预处理信号，以提高瞬态检测过滤器设置。

生成具有啁啾分量和高斯白噪声的信号。

FS = 1000;t = 0时：1 / FS：10;F1 = 300;F2 = 400;XC =啁啾（T，f1,10，F2）;X = XC + randn（1，长度（T））;

使用采样率画出kurtogramFS。

kurtogram（X，FS）

所述kurtogram显示峰度的结果为一系列窗口的长度和频率。高峰度级别对应于高电平的非平稳或非高斯行为。峰值峭度在顶部的文本提供，与它相关的窗口长度和中心频率一起。带宽是窗口长度的函数。

探索降低最高级别为5的效果。

电平= 5;kurtogram（X，FS，电平）

较低分辨率是显而易见，并且导致较低的峰值的峰度值和一个移位中心频率。

现在积峰度不指定采样率或时间。

kurtogram（x）的

所述kurtogram现在显示与归一化频率。

在图的顶部的参数提供一种用于可用于预过滤所述数据和增强非平稳成分的分化的带通滤波器的建议。你也可以拥有kurtogram返回这些值，这样他们可以输入更直接进入过滤或光谱峰度的功能。

[的kgram，F，W，FC，WC，BW] = kurtogram（X）;厕所
WC = 256
FC
FC = 2.4421
BW
BW = 0.0245

这些值匹配的最佳窗口大小，中心频率，并且所述第一情节的带宽。的kgram是实际kurtogram矩阵，并且F和w ^是伴随它的频率和窗口尺寸的向量。

开立真实脚本

该pkurtosis功能使用默认pspectrum窗口大小（时间分辨率）。您可以指定窗口的大小，而不是使用。在这个例子中，使用功能kurtogram返回一个最佳窗口大小，并使用该结果进行pkurtosis。

创建具有高斯白噪声线性调频信号。

FS = 1000;t = 0时：1 / FS：10;F1 = 300;F2 = 400;X =啁啾（T，f1,10，F2）+ randn（1，长度（T））;

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

pkurtosis（X，FS）标题（“光谱峰度与默认窗口大小”）

现在，使用计算最优窗口大小kurtogram。

kurtogram（X，FS）

所述kurtogram图还示出了300和400赫兹，和节目之间的线性调频脉冲，最佳窗口大小为256饲料W0成pkurtosis。

W0 = 256;pkurtosis（X，FS，W0）标题（“光谱峰度与256的最佳窗口大小”）

主要游览具有较高的峰度值。较高的值提高固定和非固定组件之间的差异，并增强你的提取非平稳成分为特征的能力。

输入参数

全部收缩

时间序列信号，这kurtogram返回快kurtogram，指定为矢量。

采样率或采样时间，指定为以下情况之一：

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

  • 持续时间标量 - 的连续样本之间的时间间隔X

  • 向量，持续时间阵列，或约会时间阵列 - 时间瞬间或持续时间对应于的每个元素X

对于一个示例，请参见计算非平稳信号的Kurtogram。

什么时候sampx表示时间矢量，时间样本可以是不均匀的，与所述pspectrum约束平均时间间隔，平均时间间隔必须遵守：

1 100 < 平均时间间隔 平均时间间隔 < 100。

如果您指定sampx为空，则kurtogram采用归一化频率。换句话说，它假定一个均匀间隔的频率向量跨越[0π]。

信号的时间表从kurtogram返回快kurtogram，指定为时间表包含与单个列的单个可变。XT必须包含的增加，有限的排次。如果时间表缺失或重复的时间点，你可以使用提示修复清洁时间表与缺失，重复，或不均匀时报（MATLAB）。XT可以非均匀采样，用pspectrum约束平均时间间隔，平均时间间隔必须遵守：

1 100 < 平均时间间隔 平均时间间隔 < 100。

最大kurtogram水平，其驱动情况下计算的数，指定为一个正整数。级驱动器的频率窗口大小是kurtogram使用，因此，频率分辨率。由于频率分辨率，上升，时间分辨率下降。如果任分辨率太低的光谱峰度估计会很差。当您指定水平，平衡在时间和频率分辨率的影响。

输出参数

全部收缩

快速kurtogram，与由定义的维度返回作为基质水平。的kgram有2 *水平行3 * 2水平列。该矩阵的每一行表示在频向量的每个元素的谱峰度的结果，以及用于通过该行的位置限定，相对于所述序列中的窗口大小：

[ 0 ， 1 ， 日志 2 （ 3 ） ， 2 ， 1 + 日志 2 （ 3 ） ， 3 ， 3 + 日志 2 （ 3 ） ， ... ， ñ ， ñ + 日志 2 （ 3 ） ， ... ， 升 Ë v Ë 升 ] ，

其中等效窗口大小的电平<Ëm class="varname">ñ是2ñ+1样本。

变频矢量关联的kgram时，返回作为载体。长度F等于列的数量的kgram。

窗口大小的矢量具有相关联的kgram时，返回作为载体。长度F等于列的数量的kgram。

在最大频谱的峰度值的频率的kgram，返回一个标量：

  • 在拉德/秒，如果您还没有指定sampx造成kurtogram使用归一化频率

  • 在赫兹，如果sampx被定义为

您可以使用FC作为一个最佳的带通滤波器最大化经滤波的信号的包络的峰度的中心频率。最大化信封峰度可以让你更轻松地提取形成的冲动成分为特征。

在最大频谱的峰度值的窗口大小的kgram时，返回作为采样的标量。您可以使用厕所以提供最佳的窗口大小pkurtosis。对于一个示例，请参见绘制光谱峰度使用自定义窗口大小。

建议带宽用于最佳带通滤波器，返回作为一个标量：

  • 在拉德/秒，如果您还没有指定sampx造成kurtogram使用归一化频率

  • 在赫兹，如果您指定sampx

您可以使用BW创建一个过滤器最大化经滤波的信号的包络的峰度。BW等于FX/厕所，其中FX是信号采样频率kurtogram源自于sampx。

更多关于

全部收缩

Kurtogram

该kurtogram功能提供关键信息，您可以使用，当你正在执行光谱峰度分析使用pkurtosis。kurtogram计算用于使用快速算法kurtogram多个窗口大小的频谱的峰度。随着kurtogram及其相关联的频率和窗口载体，kurtogram返回最佳窗口大小和其他滤波器调谐参数。它可以直观的计算结果。

快速kurtogram算法使用带通滤波以及简化计算近似为每个窗口的大小和频率，而不是计算的频谱的峰度的短时傅立叶变换（STFT），作为更高保真度pkurtosis确实。它还减少迭代次数的算法要求的频率窗相对于全kurtogram平面跨越[1]。

参考

[1]安东尼，J.，和R. B.兰德尔。“在Kurtogram为瞬时性故障的检测的快速计算。”<Ëm class="citetitle">机械系统和信号处理。卷。20，第1版，2007年，第108-124。

扩展功能

也可以看看

||

介绍了在R2018a

信号处理工具箱文档
金宝app
5G发展与MATLAB

5G发展与MATLAB

下载电子书