CN111089720B - 一种适用于变转速下齿轮故障诊断的正则化稀疏滤波方法 - Google Patents

一种适用于变转速下齿轮故障诊断的正则化稀疏滤波方法 Download PDF

Info

Publication number
CN111089720B
CN111089720B CN202010044852.6A CN202010044852A CN111089720B CN 111089720 B CN111089720 B CN 111089720B CN 202010044852 A CN202010044852 A CN 202010044852A CN 111089720 B CN111089720 B CN 111089720B
Authority
CN
China
Prior art keywords
sparse
training
samples
rotating speed
fault diagnosis
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.)
Active
Application number
CN202010044852.6A
Other languages
English (en)
Other versions
CN111089720A (zh
Inventor
韩宝坤
王金瑞
鲍怀谦
初振云
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University of Science and Technology
Original Assignee
Shandong University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shandong University of Science and Technology filed Critical Shandong University of Science and Technology
Priority to CN202010044852.6A priority Critical patent/CN111089720B/zh
Publication of CN111089720A publication Critical patent/CN111089720A/zh
Application granted granted Critical
Publication of CN111089720B publication Critical patent/CN111089720B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/02Gearings; Transmission mechanisms
    • G01M13/021Gearings
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/02Gearings; Transmission mechanisms
    • G01M13/028Acoustic or vibration analysis

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明公开了一种适用于变转速下齿轮故障诊断的正则化稀疏滤波方法,涉及旋转机械振动信号的故障诊断领域。本发明方法如下:首先获取齿轮箱在不同故障类型下的时变转速振动信号;然后建立基于L1/2范数正则化的稀疏滤波模型(L1/2‑稀疏滤波),将振动信号直接输入到模型中进行训练并提取特征;最后将提取的特征输入到Softmax分类器中进行智能诊断。本发明通过正则化方法提升了稀疏滤波的特征提取能力,准确可靠的实现了时变转速下齿轮箱的智能故障诊断。

Description

一种适用于变转速下齿轮故障诊断的正则化稀疏滤波方法
技术领域
本发明涉及旋转机械振动信号的故障诊断领域,具体涉及一种适用于变转速下齿轮故障诊断的正则化稀疏滤波方法。
背景技术
在航空航天、轨道交通、工程车辆等领域广泛使用的机械设备都在向着高速率、自动化、多功能和精密化方向进化,使得机械设备结构日趋复杂,产生故障的隐患也随之持续上涨。而齿轮是实现这些机械设备功率传递和运动转换的重要机构,也是机械故障诊断的主要研究对象。然而传统的故障诊断方法大都是在恒定转速下进行,对于变转速下的故障诊断研究较少。齿轮在工作过程中转速往往存在较大变化,例如:车辆行驶受道路状况的影响进行调速,风力发电机中齿轮箱的转速会随着风力和风向的改变而不断变化,在大型军舰、歼击机和直升机等军用机械中传动系统中齿轮的转速也会随着机动动作和巡航速度的改变而改变。综上所述,可以看出变转速工况在各个领域中的各种关键设备运行中均会涉及。在应对如此众多复杂机械设备齿轮的故障诊断问题成为了当前研究的一个难点。
当齿轮在恒定转速下工作时,采集的振动信号会与不同的故障类型之间产生很强的相关性,也就是说不同的信号征兆能与不同的故障类型产生相应的映射。因此,传统的故障诊断方法往往是参照设备在不同故障类型下振动信号的时域、频域或时频域信息的特征进行比对,以此来分析得出设备当前的健康状况。但是在变转速工况下,信号与故障类型间的映射不再如此单一,给机械故障诊断带来了严峻的挑战。
综上所述,变转速下齿轮产生的振动信号是一种异常复杂的非平稳信号,它不仅会使传统的故障诊断方法无法应对相同健康状况下的样本之间产生频率调制现象和统计特征参数变化,同时也降低了当下流行的许多数据驱动模型的泛化能力。因此,如何通过更加先进的信号处理技术对变转速工况引起的信号调制现象进行有效抑制,是当前故障诊断领域中的热点和难点问题。因此开展变转速工况下齿轮的智能故障诊断技术是一项急迫发展,也是十分具有挑战的技术。
发明内容
本发明的目的是针对上述不足,提出了采用L1/2-稀疏滤波网络从齿轮的原始变转速信号中自动提取故障特征,并基于该故障特征对齿轮的健康状况进行智能诊断的适用于变转速下齿轮故障诊断的正则化稀疏滤波方法。
本发明具体采用如下技术方案:
一种适用于变转速下齿轮故障诊断的正则化稀疏滤波方法,包括以下步骤:
(1)设置不同的齿轮健康状态,通过手动调节变频器控制电机转速在一定范围内随机波动,由此获得变转速工况,并采集不同故障齿轮在变转速工况下的振动信号;
(2)根据L2、L1和L1/2范数对不同特征分布条件下稀疏解的提取效果,选择采用L1/2范数约束稀疏滤波网络的目标函数,建立基于L1/2范数正则化的稀疏滤波网络模型,模型如式(1)所示,
Figure BDA0002368996110000021
其中,W为权值矩阵,M为样本个数,
Figure BDA0002368996110000022
表示列归一化特征,
Figure BDA0002368996110000023
表示行归一化特征,λ为正则化参数,w为正则化项,C为权值矩阵的行数,R为权值矩阵的列数;
(3)随机选取50%的样本训练L1/2-稀疏滤波网络模型,提取变转速振动信号的波动特征;
(4)将提取的波动特征输入到Softmax分类器中进行样本分类,实现故障诊断。
优选地,步骤(3)中,对L1/2-稀疏滤波网络模型进行训练,包括以下步骤:
(3-1)采用固定步长的卷积操作对转速波动的信号样本做逐步取段处理,组成训练样本集,以此学习时变转速信号的波动特征;
(3-2)对训练集进行Z-Score标准化,使之落入一个特定区间,为网络训练提供一个简易的初始环境,即如式(2),
Figure BDA0002368996110000024
其中,μ为均值,σ为方差。
(3-3)将标准化后的训练集输入到L1/2-稀疏滤波模型中训练出权值矩阵W;
(3-4)对样本进行无卷积均匀取段,取段的长度与卷积取段长度相同;
(3-5)采用权值矩阵对无卷积均匀取段的样本进行特征映射,并采用软阈值函数作为网络的激活函数,得到样本的局部波动特征,即如式(3),
Figure BDA0002368996110000025
(3-6)样本的波动特征fi通过使用全局平均池化来组合局部波动特征获得。
本发明具有如下有益效果:
本发明直接采用L1/2-稀疏滤波网络从齿轮的原始变转速信号中自动提取故障特征,并基于该故障特征对齿轮的健康状况进行智能诊断,消除了转速波动对特征提取的影响,精准智能地实现了变转速工况下齿轮的故障诊断。
附图说明
图1为适用于变转速下齿轮故障诊断的正则化稀疏滤波方法的流程图;
图2a为相同特征分布条件下生成稀疏解的可能性图;
图2b为特征分条件改变下生成稀疏解的可能性图。
具体实施方式
下面结合附图和具体实施例对本发明的具体实施方式做进一步说明:
结合图1,一种适用于变转速下齿轮故障诊断的正则化稀疏滤波方法,包括以下步骤:
(1)设置不同的齿轮健康状态,通过手动调节变频器控制电机转速在一定范围内随机波动,由此获得变转速工况,并采集不同故障齿轮在变转速工况下的振动信号;采用原始振动信号作为训练样本,通过卷积取段法从训练样本中采集Z段信号,这些分段组成训练集
Figure BDA0002368996110000037
,其中
Figure BDA0002368996110000036
是第j个取段包含Nin个数据点,Nin用于表示稀疏滤波的输入维数,Nout表示输出维数。
(2)参照图2a、2b,研究了基于L2、L1和L1/2范数对相同和不同批次样本提取稀疏解的学习过程,可以看出L1/2范数在样本特征分布改变较大的情况下仍能提取最具稀疏性的解,因此更适用于处理变转速工况下齿轮的波动数据,根据L2、L1和L1/2范数对不同特征分布条件下稀疏解的提取效果,选择采用L1/2范数约束稀疏滤波网络的目标函数,建立基于L1/2范数正则化的稀疏滤波网络模型,模型如式(1)所示,
Figure BDA0002368996110000031
其中,W为权值矩阵,M为样本个数,
Figure BDA0002368996110000032
表示列归一化特征,
Figure BDA0002368996110000033
表示行归一化特征,λ为正则化参数,w为正则化项,C为权值矩阵的行数,R为权值矩阵的列数;
(3)随机选取50%的样本训练L1/2-稀疏滤波网络模型,提取变转速振动信号的波动特征,对L1/2-稀疏滤波网络模型进行训练,包括以下子步骤:
(3-1)采用固定步长的卷积操作对转速波动的信号样本做逐步取段处理,组成训练样本集,以此学习时变转速信号的波动特征;
(3-2)对训练集进行Z-Score标准化,数据集
Figure BDA0002368996110000034
写作矩阵形式
Figure BDA0002368996110000035
然后进行Z-Score标准化处理,使之落入一个特定区间,为网络训练提供一个简易的初始环境,即如式(2),
Figure BDA0002368996110000041
其中,μ为均值,σ为方差。
(3-3)将标准化后的训练集输入到L1/2-稀疏滤波模型中训练出权值矩阵W;
(3-4)对样本进行无卷积均匀取段,取段的长度与卷积取段长度相同;
(3-5)训练样本xi通过无卷积取段成K段,其中K是整数且等于N/Nin。K段组成数据集
Figure BDA0002368996110000042
Figure BDA0002368996110000043
然后使用权值矩阵W来映射样本,并采用软阈值函数作为网络的激活函数,得到样本的局部波动特征,即如式(3),
Figure BDA0002368996110000044
(3-6)样本的波动特征fi通过使用全局平均池化来组合局部波动特征获得。
(4)将提取的波动特征输入到Softmax分类器中进行样本分类,实现故障诊断;
将所有样本的学习特征fi与样本标签相结合输入到Softmax分类器中进行训练,然后采用剩余样本作为测试样本来测试提出方法的准确度。假设yi∈{1,2,...,k}为含有k个类型的标签集,p(yi=j|x)为每一个样本的条件概率。假设函数用于计算归于每一个分类j的概率。因此输出hθ(xi)为一个k维向量:
Figure BDA0002368996110000045
其中θ12,…θk为模型的参数,
Figure BDA0002368996110000046
用于对分布进行标准化,以保证p(yi=j|x)的和为1。
方便起见,1{·}表示指示函数,即1{true}=1,1{false}=0。隐层回归模型的代价函数为:
Figure BDA0002368996110000047
其中m为训练样本个数,λ为权重惩罚项,用于对较大值的参数进行惩罚,并通过全局优化达到更高的精度。
首先对每一个测试样本进行无卷积取段,然后采用训练好的正则化稀疏滤波模型对测试样本的取段信息进行局部波动特征学习,并对学习到的特征进行取平均值,最后将学习到的波动特征通过Softmax回归模型与真实的故障类型标签相结合,进行分类得出测试准确率。
通过对一个特殊设计的行星齿轮箱故障试验台进行时变转速下的振动信号采集和智能诊断,对本发明作进一步描述。
共设置10种不同的齿轮健康状况:正常;三种太阳轮故障:裂纹、点蚀和磨损,分别命名为太阳轮1、太阳轮2、太阳轮3;三种行星轮故障:裂纹、点蚀和磨损,分别命名为行星轮1、行星轮2、行星轮3;三种复合齿轮故障:太阳轮磨损和行星轮磨损复合、太阳轮点蚀和行星轮裂纹复合、太阳轮点蚀和行星轮磨损复合,分别命名为复合1、复合2、复合3。齿轮的转速在700-1500r/min之间进行无规律波动,设置采样频率为12.8kHz,每个齿轮健康状况下采集500个样本,每个样本1200个数据点。
使用本发明方法诊断变转速齿轮数据集,针对该数据集,选择100作为稀疏滤波的输入和输出维数,随机选取每种健康状况下50%的样本来训练正则化的稀疏滤波模型,然后设置卷积尺寸为100、卷积步长为20,对每个样本采集50段,随后使用这些数据段组成的训练集来训练正则化的稀疏滤波模型,并获得权值矩阵W。在计算局部波动特征时,将每个样本无卷积取12段,每段100个数据点。对每一个样本来说,通过训练得到的稀疏滤波模型可以提取100维的学习特征。使用本发明方法,可以达到100%的训练准确率和99.03%的测试准确率。为验证本发明方法的有效性,将未加正则化的稀疏滤波模型和基于L2、L1范数正则化的稀疏滤波模型进行对比,这四种方法的比较结果如下表所示:
表1
Figure BDA0002368996110000051
可以看出L1-稀疏滤波方法的准确率结果与未加正则化的稀疏滤波方法几乎没有差别,在训练和测试集上均表现很差,因此L1正则化对处理时变转速工况没有任何作用,这也与图2a、2b的分析结果相一致。L2-稀疏滤波方法和L1/2-稀疏滤波方法训练准确率均为100%,在测试准确率上L1/2-稀疏滤波方法的表现更好,由此可以看出L1/2-稀疏滤波对处理齿轮变转速工况有明显效果。
通过以上变转速下齿轮故障诊断的结果分析和对比实验可知,本发明的L1/2-稀疏滤波方法能够克服时变转速工况带来的数据特征变化,学习到的波动特征向量也更有利于进行故障识别和分类,拥有处理时变转速下齿轮故障诊断问题的优越性能。
当然,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换,也应属于本发明的保护范围。

Claims (1)

1.一种适用于变转速下齿轮故障诊断的正则化稀疏滤波方法,其特征在于,包括以下步骤:
(1)设置不同的齿轮健康状态,通过手动调节变频器控制电机转速在一定范围内随机波动,由此获得变转速工况,并采集不同故障齿轮在变转速工况下的振动信号;
(2)根据L2、L1和L1/2范数对不同特征分布条件下稀疏解的提取效果,选择采用L1/2范数约束稀疏滤波网络的目标函数,建立基于L1/2范数正则化的稀疏滤波网络模型,模型如式(1)所示,
Figure FDA0003172658730000011
其中,W为权值矩阵,M为样本个数,
Figure FDA0003172658730000012
表示列归一化特征,
Figure FDA0003172658730000013
表示行归一化特征,λ为正则化参数,w为正则化项,C为权值矩阵的行数,R为权值矩阵的列数;
(3)随机选取50%的样本训练L1/2-稀疏滤波网络模型,提取变转速振动信号的波动特征;
(4)将提取的波动特征输入到Softmax分类器中进行样本分类,实现故障诊断;
步骤(3)中,对L1/2-稀疏滤波网络模型进行训练,包括以下步骤:
(3-1)采用固定步长的卷积操作对转速波动的信号样本做逐步取段处理,组成训练样本集,以此学习时变转速信号的波动特征;
(3-2)对训练集进行Z-Score标准化,使之落入一个特定区间,为网络训练提供一个简易的初始环境,即如式(2),
Figure FDA0003172658730000014
其中,μ为均值,σ为方差;
(3-3)将标准化后的训练集输入到L1/2-稀疏滤波模型中训练出权值矩阵W;
(3-4)对样本进行无卷积均匀取段,取段的长度与卷积取段长度相同;
(3-5)采用权值矩阵对无卷积均匀取段的样本进行特征映射,并采用软阈值函数作为网络的激活函数,得到样本的局部波动特征,即如式(3),
Figure FDA0003172658730000015
(3-6)样本的波动特征fk通过使用全局平均池化来组合局部波动特征获得。
CN202010044852.6A 2020-01-16 2020-01-16 一种适用于变转速下齿轮故障诊断的正则化稀疏滤波方法 Active CN111089720B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010044852.6A CN111089720B (zh) 2020-01-16 2020-01-16 一种适用于变转速下齿轮故障诊断的正则化稀疏滤波方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010044852.6A CN111089720B (zh) 2020-01-16 2020-01-16 一种适用于变转速下齿轮故障诊断的正则化稀疏滤波方法

Publications (2)

Publication Number Publication Date
CN111089720A CN111089720A (zh) 2020-05-01
CN111089720B true CN111089720B (zh) 2021-10-08

Family

ID=70400492

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010044852.6A Active CN111089720B (zh) 2020-01-16 2020-01-16 一种适用于变转速下齿轮故障诊断的正则化稀疏滤波方法

Country Status (1)

Country Link
CN (1) CN111089720B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111855208B (zh) * 2020-07-02 2021-09-21 南京航空航天大学 一种基于阶次分析和稀疏编码的旋转机械微弱故障信号提取方法
CN112362343A (zh) * 2020-10-28 2021-02-12 华南理工大学 基于调频字典的齿轮箱变转速下分布型故障特征提取方法
CN112487890B (zh) * 2020-11-17 2022-03-04 山东科技大学 一种基于并行稀疏滤波的轴承声信号故障诊断方法
CN113075546A (zh) * 2021-03-24 2021-07-06 河南中烟工业有限责任公司 电机振动信号的特征提取方法及系统
CN113076844B (zh) * 2021-03-26 2022-05-20 华中科技大学 一种旋转部件故障诊断模型构建方法及应用
CN114061947B (zh) * 2021-09-29 2023-01-24 上海交通大学 基于稀疏时频分析的齿轮箱变转速故障诊断方法及系统
CN114217256B (zh) * 2021-11-25 2022-09-16 华中科技大学 基于自适应稀疏滤波的电机驱动系统变流器故障诊断方法
JP7459857B2 (ja) 2021-12-10 2024-04-02 トヨタ自動車株式会社 異音解析装置および方法
CN115389247B (zh) * 2022-11-01 2023-02-03 青岛睿发工程咨询服务合伙企业(有限合伙) 一种基于速度自适应自编码器的旋转机械故障监测方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108509701A (zh) * 2018-03-21 2018-09-07 南京航空航天大学 一种基于振动信号的旋转机械故障直接智能诊断方法
CN108805059A (zh) * 2018-05-29 2018-11-13 东华大学 稀疏正则化滤波与自适应稀疏分解的齿轮箱故障诊断方法
CN108827605A (zh) * 2018-03-20 2018-11-16 南京航空航天大学 一种基于改进稀疏滤波的机械故障特征自动提取方法
CN109186973A (zh) * 2018-08-20 2019-01-11 东南大学 一种无监督深度学习网络的机械故障诊断方法
CN110428004A (zh) * 2019-07-31 2019-11-08 中南大学 数据失衡下基于深度学习的机械零部件故障诊断方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10367477B2 (en) * 2015-06-12 2019-07-30 Analog Devices, Inc. Sparse cascaded-integrator-comb filters
CN106017876A (zh) * 2016-05-11 2016-10-12 西安交通大学 基于等权局部特征稀疏滤波网络的轮对轴承故障诊断方法
CN106067042B (zh) * 2016-06-13 2019-02-15 西安电子科技大学 基于半监督深度稀疏滤波网络的极化sar分类方法
CN108871762A (zh) * 2018-06-29 2018-11-23 广东工业大学 一种风电机组齿轮箱故障诊断方法
CN110580471B (zh) * 2019-09-12 2021-11-02 北京航空航天大学 基于编码器信号瞬变特征的机械设备故障诊断方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108827605A (zh) * 2018-03-20 2018-11-16 南京航空航天大学 一种基于改进稀疏滤波的机械故障特征自动提取方法
CN108509701A (zh) * 2018-03-21 2018-09-07 南京航空航天大学 一种基于振动信号的旋转机械故障直接智能诊断方法
CN108805059A (zh) * 2018-05-29 2018-11-13 东华大学 稀疏正则化滤波与自适应稀疏分解的齿轮箱故障诊断方法
CN109186973A (zh) * 2018-08-20 2019-01-11 东南大学 一种无监督深度学习网络的机械故障诊断方法
CN110428004A (zh) * 2019-07-31 2019-11-08 中南大学 数据失衡下基于深度学习的机械零部件故障诊断方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于迭代经验小波变换的齿轮故障诊断方法;辛玉等;《仪器仪表学报》;20181130;第39卷(第11期);全文 *

Also Published As

Publication number Publication date
CN111089720A (zh) 2020-05-01

Similar Documents

Publication Publication Date Title
CN111089720B (zh) 一种适用于变转速下齿轮故障诊断的正则化稀疏滤波方法
Kaplan et al. An improved feature extraction method using texture analysis with LBP for bearing fault diagnosis
Wang et al. Convolutional neural network-based hidden Markov models for rolling element bearing fault identification
Li et al. Understanding and improving deep learning-based rolling bearing fault diagnosis with attention mechanism
CN110595775B (zh) 基于多分支多尺度卷积神经网络的滚动轴承故障诊断方法
CN110866365A (zh) 基于部分迁移卷积网络的机械设备智能故障诊断方法
CN110543860B (zh) 基于tjm迁移学习的机械故障诊断方法及系统
CN111397896B (zh) 一种旋转机械故障诊断方法、系统及存储介质
CN112668459B (zh) 基于有监督对比学习和卷积神经网络的滚动轴承故障诊断方法
CN112084974A (zh) 一种基于元学习的多标签滚动轴承故障诊断方法
CN112881942B (zh) 基于小波分解与经验模态分解的异常电流诊断方法和系统
Wang et al. Construction of the efficient attention prototypical net based on the time–frequency characterization of vibration signals under noisy small sample
CN113095179A (zh) 元度量学习驱动的直升机行星齿轮箱故障诊断方法
CN115219197A (zh) 一种轴承故障智能诊断方法、系统、介质、设备及终端
CN115905976B (zh) 基于注意力机制的Highway Bi-LSTM轴承故障诊断方法、系统及设备
CN115358259A (zh) 一种基于自学习的无监督跨工况轴承故障诊断方法
Cao et al. Fault diagnosis of rolling bearing based on multiscale one-dimensional hybrid binary pattern
CN114169377A (zh) 基于g-mscnn的有噪环境中滚动轴承故障诊断方法
CN111665050A (zh) 一种基于聚类k-svd算法的滚动轴承故障诊断方法
Han et al. An intelligent fault diagnosis method of variable condition gearbox based on improved DBN combined with WPEE and MPE
Wang et al. A novel data augmentation method for intelligent fault diagnosis under speed fluctuation condition
CN116070128A (zh) 一种基于对比自监督学习的时间序列设备故障诊断方法
He et al. Fault identification of rotating machinery based on dynamic feature reconstruction signal graph
CN114199570B (zh) 基于改进变分模态分解的风电齿轮箱轴承故障识别方法
CN113158769A (zh) 基于CEEMDAN和FastICA的机电设备轴承振动信号除噪方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant