kurtogram

可视化谱峰度

折叠所有页面

语法

kgram = kurtogram (x)

kgram = kurtogram (x, sampx)

kgram = kurtogram (xt)

kgram = kurtogram (＿＿＿水平)

[kgram, f, w, fc, wc, bw) = kurtogram (＿＿＿）

kurtogram (＿＿＿）

描述

例子

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

例子

kgram= kurtogram (x，sampx）快速返回kurtogram信号矢量x按速率或时间间隔采样sampx，作为矩阵。

kgram= kurtogram (xt）快速返回kurtogramkgram的时间表xt作为一个矩阵。

例子

kgram= kurtogram (＿＿＿，水平）快速返回kurtogram使用指定的水平。水平确定要使用的窗口分辨率的水平，因此需要计算多少谱峰情况。

例子

［kgram，f，w，足球俱乐部，wc，bw) = kurtogram (＿＿＿）快速返回kurtogram连同一组参数，您可以用于后续的带通滤波器设计和光谱峰度:

f-频率矢量kgram
w窗口大小矢量kgram
足球俱乐部-谱峰最大位置的频率
wc-峰度图上最大谱峰的窗口大小
bw—建议的最佳带通滤波器带宽

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

kurtogram (＿＿＿）绘制峰度图，以及关键的关键优化参数，不返回任何数据。您可以将此语法与前面语法中的任何输入参数一起使用。

例子

全部折叠

计算非平稳信号的峰度图

打开生活的脚本

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

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

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

使用采样率绘制峰度图fs。

kurtogram (x, fs)

$图中包含一个坐标轴。在第7层，标题为K_{max} = 8.5267、最优窗长= 256、中心频率= 388.6719 Hz、带宽= 3.9062 Hz的轴包含一个类型为image的对象。$

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

探索将最高等级降至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的对象。$

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

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

[kgram, f, w, fc, wc, bw) = kurtogram (x);wc

wc = 256

足球俱乐部

fc = 2.4421

bw

bw = 0.0245

这些值与第一个图的最佳窗口大小、中心频率和带宽相匹配。kgram是实际的峰度矩阵吗f和w是伴随它的频率和窗口大小向量。

使用定制的窗口大小绘制光谱峰度

打开生活的脚本

的pkurtosis函数使用默认值pspectrum窗口大小(时间分辨率)。您可以指定要使用的窗口大小。在本例中，使用该函数kurtogram返回最佳窗口大小并将该结果用于pkurtosis。

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

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

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

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

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

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

kurtogram (x, fs)

$图中包含一个坐标轴。在第7层，标题为K_{max} = 8.5267、最优窗长= 256、中心频率= 388.6719 Hz、带宽= 3.9062 Hz的轴包含一个类型为image的对象。$

峰度图还说明了在300和400 Hz之间的啁啾，并表明最佳窗口大小为256。饲料w0成pkurtosis。

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

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

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

输入参数

全部折叠

`x`- - - - - -时间序列信号
向量

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

`sampx`- - - - - -信号的采样率或采样时间
归一化频率(默认)|积极的数字标量|`持续时间`标量|数字向量(秒)|`持续时间`数组|`datetime`数组

采样率或采样时间，指定为下列之一:

正数值标量-频率单位为赫兹
持续时间的连续样本之间的时间间隔X
向量,持续时间数组,或datetime的每个元素对应的时间瞬间或持续时间x

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

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

$\frac{1}{One hundred.} < \frac{平均时间间隔}{平均时间间隔} < One hundred.$

如果您指定sampx是空的,那么kurtogram使用归一化频率。换句话说，它假设一个均匀间隔的频率向量生成[0 π]。

`xt`- - - - - -信号的时间表
时间表

信号时间表kurtogram返回指定为的快速峰度图时间表它包含单个变量和单个列。xt必须包含递增的有限行时间。如果时间表有缺失或重复的时间点，你可以使用提示来修正它清除时间表缺失，重复，或不一致的时间。xt可以非均匀采样，用pspectrum中值时间间隔和平均时间间隔必须服从的约束:

$\frac{1}{One hundred.} < \frac{平均时间间隔}{平均时间间隔} < One hundred.$

`水平`- - - - - -最大kurtogram水平
正整数

最大峰度图级别，用于驱动要计算的案例数，指定为正整数。电平驱动的频率窗口大小kurtogram使用，因此，频率分辨率。频率分辨率上升，时间分辨率下降。如果任何一个分辨率过低，谱峰度估计都会很差。当你指定水平，平衡对时间和频率分辨率的影响。

输出参数

全部折叠

`kgram`——快kurtogram
矩阵

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

$［ 0 ， 1 ， {日志}_{2} （ 3. ）， 2 ， 1 + {日志}_{2} （ 3. ）， 3. ， 3. + {日志}_{2} （ 3. ）， … ， n ， n + {日志}_{2} （ 3. ）， … ， l e v e l ］，$

一个关卡的等效窗口大小在哪里n是2^n＋１样本。

`f`——频率向量
向量

频率矢量kgram，返回为向量。的长度f等于kgram。

`w`- - - - - -窗口大小向量
向量

关联的窗口大小向量kgram，返回为向量。的长度f等于kgram。

`足球俱乐部`-最大谱峰度值的频率
标量

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

如果没有指定，单位是rad/秒sampx,导致kurtogram使用归一化频率
在赫兹,如果sampx定义

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

`wc`-最大谱峰度值的窗口大小
标量

中最大谱峰度值的窗口大小kgram，在示例中作为标量返回。您可以使用wc提供的最佳窗口大小pkurtosis。例如，请参见使用定制的窗口大小绘制光谱峰度。

`bw`-最佳带通滤波器的建议带宽
标量

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

如果没有指定，单位是rad/秒sampx,导致kurtogram使用归一化频率
赫兹，如果你已经指定sampx

您可以使用bw创建一个滤波器，使滤波信号的包络峰度最大化。bw等于外汇/wc,在那里外汇信号采样频率是这个吗kurtogram来源于sampx。

参考文献

安东尼，J, R. B.兰德尔。暂态故障检测的峰度图快速计算机械系统和信号处理。第20卷，第1期，2007,108-124页。

扩展功能

C / c++代码生成
使用MATLAB®Coder™生成C和c++代码。

使用注意事项及限制:

输入必须是双精度的。
代码生成不支持时间表。万博1manbetx

另请参阅

pentropy|pkurtosis|pspectrum

介绍了R2018a

kurtogram

语法

描述

例子

计算非平稳信号的峰度图

使用定制的窗口大小绘制光谱峰度

输入参数

`x`- - - - - -时间序列信号
向量

`sampx`- - - - - -信号的采样率或采样时间
归一化频率(默认)|积极的数字标量|`持续时间`标量|数字向量(秒)|`持续时间`数组|`datetime`数组

`xt`- - - - - -信号的时间表
时间表

`水平`- - - - - -最大kurtogram水平
正整数

输出参数

`kgram`——快kurtogram
矩阵

`f`——频率向量
向量

`w`- - - - - -窗口大小向量
向量

`足球俱乐部`-最大谱峰度值的频率
标量

`wc`-最大谱峰度值的窗口大小
标量

`bw`-最佳带通滤波器的建议带宽
标量

更多关于

Kurtogram

参考文献

扩展功能

C / c++代码生成
使用MATLAB®Coder™生成C和c++代码。

另请参阅

信号处理工具箱文档

万博1manbetx

基于MATLAB的信号处理深度学习

kurtogram

语法

描述

例子

计算非平稳信号的峰度图

使用定制的窗口大小绘制光谱峰度

输入参数

x- - - - - -时间序列信号向量

sampx- - - - - -信号的采样率或采样时间归一化频率(默认)|积极的数字标量|持续时间标量|数字向量(秒)|持续时间数组|datetime数组

xt- - - - - -信号的时间表时间表

水平- - - - - -最大kurtogram水平正整数

输出参数

kgram——快kurtogram矩阵

f——频率向量向量

w- - - - - -窗口大小向量向量

足球俱乐部-最大谱峰度值的频率标量

wc-最大谱峰度值的窗口大小标量

bw-最佳带通滤波器的建议带宽标量

更多关于

Kurtogram

参考文献

扩展功能

C / c++代码生成使用MATLAB®Coder™生成C和c++代码。

另请参阅

信号处理工具箱文档

万博1manbetx

基于MATLAB的信号处理深度学习

`x`- - - - - -时间序列信号
向量

`sampx`- - - - - -信号的采样率或采样时间
归一化频率(默认)|积极的数字标量|`持续时间`标量|数字向量(秒)|`持续时间`数组|`datetime`数组

`xt`- - - - - -信号的时间表
时间表

`水平`- - - - - -最大kurtogram水平
正整数

`kgram`——快kurtogram
矩阵

`f`——频率向量
向量

`w`- - - - - -窗口大小向量
向量

`足球俱乐部`-最大谱峰度值的频率
标量

`wc`-最大谱峰度值的窗口大小
标量

`bw`-最佳带通滤波器的建议带宽
标量

C / c++代码生成
使用MATLAB®Coder™生成C和c++代码。