CN103226649A - 机械早期及复合故障的集成噪声重构经验模式分解方法 - Google Patents
机械早期及复合故障的集成噪声重构经验模式分解方法 Download PDFInfo
- Publication number
- CN103226649A CN103226649A CN2013100975026A CN201310097502A CN103226649A CN 103226649 A CN103226649 A CN 103226649A CN 2013100975026 A CN2013100975026 A CN 2013100975026A CN 201310097502 A CN201310097502 A CN 201310097502A CN 103226649 A CN103226649 A CN 103226649A
- Authority
- CN
- China
- Prior art keywords
- noise
- component
- imf
- signal
- emd
- 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.)
- Granted
Links
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种机械早期及复合故障的集成噪声重构经验模式分解方法,它首先对信号中的噪声分量进行估算和重采样,并重构若干输入信号。利用噪声极值点的变化来改善信号不连续问题,为避免模式混淆的信号分解提供前提。然后,对这些重构输入信号进行经验模式分解(EMD)得到相应的基本模式分量(IMF)。最后,利用随机序列统计均值为零的原理,将所有对应的IMF做总体平均运算,消除信号中噪声分量的影响,从而达到信号降噪与特征增强的目的。本发明非常适合于信噪比相对较低的机械早期及复合故障信号处理及特征提取,识别精度高,运算实时性好,简单易行,便于工程实践中使用。
Description
技术领域
本发明涉及大型机械设备(航空发动机、汽轮发电机组、连铸连轧机组、风电装备等)故障诊断技术,具体涉及一种机械设备的早期故障和复合故障的特征提取与识别方法。
背景技术
机械设备一旦出现事故,将带来巨大的经济损失和人员伤亡。在工程实践中,我们期望对故障的发生与发展能够做到防微杜渐,不期望亡羊补牢般地处理事故。因此,研究开发有效的早期故障和复合故障诊断技术具有重要的工程应用价值。早期故障包含两方面含义,其一是指处于早期阶段的微弱故障或潜在故障,具有症状不明显、特征信息微弱、信噪比低等特点;其二是从物理意义上讲,某一故障是另一故障的早期阶段,并随着时间推移进一步诱发复合故障。复合故障是指关键设备同时发生两个或多个互相关联、交叉影响的故障,其动态信号一般表现为多个故障特征的相互耦合与调制。早期故障和复合故障的特征提取、故障定位与定量比较困难,已成为国内外本领域研究的难点和热点。
近年来兴起的第二代小波变换、多小波变换、双树复小波变换、谱峭度等方法在早期故障特征提取上展现出广阔的应用前景。但是这些方法中如何构造和选择与故障特征波形相匹配、且具有优良性质的合适基函数成为故障诊断的关键所在。同时,这些方法在进行抑制噪声的同时,可能会“过扼杀”信号中的重要微弱故障特征信息,造成故障特征信息的丢失。
经验模式分解(Empirical Mode Decomposition,EMD)方法可根据信号本身具有的时间尺度特征自适应地将其分解为一系列本征模式分量(intrinsic mode function,IMF)的线性和,已在设备早期故障诊断领域得到应用。但传统EMD方法对信号进行分解时,时常会出现由于某一或几个时间尺度的IMF不连续,产生模式混淆的现象,从而影响EMD方法进一步应用。为克服模式混淆问题,集成经验模式分解(Ensemble Empirical Mode Decomposition,EEMD)方法应运而生。该方法从噪声利用的角度,对信号添加白噪声来改善传统EMD方法中的模式混淆现象。然而EEMD应用效果依赖于添加的白噪声大小和总体平 均次数,目前尚缺乏可靠、通用的噪声协助准则。同时,对于信噪比相对较低的机械故障信号,EEMD中的平均过程仅能抵消所添加的白噪声,不能消除信号中的原有噪声从而提取信号中早期及复合故障特征。
发明内容
本发明的目的在于提供一种应用于机械早期及复合故障特征提取的集成噪声重构经验模式分解(ensemble noise-reconstructed EMD)的方法,该方法不仅可继承EMD方法在自适应分解非平稳信号上的优势,同时它利用原信号中固有噪声分量可实现信号降噪并改善EMD中的模式混淆问题。
为达到以上目的,本发明是采取如下技术方案予以实现的:
一种机械早期及复合故障的集成噪声重构经验模式分解方法,其特征在于,包括下述步骤:
a.首先对x(t)进行EMD分解,得到各IMF{ck(t),k=1,···,n}和残差分量rn(t),其中,k代表IMF个数;
b.计算各个{ck(t)}的能量{Ek,k=1,···n},设定第一个IMF分量c1(t)为纯白噪声,即c1(t)中的纯白噪声能量
e.将以上待处理的含噪IMF组成{cl(t),l=1,···};
4)重复步骤2)和步骤3)直至满足以下运算停止准则:
a.假设进行了m次EMD运算,则得到平均运算后的IMF分量和残差分量如下:
并重构得到降噪信号
若得到的相对噪声误差比例er小于给定的误差允许值ε,则停止运算,否则重复运算直至满足该误差,其中ε∈[0.01,0.1];
5)当满足停止准则时,对步骤3)中计算得到的所有IMF分量进行平均运算处理,得到最终平均IMF分量和平均残差分量:
上述方法中,所述对信号x(t)进行EMD分解包含以下步骤:
1)把x(t)作为待处理信号,确定其所有局部极值点,然后将所有极大值点和所有极小值点分别用三次样条曲线连接起来,得到x(t)的上、下包络线,使x(t)的所有数据点都处于这两条包络线之间;
2)从待处理信号x(t)中减去其上、下包络线均值m(t),得到h1(t)=x(t)-m(t);
3)检测h1(t)是否满足基本模式分量的两个条件,如果不满足,则把h1(t)作为待处理信号,重复上述操作,直至h1(t)是一个基本模式分量c1(t)=h1(t);
4)分解出第一个基本模式分量c1(t)之后,从x(t)中减去c1(t),得到剩余值序列r1(t)=x(t)-c1(t);
5)把r1(t)作为新的“原始”信号重复上述操作,依次可得第二、第三直至第n个基本模式分量,记为c1(t),c2(t),···,cn(t),这个处理过程在满足预先设定的停止准则后即可停止,最后剩下原始信号的余项rn(t)。
本发明继承了EMD方法在自适应分解非平稳信号上的优势,并具有以下优点:
1)与传统基于噪声抑制的降噪方法不同,本发明利用了原信号中固有噪声分量的相互抵消作用进行降噪,可避免“过扼杀”重要微弱故障特征信息的不足,从而有效识别故障信号中的早期及复合故障特征。
2)与传统EEMD方法不同,本发明不必再额外添加白噪声序列,仅利用原有噪声的重采样便可有效改善和消除模式混淆现象,避免了人为添加白噪声的干扰。
3)本发明识别精度高,运算实时性好,简单易行,非常适合于信噪比相对较低的机械早期及复合故障信号处理及特征提取。
附图说明
以下结合附图及具体实施方式对本发明作进一步的详细描述。
图1为实施本发明方法的仿真信号图。
图2为本发明方法的流程框图。
图5为图2信号中估计噪声和实际噪声的对比图。
图6为图2信号中估计噪声和实际噪声的概率密度分布。
图7为图2信号中采用本发明的分析结果。其中,(a)图为采用本发明分析得到的IMF分量;(b)图为IMF分量c2的放大图;(c)图为IMF分量c4的放大图。
图8为图2信号中采用EMD的分析结果。其中,(a)图为采用EMD分析得到的IMF分量;(b)图为IMF分量c2的放大图;(c)图为IMF分量c4的放大图。
图9为图2信号中采用EEMD的分析结果。其中,(a)为采用EEMD分 析得到的IMF分量;(b)图为IMF分量c2的放大图;(c)图为IMF分量c4的放大图。
图10为实例1齿轮早期故障中振动信号[(a)图]及其频谱[(b)图]。
图11为采用本发明对图10故障的分析结果。其中,(a)图为采用本发明分析得到的IMF分量;(b)图为IMF分量c2的放大图;(c)图为IMF分量c4的放大图。
图12为实例2齿轮复合故障中振动信号[(a)图]及其频谱[(b)图]。
图13为采用本发明对图12故障的分析结果。齿轮复合故障中采用本发明的分析结果。其中,(a)图为采用本发明分析得到的IMF分量;(b)图为IMF分量c2的放大图;(c)图为IMF分量c4的放大图。
具体实施方式
本发明具体方法包含以下内容:
首先对x(t)进行EMD分解,得到各IMF{ck(t),k=1,···,n}和残差分量rn(t),其中,k代表IMF个数;计算各个{ck(t)}的能量{Ek,k=1,···n},设定第一个IMF分量c1(t)为纯白噪声,即c1(t)中的纯白噪声能量
根据纯白噪声下IMF能量公式,分别计算各个{ck(t)}得到各IMF在置信区间95%的纯白噪声能量和置信区间99%下的纯白噪声能量 将各IMF的能量{Ek}与相应的纯白噪声能量对比:若 或 或 其中α∈[1,2],为一个给定容差值,则将ck(t)作为待处理的含噪IMF;否则,ck(t)作为估算的故障特征分量;
将以上待处理的含噪IMF组成{cl(t),l=1,···};
重复噪声分量估算和EMD分解,直至满足以下运算停止准则:
假设进行了m次EMD运算,则得到平均运算后的IMF分量和残差分量如下:
并重构得到降噪信号
若得到的相对噪声误差比例er小于给定的误差允许值ε,则停止运算,否则重复运算直至满足该误差,其中ε∈[0.01,0.1];
当满足停止准则时,对多次EMD分解得到的所有相应的IMF分量进行平均运算处理,得到最终平均IMF分量和平均残差分量:
对信号x(t)进行EMD分解包含以下步骤:
把x(t)作为待处理信号,确定其所有局部极值点,然后将所有极大值点和所有极小值点分别用三次样条曲线连接起来,得到x(t)的上、下包络线,使x(t)的所有数据点都处于这两条包络线之间;
从待处理信号x(t)中减去其上、下包络线均值m(t),得到h1(t)=x(t)-m(t);
检测h1(t)是否满足基本模式分量的两个条件,如果不满足,则把h1(t)作为待处理信号,重复上述操作,直至h1(t)是一个基本模式分量c1(t)=h1(t);
分解出第一个基本模式分量c1(t)之后,从x(t)中减去c1(t),得到剩余值序列r1(t)=x(t)-c1(t);
把r1(t)作为新的“原始”信号重复上述操作,依次可得第二、第三直至第n个 基本模式分量,记为c1(t),c2(t),···,cn(t),这个处理过程在满足预先设定的停止准则后即可停止,最后剩下原始信号的余项rn(t)。
为了验证本发明所述方法的正确性,给出两个具体实例进一步说明。
实例1
构造包含正弦分量(模拟旋转机械特征信号)、冲击分量(模拟局部故障特征信号)和噪声分量的混合信号x(t)为
x(t)=s1(t)+s2(t)+λn(t)
s1(t)=sin(160πt), f(t)=e-200|t|sin(1200πt)
其中,T=0.125s,n(t)为均值为0方差为1的高斯白噪声,且λ=0.35。对该仿真信号采用采样频率fs=5120Hz以及点数n=5120进行仿真,构造的仿真信号如图1所示。
对该仿真信号采用本发明进行特征提取,按照图2所示的流程进行分析。首先对仿真信号进行EMD分解,计算得到各IMF分量的{Ek,k=1,···n}, 和 如图3所示。根据 或者或者(其中,参数α=2),选择{cl(t),l=1,2,3,5,6,8,9}组成待处理的含噪IMF。对这些IMF进行反硬阈值处理(其中,参数C=0.5),如图4所示。图中,横线代表各IMF的阈值门限{Tl},深色曲线代表原信号{cl(t)},阈值门限内的浅色曲线代表反硬阈值处理后的信号 将以上阈值处理后的信号合并得到信号估算噪声分量如图5所示。图中,实线表示仿真信号中所估计的噪声分量虚线表示仿真实际噪声分量λn(t),可以看到,估计噪声分量和实际噪声分量的幅值大小基本一致。同时,图6为仿真信号中估计噪声分量(实线)和仿真实际噪声分量λn(t)(虚线)的概率密度分布,可以看到,两者的概论密度分布也非常相似,也证明了本发明中的噪声估计技术的有效性。对估算的噪声分量进行随机重排的重采样得到噪声分量然后重构相应的含噪信号并进行EMD分解。计算噪声重构及总体平均后的噪声分量和相对噪声误差比er。重复以上步骤直至er小于给定的误差允许值ε(参数ε=0.1)。最后平均所有IMF分量,输出的分析结果如图7所示。从图7中可以看出,IMF分量c2和c4分别对应于仿真分量s2(t)和s1(t),且c2中的周期性冲击序列非常明显。采用峭度指标ks来衡量本发明对于冲击序列的提取效果,并且采用s1(t)和c4的相 关系数e来评价正弦信号的完整性,也即评估改善模式混淆现象的效果。本发明中,c2的ks=17.9129,且c4的e=0.9699。
对该仿真信号采用EMD和EEMD方法进行分析,分别如图8和图9所示。从图8中可以看到,周期性冲击序列在c2中不易识别,且在c4和c5中发生了模式混淆现象(如图中双箭头所示)。从图9中可以看到,虽然EEMD可以改善模式混淆现象,但是c2中提取的周期性冲击特征远不如本发明中的分析结果。同时,在EMD中峭度指标ks=5.1656且相关系数e=0.9203,而在EEMD中ks=5.9622且e=0.9600,这些指标都不如本发明的指标值。因此,本实施案例中冲击序列分析说明本发明可以很好地提高信号的信噪比,提取与识别微弱特征信息,并且本案例中的正弦分量分析说明本方法可以明显改善EMD方法中的模式混淆现象。
实施案例2:
在齿轮试验台上采用齿轮试验来验证本发明对早期故障和复合故障特征提取能力。其中,大齿轮齿数223,小齿轮齿数48,其中小齿轮有一个轮齿被切除用于计数。采样频率设定为25kHz,大齿轮转速为145r/min,小齿轮转速为659r/min。
图10为齿轮早期裂纹故障的振动信号及其频谱。按照图2所示的流程,采用本发明对齿轮早期故障信号进行分析,输出的分析结果如图11所示。从图11中可以看出,在c2和c4中分别出现明显的周期性冲击特征。c4中周期性冲击是由小齿轮中缺少的一个轮齿所引起。而c2中周期性冲击间隔为0.4144s,接近于大齿轮的旋转频率,说明该大齿轮存在局部故障,如齿根裂纹或疲劳剥落,与齿轮实际故障状态相符。
实施案例3:
图12为上述齿轮试验中齿轮存在2个断齿时的振动信号及其频谱。同样按照图2所示的流程,采用本发明对齿轮复合故障信号进行分析,输出的分析结果如图13所示。其中,c4中周期性冲击特征依旧对应于小齿轮缺少的轮齿所产生的振动特征。而c2呈现为若干包含双周期性冲击特征信号,与大齿轮存在两个断齿故障相符。同时,两组周期性冲击的间隔为0.0483s,说明两个断齿在圆周上相隔26个轮齿,与实际相符。
Claims (2)
1.一种机械早期及复合故障的集成噪声重构经验模式分解方法,其特征在于,包括如下步骤:
a.首先对x(t)进行EMD分解,得到各IMF分量{ck(t),k=1,…,n}和残差分量rn(t),其中,k代表IMF个数;
e.将以上待处理的含噪IMF组成{cl(t),l=1,…};
4)重复步骤2)和步骤3)直至满足以下运算停止准则:
a.假设进行了m次EMD运算,则得到平均运算后的IMF分量和残差分量如下:
若得到的相对噪声误差比例er小于给定的误差允许值ε,则停止运算,否则重复运算直至满足该误差,其中ε∈[0.01,0.1];
5)当满足停止准则时,对步骤3)中计算得到的所有IMF分量进行平均运算处理,得到最终平均IMF分量和平均残差分量:
2.如权利要求1所述的机械早期及复合故障的集成噪声重构经验模式分解方法,其特征在于,所述对信号x(t)进行EMD分解包含以下步骤:
1)把x(t)作为待处理信号,确定其所有局部极值点,然后将所有极大值点和所有极小值点分别用三次样条曲线连接起来,得到x(t)的上、下包络线,使x(t)的所有数据点都处于这两条包络线之间;
2)从待处理信号x(t)中减去其上、下包络线均值m(t),得到h1(t)=x(t)-m(t);
3)检测h1(t)是否满足基本模式分量的两个条件,如果不满足,则把h1(t)作为待处理信号,重复上述操作,直至h1(t)是一个基本模式分量c1(t)=h1(t);
4)分解出第一个基本模式分量c1(t)之后,从x(t)中减去c1(t),得到剩余值序列r1(t)=x(t)-c1(t);
5)把r1(t)作为新的“原始”信号重复上述操作,依次可得第二、第三直至第n个基本模式分量,记为c1(t),c2(t),…,cn(t),这个处理过程在满足预先设定的停止准则后即可停止,最后剩下原始信号的余项rn(t)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310097502.6A CN103226649B (zh) | 2013-03-25 | 2013-03-25 | 机械早期及复合故障的集成噪声重构经验模式分解方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310097502.6A CN103226649B (zh) | 2013-03-25 | 2013-03-25 | 机械早期及复合故障的集成噪声重构经验模式分解方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103226649A true CN103226649A (zh) | 2013-07-31 |
CN103226649B CN103226649B (zh) | 2016-08-17 |
Family
ID=48837094
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310097502.6A Expired - Fee Related CN103226649B (zh) | 2013-03-25 | 2013-03-25 | 机械早期及复合故障的集成噪声重构经验模式分解方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103226649B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103954443A (zh) * | 2014-04-29 | 2014-07-30 | 华电电力科学研究院 | 自适应总体平均经验模式分解eemd协助噪声大小确定方法 |
CN104748961A (zh) * | 2015-03-30 | 2015-07-01 | 中国矿业大学 | 基于svd分解降噪和相关性eemd熵特征的齿轮故障诊断方法 |
CN105092239A (zh) * | 2014-05-09 | 2015-11-25 | 潍坊学院 | 一种齿轮早期故障检测方法 |
CN105698922A (zh) * | 2016-02-04 | 2016-06-22 | 国网福建省电力有限公司 | 基于改进emd和谱峭度法变压器振动故障特征提取方法 |
CN105973999A (zh) * | 2016-04-28 | 2016-09-28 | 西安交通大学 | 基于增强相位瀑布图的转子裂纹微弱分数谐波特征识别方法 |
CN106096201A (zh) * | 2016-06-29 | 2016-11-09 | 潍坊学院 | 一种旋转机械的eemd和三次样条平滑包络分析方法 |
CN106525435A (zh) * | 2016-10-25 | 2017-03-22 | 天津大学 | 一种发动机噪声源识别方法 |
CN108844617A (zh) * | 2018-05-03 | 2018-11-20 | 武汉理工大学 | 主减速器强噪声振动信号的自适应稀疏树结构降噪方法 |
CN110146292A (zh) * | 2019-06-04 | 2019-08-20 | 昆明理工大学 | 一种基于集成噪声重构的总体局部均值分解的滚动轴承故障检测方法 |
CN110716228A (zh) * | 2019-11-18 | 2020-01-21 | 上海理工大学 | 用于回转式多余物检测特征信号提取方法 |
CN113192489A (zh) * | 2021-05-16 | 2021-07-30 | 金陵科技学院 | 一种基于多尺度增强BiLSTM模型的喷漆机器人语音识别方法 |
CN117167903A (zh) * | 2023-11-03 | 2023-12-05 | 江苏中安建设集团有限公司 | 基于人工智能的暖通设备异物故障检测方法 |
CN117540328A (zh) * | 2024-01-09 | 2024-02-09 | 山西众诚安信安全科技有限公司 | 一种煤矿噪声高精度测量过程中噪声处理方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6311130B1 (en) * | 1996-08-12 | 2001-10-30 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Computer implemented empirical mode decomposition method, apparatus, and article of manufacture for two-dimensional signals |
CN102254103A (zh) * | 2011-07-27 | 2011-11-23 | 西安交通大学 | 自适应总体平均经验模式分解eemd筛选次数确定方法 |
CN102620928A (zh) * | 2012-03-02 | 2012-08-01 | 燕山大学 | 基于小波半软阈值和emd的风电齿轮箱故障诊断方法 |
CN102854015A (zh) * | 2012-10-15 | 2013-01-02 | 哈尔滨理工大学 | 一种滚动轴承故障位置及性能退化程度诊断方法 |
-
2013
- 2013-03-25 CN CN201310097502.6A patent/CN103226649B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6311130B1 (en) * | 1996-08-12 | 2001-10-30 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Computer implemented empirical mode decomposition method, apparatus, and article of manufacture for two-dimensional signals |
CN102254103A (zh) * | 2011-07-27 | 2011-11-23 | 西安交通大学 | 自适应总体平均经验模式分解eemd筛选次数确定方法 |
CN102620928A (zh) * | 2012-03-02 | 2012-08-01 | 燕山大学 | 基于小波半软阈值和emd的风电齿轮箱故障诊断方法 |
CN102854015A (zh) * | 2012-10-15 | 2013-01-02 | 哈尔滨理工大学 | 一种滚动轴承故障位置及性能退化程度诊断方法 |
Non-Patent Citations (2)
Title |
---|
陈略等: "噪声协助的EMD-1.5维谱信号抗混分解与特征提取", 《振动与冲击》 * |
陈略等: "总体平均经验模式分解与1.5维谱方法的研究", 《西安交通大学学报》 * |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103954443B (zh) * | 2014-04-29 | 2017-01-04 | 华电电力科学研究院 | 自适应总体平均经验模式分解eemd协助噪声大小确定方法 |
CN103954443A (zh) * | 2014-04-29 | 2014-07-30 | 华电电力科学研究院 | 自适应总体平均经验模式分解eemd协助噪声大小确定方法 |
CN105092239A (zh) * | 2014-05-09 | 2015-11-25 | 潍坊学院 | 一种齿轮早期故障检测方法 |
CN105092239B (zh) * | 2014-05-09 | 2017-06-20 | 潍坊学院 | 一种齿轮早期故障检测方法 |
CN104748961A (zh) * | 2015-03-30 | 2015-07-01 | 中国矿业大学 | 基于svd分解降噪和相关性eemd熵特征的齿轮故障诊断方法 |
CN105698922A (zh) * | 2016-02-04 | 2016-06-22 | 国网福建省电力有限公司 | 基于改进emd和谱峭度法变压器振动故障特征提取方法 |
CN105973999A (zh) * | 2016-04-28 | 2016-09-28 | 西安交通大学 | 基于增强相位瀑布图的转子裂纹微弱分数谐波特征识别方法 |
CN105973999B (zh) * | 2016-04-28 | 2018-10-30 | 西安交通大学 | 基于增强相位瀑布图的转子裂纹微弱分数谐波特征识别方法 |
CN106096201B (zh) * | 2016-06-29 | 2019-02-22 | 潍坊学院 | 一种旋转机械的eemd和三次样条平滑包络分析方法 |
CN106096201A (zh) * | 2016-06-29 | 2016-11-09 | 潍坊学院 | 一种旋转机械的eemd和三次样条平滑包络分析方法 |
CN106525435A (zh) * | 2016-10-25 | 2017-03-22 | 天津大学 | 一种发动机噪声源识别方法 |
CN108844617A (zh) * | 2018-05-03 | 2018-11-20 | 武汉理工大学 | 主减速器强噪声振动信号的自适应稀疏树结构降噪方法 |
CN110146292A (zh) * | 2019-06-04 | 2019-08-20 | 昆明理工大学 | 一种基于集成噪声重构的总体局部均值分解的滚动轴承故障检测方法 |
CN110146292B (zh) * | 2019-06-04 | 2021-08-31 | 昆明理工大学 | 一种基于集成噪声重构的总体局部均值分解的滚动轴承故障检测方法 |
CN110716228A (zh) * | 2019-11-18 | 2020-01-21 | 上海理工大学 | 用于回转式多余物检测特征信号提取方法 |
CN113192489A (zh) * | 2021-05-16 | 2021-07-30 | 金陵科技学院 | 一种基于多尺度增强BiLSTM模型的喷漆机器人语音识别方法 |
CN117167903A (zh) * | 2023-11-03 | 2023-12-05 | 江苏中安建设集团有限公司 | 基于人工智能的暖通设备异物故障检测方法 |
CN117167903B (zh) * | 2023-11-03 | 2024-01-30 | 江苏中安建设集团有限公司 | 基于人工智能的暖通设备异物故障检测方法 |
CN117540328A (zh) * | 2024-01-09 | 2024-02-09 | 山西众诚安信安全科技有限公司 | 一种煤矿噪声高精度测量过程中噪声处理方法 |
CN117540328B (zh) * | 2024-01-09 | 2024-04-02 | 山西众诚安信安全科技有限公司 | 一种煤矿噪声高精度测量过程中噪声处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103226649B (zh) | 2016-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103226649A (zh) | 机械早期及复合故障的集成噪声重构经验模式分解方法 | |
CN102721545B (zh) | 一种基于多特征参量的滚动轴承故障诊断方法 | |
CN102866027A (zh) | 基于lmd和局域时频熵的旋转机械故障特征提取方法 | |
Zhang et al. | An improved filtering method based on EEMD and wavelet-threshold for modal parameter identification of hydraulic structure | |
CN102636577B (zh) | 一种基于小波变换提取轮辋超声波探伤信号的方法 | |
CN101900789B (zh) | 基于小波变换与分形维数的容差模拟电路故障诊断方法 | |
CN105928702B (zh) | 基于形态分量分析的变工况齿轮箱轴承故障诊断方法 | |
CN104748961A (zh) | 基于svd分解降噪和相关性eemd熵特征的齿轮故障诊断方法 | |
CN102288843B (zh) | 一种电能质量扰动信号检测方法 | |
CN102620928A (zh) | 基于小波半软阈值和emd的风电齿轮箱故障诊断方法 | |
CN104729853A (zh) | 一种滚动轴承性能退化评估装置及方法 | |
CN104390781A (zh) | 一种基于lmd和bp神经网络的齿轮故障诊断方法 | |
CN104236905A (zh) | 一种轴承故障诊断方法 | |
CN105424359A (zh) | 一种基于稀疏分解的齿轮和轴承混合故障特征提取方法 | |
CN105469611A (zh) | 一种短时交通流预测模型方法 | |
CN102436586A (zh) | 一种基于小波阈值降噪和经验模态分解的高光谱图像分类方法 | |
CN105043767A (zh) | 一种滚动球轴承外圈剥落故障双冲击特征提取方法及系统 | |
CN104215459A (zh) | 一种轴承故障的诊断方法 | |
CN104949840A (zh) | 一种基于振动分析的柴油机故障诊断方法 | |
CN113390631A (zh) | 一种柴油机齿轮箱故障诊断方法 | |
CN106771598A (zh) | 一种自适应谱峭度信号处理方法 | |
Lin et al. | Application of the EEMD method to multiple faults diagnosis of gearbox | |
Lv et al. | Generalized synchroextracting-based stepwise demodulation transform and its application to fault diagnosis of rotating machinery | |
CN104089778B (zh) | 一种水轮机振动故障诊断方法 | |
Yu et al. | Sparse time–frequency representation for the transient signal based on low-rank and sparse decomposition |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160817 Termination date: 20190325 |