MATLAB R2018a中信号处理相关的新函数 : pkurtosis与kurtogram

最新推荐文章于 2025-10-19 13:46:24 发布
转载 最新推荐文章于 2025-10-19 13:46:24 发布 · 3.2k 阅读 · 14 ·
CC 4.0 BY-SA版权
原文链接：https://zhuanlan.zhihu.com/p/34582059

本文介绍了一种基于频域峰度的信号分析方法，通过pkurtosis和kurtogram函数，能够有效检测信号的非平稳性和非白噪声特性。特别地，文章详细解释了如何利用这些工具确定最佳窗口长度，以实现频域峰度的最大化。

pkurtosis: Spectral kurtosis from signal or spectrogram

用来判断是否为平稳的(stationary)白噪信号, 如果是的话, 取值在各个频率上都是较低的(在置信区间内).

如果某些频率上在置信区间外面的话, 说明在对应的频率下, 存在非平稳信号或者是非白噪信号.

基于STFT的结果, 进一步计算频域的峰度.

 

 

一个300Hz到400Hz的非平稳信号 + 白噪声做实验.

 

很清晰看到, 检测到300Hz到400Hz是非平稳或非白噪

 

kurtogram: Visualize spectral kurtosis

相关函数为pkurtosis.

 

pkurtosis使用了一个窗口长度来计算STFT(进一步计算频域峰度) , 而kurtogram使用了很多个窗口长度来计算STFT(进一步计算频域峰度)

 

kurtogram可以用来观察不同窗口长度下, 频域峰度的变化, 还可以用来确定最佳窗口长度.

最佳指的是频域峰度最大化.

 

 

使用最佳的窗口长度来计算pkurtosis

 

可以看出, 最佳窗口长度下, 频域峰度最大到了10, 而使用默认的窗口长度下, 频域峰度最大才到6. 说明默认的窗口长度不是最佳窗口长度.

jacksong2021
关注
R语言使用moments包计算偏度（Skewness）和峰度（Kurtosis）实战：计算偏度（Skewness）和峰度（Kurtosis）、确定样本数据是否具有与正态分布匹配的偏度和峰度（假设检验）
data+scenario+science+insight
01-04 9400
R语言使用moments包计算偏度（Skewness）和峰度（Kurtosis）实战：计算偏度（Skewness）和峰度（Kurtosis）、确定样本数据是否具有与正态分布匹配的偏度和峰度（假设检验） 目录 R语言使用moments包计算偏度（Skewness）和峰度（Kurtosis）实战：计算偏度（Skewness）和峰度（Kurtosis）、确定样本数据是否具有与正态分布匹配的偏度和峰度（假设检验） #偏度和峰度 #仿真数据及可视化 #计算偏度（Skewness）和峰度（Kurtosis
轴承故障诊断之快速谱峭度法 Matlab 代码说明
cocofisher的博客
06-17 1万+
文章目录文章介绍代码说明1、主程序2、Fast_Kurtogram（快速谱峭度主函数）3、kurt（计算峭度）4、max_IJ（搜索最大峭度对应的信号子带）5、K_wpQ（计算峭度矩阵）6、DBFB（信号通过1/2树状滤波器）7、TBFB（信号通过1/3树状滤波器）8、K_wpQ_local（计算每层的峭度）level=3----xlevel=2----a_level3level=1----a_level2_from_a_level3level=1----d_level2_from_a_level3leve
matlab求峰度(kurtosis）、偏度(skewness)得函数及python对应的峰度(kurtosis）、偏度(skewness)函数
热门推荐
雾隐雾现的至渝博客
11-07 4万+
                    正在做一个把matlab程序转python的工作，记录下遇到的问题与解决方案           定义 峰度kurtosis：用于度量x偏离某分布的程度。 正太分布的峰度是3； 当时间序列的曲线峰值比正太分布的高时，峰度大于3； 当比正太分布的低时，峰度小于3。      
小波包matlab代码-Subband-averaging-kurtogram-with-dual-tree-complex-wavelet-
05-24
小波包matlab代码双树复数小波包的子带平均峰度图 该存储库提供了带有双树复数小波包变换的纸带平均峰度图的代码，用于旋转机械故障诊断 Matlab SAK的实现。 本文提出了一种称为子带平均峰度图（SAK）的方法，并结合了双树复数小波包变换（DTCWPT），以提高快速峰度图（FK）在旋转机械故障诊断中的性能。 所提出的方法首先通过滑动窗口将信号分割成M个子信号，然后计算通过每个子信号的DTCWPT获得的子带的峰度。 最后，计算相应子带的平均峰度，以获得SAK，该SAK表示包络分析的最佳频带。 FK很容易被非高斯噪声（例如，偶发的脉冲干扰）误导，而SAK可以克服此问题。 此外，DTCWPT同时细分了高频和低频频段，提供了所需的近似位移不变性，同时保持了较低的计算量。 当原始DTCWPT在高通通道上迭代滤波器组时，所获得的信号子带不会以中心频率的单调顺序排列。 可以通过根据带通特性交换反向滤波器组来解决此问题。 所提出的方法与FK相比提供了改进的性能，特别是用于从包含各种干扰的噪声信号中提取周期性瞬变。
线性偏度与峰度分析实战：非高斯信号检测与数据分布特征提取
最新发布
weixin_33557333的博客
10-19 1021
设随机变量 $ X $ 具有期望 $ \mu = E[X] $ 和标准差 $ \sigma = \sqrt{E[(X - \mu)^2]} $，则其三阶中心矩定义为：这一数学表达捕捉的是偏差立方后的平均值。由于立方运算保留符号信息（即 $(x - \mu)^3$ 在 $ x < \mu $ 时为负，在 $ x > \mu $ 时为正），因此 $\mu_3$ 的符号直接反映了分布偏离对称轴的方向。
7.带你入门matlab偏斜度和峰度（matlab程序）
m0_57943157的博客
07-14 4058
峰度衡量数据分布的平坦度（flatness），即数据取值分布形态陡缓程度的统计量。尾部大的数据分布，其峰度值较大。正态分布的峰度值为3。图 1所示，黑线服从尖峰（leptokurtic）、厚尾（thick-tailed）分布的峰度值大于3。式中， K表示峰度（无量纲）；峰度（Kurtosis）与偏态（Skewness）就是量测数据正态分布特性的两个指标。如图2所示，右偏态为正，表明平均值大于中位数。0说明是最完美的对称性，正态分布的偏态就是0。峰度kurtosis：用于度量x偏离某分布的程度。
181225 Matlab图解峰度kurtosis与偏度skewness
专注机器学习之路
12-25 1万+
峰度kurtosis：用于度量x偏离某分布的程度。 正太分布的峰度是3； 当时间序列的曲线峰值比正太分布的高时，峰度大于3； 当比正太分布的低时，峰度小于3。 偏度skewness：用于衡量x的对称性。 对于正太分布，偏度为0; 若偏度为正，则x均值左侧的离散度比右侧弱； 若偏度为负，则x均值左侧的离散度比右侧强。 正态分布的偏度系数和峰度系数 ...
Fast_Kurtogram.m
05-27
matlab快速峭度谱，故障诊断专用，希望大家能够使用的到，不求积分，只求帮助大家。
基于频率滑动广义互相关算法的信号时延精确估计方法在MATLAB R2018A中的应用与研究 时延估计
09-03
基于频率滑动广义互相关的信号时延估计方法，探讨了其在MATLAB R2018A环境下的具体实现。时延估计作为信号处理领域的核心技术，对于确定信源距离、方位、速度等关键信息至关重要。文中首先概述了时延估计的重要性和...
MATLAB R2018a中的变分模态分解与Cramer von Mises检验结合的信号降噪方法
08-24
内容概要：本文探讨了一种创新的信号降噪方法，该方法结合了变分模态分解（VMD）和Cramer von Mises检验，在MATLAB R2018a环境中实现了对一维时间序列信号的有效降噪。文中详细介绍了算法的工作流程，从数据准备到...
MATLAB R2018A中分数阶全变分泊松噪声处理：非盲与盲反卷积算法研究
07-29
在MATLAB R2018A环境中，针对分数阶全变分泊松噪声的非盲和盲反卷积方法的研究与实践。反卷积作为一种重要图像处理技术，旨在恢复因模糊和噪声导致的质量下降的图像。文中首先解释了非盲反卷积和盲反卷积的区别及其...
基于频率滑动广义互相关算法的信号时延精确估计方法在MATLAB R2018A中的应用与研究,基于MATLAB R2018A的频率滑动广义互相关算法：时延估计及其在信号处理领域的广泛应用,基于频率滑动广
02-14
基于频率滑动广义互相关算法的信号时延精确估计方法在MATLAB R2018A中的应用与研究,基于MATLAB R2018A的频率滑动广义互相关算法：时延估计及其在信号处理领域的广泛应用,基于频率滑动广义互相关的信号时延估计方法...
matlab中统计工具箱函数大全
RickyWong_的专栏
03-03 4153
转自：http://blog.sina.com.cn/s/blog_7671b3eb01013vma.html MATLAB统计工具箱包括概率分布、方差分析、假设检验、分布检验、非参数检验、回归分析、判别分析、主成分分析、因子分析、系统聚类分析、K均值聚类分析、试验设计、决策树、多元方差分析、统计过程控制和统计图形绘制等。优化工具箱包括无约束最优化、有约束最优化、二次规划、多目标规划、最大最
图像峰度计算函数（matlab）
05-24
基于matlab的图像峰度计算函数，可运行。基于matlab的图像峰度计算函数，可运行。
Fast- Kurtogram_齿轮故障诊断_快速谱峭度源代码_包络滤波_齿轮故障_包络峭度
09-11
计算信号的峭度，寻找共振频带，并对共振频带进行滤波后包络解调，适合于轴承、齿轮等故障诊断的信号处理中。
Kuramoto 的同步振荡器模型：Kuramoto 模型是耦合振荡器的非线性动态系统，最初具有随机的固有频率和相位。-matlab开发
05-29
Kuramoto 模型是耦合振荡器的非线性动态系统，最初具有随机的固有频率和相位。 如果耦合足够强，系统将演变为所有振荡器同相的系统。 见克利夫角， https: //blogs.mathworks.com/cleve/2019/10/30/stability-of-kuramoto-oscillators/
【故障诊断】基于负熵诱导灰狼优化算法的多目标信息频带选择用于滚动轴承故障诊断（Matlab代码实现）
weixin_46039719的博客
04-03 818
在SK的上述改进中，这些指标都用于定量评估带通滤波信号或其解析包络，以检测脉冲可以归类为时域估计器，其最大的缺点是容易受到异常值的影响，如上所述。因此，时域和频域的稀疏性应该达到相对平衡，即IFB选择的折衷。]被视为IFB选择的一个里程碑，其中具有不同中心频率和带宽的准分析滤波器组输出处的系数峰度在频率/频率分辨率平面中用于表示和检测振动信号中的非平稳性。[2] 顾晓， 杨淑， 刘彦， 等. 基于负熵诱导灰狼优化器的多目标信息频带选择用于滚动轴承故障诊断.传感器 2020， 20（7）， 1845.
matlab峰度和偏度,峰度和偏度
weixin_42365065的博客
03-20 9535
峰度和偏度峰度(kurtosis)和偏度(skewness)张贴者： EasyTrader这是一个应用统计资料分佈的峰度(kurtosis)和偏度(skewness)所作的交易策略。当资料呈现趋势性，并且潜在趋势为正时，我们做多。当资料呈现趋势性，并且潜在趋势为负时，我们做空。当趋势发生反转后，我们平仓。那么，我们如何确定趋势和趋势的强度呢？让我们先来复习一下峰度和偏度的定义。首先来看常态分佈 (...
MATLAB常见函数的使用
让一让的博客
01-22 1819
基本统计量 对于给定的随机变量X,可通过如下函数求得其常见的统计量： 均值：mean(X) 中位数：median(X) 标准差：std(X) 方差：var(X) 偏度：skewness(X) 蜂度：kurtosis(X) 常见概率分布 MATLAB统计工具箱中有20多种概率分布，其中较为常见到概率分布命令的字符如下： 卡方分布：chi2 t 分布：t F分布：F β 分布...
Ubuntu 18.04.4 安装MATLAB R2018a：权限与教程
本文档主要介绍了在Ubuntu 18.04.4操作系统环境下安装MATLAB R2018a的具体步骤，以及作者对Ubuntu安装软件特点的理解，特别是关于权限管理的问题。以下是详细的步骤总结： 1. **下载与授权**: - 从指定的百度网盘...
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
表情包 代码片
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值