基于DMD降噪的滚动轴承故障诊断方法.pdf

第4 6卷第5期武汉科技大学学报V o l.4 6,N o.52 0 2 3年1 0月J o u r n a l o fW u h a nU n i v e r s i t yo fS c i e n c ea n dT e c h n o l o g yO c t.2 0 2 3 收稿日期:2 0 2 2-1 0-0 3 基金项目:国家自然科学基金资助项目(5 1 9 7 5 4 2 7).作者简介:涂福泉(1 9 7 0-),男,武汉科技大学教授,博士.E-m a i l:t u f u q u a n w u s t.e d u.c nD O I:1 0.3 9 6 9/j.i s s n.1 6 7 4-3 6 4 4.2 0 2 3.0 5.0 0 8基于DMD降噪的滚动轴承故障诊断方法涂福泉,杨家瑜,陈 超,罗迎九,吴维崧(武汉科技大学冶金装备及其控制教育部重点实验室,湖北 武汉,4 3 0 0 8 1)摘要:针对轴承振动信号难以剔除噪声的问题,提出一种将H i l b e r t变换、动力学模态分解(DMD)和奇异值差分谱相结合的滚动轴承故障诊断方法。首先将原始信号进行H i l b e r t变换得到包络信号,由包络信号构造H a n k e l矩阵进行动力学模态分解,利用奇异值差分谱确定合适的截断秩后进行信号重构,最后通过频谱分析来提取故障特征。采用该方法对滚动轴承故障仿真信号和实验数据进行分析,结果表明其降噪效果显著,能有效获取轴承故障特征频率。关键词:故障诊断;滚动轴承;动力学模态分解;H i l b e r t变换;奇异值差分谱;信号降噪中图分类号:TH 1 3 3.3 3 文献标志码:A 文章编号:1 6 7 4-3 6 4 4(2 0 2 3)0 5-0 3 7 6-0 8 滚动轴承是机械传动中的重要组成部件,广泛应用于各种机械设备,如风电机组、高速列车、航空发动机等。由于工作环境的复杂性,滚动轴承容易出现点蚀、断裂、磨损等故障,从而影响设备运行的稳定性和可靠性,因此对滚动轴承进行故障诊断具有重要意义。在设备故障诊断时,通常是利用传感器收集的振动信号进行分析。这些振动信号中有大量的噪声,导致其故障特征不明显,难以直接应用,需要将原始信号进行降噪处理,再进行下一步的分析1。经典的信号降噪和特征提取方法包括小波变换(WT)、经验模态分解(EMD)、变分模态分解(VMD)、奇异值分解(S V D)等,然而这些方法都有一些不足之处。例如,WT需要选择合适的母小波及分解层数;EMD公式缺乏严格的理论推导,并且存在模态混叠和端点效应等问题;VMD的参数选择对其分解效果影响较大;S V D对强噪声背景下的信号特征提取效果欠佳。因此,研究人员以上述方法为基础提出了一系列的改进措施。Y a n g等2针对传统的采用单一准则优化的小波参数无法完全描述故障特征的问题,设计了一种多目标优化小波滤波器用于提取轴承故障特征,在强噪声背景下的应用效果较好。滚动轴承早期故障信号存在大量噪声使得故障特征提取十分困难,为此曹玲玲等3提出一种将互补集合经验模态分解(C E EMD)和改进小波阈值算法相结合的方法,可以提高故障信号的信噪比,降噪效果明显。Y i n等4提出一种滚动轴承弱故障特征提取方法,能选择出具有故障特征的内模函数,并且通过新构造的自适应阈值进行信号去噪。赵玮5采用基于VMD和奇异值差分谱的故障特征提取方法,可 以 有 效 降 低 噪 声 的 影 响。鉴 于 传 统VMD方法对高速列车轮对轴承故障特征提取不够准确,李翠省等6设计一种结合集成经验模态分解(E EMD)和改进参数自适应VMD的方法,其在强噪声干扰下仍具有优良的故障特征提取效果。李华等7为了提升S V D方法的降噪性能,提出一种基于相关奇异值比的S V D方法并成功应用于轴承故障诊断。C h e n g等8提出自适应周期模态分解方法,以解决现有滚动轴承故障诊断方法对周期性冲击故障识别能力差的问题。动力学模态分解(d y n a m i cm o d ed e c o m p o s i-t i o n,DMD)是基于K o o p m a n算子的无方程数据驱动分析方法,它不仅有S V D算法的优点,还考虑了机械系统演化过程中的时空特性,具有时空耦合建模的独特优势,在信号降噪和特征提取方面已得到有效应用9。DMD本质上是一种融合奇异值分解和模态分解理论的空间降维方法,因此它也不可避免地继承了两者的缺点,例如在截断秩的选取和含调制成分信号的处理方面存在困难。针对上述问题,本文提出一种H i l b e r t变换和DMD相结合并通过奇异值差分谱确定截断秩的信号降噪和特征提取方法,用于滚动轴承故障诊断。首先使用H i l b e r t变换对原始信号进行包2 0 2 3年第5期涂福泉,等:基于DMD降噪的滚动轴承故障诊断方法络解调,消除原信号中调制成分对后续分析的影响,然后进行动力学模态分解,通过奇异值差分谱选择合适的截断秩进行信号重构,从而减少噪声的干扰,以达到较好的特征提取效果。本文最后利用仿真信号和美国凯斯西储大学公开的轴承故障数据集进行方法验证。1 基本理论1.1 H i l b e r t变换对一个连续信号x(t)进行H i l b e r t变换可表示为1 0 x(t)=Hx(t)=1n-x()t-d=x(t)1t(1)式中:Hx(t)表示时域中的H i l b e r t变换,为卷积算子。原信号的解析信号为y(t)=x(t)+jx(t)(2)从解析信号中可以得到原信号x(t)的包络信号为B(t)=x2(t)+x2(t)(3)1.2 DMD算法假设两个连续序列X、X(用矩阵形式表示)满足与最优线性算子的映射关系,即X=A X,DMD算法利用矩阵A的低秩表达式来捕捉连续序列潜在的动态特性1 1-1 2。由于DMD最初是应用于流体力学领域,所收集的信号通常是二维矩阵形式,而轴承振动信号是一维的,因此需要将一维振动信号转换为矩阵形 式,这 里 构 造H a n k e l矩 阵 以 满 足 上 述 要求1 3。对于振动信号X(t)=x(1),x(2),x(3),x(N),其H a n k e l矩阵形式为X=x(1)x(n-1)x(m)x(N-1),X=x(2)x(n)x(m+1)x(N)(4)式中:m=N/2,1nN且n=N-m+1。DMD算法步骤为:1)对X进行奇异值分解X=U V*(5)式中:U、V为酉矩阵,为奇异值对角矩阵。U的每个列 向 量 相 互 正 交,是X的 本 征 正 交 分 解(P O D)模态。取U的前r阶来截断数据,通过保留主要模态来降维。2)利用X的伪逆来求AA=X XT=X V-1U*(6)通常情况下,A包含的数据量很大,故通过合适的投影矩阵将其降维为r阶的矩阵A-:A-=U*rA Ur=U*rX V-1(7)3)计算矩阵A-的特征值和特征向量A-W=W(8)式中:W是矩阵A-的特征向量矩阵,是由矩阵A-的特征值i组成的对角矩阵。4)计算重构矩阵定义矩阵:=X V-1W(9)中每一个列向量i称为DMD的模态,并且每个模态对应着一个特征值i。最后得到由若干个模态组成的DMD重构矩阵:XDMD=ri=1ie x p(wit)bi(1 0)式中:wi=l n(i)/t,bi是每个模态的幅值。1.3 奇异值差分谱奇异值差分谱可以表征奇异值序列的变化,差分谱计算公式如下1 4:di=i-i+1ni=1(i-i+1),i=1,2,n(1 1)式中:i为奇异值。所有的di构成奇异值差分谱,差分谱的峰值与信号特征相关,峰值越大,特征越明显,当峰值降低到一定程度便可截断。2 本文方法本文方法的基本流程如图1所示。首先将原始信号进行H i l b e r t变换得到包络信号,再根据包络信号构造H a n k e l矩阵进行动力学模态分解,然后通过奇异值差分谱确定DMD的截断秩,使用该截断秩进行重构得到重构矩阵,其第一行图1 本文方法的基本流程F i g.1F l o w c h a r t o f t h ep r o p o s e dm e t h o d773武汉科技大学学报2 0 2 3年第5期为重构信号,根据重构信号进行频谱分析,提取出准确的信号特征,从而实现滚动轴承故障的精确诊断。3 仿真信号分析下面采用滚动轴承内圈故障仿真信号来验证本文方法的有效性。仿真信号表达式如下1 5:x(t)=s(t)+n(t)=iA(t)h(t-i T)+n(t)h(t)=e x p(-C t)c o s(2 fnt)A(t)=1+A0c o s(2 frt)(1 2)式中:fr为转频,fn为轴承的共振频率,T为故障周期,C为衰减系数,A0为幅值,n(t)为高斯白噪声。设定采样点数N=40 9 6,高斯白噪声的信噪比为-5d B,仿真信号其余参数的具体取值如表1所示,其中fs为系统的采样频率。表1 滚动轴承内圈故障仿真信号的参数T a b l e 1P a r a m e t e r s o f t h e s i m u l a t e d s i g n a l o f r o l l i n gb e a r i n gw i t h i n n e r-r a c e f a u l t参数A0fr/H zCT/sfs/k H zfn/k H z取值0.3 52 57 5 01/1 2 01 23 图2为仿真信号的时域波形图。由图2(b)可见,时域信号的故障特征已经被噪声淹没,从中无法获得轴承故障频率。为了分析信号的故障特征,对加噪声后的仿真信号进行包络谱分析,结果如图3所示。在图3中可以找到信号的特征频率和部分多倍频,但是包含大量的噪声信息。(a)未加噪声的仿真信号(b)加噪声的仿真信号图2 仿真信号的时域波形图F i g.2T i m e-d o m a i nw a v e f o r m so f t h e s i m u l a t e ds i g n a l s图3 加噪声仿真信号的包络谱F i g.3E n v e l o p e s p e c t r u mo f t h e s i m u l a t e ds i g n a lw i t hn o i s e 为了更加有效地提取轴承的故障特征,使用本文方法进行信号分析。首先将原始信号通过H i l b e r t变换构造解析信号,对解析信号取模得到包络信号,如图4所示。8732 0 2 3年第5期涂福泉,等:基于DMD降噪的滚动轴承故障诊断方法图4 本文方法得到的包络信号F i g.4E n v e l o p e s i g n a l b y t h ep r o p o s e dm e t h o d 根据包络信号构造H a n k e l矩阵,进行动力学模态分解,通过奇异值差分谱确定截断秩的阶数,如图5所示。需要在奇异值分布曲线上选取一个合适的截断秩阶数,但是仅根据图5(a)难以确定最优点,因此可以根据奇异值差分谱进行分析。在理想状态下,差分谱峰值的大小与信号特征相关,峰值越大,特征越明显,当差分谱峰值降低到一定程度而难以区分特征信息时便可截断。如图5(b)所示,通过奇异值差分谱分析,设定奇异值差分阈值为0.0 0 1,取最后一个峰值大于阈值的差分阶数作为截断秩阶数,最终确定的截断秩阶数为3 1。(a)奇异值分布曲线(b)奇异值差分谱图5 仿真信号的奇异值分布及奇异值差分谱F i g.5S i n g u l a rv a l u ed i s t r i b u t i o na n ds i n g u l a rv a l u ed i f f e r e n c e s p e c t r u mo f t h e s i m u l a t e ds i g n a l 根据设定的截断秩阶数,可以得到3 1个精确DMD模态,然后构造DMD重构矩阵,选取其中第一行为重构信号。通过DMD重构后的信号如图6(a)所示,从中可以观察到故障特征。利用快速傅里叶变换(F F T)对重构信号进行频谱分析,结果如图6(b)所示,可以观察到特征频率及其多倍频,谱线突出,特征频率附近存在一对边频,4倍频附近存在一个上边频。上述结果表明,本文方法能有效降低轴承振动信号中的噪声干扰,可提取出轴承故障特征频率。(a)时域波形(b)频率谱图6 本文方法对仿真信号的分析结果F i g.6A n a l y s i s r e s u l t so f t h e s i m u l a t e ds i g n a l b y t h ep r o p o s e dm e t h o d 图7所示为对比方法(包括S V D、自适应DMD、VMD、E EMD)对仿真 信 号 的 分 析 结 果。S V D重构信号的包络谱中特征频率和部分多倍频谱线明显,但是存在部分噪声;自适应DMD重973武汉科技大学学报2 0 2 3年第5期构信号的包络谱中特征频率和部分多倍频谱线突出,但高于5倍频的特征不明显,部分特征信息丢失;选 取VMD重 构 信 号 中 合 适 的 模 态 分 量(I MF 3)进行分析,其包络谱中故障频率、2倍频和3倍频的谱线较明显,但伴随有部分噪声;选取E EMD重构信号中合适的模态分量(I MF 1)进行分析,其包络谱中故障频率、部分多倍频的谱线较明显,但噪声较多。(a)S V D重构信号的包络谱(b)自适应DMD重构信号的包络谱(c)VMD重构信号I MF 3分量的包络谱(d)E EMD重构信号I MF 1分量的包络谱图7 其他方法对仿真信号的分析结果F i g.7A n a l y s i s r e s u l t so f t h e s i m u l a t e ds i g n a l b yo t h e rm e t h o d s4 实验验证为了进一步验证本文方法在实际应用中的有效性,以美国凯斯西储大学公开的轴承故障数据集中的内圈故障数据为案例进行分析。轴承故障直径为0.0 0 7英寸(0.0 1 7 7 8c m),转速为1 7 9 7r/m i n,采样频率为1 2k H z。通过理论计算,轴承内圈故障特征频率为1 6 2.1 9H z。由轴承内圈故障信号的时域波形(图8)可知,故障特征被噪声淹没,无法提取出有效的特征信息。该信号的包络谱分析结果见图9,从包络谱中可以提取轴承故障特征,但是还存在噪声干扰。为了更好地提取故障特征信息,使用本文方法对原始信号进行降噪处理。首先将原始信号通过H i l b e r t变换得到包络信号,如图1 0所示。根据包络信号构造H a n k e l矩阵后进行动力学模态分解,并通过奇异值差分谱确定截断秩的阶数。轴承内圈故障信号的奇异值分布及奇异值差分谱如图1 1所示,据此选取合适的截断秩阶数。从图1 1(a)来看,应选取曲线变化趋于稳定的一个点。设定奇异值差分阈值为0.0 0 1,结合图1 1(b),确定截断秩阶数为1 0 5。图8 轴承内圈故障信号的时域波形图F i g.8T i m e-d o m a i nw a v e f o r mo fb e a r i n g i n n e r-r a c e f a u l t s i g n a l0832 0 2 3年第5期涂福泉,等:基于DMD降噪的滚动轴承故障诊断方法图9 轴承内圈故障信号的包络谱F i g.9E n v e l o p e s p e c t r u mo fb e a r i n g i n n e r-r a c e f a u l t s i g n a l图1 0 轴承内圈故障包络信号的时域波形图F i g.1 0T i m e-d o m a i nw a v e f o r mo f t h ee n v e l o p e s i g n a l o fb e a r i n g i n n e r-r a c e f a u l t(a)奇异值分布曲线(b)奇异值差分谱图1 1 轴承内圈故障信号的奇异值分布及奇异值差分谱F i g.1 1S i n g u l a rv a l u ed i s t r i b u t i o na n ds i n g u l a rv a l u ed i f f e r e n c e s p e c t r u mo fb e a r i n g i n n e r-r a c e f a u l t s i g n a l 截断秩为1 0 5时的DMD重构信号如图1 2(a)所示,通过F F T得到的频率谱如图1 2(b)所示。将图1 2与图8图9对比可知,尽管从原始信号的包络谱中可以找到特征频率和部分多倍频,但其中伴随着大量的噪声,而经过H i l b e r t变换和DMD处理后,重构信号的频率谱中噪声明(a)时域波形图(b)频率谱图1 2 本文方法对轴承内圈故障信号的分析结果F i g.1 2A n a l y s i s r e s u l t so fb e a r i n g i n n e r-r a c e f a u l t s i g n a l b y t h ep r o p o s e dm e t h o d183武汉科技大学学报2 0 2 3年第5期显减少,比较容易找到特征频率和多倍频。由此可见,采用本文提出的将H i l b e r t变换与DMD相结合、再通过奇异值差分谱确定截断秩的方法可有效降低振动信号噪声,便于提取轴承故障特征。另外,在采用i 5-8 3 0 0 H处理器、1 6G B内存、G T X1 0 5 0T i显卡的实验配置下,本文方法对上述轴承内圈故障信号的计算时间为9.6 6 3 1 2s,计算时间较短。因此,当设备性能较好、采样点数量较少时可通过该方法进行轴承故障信号实时分析。图1 3所示为对比方法对轴承内圈故障信号的分析结果。从S V D重构信号的包络谱中可以分辨出故障频率和部分多倍频,但是包络谱中噪声较多;从自适应DMD重构信号的包络谱中可以看出部分多倍频的谱线较为突出,但是噪声干扰较大,部分倍频被噪声淹没;选取VMD重构信号合适的模态分量(I MF 3)进行分析,从包络谱中只能观察到特征频率和2倍频,并伴随有噪声影响;对E EMD重构信号的I MF 1分量进行分析,其包络谱中故障频率比较明显,可以观察到部分多 倍频,但是噪声也 较多。相对 而言,VMD和E EMD重构信号中特征丢失和噪声剔除不理想的情况更严重。(a)S V D重构信号的包络谱(b)自适应DMD重构信号的包络谱(c)VMD重构信号I MF 3分量的包络谱(d)E EMD重构信号I MF 1分量的包络谱图1 3 其他方法对轴承内圈故障信号的分析结果F i g.1 3A n a l y s i s r e s u l t so fb e a r i n g i n n e r-r a c e f a u l t s i g n a l b yo t h e rm e t h o d s 从计算效率角度分析,在实验配置不变的条件下,S V D方法的计算时间为2 8.5 3 6 4 8s,计算时间 较 长;自 适 应DMD方 法 的 计 算 时 间 为1 1.6 3 3 3 8s,计算时间稍长。与S V D和自 适应DMD相比,本文方法计算速度更快、降噪性能更优秀。5 结语针对轴承故障诊断中提取故障特征时难以剔除信号噪声的问题,本文提出了一种将H i l b e r t变换、DMD和奇异值差分谱相结合的方法。首先对原始信号进行H i l b e r t变换求得包络信号,以分离出原信号中的调幅和载波信息,去除原信号中高频调制部分对后续分析的影响。根据包络信号构造H a n k e l矩阵,再对该矩阵进行DMD处理。鉴于DMD算法的截断秩选取较为困难,本文通过奇异值差分谱来确定截断秩,当差分值降低到一定程度难以区分特征信息时便可以截断。取DMD重构矩阵的第一行 为重构信号,通过F F T得到重构信号的频域谱,从而提取出轴承故障特征。本文方法在对仿真信号和实际的滚动轴承故障信号的分析中均取得了较好的降噪效果,为滚动轴承故障诊断提供了一个有效手段。参考文献1 董绍江,裴雪武,吴文亮,等.基于多层降噪技术及改进卷积神经网络的滚动轴承故障诊断方法J.机械工程学报,2 0 2 1,5 7(1):1 4 8-1 5 6.2832 0 2 3年第5期涂福泉,等:基于DMD降噪的滚动轴承故障诊断方法2 Y a n gSP,G uX H,L i uY Q,e ta l.Ag e n e r a l m u l t i-o b j e c t i v eo p t i m i z e dw a v e l e t f i l t e ra n d i t sa p p-l i c a t i o n s i nf a u l td i a g n o s i so fw h e e l s e tb e a r i n g sJ.M e c h a n i c a lS y s t e m sa n dS i g n a lP r o c e s s i n g,2 0 2 0,1 4 5:1 0 6 9 1 4.3 曹玲玲,李晶,彭镇,等.基于改进小波阈值降噪的滚动轴承故障诊断方法J.振动工程学报,2 0 2 2,3 5(2):4 5 4-4 6 3.4 Y i nC,W a n gYL,M aGC,e ta l.W e a kf a u l t f e a-t u r ee x t r a c t i o n o fr o l l i n g b e a r i n g sb a s e d o ni m-p r o v e de n s e m b l en o i s e-r e c o n s t r u c t e d EMD a n da-d a p t i v e t h r e s h o l dd e n o i s i n gJ.M e c h a n i c a lS y s t e m sa n dS i g n a lP r o c e s s i n g,2 0 2 2,1 7 1:1 0 8 8 3 4.5 赵玮.基于VMD和奇异差分谱的滚动轴承早期故障诊断J.机械设计与制造,2 0 1 9(1 2):2 3 0-2 3 4.6 李翠省,廖英英,刘永强.基于E EMD和参数自适应VMD的高速列车轮对轴承故障诊断J.振动与冲击,2 0 2 2,4 1(1):6 8-7 7.7 李华,刘韬,伍星,等.相关奇异值比的S V D在轴承故障诊断中的应用J.机械工程学报,2 0 2 1,5 7(2 1):1 3 8-1 4 9.8 C h e n gJ,Y a n gY,L iX,e ta l.A d a p t i v ep e r i o d i cm o d ed e c o m p o s i t i o na n di t sa p p l i c a t i o ni nr o l l i n gb e a r i n gf a u l td i a g n o s i sJ.M e c h a n i c a lS y s t e m sa n dS i g n a lP r o c e s s i n g,2 0 2 1,1 6 1:1 0 7 9 4 3.9 D a n gZ,L vY,L iY Y,e ta l.O p t i m i z e dd y n a m i cm o d ed e c o m p o s i t i o nv i an o n-c o n v e xr e g u l a r i z a t i o na n d m u l t i s c a l ep e r m u t a t i o ne n t r o p yJ.E n t r o p y,2 0 1 8,2 0(3):1 5 2.1 0郭田雨,严荣国,方旭晨,等.基于希尔伯特变换和自适应阈值的R波检测算法J.计算机与现代化,2 0 2 2(2):1 1 4-1 1 9.1 1K u t zJN,B r u n t o nSL,B r u n t o nB W,e ta l.D y-n a m i cm o d ed e c o m p o s i t i o n:d a t a-d r i v e nm o d e l i n go fc o m p l e xs y s t e m sM.P h i l a d e l p h i a:S o c i e t yf o rI n-d u s t r i a l a n dA p p l i e dM a t h e m a t i c s,2 0 1 6.1 2A k a nA,C u r aOK.T i m e-f r e q u e n c ys i g n a l p r o c e s s-i n g:t o d a ya n d f u t u r eJ.D i g i t a l S i g n a lP r o c e s s i n g,2 0 2 1,1 1 9:1 0 3 2 1 6.1 3刘佳音,于晓光,王琦,等.基于H a n k e l矩阵与奇异值分解降噪方法的齿轮故障诊断研究J.机床与液压,2 0 1 8,4 6(1):1 5 8-1 6 2,1 7 2.1 4赵洪山,郭双伟,高夺.基于奇异值分解和变分模态分解 的 轴 承 故 障 特 征 提 取 J.振 动 与 冲 击,2 0 1 6,3 5(2 2):1 8 3-1 8 8.1 5唐贵基,王晓龙.I VMD融合奇异值差分谱的滚动轴承早期故障诊断J.振动、测试与诊断,2 0 1 6,3 6(4):7 0 0-7 0 7.F a u l td i a g n o s i s f o r r o l l i n gb e a r i n g sb a s e do nDMDd e n o i s i n gT uF u q u a n,Y a n gJ i a y u,C h e nC h a o,L u oY i n g j i u,W uW e i s o n g(K e yL a b o r a t o r yo fM e t a l l u r g i c a lE q u i p m e n t a n dC o n t r o l o fM i n i s t r yo fE d u c a t i o n,W u h a nU n i v e r s i t yo fS c i e n c ea n dT e c h n o l o g y,W u h a n4 3 0 0 8 1,C h i n a)A b s t r a c t:Ar o l l i n gb e a r i n gf a u l td i a g n o s i sm e t h o dc o m b i n i n gH i l b e r tt r a n s f o r m,d y n a m i cm o d ed e-c o m p o s i t i o n(DMD)a n ds i n g u l a rv a l u ed i f f e r e n c es p e c t r u m w a sp r o p o s e dt oe l i m i n a t et h en o i s ei nb e a r i n gv i b r a t i o ns i g n a l s.F i r s t l y,t h ee n v e l o p es i g n a lw a so b t a i n e db yH i l b e r t t r a n s f o r mo f t h eo r i g i-n a l s i g n a l,a n d t h eH a n k e lm a t r i xc o n s t r u c t e d f r o mt h ee n v e l o p e s i g n a lw a sp r o c e s s e db yDMD.T h e nt h es i g n a lw a sr e c o n s t r u c t e da f t e rd e t e r m i n i n gt h ea p p r o p r i a t et r u n c a t i o nr a n kv i as i n g u l a rv a l u ed i f f e r e n c es p e c t r u m.F i n a l l y,t h e f a u l t f e a t u r e sw e r e e x t r a c t e db ys p e c t r u ma n a l y s i s.T h i sm e t h o dw a su s e dt oa n a l y z et h es i m u l a t e ds i g n a la n dt h ee x p e r i m e n t a ld a t ao fr o l l i n gb e a r i n gf a u l t s.T h er e s u l t ss h o wt h a t t h ed e n o i s i n ge f f e c t i s s i g n i f i c a n t a n d t h eb e a r i n g f a u l t f e a t u r e f r e q u e n c i e s c a nb ee f f e c t i v e l ye x t r a c t e d.K e yw o r d s:f a u l td i a g n o s i s;r o l l i n gb e a r i n g;DMD;H i l b e r t t r a n s f o r m;s i n g u l a rv a l u ed i f f e r e n c es p e c-t r u m;s i g n a l d e n o i s i n g 责任编辑 尚 晶383