CN109948286A - Based on the signal decomposition method for improving experience wavelet decomposition - Google Patents
Based on the signal decomposition method for improving experience wavelet decomposition Download PDFInfo
- Publication number
- CN109948286A CN109948286A CN201910253086.1A CN201910253086A CN109948286A CN 109948286 A CN109948286 A CN 109948286A CN 201910253086 A CN201910253086 A CN 201910253086A CN 109948286 A CN109948286 A CN 109948286A
- Authority
- CN
- China
- Prior art keywords
- decomposition
- signal
- fault
- result
- spectrum
- 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
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Complex Calculations (AREA)
Abstract
The present invention relates to a kind of based on the signal decomposition method for improving experience wavelet decomposition, it include: to carry out Fourier power spectral density spectral line to fault-signal to seek, it chooses several different size threshold values and rejecting processing is carried out to above-mentioned interference component spectrum, and fault-signal is subjected to resolution process based on New cultivars, utilize fault signature energy ratio information screening optimal threshold, optimal Decomposition result and modal components containing abundant fault characteristic information.The present invention is under jamming pattern, effectively optimal Decomposition is remained to go out the modal components for containing abundant fault characteristic information, eliminate interference, inhibit modal overlap and excessive decomposition phenomenon, desirable broken-down is obtained as a result, compensating for the deficiency of EWT, enriches the theoretical method of mode decomposition.
Description
Technical field
It is specifically a kind of based on improving experience wavelet decomposition the present invention relates to the fault-signal processing method of rotating machinery
Signal decomposition method.
Background technique
The rotating machineries such as hydraulic pump, hydraulic motor, motor, bearing, gear and rotor are widely used in each important
Industrial circle, the working environment that they are faced also tend to the bad working environments such as high temperature, high pressure, high speed, heavy duty, accelerate rotation
The deterioration for tool health status of making a connection, therefore the intelligent trouble diagnosis of rotating machinery is had a very important significance.Rotating machinery
The vibration mode and transmission path of failure are sufficiently complex, fault-signal have the characteristics that it is non-linear and non-stationary, and be easy quilt
Equipment introduces the interference such as noise, ambient noise and electromagnetism.Therefore, how to efficiently extract containing fault characteristic information mould abundant
State component inhibits interference just to become a critical issue.
The shortcomings that wavelet transformation are as follows: none principle that can be followed and criterion come select suitable wavelet basis function with
Signal aspect feature carries out accurate match;Once it is determined that basic function, during decomposition it will not be according to the shape of signal
State feature and change, result in it and do not have adaptivity truly;Again, although small echo has " school microscop "
Function, once it is determined that its resolution ratio also just determines therewith after basic function and scale factor;Finally, with Fourier in short-term
The limitation of transformation is the same, although wavelet transformation has multiple dimensioned more resolution capabilities, requires the signal in small echo window necessary
It is approximate steady (pseudo- steady).
Traditional empirical mode decomposition disadvantage are as follows: such as modal overlap, it cannot effectively be based on signal aspect
Feature separates the mode function of specific time scale, so that different modalities component appears in same decomposition result, either
The same modal components are decomposed in multiple decomposition results;For example end effect, the data edges of modal components will appear hair
Phenomenon is dissipated, and in an iterative process, this edge diverging gradually can inwardly pollute, and with the increase of the number of iterations, can make
Data sequence serious distortion is obtained, modal overlap and chaff component occurs.
In view of the above-mentioned problems, Gilles proposed in paper " Empirical Wavelet Transform " it is a kind of non-thread
The New Method for Processing of property and non-stationary signal, i.e. experience wavelet decomposition (Empirical Wavelet Transform, EWT), together
When its advantages of having merged both wavelet transformation and empirical mode decomposition.EWT passes through segmentation Fourier amplitude spectrum, and each
Wavelet orthogonal basis is established in segmentation section, can be several amplitude modulation-with compactly support frequency spectrum by a multi-modal signal decomposition
The sum of FM signal.But under jamming pattern, it can not effectively be decomposited containing abundant fault characteristic information modal components, pole
It is possible that there is modal overlap and excessive decomposition phenomenon.
Summary of the invention
The present invention is directed under jamming pattern, can not effectively be decomposed through EWT containing abundant fault characteristic information mode
Problem provides a kind of based on improvement experience wavelet decomposition (Improved Empirical Wavelet Transform, IEWT)
Signal decomposition method, this method use condition is broader, and discomposing effect becomes apparent from.
The technical solution used to solve the technical problems of the present invention is that:
A kind of signal decomposition method based on improvement experience wavelet decomposition, includes the following steps:
(1) power density spectrum is sought
The calculating of power density spectrum is carried out to fault-signal;
(2) threshold value based on power density spectrum is rejected
The spectrum for being less than the threshold value in the power density spectrum of fault-signal is rejected using different threshold values, obtains new spectrum sequence
Column;
(3) fault-signal resolution process
N number of continuum is divided into [0 π] range to power density spectrum, establishes wavelet orthogonal basis in each section, it will
Fault-signal is decomposed into the sum of N number of modal components;
(4) optimal Decomposition result is screened
The characteristic energy ratio of each component in decomposition result corresponding to each threshold value is solved, and will be in each decomposition result
Maximum characteristic energy ratio compared and analyzed with time big characteristic energy ratio, and the maximum contrast value based on all decomposition results
In, maximum contrast is filtered out as a result, decomposition result corresponding to the maximum contrast value is optimal Decomposition as a result, corresponding threshold
Value is then known as optimal Decomposition threshold value, and maximum characteristic energy contains most abundant fault characteristic information than corresponding modal components.
The present invention by adopting the above technical scheme, compared with prior art EWT, beneficial effect is:
Under jamming pattern, remains to effectively optimal Decomposition and go out the modal components for containing abundant fault characteristic information, eliminate
Interference inhibits modal overlap and excessive decomposition phenomenon, obtains desirable broken-down result.The deficiency for compensating for EWT enriches mode point
The theoretical method of solution.
Detailed description of the invention
Fig. 1 is piston shoes wear-out failure signal time-domain diagram;
Fig. 2 is in the piston shoes wear-out failure signal decomposition result obtained based on IEWT method containing most abundant fault signature letter
Cease component time-domain diagram;
Fig. 3 is in the piston shoes wear-out failure signal decomposition result obtained based on IEWT method containing most abundant fault signature letter
Cease component power spectrum density figure;
Fig. 4 is in the piston shoes wear-out failure signal decomposition result obtained based on EWT method containing most abundant fault signature letter
Cease component time-domain diagram;
Fig. 5 is in the piston shoes wear-out failure signal decomposition result obtained based on EWT method containing most abundant fault signature letter
Cease component power spectrum density figure;
Fig. 6 is flow chart of the invention.
Specific embodiment
The present invention is further described with reference to the accompanying drawings and embodiments.
The present embodiment illustrates by taking piston shoes wear-out failure as an example.
Referring to Fig. 6, one kind being based on improvement experience wavelet decomposition (Improved Empirical Wavelet
Transform, IEWT) signal decomposition method, the specific steps are as follows:
Step (1), power density spectrum is sought
The calculating of power density spectrum is carried out to collected discrete state fault-signal f (n), spectrum distribution series are denoted as
P。
Step (2), the threshold value based on power density spectrum are rejected
The spectrum for being less than the threshold value in P is picked respectively using L different size of threshold value coefficient × mean (P)
It removes, obtains L new spectrum distribution series Pcoefficient, wherein coefficient=1, integer, the mean (P) of 2 ..., L is flat
Equal spectrum.Since interference spectrum is generally little, so the step largely eliminates interference component spectrum information to boundary
The adverse effect of adaptivenon-uniform sampling, to avoid modal overlap and excessive decomposition phenomenon on as far as possible.
Step (3), fault-signal resolution process
For PcoefficientIt is handled as follows: N number of continuum is divided into [0 π] range to power spectral density spectrum,
Thus N+1 boundary line is generated, it is assumed that signal has N number of mode to constitute, the maximum in spectrogram is screened, by M maximum
Descending arrangement is carried out, there can be two following situations:
1) M >=N illustrates that this method has filtered out sufficient maximum, and retains top n maximum;
2) M < N illustrates that the single mode quantity in signal is less than expected number of components N, retains all maximum, lay equal stress on
Set N.
Either segment is represented by Λn=[ωn-1,ωn],Wherein ωnIt is two
Midpoint (the ω of a continuous maximumn=0, ωN=π).Therefore by center ωn, width Tn=2 τnForm a changeover portion.
To ΛnWindowing process is carried out, is based on Meyer wavelet method, experience scaling function and experience wavelet function are respectively such as
Shown in lower.
Wherein, γ < minn[ωn+1-ωn/ωn+1+ωn], β (x)=35x4-8x5+70x6-20x7。
According to wavelet theory, formula (1) and (2), the wavelet transformation detail coefficients W of signalx(n, t) and Coefficients of Approximation Wx(0,
T) formula (3) and formula (4) be may be defined as:
Signal is reconstructed, as a result as shown in formula (5).
Modal components f is obtained by above-mentionedk, can be indicated by following formula (6) and (7)
Step (4), the screening of optimal Decomposition result
It is obtained in L decomposition result based on above-mentioned, each and every one several modal components is contained in every kind of decomposition result, to every kind point
Each modal components of solution result seek feature energy energy ratio FER (energy and and frequency band at fault characteristic frequency and its frequency multiplication
The ratio of gross energy), and by the maximum characteristic energy ratio FER of each decomposition resultcoefficient,maxWith secondary big characteristic energy ratio
FERcoefficient,secondmaxIt compares and analyzes, i.e. Acoefficient=(FERcoefficient,max-FERcoefficient,secondmax)/
FERcoefficient,secondmax, in certain larger FERcoefficient,maxMaximum A is found out in rangecoefficient, it is denoted as Amax, Amax
Corresponding decomposition result is optimal Decomposition as a result, corresponding threshold value is known as optimal threshold, the FER of the resultcoefficient,max
Corresponding component is the component for containing abundant fault characteristic information.
Referring to fig. 2, Fig. 3, Fig. 4, Fig. 5, the present embodiment are based on IEWT method and prior art EWT method and grind respectively to piston shoes
Damage fault-signal is decomposed, and 8 and 53 modal components are respectively obtained, and is chosen in optimal Decomposition result and is contained most abundant event
Hinder characteristic information modal components.Clearly it is in based on the obtained Fig. 2 containing most abundant fault characteristic information modal components of IEWT
Periodic impulse characteristic information is showed, although also clear based on Fig. 4 containing most abundant fault characteristic information modal components that EWT is obtained
Periodic impulse characteristic information is presented clearly, but its spectrum is much smaller compared with obtaining based on IEWT;Contained based on what IEWT was obtained
There is the piston shoes wear-out failure original signal in the information from objective pattern and Fig. 1 of Fig. 2 of most abundant fault characteristic information modal components
Information from objective pattern has very high similitude, and contains most abundant fault characteristic information modal components compared with what is obtained based on EWT
Fig. 4 information from objective pattern and Fig. 1 in piston shoes wear-out failure original signal information from objective pattern similitude it is much higher.Base
In the modal components result figure 3 that IEWT optimal Decomposition obtains clearly present piston shoes wear-out failure characteristic frequency 171.5Hz and
Spectral line at its frequency multiplication, although being also clearly in based on the obtained Fig. 5 containing most abundant fault characteristic information modal components of EWT
The spectral line at piston shoes wear-out failure characteristic frequency 171.5Hz and its frequency multiplication is showed, but its spectrum is small compared with obtaining based on IEWT
Very much.
The above is only a specific embodiments of the invention, it is noted that although making referring to the embodiment to the present invention
It is described in detail, for those of ordinary skill in the art, several improvement and profit can be carried out to technical solution of the present invention
Decorations, but the spirit and scope of technical solution of the present invention is not departed from, these improvements and modifications are also considered as the protection model of the invention patent
It encloses.
Claims (1)
1. a kind of based on the signal decomposition method for improving experience wavelet decomposition, which comprises the steps of:
(1) power density spectrum is sought
The calculating of power density spectrum is carried out to fault-signal;
(2) threshold value based on power density spectrum is rejected
The spectrum for being less than the threshold value in the power density spectrum of fault-signal is rejected using different threshold values, obtains new spectrum sequence;
(3) fault-signal resolution process
N number of continuum is divided into [0 π] range to power density spectrum, wavelet orthogonal basis is established in each section, by failure
Signal decomposition is the sum of N number of modal components;
(4) optimal Decomposition result is screened
Solve the characteristic energy ratio of each component in decomposition result corresponding to each threshold value, and by each decomposition result most
Big characteristic energy ratio is compared and analyzed with time big characteristic energy ratio, and in the maximum contrast value based on all decomposition results, sieve
Select maximum contrast as a result, decomposition result corresponding to the maximum contrast value be optimal Decomposition as a result, corresponding threshold value then
Referred to as optimal Decomposition threshold value, maximum characteristic energy contain most abundant fault characteristic information than corresponding modal components.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910253086.1A CN109948286B (en) | 2019-03-29 | 2019-03-29 | Signal decomposition method based on improved empirical wavelet decomposition |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910253086.1A CN109948286B (en) | 2019-03-29 | 2019-03-29 | Signal decomposition method based on improved empirical wavelet decomposition |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109948286A true CN109948286A (en) | 2019-06-28 |
CN109948286B CN109948286B (en) | 2023-10-03 |
Family
ID=67013155
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910253086.1A Active CN109948286B (en) | 2019-03-29 | 2019-03-29 | Signal decomposition method based on improved empirical wavelet decomposition |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109948286B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110441654A (en) * | 2019-07-29 | 2019-11-12 | 三峡大学 | Based on the distribution network electric energy quality disturbance detecting method for improving EWT and CMPE |
CN110795843A (en) * | 2019-10-24 | 2020-02-14 | 北京建筑大学 | Method and device for identifying faults of rolling bearing |
CN113111618A (en) * | 2021-03-09 | 2021-07-13 | 电子科技大学 | Analog circuit fault diagnosis method based on improved empirical wavelet transform |
CN113554103A (en) * | 2021-07-28 | 2021-10-26 | 大连海天兴业科技有限公司 | Fault diagnosis algorithm for rolling bearing of train running gear |
CN115146672A (en) * | 2022-06-20 | 2022-10-04 | 中国人民解放军96963部队 | Dense-frequency modal separation reconstruction method and device |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104776908A (en) * | 2015-04-17 | 2015-07-15 | 南京理工大学 | EMD generalized energy-based wheeltrack vibration signal fault feature extraction method |
CN105004923A (en) * | 2015-07-10 | 2015-10-28 | 湘潭大学 | Magnetic control submerged-arc welding seam tracking signal analyzing method based on experience wavelet transformation |
CN105044478A (en) * | 2015-07-23 | 2015-11-11 | 国家电网公司 | Transmission line audible noise multi-channel signal extraction method |
WO2017144007A1 (en) * | 2016-02-25 | 2017-08-31 | 深圳创维数字技术有限公司 | Method and system for audio recognition based on empirical mode decomposition |
CN108375472A (en) * | 2018-02-12 | 2018-08-07 | 武汉科技大学 | Based on the Method for Bearing Fault Diagnosis and system and device for improving experience wavelet transformation |
CN108458871A (en) * | 2017-12-22 | 2018-08-28 | 浙江工业大学 | A kind of gearbox fault recognition methods based on improvement experience wavelet transformation |
CN108760327A (en) * | 2018-08-02 | 2018-11-06 | 南昌航空大学 | A kind of diagnostic method of aeroengine rotor failure |
CN109211568A (en) * | 2018-09-19 | 2019-01-15 | 四川大学 | Fault Diagnosis of Roller Bearings based on condition experience wavelet transformation |
-
2019
- 2019-03-29 CN CN201910253086.1A patent/CN109948286B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104776908A (en) * | 2015-04-17 | 2015-07-15 | 南京理工大学 | EMD generalized energy-based wheeltrack vibration signal fault feature extraction method |
CN105004923A (en) * | 2015-07-10 | 2015-10-28 | 湘潭大学 | Magnetic control submerged-arc welding seam tracking signal analyzing method based on experience wavelet transformation |
CN105044478A (en) * | 2015-07-23 | 2015-11-11 | 国家电网公司 | Transmission line audible noise multi-channel signal extraction method |
WO2017144007A1 (en) * | 2016-02-25 | 2017-08-31 | 深圳创维数字技术有限公司 | Method and system for audio recognition based on empirical mode decomposition |
CN108458871A (en) * | 2017-12-22 | 2018-08-28 | 浙江工业大学 | A kind of gearbox fault recognition methods based on improvement experience wavelet transformation |
CN108375472A (en) * | 2018-02-12 | 2018-08-07 | 武汉科技大学 | Based on the Method for Bearing Fault Diagnosis and system and device for improving experience wavelet transformation |
CN108760327A (en) * | 2018-08-02 | 2018-11-06 | 南昌航空大学 | A kind of diagnostic method of aeroengine rotor failure |
CN109211568A (en) * | 2018-09-19 | 2019-01-15 | 四川大学 | Fault Diagnosis of Roller Bearings based on condition experience wavelet transformation |
Non-Patent Citations (3)
Title |
---|
孟宗;李姗姗: "小波改进阈值去噪和经验模态分解相结合的旋转机械故障特征提取", vol. 36, no. 01, pages 24 - 29 * |
李力 等: "基于改进小波去噪预处理和EEMD的采煤机齿轮箱故障诊断", vol. 47, no. 10, pages 3394 - 3400 * |
贾广飞;孙师泽;武哲;: "经验模态分解在切削振动信号分析中的应用", 河北工业科技, no. 03, pages 67 - 71 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110441654A (en) * | 2019-07-29 | 2019-11-12 | 三峡大学 | Based on the distribution network electric energy quality disturbance detecting method for improving EWT and CMPE |
CN110441654B (en) * | 2019-07-29 | 2021-05-04 | 三峡大学 | Power quality disturbance detection method for power distribution network based on improved EWT and CMPE |
CN110795843A (en) * | 2019-10-24 | 2020-02-14 | 北京建筑大学 | Method and device for identifying faults of rolling bearing |
CN110795843B (en) * | 2019-10-24 | 2024-03-29 | 北京建筑大学 | Method and device for identifying faults of rolling bearing |
CN113111618A (en) * | 2021-03-09 | 2021-07-13 | 电子科技大学 | Analog circuit fault diagnosis method based on improved empirical wavelet transform |
CN113111618B (en) * | 2021-03-09 | 2022-10-18 | 电子科技大学 | Analog circuit fault diagnosis method based on improved empirical wavelet transform |
CN113554103A (en) * | 2021-07-28 | 2021-10-26 | 大连海天兴业科技有限公司 | Fault diagnosis algorithm for rolling bearing of train running gear |
CN115146672A (en) * | 2022-06-20 | 2022-10-04 | 中国人民解放军96963部队 | Dense-frequency modal separation reconstruction method and device |
Also Published As
Publication number | Publication date |
---|---|
CN109948286B (en) | 2023-10-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109948286A (en) | Based on the signal decomposition method for improving experience wavelet decomposition | |
CN111623982B (en) | Planetary gearbox early fault diagnosis method based on APEWT and IMOMEDA | |
Lei et al. | A review on empirical mode decomposition in fault diagnosis of rotating machinery | |
Wang et al. | A new image denoising method based on Gaussian filter | |
Dong et al. | Non-iterative denoising algorithm for mechanical vibration signal using spectral graph wavelet transform and detrended fluctuation analysis | |
CN110276726B (en) | Image deblurring method based on multichannel network prior information guidance | |
US8594456B2 (en) | Image denoising method | |
CN110081967A (en) | Mechanical oscillation signal Threshold Noise Reduction Methods based on spectrogram wavelet transformation | |
CN105651376A (en) | Mechanical equipment offline detection system vibration signal spectrum analysis and alarm method | |
Selesnick | TQWT toolbox guide | |
CN114970602A (en) | Signal denoising method and system based on improved empirical mode decomposition and wavelet threshold function | |
CN110084208B (en) | Calculation order tracking method capable of adaptively reducing noise and avoiding order aliasing | |
Zhang et al. | Bearing fault diagnosis under various operation conditions using synchrosqueezing transform and improved two-dimensional convolutional neural network | |
CN114325598A (en) | Signal denoising method and device | |
CN109359633B (en) | Signal joint classification method based on Hilbert-Huang transform and wavelet ridge line | |
CN112487882A (en) | Method for generating non-sparse index-guided enhanced envelope spectrum based on spectrum coherence | |
CN111881848A (en) | Motor fault signal extraction method based on variational modal decomposition and improved particle swarm | |
CN112485028A (en) | Vibration signal characteristic frequency spectrum extraction method and mechanical fault diagnosis analysis method | |
CN115115864B (en) | Dynamic balance testing method and system for pump precision rotor shaft | |
CN116725542A (en) | Composite filtering method of NV color center magnetic field measurement signal | |
CN111241902A (en) | High-precision multiple synchronous compression generalized S-transform time-frequency analysis method | |
CN113657268B (en) | Signal automatic decomposition method applied to wind turbine generator gearbox fault diagnosis | |
Neelamani et al. | Wavelet-domain regularized deconvolution for ill-conditioned systems | |
CN114894478A (en) | Method for extracting weak fault features of rolling bearing | |
Wu et al. | Faults classification of a scooter engine platform using wavelet transform and artificial neural network |
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 |