CN106289508A - 一种面向机械故障诊断的振动信号重构方法 - Google Patents
一种面向机械故障诊断的振动信号重构方法 Download PDFInfo
- Publication number
- CN106289508A CN106289508A CN201610569483.6A CN201610569483A CN106289508A CN 106289508 A CN106289508 A CN 106289508A CN 201610569483 A CN201610569483 A CN 201610569483A CN 106289508 A CN106289508 A CN 106289508A
- Authority
- CN
- China
- Prior art keywords
- sparse coding
- vibration signal
- negative
- coding matrix
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M99/00—Subject matter not provided for in other groups of this subclass
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明提供了一种面向机械故障诊断的振动信号重构方法,在非负条件下,对原始振动信号进行更为稀疏的编码,提高信号频谱检测的有效性,属于振动信号处理技术领域。1、对机械系统振动信号采样,将原始振动信号重叠分割得到列向量构成训练样本。2、初始化学习字典,基于非负基追踪算法对训练样本进行稀疏编码得到初始化稀疏编码矩阵。3、固定稀疏编码矩阵,基于非负K奇异值分解算法优化更新学习字典。4、固定学习字典,基于非负基追踪算法进行稀疏编码更新对应的稀疏编码矩阵。5、循环前两步,交替迭代更新得到最终学习字典和最终稀疏编码矩阵。6、计算更新样本,逆叠加重构信号。7、提取重构信号的包络谱实现故障诊断。用于机械故障诊断。
Description
技术领域
本发明属于振动信号处理技术领域,主要涉及机械设备振动信号的处理技术。
背景技术
机械设备运行状态信息的获取是系统故障诊断的前提。由于振动信号自身实时性、周期性以及加速度传感器体积小,频带宽,易安装,高频响应稳定等特性,振动信号已经成为当前机械系统故障诊断中应用最广泛的信息载体。针对旋转机械关键零部件,基于振动信号频谱检测实现故障诊断的方法应用极为广泛。诊断过程中,振动信号滤波去噪(预处理)是关键,也是信号处理领域的研究热点和难点,其有效性直接决定了方法的可行性。常用的去噪方法通常基于小波变换,经验模态分解,谱峭度,盲源分离,自回归模型等算法,均取得不错的效果,但也存在过度依赖于经验知识和人工确定参数等局限性。
近年来,自压缩采样理论被论证,根据数据的可压缩性,通过低维欠奈奎斯特数据实现高维采样,极大地推动了字典学习在数据分类,信息压缩,图像处理等领域的发展。为了提高故障诊断的有效性,本发明更注重机械系统结构和信号本质,区别于传统去噪方法单纯减弱或消除振动信号噪声成分,采用基于非负K奇异值分解(NN-K-SVD)算法改进的字典学习算法重构振动信号,得到一种面向机械故障诊断的振动信号重构方法。
发明内容
本发明的目的是提供一种面向机械故障诊断的振动信号重构方法,它能有效地重构振动信号并解决机械系统故障诊断问题。
本发明的目的是通过以下技术方案来实现的:一种面向机械故障诊断的振动信号重构方法,在非负条件下,对振动信号进行更为稀疏的编码,提高振动信号频谱检测的有效性。首先,由加速度传感器采集机械系统一维原始振动信号,重叠、分割得到列向量,并构成训练样本M。初始化学习字典D,基于非负基追踪(NN-BP)算法和非负K奇异值分解(NN-K-SVD)算法对学习字典D进行改进,通过更新稀疏编码矩阵C和学习字典D,交替迭代生成最终学习字典和对应的最终稀疏编码矩阵计算得到更新样本逆叠加重构信号。原始振动信号通过改进的字典学习算法更逼近信号本质结构,达到了较优的预处理效果,有利于通过提取重构信号的包络谱实现故障诊断。
具体实现步骤包括:
第一步,对机械系统连续的一维原始振动信号进行采样,将原始振动信号重叠、分割得到列向量,构成训练样本M:
M=M0+δ
式中,M0表征理论上无噪声成分,δ表征噪声成分;
第二步,初始化学习字典D=[d1,d2,...],每一列对应一原子d,基于非负基追踪算法对训练样本M进行稀疏编码,得到初始化稀疏编码矩阵C:
式中,学习字典D、训练样本M、稀疏编码矩阵C的任意元素均大于0,对于任意原子d,||d||=1,常数ω>0,表征振动信号重构精度和稀疏性的相对数值关系;第一项求M-DC的l2范数重构误差,第二项是稀疏编码惩罚项,约束稀疏分布,确保能量函数逐步减小。
第三步,固定稀疏编码矩阵C,基于非负K奇异值分解算法优化更新学习字典D:
Dt=Dt-1-υ(Dt-1Ct-1-M)(Ct-1)'
式中,t是迭代次数,(.)'是矩阵转置,迭代步长υ>0,原子d<0时归零,确保学习字典D非负。
第四步,固定学习字典D,基于非负基追踪算法进行稀疏编码,更新对应的稀疏编码矩阵C:
Ct=Ct-1.×(Dt'M)./(Dt'DtCt-1+ω)
式中,.×和./是矩阵对应位置元素相乘除,稀疏编码矩阵C非负,确保达到要求精度下全局最小值。
第五步,循环第三步和第四步,交替迭代更新得到最终学习字典和最终稀疏编码矩阵
第六步,计算更新样本逆叠加重构信号:
式中,求得的更新样本相较于训练样本M少了部分噪声成分。
第七步,逆叠加重构信号,提取重构信号的包络谱实现故障诊断。
本发明与现有技术相比的优点和效果:
(1)本发明主要思想是通过原始振动信号自身驱动构成训练样本进行信号预处理,区别于传统去噪方法单纯减弱或消除振动信号噪声成分,重构信号更符合振动信号本质和内在规律。
(2)本发明对原始振动信号在学习字典域进行稀疏表示,降低了数据在存储和传输的成本。
(3)本发明采用改进的字典学习算法重构振动信号,对K-SVD算法和NN-K-SVD算法进行改进,添加了非负约束,实现信号局部特征提取,得到的编码更稀疏,有益于数据压缩,优化和简化了算法。
附图说明
图1为本发明的面向机械系统故障诊断的振动信号重构方法框图
图2为本发明实施例原始振动信号
图3为本发明实施例不含噪声的振动信号
图4为本发明实施例改进的字典学习算法重构振动信号
具体实施方式
下面结合附图,对本发明作进一步的描述,方法框图如图1所示。
(1)仿真采集得到机械系统振动信号如图2,其中包含的无噪成分如图3,随机截取2048个数据点,将原始信号重叠72个数据点分割得到含有80个数据点的列向量,构成训练样本M。
M=M0+δ
式中,M0表征理论上无噪声成分,δ表征噪声成分。
(2)初始化学习字典D=[d1,d2,...],每一列对应一原子d,尺寸80×10。基于非负基追踪算法对训练样本M进行稀疏编码,得到尺寸10×247的初始化稀疏编码矩阵C:
式中,学习字典D、训练样本M、稀疏编码矩阵C的任意元素均大于0,对于任意原子d,||d||=1,常数ω>0,表征振动信号重构精度和稀疏性的相对数值关系。第一项求M-DC的l2范数重构误差,第二项是稀疏编码惩罚项,约束稀疏分布,确保能量函数逐步减小。
(3)固定稀疏编码矩阵C,基于非负K奇异值分解算法优化更新学习字典D:
Dt=Dt-1-υ(Dt-1Ct-1-M)(Ct-1)'
式中,t是迭代次数,(.)'是矩阵转置,迭代步长υ>0,原子d<0时归零,确保学习字典D非负。
(4)固定学习字典D,取稀疏阀值3,基于非负基追踪算法进行稀疏编码,更新对应的稀疏编码矩阵C:
Ct=Ct-1.×(Dt'M)./(Dt'DtCt-1+ω)
式中,.×和./是矩阵对应位置元素相乘除,稀疏编码矩阵C非负,确保达到要求精度下全局最小值。
(5)循环(3)和(4),交替迭代20次,得到最终学习字典和最终稀疏编码矩阵
(6)计算更新样本逆叠加得到重构信号如图4。
式中,求得的更新样本相较于训练样本M少了部分噪声成分。
(7)逆叠加重构信号,提取重构信号的包络谱实现故障诊断。
本实施例中核心算法由C语言编写完成,人机交互界面和逻辑操作程序由Python语言编写完成。结果可知,原始振动信号经本发明提供的方法得到的重构信号,有效地保留了原始振动信号本质和内在规律成分,获得预处理的有益效果,达到压缩采样目的,提高振动信号频谱检测的有效性,实现机械系统故障诊断。
Claims (1)
1.一种面向机械故障诊断的振动信号重构方法,具体实施步骤包括:
第一步、对机械系统连续的原始振动信号进行采样,将原始振动信号重叠、分割得到列向量,构成训练样本M:
M=M0+δ
式中,M0表征理论上无噪声成分,δ表征噪声成分;
第二步、初始化学习字典D=[d1,d2,...],每一列对应一原子d,基于非负基追踪算法对训练样本M进行稀疏编码,得到初始化稀疏编码矩阵C:
式中,学习字典D、训练样本M、稀疏编码矩阵C的任意元素均大于0,对于任意原子d,||d||=1,常数ω>0,表征振动信号重构精度和稀疏性的相对数值关系;第一项求M-DC的l2范数重构误差,第二项是稀疏编码惩罚项,约束稀疏分布,确保能量函数逐步减小;
第三步、固定稀疏编码矩阵C,基于非负K奇异值分解算法优化更新学习字典D:
Dt=Dt-1-υ(Dt-1Ct-1-M)(Ct-1)'
式中,t是迭代次数,(.)'是矩阵转置,迭代步长υ>0,原子d<0时归零,确保学习字典D非负;
第四步、固定学习字典D,基于非负基追踪算法进行稀疏编码,更新对应的稀疏编码矩阵C:
Ct=Ct-1.×(Dt'M)./(Dt'DtCt-1+ω)
式中,.×和./是矩阵对应位置元素相乘除,稀疏编码矩阵C非负,确保达到要求精度下全局最小值;
第五步、循环第三步和第四步,交替迭代更新得到最终学习字典和最终稀疏编码矩阵
第六步、计算更新样本逆叠加重构信号;
式中,求得的更新样本相较于训练样本M少了部分噪声成分;
第七步、逆叠加重构信号,提取重构信号的包络谱实现故障诊断。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610569483.6A CN106289508A (zh) | 2016-07-19 | 2016-07-19 | 一种面向机械故障诊断的振动信号重构方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610569483.6A CN106289508A (zh) | 2016-07-19 | 2016-07-19 | 一种面向机械故障诊断的振动信号重构方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106289508A true CN106289508A (zh) | 2017-01-04 |
Family
ID=57651657
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610569483.6A Pending CN106289508A (zh) | 2016-07-19 | 2016-07-19 | 一种面向机械故障诊断的振动信号重构方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106289508A (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111126425A (zh) * | 2018-11-01 | 2020-05-08 | 中南大学 | 复杂工业过程工况监测和分类方法 |
WO2020113575A1 (zh) * | 2018-12-07 | 2020-06-11 | 广东省智能制造研究所 | 基于约束半非负矩阵分解的声音分类方法、装置及介质 |
CN111582137A (zh) * | 2020-04-30 | 2020-08-25 | 燕山大学 | 一种滚动轴承信号重构方法及系统 |
CN112183300A (zh) * | 2020-09-23 | 2021-01-05 | 厦门大学 | 一种基于多层次稀疏表示的ais辐射源识别方法及系统 |
CN112304419A (zh) * | 2020-10-25 | 2021-02-02 | 广东石油化工学院 | 一种利用泛化稀疏编码的振声检测信号重构方法和系统 |
CN113449683A (zh) * | 2021-07-15 | 2021-09-28 | 江南大学 | 基于k-svd训练局部字典的高频超声稀疏去噪方法及系统 |
CN113743240A (zh) * | 2021-08-12 | 2021-12-03 | 武汉钢铁有限公司 | 故障检测方法、装置及计算机可读存储介质 |
CN114970645A (zh) * | 2022-07-27 | 2022-08-30 | 中国工业互联网研究院 | 基于5g边云协同的民用航空发动机故障诊断系统及方法 |
CN115524150A (zh) * | 2022-09-13 | 2022-12-27 | 西安交通大学 | 基于稀疏时域同步平均的旋转机械故障检测方法 |
CN116070091A (zh) * | 2023-03-01 | 2023-05-05 | 北京大学 | 基于在线字典学习的多模态机械加工中心信号的重构方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20010037180A1 (en) * | 2000-04-20 | 2001-11-01 | Hidemichi Komura | Fault diagnosis method and apparatus |
JP5428994B2 (ja) * | 2010-03-26 | 2014-02-26 | 新日鐵住金株式会社 | 機械設備の異常診断方法、装置及びプログラム |
CN103679662A (zh) * | 2013-12-25 | 2014-03-26 | 苏州市职业大学 | 基于类别先验非负稀疏编码字典对的超分辨率图像恢复方法 |
CN104935349A (zh) * | 2015-06-04 | 2015-09-23 | 西南交通大学 | 一种振动信号压缩采样方法 |
CN105258947A (zh) * | 2015-11-06 | 2016-01-20 | 北京航空航天大学 | 一种基于压缩感知的工况扰动条件下的滚动轴承故障诊断方法 |
-
2016
- 2016-07-19 CN CN201610569483.6A patent/CN106289508A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20010037180A1 (en) * | 2000-04-20 | 2001-11-01 | Hidemichi Komura | Fault diagnosis method and apparatus |
JP5428994B2 (ja) * | 2010-03-26 | 2014-02-26 | 新日鐵住金株式会社 | 機械設備の異常診断方法、装置及びプログラム |
CN103679662A (zh) * | 2013-12-25 | 2014-03-26 | 苏州市职业大学 | 基于类别先验非负稀疏编码字典对的超分辨率图像恢复方法 |
CN104935349A (zh) * | 2015-06-04 | 2015-09-23 | 西南交通大学 | 一种振动信号压缩采样方法 |
CN105258947A (zh) * | 2015-11-06 | 2016-01-20 | 北京航空航天大学 | 一种基于压缩感知的工况扰动条件下的滚动轴承故障诊断方法 |
Non-Patent Citations (2)
Title |
---|
PATRIK O. HOYER, ET AL.: "NON-NEGATIVE SPARSE CODING", 《NETWORKS FOR SIGNAL PROCESSING》 * |
周亚同 等: "缩感知框架下基于K-奇异值分解字典学习的地震数据重建", 《石油地球物理勘探》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111126425B (zh) * | 2018-11-01 | 2023-04-25 | 中南大学 | 复杂工业过程工况监测和分类方法 |
CN111126425A (zh) * | 2018-11-01 | 2020-05-08 | 中南大学 | 复杂工业过程工况监测和分类方法 |
WO2020113575A1 (zh) * | 2018-12-07 | 2020-06-11 | 广东省智能制造研究所 | 基于约束半非负矩阵分解的声音分类方法、装置及介质 |
CN111582137A (zh) * | 2020-04-30 | 2020-08-25 | 燕山大学 | 一种滚动轴承信号重构方法及系统 |
CN112183300A (zh) * | 2020-09-23 | 2021-01-05 | 厦门大学 | 一种基于多层次稀疏表示的ais辐射源识别方法及系统 |
CN112183300B (zh) * | 2020-09-23 | 2024-03-22 | 厦门大学 | 一种基于多层次稀疏表示的ais辐射源识别方法及系统 |
CN112304419A (zh) * | 2020-10-25 | 2021-02-02 | 广东石油化工学院 | 一种利用泛化稀疏编码的振声检测信号重构方法和系统 |
CN113449683A (zh) * | 2021-07-15 | 2021-09-28 | 江南大学 | 基于k-svd训练局部字典的高频超声稀疏去噪方法及系统 |
CN113743240A (zh) * | 2021-08-12 | 2021-12-03 | 武汉钢铁有限公司 | 故障检测方法、装置及计算机可读存储介质 |
CN114970645A (zh) * | 2022-07-27 | 2022-08-30 | 中国工业互联网研究院 | 基于5g边云协同的民用航空发动机故障诊断系统及方法 |
CN115524150A (zh) * | 2022-09-13 | 2022-12-27 | 西安交通大学 | 基于稀疏时域同步平均的旋转机械故障检测方法 |
CN115524150B (zh) * | 2022-09-13 | 2024-04-05 | 西安交通大学 | 基于稀疏时域同步平均的旋转机械故障检测方法 |
CN116070091A (zh) * | 2023-03-01 | 2023-05-05 | 北京大学 | 基于在线字典学习的多模态机械加工中心信号的重构方法 |
CN116070091B (zh) * | 2023-03-01 | 2024-05-24 | 北京大学 | 基于在线字典学习的多模态机械加工中心信号的重构方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106289508A (zh) | 一种面向机械故障诊断的振动信号重构方法 | |
Guo et al. | A review of wavelet analysis and its applications: Challenges and opportunities | |
Allen et al. | Signal analysis: time, frequency, scale, and structure | |
CN105827250A (zh) | 一种基于自适应字典学习的电能质量数据压缩重构方法 | |
Hu et al. | A new method of wind turbine bearing fault diagnosis based on multi-masking empirical mode decomposition and fuzzy C-means clustering | |
CN113865866B (zh) | 一种基于改进局部非负矩阵分解的轴承复合故障诊断方法 | |
CN110164465A (zh) | 一种基于深层循环神经网络的语音增强方法及装置 | |
CN104458250A (zh) | 智能化齿轮箱故障诊断方法 | |
CN106295690A (zh) | 基于非负矩阵分解的时间序列数据聚类方法及系统 | |
CN109839271A (zh) | 一种基于匹配追踪稀疏算法的轴承故障特征提取方法 | |
CN105869126A (zh) | 高光谱图像的压缩感知去噪方法 | |
CN104915931A (zh) | 基于压缩感知理论的自适应光学波前信号去噪方法 | |
CN112990139A (zh) | 基于变模态分解加权重构信号结合小波阈值去噪方法 | |
CN109766798A (zh) | 基于经验小波的手势数据处理方法、服务器及感知设备 | |
Bright et al. | Classification of spatiotemporal data via asynchronous sparse sampling: Application to flow around a cylinder | |
CN103093431A (zh) | 基于pca字典和结构先验信息的压缩感知重构方法 | |
CN112884062B (zh) | 一种基于cnn分类模型和生成对抗网络的运动想象分类方法及系统 | |
Zhang et al. | Variable spectral segmentation empirical wavelet transform for noisy signal processing | |
CN104915935A (zh) | 基于非线性压缩感知与字典学习的压缩光谱成像方法 | |
Chen et al. | Physics-enhanced pca for data compression in edge devices | |
CN104391325A (zh) | 不连续非均质地质体检测方法和装置 | |
Yu et al. | Sparse time–frequency representation for the transient signal based on low-rank and sparse decomposition | |
Gao et al. | Application of multi-layer denoising based on ensemble empirical mode decomposition in extraction of fault feature of rotating machinery | |
Liu et al. | Sparse coefficient fast solution algorithm based on the circulant structure of a shift-invariant dictionary and its applications for machine fault diagnosis | |
CN116046396A (zh) | 一种轴承故障诊断方法、装置、设备及介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170104 |