CN106096198B - 一种基于变分模式分解和谱峭度的包络分析方法 - Google Patents
一种基于变分模式分解和谱峭度的包络分析方法 Download PDFInfo
- Publication number
- CN106096198B CN106096198B CN201610492074.0A CN201610492074A CN106096198B CN 106096198 B CN106096198 B CN 106096198B CN 201610492074 A CN201610492074 A CN 201610492074A CN 106096198 B CN106096198 B CN 106096198B
- Authority
- CN
- China
- Prior art keywords
- signal
- envelope
- spectrum
- data
- frequency
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于变分模式分解和谱峭度的包络分析方法,该方法首先利用变分模式分解方法对原始信号进行分解,然后利用数据的重排和替代操作排除分解结果中的噪声分量和趋势项,接着再采用谱峭度方法对第一次滤波后的信号进行分析,得到最优滤波器的中心频率和带宽,然后利用该滤波器对第一次滤波后的信号再进行第二次滤波,然后采用三次样条迭代平滑包络分析方法对第二次滤波后的信号进行包络分析,最后根据包络谱确定旋转机械的故障类型。本发明适合于处理复杂的旋转机械故障信号,能够准确地判定出旋转机械的故障类型,具有良好的抗噪性和鲁棒性,便于工程应用。
Description
技术领域
本发明涉及旋转机械状态监测与故障诊断领域,具体涉及一种基于变分模式分解和谱峭度的包络分析方法。
背景技术
包络分析技术广泛应用于齿轮和滚动轴承的故障诊断中。现有的包络分析技术有下面三个缺陷:①现有的包络分析技术或者是直接对原始信号进行分析,或者是仅对原始信号进行简单的滤波后再进行分析,因此现有的方法容易受到噪声、趋势及其它成分的干扰,从而导致现有技术的分析精度较低;②现有的包络分析技术是以Hilbert变换为基础,而Hilbert变换要求被分析的信号必须是单分量的窄带信号,否则信号的频率调制部分将要污染信号的幅值包络分析结果,但是目前待分析的信号都不严格满足单分量且窄带的条件,这样就会导致现有技术因精度不高而容易出现误判问题;③由传统方法得到的包络谱存在着端点效应。
发明内容
本发明要解决的问题是针对以上不足,提出一种基于变分模式分解和谱峭度的包络分析方法,采用本发明的包络分析方法后,具有分析结果准确度和精确度高,并能准确地检测出旋转机械故障类型的优点。
为解决以上技术问题,本发明采取的技术方案如下:一种基于变分模式分解和谱峭度的包络分析方法,其特征在于,包括以下步骤:
步骤1:利用加速度传感器以采样频率fs测取旋转机械的振动信号x(k), (k=1,2, …,N),N为采样信号的长度;
步骤2:采用变分模式分解(Variational Mode Decomposition)算法将信号x(k)分解成n个分量之和,即,其中,ci(k)代表由变分模式分解算法得到的第i个分量,变分模式分解已公知,见文献
Konstantin Dragomiretskiy, Dominique Zosso. Variational ModeDecomposition, IEEE TRANSACTIONS ON SIGNAL PROCESSING, 2014, 62(3) : 531-544;本例中,设置模态数目为10;
步骤3:对ci(k)执行重排操作和替代操作,经重排操作得到的数据用ci shuffle(k)表示,替代操作后得到数据用ci FTran(k)表示;
步骤4:对ci(k)、ci shuffle(k)和ci FTran(k)分别执行多重分形去趋势波动分析(Multifractal Detrended Fluctuation Analysis, MFDFA),得到广义Hurst指数曲线,ci(k)的广义Hurst指数曲线用Hi(q)表示;ci shuffle(k)的广义Hurst指数曲线用Hi shuffle(q)表示;ci FTran(k)的广义Hurst指数曲线用Hi FTran(q)表示;
步骤5:如果Hi(q) 与Hi shuffle(q)或Hi(q) 与Hi FTran(q)之间的相对误差小于5%,或者Hi(q) 、Hi shuffle(q) 和Hi FTran(q)三者都不随q而变化,则抛弃对应的ci(k)分量;
步骤6:对剩余的ci(k)分量求和,将该和记为信号经重排和替代滤波后的结果xf1(k);
步骤7:对xf1(k)执行谱峭度分析,求出信号峭度最大处所对应的中心频率f0和带宽B;
步骤8: 根据中心频率f0和带宽B对xf1(k)进行带通滤波,得到xf2(k);
步骤9:对信号xf2(k)执行三次样条迭代平滑包络分析,得到信号包络eov(k);
步骤10:对得到的信号包络eov(k)执行离散傅里叶变换得到包络谱,根据包络谱特征频率判断机器的故障类型。
一种优化方案,所述步骤3中数据重排操作包括以下步骤:
随机打乱分量ci(k)的排列顺序。
进一步地,所述步骤3中数据替代操作包括以下步骤:
1) 对分量ci(k)执行离散傅里叶变换,获得分量ci(k)的相位;
2) 用一组位于(-π,π)区间内的伪独立同分布数来代替分量ci(k)的原始相位;
3) 对经过相位替代后的频域数据执行离散傅里叶逆变换得到数据ci IFFT(k),求取数据ci IFFT(k)的实部。
进一步地,所述步骤4中MFDFA方法包括以下步骤:
1)构造x(k)(k=1,2,…,N)的轮廓Y(i):
x(k)代表所述步骤4中的ci(k)或ci shuffle(k)或ci FTran(k);
2)将信号轮廓Y(i)分成不重叠的NS段长度为s的数据,由于数据长度N通常不能整除s,所以会剩余一段数据不能利用;
为了充分利用数据的长度,再从数据的反方向以相同的长度分段,这样一共得到2NS段数据;
3)利用最小二乘法拟合每段数据的多项式趋势,然后计算每段数据的方差:
yv(i)为拟合的第v段数据的趋势,若拟合的多项式趋势为m阶,则记该去趋势过程为(MF-)DFAm;本例中,m=1;
4) 计算第q阶波动函数的平均值:
;
5)如果x(k)存在自相似特征,则第q阶波动函数的平均值Fq(s)和时间尺度s之间存在幂律关系:
当q=0时,步骤4)中的公式发散,这时H(0)通过下式所定义的对数平均过程来确定:
6)对步骤5)中的公式两边取对数可得ln[Fq(s)]=H(q)ln(s)+c(c为常数),由此可以获得直线的斜率H(q)。
进一步地,所述步骤7中的谱峭度方法包括以下步骤:
1)构造一个截止频率为fc=0.125+ε的低通滤波器h(n);ε>0,本例中fc=0.3;
2)基于h(n)构造通频带为[0, 0.25]的准低通滤波器h0(n)和通频带为[0.25,0.5]的准高通滤波器h1(n),
;
3)信号ci k(n)经 h0(n)、 h1(n)滤波并降采样后分解成低频部分c2i k+1(n)和高频部分c2i+1 k+1(n),降采样的因子为2,再经多次迭代滤波后形成滤波器树,第k层有2k个频带,其中ci k(n)表示滤波器树中第k层上的第i个滤波器的输出信号,i=0,…, 2k-1,0≤k≤K-1,本例中K=8;c0 (n)代表所述步骤7中xf1(k);
4)分解树中第k层上的第i个滤波器的中心频率fki和带宽Bk分别为
;
5)计算每一个滤波器结果ci k(n)( i=0,…, 2k-1) 的峭度;
6)将所有的谱峭度汇总,得到信号总的谱峭度。
进一步地,所述步骤9中的三次样条迭代平滑包络分析方法包括以下步骤:
1)计算信号z(k)的绝对值∣z(k)∣的局部极值;在第1次迭代中,z(k)代表所述步骤9中xf2(k);
2)采用三次样条曲线拟合局部极值点得到包络线eov1(k);
3)对z(k)进行归一化处理得到;
4)第2次迭代:把z 1(k)重新作为新数据,重复执行上述步骤1)~3),得到;
5)第i次迭代:把z i-1(k) 重新作为新数据,重复执行上述步骤1)~3),得到;
6) 如果第n次迭代得到的z n (k)的幅值小于或等于1,则迭代过程停止,最后得到信号z(k)的包络为。
本发明采用以上技术方案,与现有技术相比,本发明具有以下优点:
1)利用变分模式分解对原始信号进行分解,然后利用数据的重排和替代操作排除其中的噪声和趋势分量,仅仅保留信号分量中的有用成分,从而避免了噪声和趋势分量对包络分析结果的影响,分析结果准确度和精确度高。
2)利用三次样条迭代平滑包络分析方法将信号包络与频率调制部分完全分离,能够避免频率调制部分对信号包络分析结果的影响,从而提高包络分析的精度。
3) 能够准确地检测出旋转机械的故障类型。
4) 由传统方法得到的包络谱存在端点效应,而由本发明得到的包络谱能够避免端点效应。
下面结合附图和实施例对本发明做进一步说明。
附图说明
附图1为本发明实施例中本发明方法的流程图。
附图2为本发明实施例中采用低通滤波器和高通滤波器对信号进行初步分解的示意图。
附图3为本发明实施例中采用树状滤波器结构快速计算谱峭度的示意图。
附图4为本发明实施例中具有内圈故障的滚动轴承振动信号。
附图5为本发明实施例中采用传统包络分析方法对内圈故障滚动轴承振动信号的分析结果。
附图6为本发明实施例中本发明对内圈故障滚动轴承振动信号的分析结果。
附图7为本发明实施例中具有外圈故障的滚动轴承振动信号。
附图8为本发明实施例中采用传统包络分析方法对外圈故障滚动轴承振动信号的分析结果。
附图9为本发明实施例中本发明对外圈故障滚动轴承振动信号的分析结果。
具体实施方式
实施例,如图1、图2、图3所示,一种基于变分模式分解和谱峭度的包络分析方法,包括以下步骤:
步骤1:利用加速度传感器以采样频率fs测取旋转机械的振动信号x(k), (k=1,2, …,N),N为采样信号的长度;
步骤2:采用变分模式分解(Variational Mode Decomposition)算法将信号x(k)分解成n个分量之和,即 ,其中,ci(k)代表由变分模式分解算法得到的第i个分量,变分模式分解已公知,见文献
Konstantin Dragomiretskiy, Dominique Zosso. Variational ModeDecomposition, IEEE TRANSACTIONS ON SIGNAL PROCESSING, 2014, 62(3) : 531-544;本例中,设置模态数目为10;
步骤3:对ci(k)执行重排操作和替代操作,经重排操作得到的数据用ci shuffle(k)表示,替代操作后得到数据用ci FTran(k)表示;
步骤4:对ci(k)、ci shuffle(k)和ci FTran(k)分别执行多重分形去趋势波动分析(Multifractal Detrended Fluctuation Analysis, MFDFA),得到广义Hurst指数曲线,ci(k)的广义Hurst指数曲线用Hi(q)表示;ci shuffle(k)的广义Hurst指数曲线用Hi shuffle(q)表示;ci FTran(k)的广义Hurst指数曲线用Hi FTran(q)表示;
步骤5:如果Hi(q) 与Hi shuffle(q)或Hi(q) 与Hi FTran(q)之间的相对误差小于5%,或者Hi(q) 、Hi shuffle(q) 和Hi FTran(q)三者都不随q而变化,则抛弃对应的ci(k)分量;
步骤6:对剩余的ci(k)分量求和,将该和记为信号经重排和替代滤波后的结果xf1(k);
步骤7:对xf1(k)执行谱峭度分析,求出信号峭度最大处所对应的中心频率f0和带宽B;
步骤8: 根据中心频率f0和带宽B对xf1(k)进行带通滤波,得到xf2(k);
步骤9:对信号xf2(k)执行三次样条迭代平滑包络分析,得到信号包络eov(k);
步骤10:对得到的信号包络eov(k)执行离散傅里叶变换得到包络谱,根据包络谱特征频率判断机器的故障类型。
步骤3中数据重排操作包括以下步骤:
随机打乱分量ci(k)的排列顺序。
步骤3中数据替代操作包括以下步骤:
1) 对分量ci(k)执行离散傅里叶变换,获得分量ci(k)的相位;
2) 用一组位于(-π,π)区间内的伪独立同分布数来代替分量ci(k)的原始相位;
3) 对经过相位替代后的频域数据执行离散傅里叶逆变换得到数据ci IFFT(k),求取数据ci IFFT(k)的实部。
步骤4中MFDFA方法包括以下步骤:
1)构造x(k)(k=1,2,…,N)的轮廓Y(i):
x(k)代表所述步骤4中的ci(k)或ci shuffle(k)或ci FTran(k);
2)将信号轮廓Y(i)分成不重叠的NS段长度为s的数据,由于数据长度N通常不能整除s,所以会剩余一段数据不能利用;
为了充分利用数据的长度,再从数据的反方向以相同的长度分段,这样一共得到2NS段数据;
3)利用最小二乘法拟合每段数据的多项式趋势,然后计算每段数据的方差:
yv(i)为拟合的第v段数据的趋势,若拟合的多项式趋势为m阶,则记该去趋势过程为(MF-)DFAm;本例中,m=1;
4) 计算第q阶波动函数的平均值:
;
5)如果x(k)存在自相似特征,则第q阶波动函数的平均值Fq(s)和时间尺度s之间存在幂律关系:
当q=0时,步骤4)中的公式发散,这时H(0)通过下式所定义的对数平均过程来确定:
6)对步骤5)中的公式两边取对数可得ln[Fq(s)]=H(q)ln(s)+c(c为常数),由此可以获得直线的斜率H(q)。
步骤7中的谱峭度方法包括以下步骤:
1)构造一个截止频率为fc=0.125+ε的低通滤波器h(n);ε>0,本例中fc=0.3;
2)基于h(n)构造通频带为[0, 0.25]的准低通滤波器h0(n)和通频带为[0.25,0.5]的准高通滤波器h1(n),
;
3)信号ci k(n)经 h0(n)、 h1(n)滤波并降采样后分解成低频部分c2i k+1(n)和高频部分c2i+1 k+1(n),降采样的因子为2,再经多次迭代滤波后形成滤波器树,第k层有2k个频带,其中ci k(n)表示滤波器树中第k层上的第i个滤波器的输出信号,i=0,…, 2k-1,0≤k≤K-1,本例中K=8;c0 (n)代表所述步骤7中xf1(k);
4)分解树中第k层上的第i个滤波器的中心频率fki和带宽Bk分别为
;
5)计算每一个滤波器结果ci k(n)( i=0,…, 2k-1) 的峭度;
6)将所有的谱峭度汇总,得到信号总的谱峭度。
步骤9中的三次样条迭代平滑包络分析方法包括以下步骤:
1)计算信号z(k)的绝对值∣z(k)∣的局部极值;在第1次迭代中,z(k)代表所述步骤9中xf2(k);
2)采用三次样条曲线拟合局部极值点得到包络线eov1(k);
3)对z(k)进行归一化处理得到;
4)第2次迭代:把z 1(k)重新作为新数据,重复执行上述步骤1)~3),得到;
5)第i次迭代:把z i-1(k) 重新作为新数据,重复执行上述步骤1)~3),得到;
6) 如果第n次迭代得到的z n (k)的幅值小于或等于1,则迭代过程停止,最后得到信号z(k)的包络为。
试验1,利用具有内圈故障的滚动轴承振动数据对本发明所述算法的性能进行验证。
实验所用轴承为6205-2RS JEM SKF,利用电火花加工方法在轴承内圈上加工深度为0.2794mm、宽度为0.3556mm的凹槽来模拟轴承内圈故障,本实验负载约为0.7457kW,驱动电机转频约为29.5Hz,轴承内圈故障特征频率约为160Hz,采样频率为4.8KHz,信号采样时长为1s。
采集到的内圈故障信号如图4所示。
首先采用传统的包络分析方法对图4所示的信号进行分析,得到的分析结果如图5所示。从图5可以看出,轴承的故障特征完全被掩盖,因此传统的包络分析方法不能有效地提取轴承的故障特征;此外,从图5可以看出,包络谱的左端点存在着异常高值,这说明由传统方法得到的包络谱存在着端点效应。
采用本发明所提出的方法对图4所示的信号进行分析,得到的分析结果如图6所示。从图6可以看出,160Hz和320Hz所对应的谱线明显高于其它谱线,这两个频率分别对应轴承内圈故障特征频率的1倍频和2倍频,据此可以判断轴承具有内圈故障;从图6可以看出,由本发明得到的包络谱没有端点效应。
经多次实验表明,在负载和故障尺寸深度不变的情况下,本发明能够可靠识别的最小内圈故障尺寸宽度约为0.24 mm,而传统方法能够可靠识别的最小内圈故障尺寸宽度约为0.53mm,精度提高54.7%。
试验2,利用具有外圈故障的滚动轴承振动数据对本发明所述算法的性能进行验证。
实验所用轴承为6205-2RS JEM SKF,利用电火花加工方法在轴承外圈上加工深度为0.2794mm、宽度为0.5334mm的凹槽来模拟轴承外圈故障,本实验负载约为2.237 kW,驱动电机转频约为28.7Hz,轴承外圈故障特征频率约为103Hz,采样频率为4.8KHz,信号采样时长为1s。
采集到的外圈故障信号如图7所示。
首先采用传统的包络分析方法对图7所示的信号进行分析,得到的分析结果如图8所示。从图8可以看出,轴承的故障特征完全被掩盖,因此传统的包络分析方法不能有效地提取轴承的故障特征;此外,从图8可以看出,包络谱的左端点存在着异常高值,这说明由传统方法得到的包络谱存在着端点效应。
采用本发明所提出的方法对图7所示的信号进行分析,得到的分析结果如图9所示。从图9可以看出,103Hz和206Hz所对应的谱线明显高于其它谱线,这两个频率分别对应轴承外圈故障特征频率的1倍频和2倍频,据此可以判断轴承具有外圈故障;从图9可以看出,由本发明得到的包络谱没有端点效应。
经多次实验表明,在负载和故障尺寸深度不变的情况下,本发明能够可靠识别的最小外圈故障尺寸宽度约为0.35mm,而传统方法能够可靠识别的最小外圈故障尺寸宽度约为0.68mm,精度提高48.5%。
根据试验结果,分析后认为:
1) 传统的包络分析方法直接对原始信号进行包络分析,或者对仅经过简单处理后的原始信号进行包络分析,与传统的包络分析方法不同,本发明首先利用变分模式分解对原始信号进行分解,然后利用数据的重排和替代操作排除其中的噪声和趋势分量,仅仅保留信号分量中的有用成分,从而避免了噪声和趋势分量对包络分析结果的影响,提高了准确度和精确度。
2) 传统的包络分析方法以Hilbert变换为基础,而Hilbert变换要求被分析的信号必须是单分量的窄带信号,否则信号的频率调制部分将要污染信号的包络分析结果,但是目前待分析的信号都不严格满足单分量且窄带的条件,这样就会导致现有技术因精度不高而容易出现误判问题,与传统包络分析方法不同,本发明利用三次样条迭代平滑包络分析方法将信号包络与频率调制部分完全分离,能够避免频率调制部分对信号包络分析结果的影响,从而提高包络分析的精度。
3)能够准确地检测出旋转机械的故障类型。
4) 由传统方法得到的包络谱存在端点效应,而由本发明得到的包络谱能够避免端点效应。
5)各步骤作用:
第1)步:采集振动信号;
第2)步:将原始信号分解成不同分量和的形式,其中有些分量对应噪声和趋势项,有些分量对应有用信号;
第3)~5)步:对上述分解得到的信号执行重排操作和替代操作,剔除其中的噪声分量和趋势项,只保留有用信号;
第6)步:将剩余的有用信号求和,将该和作为信号经重排和替代滤波后的结果xf1(k);
第7)步:对滤波后的信号xf1(k)执行谱峭度分析,求出信号最大峭度处对应的中心频率f0和带宽B;
第8)步:根据中心频率f0和带宽B对xf1(k)进行带通滤波,得到信号xf2(k);
第9)步:计算信号xf2(k)的包络eov(k);
第10)步:对eov(k)执行离散傅里叶变换得到包络谱,根据包络谱判断轴承的故障类型。
本领域技术人员应该认识到,上述的具体实施方式只是示例性的,是为了使本领域技术人员能够更好的理解本发明内容,不应理解为是对本发明保护范围的限制,只要是根据本发明技术方案所作的改进,均落入本发明的保护范围。
Claims (5)
1.一种基于变分模式分解和谱峭度的包络分析方法,其特征在于,包括以下步骤:
步骤1:利用加速度传感器以采样频率fs测取旋转机械的振动信号x(k), k=1, 2, …,N,N为采样信号的长度;
步骤2:采用变分模式分解(Variational Mode Decomposition)算法将信号x(k)分解成n个分量之和,即,其中,ci(k)代表由变分模式分解算法得到的第i个分量;
步骤3:对ci(k)执行重排操作和替代操作,经重排操作得到的数据用ci shuffle(k)表示,替代操作后得到数据用ci FTran(k)表示;
步骤4:对ci(k)、ci shuffle(k)和ci FTran(k)分别执行多重分形去趋势波动分析(Multifractal Detrended Fluctuation Analysis, MFDFA),得到广义Hurst指数曲线,ci(k)的广义Hurst指数曲线用Hi(q)表示;ci shuffle(k)的广义Hurst指数曲线用Hi shuffle(q)表示;ci FTran(k)的广义Hurst指数曲线用Hi FTran(q)表示;
步骤5:如果Hi(q) 与Hi shuffle(q)或Hi(q) 与Hi FTran(q)之间的相对误差小于5%,或者Hi(q) 、Hi shuffle(q) 和Hi FTran(q)三者都不随q而变化,则抛弃对应的ci(k)分量;
步骤6:对剩余的ci(k)分量求和,将该和记为信号经重排和替代滤波后的结果xf1(k);
步骤7:对xf1(k)执行谱峭度分析,求出信号峭度最大处所对应的中心频率f0和带宽B;
步骤8: 根据中心频率f0和带宽B对xf1(k)进行带通滤波,得到xf2(k);
步骤9:对信号xf2(k)执行三次样条迭代平滑包络分析,得到信号包络eov(k);
步骤10:对得到的信号包络eov(k)执行离散傅里叶变换得到包络谱,根据包络谱特征频率判断机器的故障类型。
2.根据权利要求1所述的一种基于变分模式分解和谱峭度的包络分析方法,其特征在于,所述步骤3中数据重排操作包括以下步骤:
随机打乱分量ci(k)的排列顺序。
3.根据权利要求1所述的一种基于变分模式分解和谱峭度的包络分析方法,其特征在于:所述步骤3中数据替代操作包括以下步骤:
1) 对分量ci(k)执行离散傅里叶变换,获得分量ci(k)的相位;
2) 用一组位于(-π,π)区间内的伪独立同分布数来代替分量ci(k)的原始相位;
3) 对经过相位替代后的频域数据执行离散傅里叶逆变换得到数据ci IFFT(k),求取数据ci IFFT(k)的实部。
4.根据权利要求1所述的一种基于变分模式分解和谱峭度的包络分析方法,其特征在于:所述步骤7中的谱峭度方法包括以下步骤:
1)构造一个截止频率为fc=0.125+ε的低通滤波器h(n);ε>0,ε=0.175,本例中fc=0.3;
2)基于h(n)构造通频带为[0, 0.25]的准低通滤波器h0(n)和通频带为[0.25, 0.5]的准高通滤波器h1(n),
;
3)信号ci k(n)经 h0(n)、 h1(n)滤波并降采样后分解成低频部分c2i k+1(n)和高频部分c2i +1 k+1(n),降采样的因子为2,再经多次迭代滤波后形成滤波器树,第k层有2k个频带,其中ci k(n)表示滤波器树中第k层上的第i个滤波器的输出信号,i=0,…, 2k-1,K代表滤波器树的层数,0≤k≤K-1,本例中K=8;c0 (n)代表所述步骤7中xf1(k);
4)分解树中第k层上的第i个滤波器的中心频率fki和带宽Bk分别为
;
5)计算每一个滤波器结果ci k(n) i=0,…, 2k-1 的峭度;
6)将所有的谱峭度汇总,得到信号总的谱峭度。
5.根据权利要求1所述的一种基于变分模式分解和谱峭度的包络分析方法,其特征在于,所述步骤9中的三次样条迭代平滑包络分析方法包括以下步骤:
1)计算信号z(k)的绝对值∣z(k)∣的局部极值;在第1次迭代中,z(k)代表所述步骤9中xf2(k);
2)采用三次样条曲线拟合局部极值点得到包络线eov1(k);
3)对z(k)进行归一化处理得到;
4)第2次迭代:把z 1(k)重新作为新数据,重复执行上述步骤1)~3),得到;
5)第i次迭代:把z i-1(k) 重新作为新数据,重复执行上述步骤1)~3),得到;
6) 如果第n次迭代得到的z n (k)的幅值小于或等于1,则迭代过程停止,最后得到信号包络为。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610492074.0A CN106096198B (zh) | 2016-06-29 | 2016-06-29 | 一种基于变分模式分解和谱峭度的包络分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610492074.0A CN106096198B (zh) | 2016-06-29 | 2016-06-29 | 一种基于变分模式分解和谱峭度的包络分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106096198A CN106096198A (zh) | 2016-11-09 |
CN106096198B true CN106096198B (zh) | 2019-02-22 |
Family
ID=57214856
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610492074.0A Expired - Fee Related CN106096198B (zh) | 2016-06-29 | 2016-06-29 | 一种基于变分模式分解和谱峭度的包络分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106096198B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106680876B (zh) * | 2017-01-22 | 2019-04-12 | 中国石油大学(华东) | 一种地震数据联合去噪方法 |
CN107944199B (zh) * | 2017-12-22 | 2020-12-01 | 浙江工业大学 | 一种基于频谱趋势和变分模态分解的齿轮箱故障识别方法 |
CN108692936B (zh) * | 2018-03-27 | 2020-03-13 | 四川大学 | 基于参数自适应vmd的机械故障诊断方法 |
CN109253244B (zh) * | 2018-11-22 | 2021-03-23 | 常州信息职业技术学院 | 一种多齿轮系统大型机械齿轮箱故障检测方法 |
CN109633268B (zh) * | 2018-12-20 | 2020-03-10 | 北京航空航天大学 | 一种基于b样条和直方图的方波基频辨识方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102620928A (zh) * | 2012-03-02 | 2012-08-01 | 燕山大学 | 基于小波半软阈值和emd的风电齿轮箱故障诊断方法 |
CN103837345A (zh) * | 2014-03-25 | 2014-06-04 | 上海电机学院 | 齿轮箱故障诊断方法及装置 |
CN104019831A (zh) * | 2014-06-20 | 2014-09-03 | 哈尔滨工业大学 | 基于emd和熵权的陀螺仪故障诊断方法 |
CN104729853A (zh) * | 2015-04-10 | 2015-06-24 | 华东交通大学 | 一种滚动轴承性能退化评估装置及方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130096848A1 (en) * | 2011-10-13 | 2013-04-18 | Charles Terrance Hatch | Methods and systems for automatic rolling-element bearing fault detection |
-
2016
- 2016-06-29 CN CN201610492074.0A patent/CN106096198B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102620928A (zh) * | 2012-03-02 | 2012-08-01 | 燕山大学 | 基于小波半软阈值和emd的风电齿轮箱故障诊断方法 |
CN103837345A (zh) * | 2014-03-25 | 2014-06-04 | 上海电机学院 | 齿轮箱故障诊断方法及装置 |
CN104019831A (zh) * | 2014-06-20 | 2014-09-03 | 哈尔滨工业大学 | 基于emd和熵权的陀螺仪故障诊断方法 |
CN104729853A (zh) * | 2015-04-10 | 2015-06-24 | 华东交通大学 | 一种滚动轴承性能退化评估装置及方法 |
Non-Patent Citations (2)
Title |
---|
Fault diagnosis of rolling bearings based on multifractal detrended fluctuation analysis and Mahalanobis distance criterion;Jinshan Lin 等;《Machanical Systems and Singal Processing》;20130306;第38卷(第2期);第515-533页 |
基于多重分形去趋势互相关分析的齿轮箱故障诊断;林近山 等;《机械传动》;20160115;第40卷(第1期);第91-94、105页 |
Also Published As
Publication number | Publication date |
---|---|
CN106096198A (zh) | 2016-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106096313B (zh) | 一种基于奇异谱分解和谱峭度的包络分析方法 | |
CN106198015B (zh) | 一种滚动轴承的vmd、谱峭度和平滑迭代包络分析方法 | |
CN106096200B (zh) | 一种基于小波分解和谱峭度的包络分析方法 | |
CN106096198B (zh) | 一种基于变分模式分解和谱峭度的包络分析方法 | |
CN106153339B (zh) | 一种基于变分模式分解滤波的包络分析方法 | |
CN106198013B (zh) | 一种基于经验模式分解滤波的包络分析方法 | |
CN106096199B (zh) | 一种滚动轴承的wt、谱峭度和平滑迭代包络分析方法 | |
CN105954031B (zh) | 一种基于奇异谱分解滤波的包络分析方法 | |
CN106053069B (zh) | 一种滚动轴承的ssd、谱峭度和平滑迭代包络分析方法 | |
CN106198009B (zh) | 一种滚动轴承的emd、谱峭度和平滑迭代包络分析方法 | |
CN106198010B (zh) | 一种基于局部均值分解滤波的包络分析方法 | |
CN106198012B (zh) | 一种基于局部均值分解和谱峭度的包络分析方法 | |
CN106053060B (zh) | 一种基于非线性模式分解滤波的包络分析方法 | |
CN106153333B (zh) | 一种基于小波分解滤波的包络分析方法 | |
CN106096201B (zh) | 一种旋转机械的eemd和三次样条平滑包络分析方法 | |
CN106053061B (zh) | 一种基于非线性模式分解和谱峭度的包络分析方法 | |
CN106198014B (zh) | 一种基于经验模式分解和谱峭度的包络分析方法 | |
CN106198017B (zh) | 一种滚动轴承的lmd、谱峭度和平滑迭代包络分析方法 | |
CN105954030B (zh) | 一种基于内秉时间尺度分解和谱峭度的包络分析方法 | |
CN106198018B (zh) | 一种旋转机械的eemd和平滑迭代包络分析方法 | |
CN106053059B (zh) | 一种基于内秉时间尺度分解滤波的包络分析方法 | |
CN105973603B (zh) | 一种旋转机械的eemd和有理样条平滑包络分析方法 | |
CN106124200B (zh) | 一种旋转机械的elmd和平滑迭代包络分析方法 | |
CN106198016B (zh) | 一种滚动轴承的nmd、谱峭度和平滑迭代包络分析方法 | |
CN106153338B (zh) | 一种旋转机械的elmd和有理样条平滑包络分析方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190222 Termination date: 20210629 |
|
CF01 | Termination of patent right due to non-payment of annual fee |