CN114254679A - Filter-based feature enhancement method - Google Patents
Filter-based feature enhancement method Download PDFInfo
- Publication number
- CN114254679A CN114254679A CN202111624531.4A CN202111624531A CN114254679A CN 114254679 A CN114254679 A CN 114254679A CN 202111624531 A CN202111624531 A CN 202111624531A CN 114254679 A CN114254679 A CN 114254679A
- Authority
- CN
- China
- Prior art keywords
- signal
- points
- noise
- valley
- peak
- 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
- 238000000034 method Methods 0.000 title claims abstract description 20
- 238000001914 filtration Methods 0.000 claims description 14
- 238000000605 extraction Methods 0.000 claims description 4
- 230000000877 morphologic effect Effects 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000013016 damping Methods 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 claims 1
- 238000003745 diagnosis Methods 0.000 abstract description 13
- 230000009467 reduction Effects 0.000 description 3
- 230000002708 enhancing effect Effects 0.000 description 2
- 238000005096 rolling process Methods 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/04—Bearings
- G01M13/045—Acoustic or vibration analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
- G06F2218/06—Denoising by applying a scale-space analysis, e.g. using wavelet analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
- G06F2218/10—Feature extraction by analysing the shape of a waveform, e.g. extracting parameters relating to peaks
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
The invention belongs to the technical field of fault diagnosis, and particularly relates to a filter-based feature enhancement method. The characteristic enhancement method based on the filter provides an effective means for fault diagnosis and signal processing of the bearing, avoids safety accidents caused by bearing faults, and has certain practicability and engineering value.
Description
Technical Field
The invention belongs to the technical field of signal processing, and particularly relates to a characteristic enhancement method based on a filter.
Background
Rolling bearings are one of the most widely used components in mechanical equipment. Detection of its state and fault diagnosis have long been the technical focus of scientific research and industry. Generally, vibration signals of bearings often contain a large amount of noise, and the existence of the noise can interfere with effective identification of fault impact, which brings great difficulty to fault diagnosis of the bearings. Signals in engineering practice are often signals containing a large amount of noise, so that the research on algorithms for signal noise reduction and bearing characteristic enhancement has important research significance for engineering application.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: in the prior art, vibration signals of the bearing often contain a large amount of noise, and the existence of the noise can interfere effective identification of fault impact, so that great difficulty is brought to fault diagnosis of the bearing. The invention aims to provide a filter-based characteristic enhancement method, which aims to realize the purposes of enhancing the fault impact characteristic and weakening the noise in a vibration signal of a rolling bearing and provides basic support for the fault prediction and the health management of the bearing.
The technical scheme adopted by the invention for solving the technical problems is as follows: a method for filter-based feature enhancement, comprising the steps of:
s1, preprocessing the acceleration signal containing noise by adopting lifting wavelet transform:
acquiring a signal s (t) with a data point n by using an exponential decay function containing impact characteristics:
wherein A is a displacement constant, xi is a damping ratio, wnT is n points at a natural frequency of vibration and f is a sampling frequencysThe relationship between t and f is (1: n)/fs(ii) a Then adding noise into the s (n) signal, wherein the overall amplitude of the noise is larger than that of the impact characteristics, so that the impact characteristics s (n) in the signal can be submerged, and the characteristics are not obvious; primarily processing a signal to be processed containing noise by adopting lifting wavelet transform to obtain a new signal f (n); lifting wavelet transform is a faster and more efficient wavelet transform method, also known as second generation wavelet transform. The method does not depend on Fourier transform, and inherits the multi-resolution characteristics of the first generation wavelet; the whole lifting wavelet transformation process comprises decomposition and reconstruction;
s2, searching all peak points and all valley points in the new signal f (n):
finding all peak points in the new signal f (n), and recording the found peak points as A ═ A { (A)pFind all valley points in the new signal f (n), and mark it as B ═ B { (B) }hH is 1,2, …, and then the interval I between each adjacent peak points is calculatedA={iA│iA=Ap+1,x-Ap,xCalculating the interval G between each adjacent valley pointB={gB│gB=Bh+1,x-Bh,xIn which A isp,xFor the peak value A in the signalpCorresponding horizontal coordinate value, Bh,xAs the valley B in the signalhCorresponding horizontal coordinate values;
s3, constructing a signal library:
calculating the peak point interval IA={iA│iA=Ap+1,x-Ap,xAnd valley point spacing GB={gB│gB=Bh+1,x-Bh,xMaximum value of max (I)A,GB) Selecting the size of C1(C1≥max(IA,GB) A linear type structural element difference filtering operator carries out filtering processing on the vibration signals f (n), and the processed result is recorded as fC1(n); then using the scale C2(C2≥max(IA,GB) A linear type structural element difference filtering operator carries out filtering processing on the vibration signals f (n), and the processed result is recorded as fC2(n); and analogizing in sequence, after m times of processing, establishing all processing results into a signal library { fC1(n),fC2(n),fC3(n),...,fCm(n)};
S4, extracting coincident impact characteristics from the constructed signal library:
in signal library { fC1(n),fC2(n),fC3(n),...,fCm(n) } extraction of the coincident impact portion z (w):
z(w)=fC1{find(fC1(n)=fC2(n)=...=fCm(n))} (2)。
in step S3, the linear-type difference filter operator is an operation form of morphological filtering.
The characteristic enhancement method based on the filter has the advantages that the characteristic enhancement method based on the filter is provided based on the mathematical principle of the morphological algorithm, the method can not only realize noise reduction, but also well enhance the impact characteristic of the fault, provides good basic support for fault diagnosis of the bearing, including qualitative diagnosis, quantitative diagnosis, outer ring positioning diagnosis and the like, prevents major accidents caused by the fault of the bearing, and has important practicability and engineering application value.
Drawings
The invention is further illustrated with reference to the following figures and examples.
FIG. 1 is a flow chart of a filter-based feature enhancement method of the present invention;
FIG. 2 is a diagram of an original signal to be processed containing noise interference according to the present invention;
FIG. 3 is a signal after processing with lifting wavelets according to the present invention;
fig. 4 shows the final processed signal of the filter-based feature enhancement method of the present invention.
Detailed Description
The present invention will now be described in further detail with reference to the accompanying drawings. These drawings are simplified schematic views illustrating only the basic structure of the present invention in a schematic manner, and thus show only the constitution related to the present invention.
As shown in fig. 1, the method for enhancing features based on a filter provided by the present invention specifically includes the following steps:
s1, preprocessing the acceleration signal containing noise by adopting lifting wavelet transform:
obtaining a signal s (t) with n data points by using an exponential decay function containing impact characteristics:
in the present embodiment, the displacement constant a is 208.8, the damping ratio ξ is 0.055, and the natural frequency w is set to 208.8n4712, total points n 78644 at sampling frequency fsTotal time t of 65536 ═ 1.2 s; then adding a large amount of noise into the s (n) signal, wherein the noise submerges the impact characteristics s (n) in the signal, so that the impact characteristics are not clearly visible any more, and the signal after noise addition is shown in FIG. 2; processing a signal to be processed containing noise by lifting wavelet transform, wherein a threshold value is set as a soft threshold value, a threshold value coefficient is set as 2.5, and an update length of a predictor is set as 20, so as to obtain a new signal f (n), as shown in fig. 3; comparing fig. 2 and 3, it can be seen that fig. 3 has revealed four impulses, but the impulses are not obvious because there is still much noise interference in the signal; the purpose of this step is to remove noise from the impulse portion of the signal.
S2, searching all peak points and all valley points in the new signal f (n):
finding all peak points in the new signal f (n), and recording the found peak points as A ═ A { (A)pFind all valley points in the new signal f (n), and mark it as B ═ B { (B) }hH is 1,2, …, and then the interval I between each adjacent peak points is calculatedA={iA│iA=Ap+1,x-Ap,xCalculating the interval G between each adjacent valley pointB={gB│gB=Bh+1,x-Bh,xIn which A isp,xFor the peak value A in the signalpPoint correspondenceHorizontal coordinate value of (B)h,xAs the valley B in the signalhThe corresponding horizontal coordinate value.
S3, constructing a signal library:
calculating the peak point interval IA={iA│iA=Ap+1,x-Ap,xAnd valley point spacing GB={gB│gB=Bh+1,x-Bh,xMaximum value of max (I)A,GB) Selecting the size of C1(C1≥max(IA,GB) A linear type structural element difference filtering operator carries out filtering processing on the vibration signals f (n), and the processed result is recorded as fC1(n); then using the scale C2(C2≥max(IA,GB) A linear type structural element difference filtering operator carries out filtering processing on the vibration signals f (n), and the processed result is recorded as fC2(n); and analogizing in sequence, after m times of processing, establishing all processing results into a signal library { fC1(n),fC2(n),fC3(n),...,fCm(n)}。
S4, extracting coincident impact characteristics from the constructed signal library:
in signal library { fC1(n),fC2(n),fC3(n),...,fCm(n) } extraction of the coincident impact portion z (w):
z(w)=fC1{find(fC1(n)=fC2(n)=...=fCm(n))} (2)
the impact after extraction is shown in fig. 4. As can be seen by comparing fig. 2, 3 and 4, 4 very distinct impact signals are directly extracted from fig. 4, while no impact signal is visible in fig. 2 and no impact signal is evident in fig. 3; in addition, the amplitude in fig. 4 reaches about 3, and the amplitude in fig. 3 is only about 1, so that the impact amplitude is obviously increased.
The characteristic enhancement method based on the filter is based on the mathematical principle of the morphological algorithm, compared with the original signal and the noise-containing ratio processed by adopting the lifting wavelet transform, the method can not only realize noise reduction, but also well enhance the impact characteristic of the fault, provide good basic support for fault diagnosis of the bearing, including qualitative diagnosis, quantitative diagnosis, outer ring positioning diagnosis and the like, prevent major accidents caused by the fault of the bearing, and has important practicability and engineering application value.
In light of the foregoing description of the preferred embodiment of the present invention, many modifications and variations will be apparent to those skilled in the art without departing from the spirit and scope of the invention. The technical scope of the present invention is not limited to the content of the specification, and must be determined according to the scope of the claims.
Claims (2)
1. A method for filter-based feature enhancement, comprising the steps of:
s1, preprocessing the acceleration signal containing noise by adopting lifting wavelet transform:
acquiring a signal s (t) with a data point n by using an exponential decay function containing impact characteristics:
wherein A is a displacement constant, xi is a damping ratio, wnT is n points at a natural frequency of vibration and f is a sampling frequencysThe relationship between t and f is (1: n)/fs(ii) a Adding noise into the s (n) signal, wherein the total amplitude of the noise is larger than that of the impact characteristic, and then preprocessing the signal to be processed containing the noise by adopting a lifting wavelet transform technology to obtain a new signal f (n);
s2, searching all peak points and all valley points in the new signal f (n):
finding all peak points in the new signal f (n), and recording the found peak points as A ═ A { (A)pFind all valley points in the new signal f (n), and mark it as B ═ B { (B) }hH is 1,2, …, and then the interval I between each adjacent peak points is calculatedA={iA│iA=Ap+1,x-Ap,xCalculating the interval G between each adjacent valley pointB={gB│gB=Bh+1,x-Bh,xIn which A isp,xFor the peak value A in the signalpCorresponding horizontal coordinate value, Bh,xAs the valley B in the signalhCorresponding horizontal coordinate values;
s3, constructing a signal library:
calculating the peak point interval IA={iA│iA=Ap+1,x-Ap,xAnd valley point spacing GB={gB│gB=Bh+1,x-Bh,xMaximum value of max (I)A,GB) Selecting the size of C1(C1≥max(IA,GB) A linear type structural element difference filtering operator carries out filtering processing on the vibration signals f (n), and the processed result is recorded as fC1(n); then using the scale C2(C2≥max(IA,GB) A linear type structural element difference filtering operator carries out filtering processing on the vibration signals f (n), and the processed result is recorded as fC2(n); building all processing results into a signal library { f) until m times of processingC1(n),fC2(n),fC3(n),...,fCm(n)};
S4, extracting coincident impact characteristics from the constructed signal library:
in signal library { fC1(n),fC2(n),fC3(n),...,fCm(n) } extraction of the coincident impact portion z (w):
z(w)=fC1{find(fC1(n)=fC2(n)=...=fCm(n))} (2)。
2. the method of claim 1, wherein in step S3, the linear-type structure element difference filter operator is an operation form of morphological filtering.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111624531.4A CN114254679B (en) | 2021-12-28 | 2021-12-28 | Feature enhancement method based on filter |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111624531.4A CN114254679B (en) | 2021-12-28 | 2021-12-28 | Feature enhancement method based on filter |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114254679A true CN114254679A (en) | 2022-03-29 |
CN114254679B CN114254679B (en) | 2024-06-04 |
Family
ID=80798452
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111624531.4A Active CN114254679B (en) | 2021-12-28 | 2021-12-28 | Feature enhancement method based on filter |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114254679B (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102269644A (en) * | 2010-06-07 | 2011-12-07 | 北京化工大学 | Diagnosis method for impact type failure between rolling bearing and gear based on optimal self-adaptive wavelet filter |
US20120265534A1 (en) * | 2009-09-04 | 2012-10-18 | Svox Ag | Speech Enhancement Techniques on the Power Spectrum |
CN111623968A (en) * | 2020-05-08 | 2020-09-04 | 安徽智寰科技有限公司 | Fault feature extraction method based on adaptive morphological filtering |
CN113269048A (en) * | 2021-04-29 | 2021-08-17 | 北京工业大学 | Motor imagery electroencephalogram signal classification method based on deep learning and mixed noise data enhancement |
-
2021
- 2021-12-28 CN CN202111624531.4A patent/CN114254679B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120265534A1 (en) * | 2009-09-04 | 2012-10-18 | Svox Ag | Speech Enhancement Techniques on the Power Spectrum |
CN102269644A (en) * | 2010-06-07 | 2011-12-07 | 北京化工大学 | Diagnosis method for impact type failure between rolling bearing and gear based on optimal self-adaptive wavelet filter |
CN111623968A (en) * | 2020-05-08 | 2020-09-04 | 安徽智寰科技有限公司 | Fault feature extraction method based on adaptive morphological filtering |
CN113269048A (en) * | 2021-04-29 | 2021-08-17 | 北京工业大学 | Motor imagery electroencephalogram signal classification method based on deep learning and mixed noise data enhancement |
Also Published As
Publication number | Publication date |
---|---|
CN114254679B (en) | 2024-06-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8073634B2 (en) | Method to extract target signals of a known type from raw data containing an unknown number of target signals, interference, and noise | |
Hao et al. | Morphological undecimated wavelet decomposition for fault diagnostics of rolling element bearings | |
CN108932699B (en) | Three-dimensional matching harmonic filtering image denoising method based on transform domain | |
CN108645620B (en) | Rolling bearing early fault diagnosis method based on information entropy and multi-scale morphology | |
CN110646851B (en) | Adaptive threshold seismic random noise suppression method based on Shearlet transformation | |
Meng et al. | Self-supervised learning for seismic data reconstruction and denoising | |
CN108961181B (en) | Shearlet transform-based ground penetrating radar image denoising method | |
CN109633761B (en) | Magnetic resonance signal power frequency noise reduction method based on wavelet transformation modulus maximum value method | |
CN115730199A (en) | Method and system for noise reduction and fault feature extraction of vibration signal of rolling bearing | |
CN116153329A (en) | CWT-LBP-based sound signal time-frequency texture feature extraction method | |
CN112037151A (en) | Image denoising method based on wavelet analysis | |
CN116992217A (en) | Mechanical signal noise reduction method based on multi-scale dynamic weighting multi-dimensional residual convolution | |
CN111678691A (en) | Gear fault detection method based on improved sparse decomposition algorithm | |
CN114705427A (en) | Rolling bearing intelligent fault diagnosis method under strong noise environment | |
CN114595728A (en) | Signal denoising method based on self-supervision learning | |
CN110657807A (en) | Indoor positioning displacement measurement method for detecting discontinuity based on wavelet transformation | |
CN103455986B (en) | Random noise point detecting method based on fractional order differential gradient | |
CN114254679B (en) | Feature enhancement method based on filter | |
CN116594061A (en) | Seismic data denoising method based on multi-scale U-shaped attention network | |
CN116296243A (en) | Pneumatic identification method based on large-size nuclear dense blocks | |
CN107341519B (en) | Support vector machine identification optimization method based on multi-resolution analysis | |
CN115828069A (en) | End-to-end magnetic anomaly signal noise reduction method based on deep learning | |
CN115561817A (en) | Desert earthquake denoising method based on multiple attention mechanism | |
Pang | Improved image denoising based on Haar wavelet transform | |
CN114858430A (en) | Mechanical seal acoustic emission signal noise reduction method based on LCD-wavelet new threshold |
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 |