CN105488341B - A kind of denoising method based on mixing empirical mode decomposition - Google Patents
A kind of denoising method based on mixing empirical mode decomposition Download PDFInfo
- Publication number
- CN105488341B CN105488341B CN201510844356.8A CN201510844356A CN105488341B CN 105488341 B CN105488341 B CN 105488341B CN 201510844356 A CN201510844356 A CN 201510844356A CN 105488341 B CN105488341 B CN 105488341B
- Authority
- CN
- China
- Prior art keywords
- imf
- noise
- signal
- amplitude
- value
- 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
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Soundproofing, Sound Blocking, And Sound Damping (AREA)
Abstract
The invention discloses one kind based on mixing empirical mode decomposition(Empirical Mode Decomposition, abbreviation EMD)Denoising method, belong to digital signal processing technique field.EEMD and MEMD algorithms are combined by the present invention, and carry out parameter optimization for the mask signal frequency used in MEMD algorithms, can more effectively eliminate modal overlap, denoising effect is more preferable;The present invention further in EEMD addition noise amplitude, the filtering threshold of noise correlation IMF, noise IMF filtering threshold targetedly optimized, so as to further improve global de-noising effect.Compared with prior art, the method for the present invention can more effectively eliminate modal overlap, and denoising effect is more preferable.
Description
Technical field
The present invention relates to a kind of denoising method more particularly to one kind based on mixing empirical mode decomposition (Empirical
Mode Decomposition, abbreviation EMD) denoising method, belong to digital signal processing technique field.
Background technology
(Ensemble Empirical Mode Decomposition gather empirical modal by EMD and the EEMD of derivation
Decompose) and masking signal EMD (mask signal method empirical mode decomposition, behind abbreviation MEMD in text) be in recent years
The signal processing method to grow up.It is decomposed based on signal time scale, is believed suitable for non-linear and unstable state
Number processing, due to not needing to determine the subjective experiences parameter setting such as basic function and Decomposition order, in some cases, compares wavelet transformation
With better discomposing effect.Nevertheless, EMD also has some problems and defect of oneself., mainly there is endpoint effect in the problem of EMD
It answers, modal overlap.For EEMD although EMD can be solved the problems, such as to a certain extent, it is uncertain to add in noise amplitude.EEMD
In practical applications, especially when signal-to-noise ratio is relatively low, modal overlap problem still remains, and the noise that EEMD is added at this time can be letter
Number noise decomposed as useful signal, so as to bring modal overlap and more falseness IMF (Spurious
modes).In addition, the amplitude selection of the ensemble noise of EEMD is related with the signal-to-noise ratio of signal, noise amplitude is accurately selected
There is certain difficulty, many experiments show when signal-to-noise ratio is low, and EMD decomposition will be advantageously.MEMD is determined by adding
High frequency mask signal solve the problems, such as extreme's distribution, therefore for modal overlap, MEMD has apparent advantage and meter
Calculate efficiency.Amplitude and the frequency parameter selection of the mask signal of MEMD are a committed steps, although giving reference value in text,
But in practical applications, reference frequency is often below the required requirements of masking signal, thus cannot reach elimination mould
The ideal effect of state aliasing, therefore, for different signals, the result of the parameter calculating of recommendation is simultaneously non-optimal.IMF filtering is main
There are wavelet filter, EMD directly to filter (EMD-DT) and EMD sections filtering (EMD-IT) and related derivative filtering method,
Middle small echo threshold coefficient error is larger, is not easy to find accurately in practice;EMD-DT cannot be guaranteed the continuity of signal, for more
The signal filtering signal of a amplitude range is bad;Although EMD-IT can guarantee the continuity of signal, but cannot guarantee that low order
The filtering performance of several value signals of IMF.To sum up, feature there are two EMD-IT filtering mainly, first, find relative
After (noise is related) IMF, " noise " IMF of low order is directly abandoned, does not account for the active ingredient wherein carried;Second, filter
Wave threshold value is not adaptive, and there are large errors for the filtering of several value signals for single threshold strategies.
Invention content
The technical problem to be solved by the present invention is to overcome the deficiencies of the prior art and to provide one kind based on mixing empirical modal
EEMD and MEMD algorithms are combined by the denoising method of decomposition, this method, and for the mask signal used in MEMD algorithms
Frequency carries out parameter optimization, can more effectively eliminate modal overlap, denoising effect is more preferable.
Technical solution of the present invention is specific as follows:
A kind of denoising method based on mixing empirical mode decomposition, includes the following steps:
Step 1 decomposes original signal using gathering empirical mode decomposition EEMD methods, obtains a series of natural modes
State component IMF and discrepance, and noise IMF, noise correlation IMF, useful IMF are determined from all IMF;
Step 2, to noise IMF, IMF related to noise carries out threshold filter denoising, and using mask signal method to useful
The signal that IMF is formed carries out the elimination of modal overlap;Wherein, mask signal frequency used in mask signal method is according to following
Method determines:The signal formed first to useful IMF carries out first derivation, obtains first derivative signal;Then Xi Er is utilized
Bert transform method obtains the instantaneous frequency of the first derivative signal, then asks for the mean value of instantaneous frequency;Finally to described
The mean value of instantaneous frequency is multiplied by the weights that a value range is (1,2), and gained product is mask signal frequency;
Step 3, the noise IMF to after threshold filter denoising, noise correlation IMF eliminate the useful IMF institutes structure of modal overlap
Into signal and discrepance carry out signal reconstruction.
In order to improve the performance of EEMD, present invention further proposes following improvement project, to the addition noise in EEMD
Amplitude optimizes:
The amplitude of addition noise used in the EEMD methods predefines by the following method:First to original letter
It number is filtered, obtains filtered signal, and add the noise of different amplitudes respectively to original signal, obtain one group and add letter of making an uproar
Number;Then following handle is carried out respectively to each noisy signal:Empirical mode decomposition is carried out, and from experience to the noisy signal
Useful IMF is determined in the obtained IMF of mode decomposition;The related root mean square calculated between each useful IMF and filtered signal misses
The summation of related coefficient between the summation of difference and every two neighboring useful IMF;Finally made an uproar with minimum the adding of the sum of two summations
Amplitude of the amplitude of addition noise corresponding to signal as the addition noise used in EEMD methods.
In order to improve the threshold filter effect to noise correlation IMF, present invention further proposes following improvement projects:
Filtering threshold used in carrying out threshold filter denoising to noise correlation IMF determines in accordance with the following methods:First will
Each noise correlation IMF carries out decorrelative transformation with noise IMF respectively;According to amplitude is positive and negative will be after each decorrelative transformation
IMF points of noise correlation is two parts, obtains its amplitude mean value and energy filter threshold value respectively to each section, and the two is added
Power is summed to get the filtering threshold of the part.
In order to more effectively extract the useful information in noise IMF, present invention further propose that following improvement project:
Filtering threshold used in carrying out threshold filter denoising to noise IMF determines in accordance with the following methods:Calculate noise IMF
Amplitude mean value, and in noise IMF be more than amplitude mean value each amplitude carry out time-sequencing, obtain the time of the amplitude
Sequence;Periodic authentication is carried out respectively for the time series of obtained each amplitude, and therefrom selects periodically higher portion
Divide time series;Finally using the minimum value of amplitude corresponding to selected part-time sequence as the filtering threshold.
Compared with prior art, it the present invention and its is further improved technical solution and has the advantages that:
EEMD and MEMD algorithms of the present invention are combined, and are joined for the mask signal frequency used in MEMD algorithms
Number optimization, can more effectively eliminate modal overlap, denoising effect is more preferable;
The present invention is optimized further directed to the addition noise amplitude in EEMD, the Relative RMSE of useful IMF
The minimum selection criteria as EEMD addition noise amplitudes of the sum of value and related coefficient, so as to more effectively solve end effect
Caused problem;
The present invention is optimized further directed to the filtering threshold of noise correlation IMF, and plan is separated using positive and negative signal amplitude
Slightly, by decorrelation, take mean value and last filtering threshold determined according to the aspect of energy three, the filter to noise correlation IMF
Wave effect is more preferable;
The present invention further using the useful information in the method for periodic authentication extraction noise IMF, is avoided and is directly lost
The loss of useful information, further improves denoising effect caused by abandoning.
Description of the drawings
Fig. 1 is the flow diagram for determining EEMD addition noise amplitude methods;
Fig. 2 is the flow diagram of MEMD mask signal frequency determination methods;
The filtering threshold that Fig. 3 is noise correlation IMF determines method flow schematic diagram.
Specific embodiment
Technical scheme of the present invention is described in detail below in conjunction with the accompanying drawings:
The basic ideas of the present invention are to be combined EEMD and MEMD algorithms, and for the mask used in MEMD algorithms
Signal frequency carries out parameter optimization, can more effectively eliminate modal overlap, denoising effect is more preferable;The present invention is further in EEMD
Addition noise amplitude, the filtering threshold of noise correlation IMF, noise IMF filtering threshold targetedly optimized, thus into
One step improves global de-noising effect.
For the ease of public understanding, technical solution of the present invention is described in detail with a preferred embodiment below.
The present invention is based on the denoising methods of mixing empirical mode decomposition, specifically include following steps:
Step 1 decomposes original signal using gathering empirical mode decomposition EEMD methods, obtains a series of natural modes
State component IMF and discrepance, and noise IMF, noise correlation IMF, useful IMF are determined from all IMF.
Traditional EEMD methods can be used in this step to the decomposition of original signal or improve EEMD methods.EEMD is existing
Technology, detailed content can be found in document [Wu.Zhaohua, Huang N E, Ensemble empirical mode
decomposition:a noise assisted data analysis method,Adv.Adapt.DataAnal.1(1)
(2009)1–41.]。
As a kind of innovatory algorithm of EMD, EEMD although EMD can be solved the problems, such as to a certain extent, make an uproar by addition
Acoustic amplitude is uncertain.In practical applications, especially when signal-to-noise ratio is relatively low, modal overlap problem still remains EEMD, this
When the noises that add in of EEMD the noise of signal can be decomposed as useful signal, so as to bring modal overlap and more
False IMF (Spurious modes).In addition, the amplitude selection of the ensemble noise of EEMD and the signal-to-noise ratio of signal have
It closes, it is accurate that noise amplitude is selected to have certain difficulty.For this purpose, the present invention improves existing EEMD methods, make an uproar to addition
Acoustic amplitude optimizes, and the sum of the Relative RMSE values of useful IMF and related coefficient is minimum as EEMD addition noise width
The selection criteria of value, so as to more effectively solve the problems, such as caused by end effect.Specifically, as shown in Figure 1, the EEMD
The amplitude of addition noise used in method predefines by the following method:
Rough filtering (for the sake of simplicity, it is preferred to use median filter method) is carried out to original signal first, is filtered
Signal after wave, and add the noise of different amplitudes respectively to original signal, obtain one group of noisy signal;Preferably, it adds and makes an uproar
The amplitude section of sound is 0.001-0.4 times of original signal standard deviation;
Then following handle is carried out respectively to each noisy signal:EMD decomposition is carried out, and from EMD to the noisy signal
It decomposes and useful IMF is determined in obtained IMF;Calculate the related root-mean-square error between each useful IMF and filtered signal
The summation of related coefficient between the summation of Relative RMSE and every two neighboring useful IMF;Finally with two summations it
The amplitude of addition noise corresponding to the noisy signal of (representing the parameter by the use of Rel-Corr RMSE) minimum is as EEMD side
The amplitude of addition noise used in method.The mathematical expression of Rel-Corr RMSE is specific as follows:
Wherein,It is signal of the original signal after mean filter, ci,rel(k) it is i-th of useful IMF, corr ()
Expression seeks related coefficient to two signals in bracket.
After the amplitude for determining addition noise, you can original signal is decomposed using EEMD, and determines noise IMF, make an uproar
Acoustic correlation IMF and useful IMF.
Step 2, to noise IMF, IMF related to noise carries out threshold filter denoising, and using mask signal method to useful
The signal that IMF is formed carries out the elimination of modal overlap;Wherein, mask signal frequency used in mask signal method is according to following
Method determines:The signal formed first to useful IMF carries out first derivation, obtains first derivative signal;Then Xi Er is utilized
Bert transform method obtains the instantaneous frequency of the first derivative signal, then asks for the mean value of instantaneous frequency;Finally to described
The mean value of instantaneous frequency is multiplied by the weights that a value range is (1,2), and gained product is mask signal frequency.
The signal formed in the present embodiment using improved MEMD algorithms to useful IMF carries out the elimination of modal overlap,
And to noise IMF, IMF related to noise carries out threshold filter respectively using the filtering threshold optimized.It illustrates separately below.
(1) optimization of mask signal frequency
In practical applications, reference frequency is often below the required requirements of masking signal to existing MEMD algorithms, because
And the ideal effect for eliminating modal overlap cannot be reached, therefore, for different signals, the result that the parameter of recommendation calculates is not
It is optimal.For this purpose, the present invention improves existing MEMD algorithms, mixed by optimizing mask signal frequency to promote elimination mode
Folded effect.Specifically, as shown in Fig. 2, mask signal frequency determines in accordance with the following methods used in mask signal method:
The signal formed first to useful IMF carries out first derivation, obtains first derivative signal;Assuming that useful IMF institutes
The signal of composition is x1=x2+x3, x in ingredient2、x3Corresponding frequency is f1、f2, and assume f2>f1, then have:
x1'=x2′+x3′ (2)
First derivative is sought in apostrophe representative above.
Then the instantaneous frequency of the first derivative signal is obtained using Hilbert transform methods, and asks for instantaneous frequency
Mean value;Its mathematical expression is as follows:
Wherein, H [(x1')] it is exactly to ask for Hilbert transformation,Represent the equal of the instantaneous frequency that Hilbert transformation acquires
Value, f1' (i) represents instantaneous frequency, f2Represent the frequency values of the highest ingredient of frequency in aliasing signal.
The weights that a value range is (1,2) are finally multiplied by the mean value f of instantaneous frequency, gained product is mask letter
Number frequency fmask, i.e.,:
(2) the filtering threshold optimization of noise correlation IMF
Problem is determined to the threshold value of noise correlation IMF, it is contemplated that asymmetry of the signal in noise IMF is needed to original
The positive and negative separate computations of signal framing value.First noise correlation IMF (IMFk) and noise IMF (IMF k=2,3...1) carry out
Decorrelation operates, removal partial noise interference, as shown in formula (6).Then the mean value of signal amplitude is asked for, such as formula (7) institute
Show, the comprehensive measurement as noise and useful signal;Further according to energy threshold formulaDetermine corresponding energy
Filtering threshold is measured, processing finally is weighted to the mean value of signal amplitude and energy filter threshold value, obtains final noise threshold,
Such as formula (7), (8).
Wherein, IMFKFor noise correlation IMF, IMF1For noise IMF, | ci| for the absolute value of IMF each point amplitudes, TmeanFor
The mean value of signal amplitude, N are signal number.For IMF1Energy square,For the energy estimators of k-th of IMF, β, ρ, C
It is constant, β=0.719, ρ=2.01, C ∈ [0.1,1.4], TuniFor the filtering threshold finally used.Two during weighted sum
A weights a1、a2Usually all it is 0.5.
(3) the filtering threshold optimization of noise IMF
The prior art is typically directly to abandon noise IMF, it is possible to lose useful information therein.The present embodiment thus
It is middle to determine filtering threshold using the method for periodic authentication, fully to extract the useful information in noise IMF.Specifically, to making an uproar
Filtering threshold used in sound IMF progress threshold filter denoisings determines in accordance with the following methods:The amplitude mean value of noise IMF is calculated,
And time-sequencing is carried out to each amplitude for being more than amplitude mean value in noise IMF, obtain the time series of the amplitude;For institute
The obtained time series of each amplitude carries out periodic authentication respectively, and therefrom selects periodically higher part-time sequence;
Finally using the minimum value of amplitude corresponding to selected part-time sequence as the filtering threshold.
Periodic authentication method in the present embodiment is specific as follows:To the time series of the amplitude, first adjacent amplitude
Time difference is in chronological sequence ranked sequentially, and obtains first time difference sequence, and calculates the mean value k of these time differences;Then to
Adjacent time difference makes the difference and is in chronological sequence ranked sequentially difference again in one time difference sequence, obtains the second time difference sequence;
To each data in the second time difference sequence, the quotient obtained by itself divided by mean value k and predetermined threshold value α are compared, less than α,
1 is assigned a value of, otherwise, is assigned a value of 0, so as to obtain a binary sequence;The Lempel-Ziv for calculating the binary sequence is complicated
Degree, complexity is lower, then the time series of the amplitude corresponding to the binary sequence is periodically higher.
Its algorithm is expressed as follows:Assuming that the time series of similar magnitude through making the difference obtained second time difference sequence twice
For { c (n1),c(n2),…,c(nm), binary conversion treatment is carried out to the second time difference sequence by following formula, obtains one new two
System sequence { Bc (n1),Bc(n2),…,Bc(nm)}:
Then binary sequence { Bc (n are calculated1),Bc(n2),…,Bc(nm) Lempel-Ziv complexities CLZN, and root
Whether the time series for judging the amplitude according to Lempel-Ziv complexities has periodically:
Wherein, c (ni) for i-th of element in the second time difference sequence;Bc(ni) for i-th yuan in binary sequence
Element;CLZNFor Lempel-Ziv complexities;β is the threshold value of Lempel-Ziv complexities, is to judge whether signal has periodically
Final step generally takes 0.5 the following value (the smaller the better).
Step 3, the noise IMF to after threshold filter denoising, noise correlation IMF eliminate the useful IMF institutes structure of modal overlap
Into signal and discrepance carry out signal reconstruction.
Claims (6)
1. a kind of denoising method based on mixing empirical mode decomposition, which is characterized in that include the following steps:
Step 1 decomposes original signal using gathering empirical mode decomposition EEMD methods, obtains a series of natural mode of vibration point
IMF and discrepance are measured, and noise IMF, noise correlation IMF, useful IMF are determined from all IMF;Made in the EEMD methods
The amplitude of addition noise predefines by the following method:Original signal is filtered first, obtains filtered signal,
And add the noise of different amplitudes respectively to original signal, obtain one group of noisy signal;Then each noisy signal is distinguished
Carry out following handle:Empirical mode decomposition is carried out, and determined from the obtained IMF of empirical mode decomposition to the noisy signal
Use IMF;Calculate the summation of the related root-mean-square error between each useful IMF and filtered signal and per two neighboring useful
The summation of related coefficient between IMF;Finally with the amplitude of the addition noise corresponding to the noisy signal of the sum of two summations minimum
Amplitude as the addition noise used in EEMD methods;
Step 2, to noise IMF, IMF related to noise carries out threshold filter denoising, and using mask signal method to useful IMF institutes
The signal of composition carries out the elimination of modal overlap;Wherein, mask signal frequency used in mask signal method is in accordance with the following methods
It determines:The signal formed first to useful IMF carries out first derivation, obtains first derivative signal;Then Hilbert is utilized
Transform method obtains the instantaneous frequency of the first derivative signal, then asks for the mean value of instantaneous frequency;Finally to described instantaneous
The mean value of frequency is multiplied by the weights that a value range is (1,2), and gained product is mask signal frequency;
Step 3, the noise IMF to after threshold filter denoising, noise correlation IMF eliminate what the useful IMF of modal overlap was formed
Signal and discrepance carry out signal reconstruction.
2. denoising method as described in claim 1, which is characterized in that the width of the noise of different amplitudes added to original signal
The original signal standard deviation of ranging from 0.001-0.4 times of value.
3. denoising method as described in claim 1, which is characterized in that threshold filter denoising is carried out to noise correlation IMF and is used
Filtering threshold determine in accordance with the following methods:Each noise correlation IMF is subjected to decorrelative transformation with noise IMF respectively first;Root
It is two parts according to positive and negative IMF points of the noise correlation by after each decorrelative transformation of amplitude, obtains its width respectively to each section
It is worth mean value and energy filter threshold value, and by the two weighted sum to get the filtering threshold of the noise correlation IMF after decorrelative transformation
Value.
4. denoising method as claimed in claim 3, which is characterized in that summation is weighted to amplitude mean value and energy filter threshold value
When, the weights of the two are 0.5.
5. denoising method as described in claim 1, which is characterized in that filter used in threshold filter denoising is carried out to noise IMF
Wave threshold value determines in accordance with the following methods:The amplitude mean value of noise IMF is calculated, and to being more than each of amplitude mean value in noise IMF
A amplitude carries out time-sequencing, obtains the time series of the amplitude;The time series of obtained each amplitude is carried out respectively
Periodic authentication, and therefrom select periodically higher part-time sequence;It is finally right with selected part-time sequence institute
The minimum value of amplitude is answered as the filtering threshold.
6. denoising method as claimed in claim 5, which is characterized in that for the time series of any one amplitude, the period
Property verification it is specific as follows:To the time series of the amplitude, first the time difference of adjacent amplitude is in chronological sequence ranked sequentially, is obtained
First time difference sequence, and calculate the mean value k of these time differences;Then adjacent time difference in first time difference sequence is made the difference again
And be in chronological sequence ranked sequentially difference, obtain the second time difference sequence;It, will to each data in the second time difference sequence
Quotient obtained by itself divided by mean value k is compared with predetermined threshold value α, less than α, is assigned a value of 1, otherwise, is assigned a value of 0, so as to obtain one
A binary sequence;The Lempel-Ziv complexities of the binary sequence are calculated, complexity is lower, then binary sequence institute is right
The periodicity of the time series for the amplitude answered is higher.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510844356.8A CN105488341B (en) | 2015-11-27 | 2015-11-27 | A kind of denoising method based on mixing empirical mode decomposition |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510844356.8A CN105488341B (en) | 2015-11-27 | 2015-11-27 | A kind of denoising method based on mixing empirical mode decomposition |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105488341A CN105488341A (en) | 2016-04-13 |
CN105488341B true CN105488341B (en) | 2018-06-29 |
Family
ID=55675315
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510844356.8A Active CN105488341B (en) | 2015-11-27 | 2015-11-27 | A kind of denoising method based on mixing empirical mode decomposition |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105488341B (en) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106154344A (en) * | 2016-08-01 | 2016-11-23 | 湖南文理学院 | A kind of Magnetotelluric signal denoising method based on combined filter |
CN106706122B (en) * | 2017-01-24 | 2019-07-12 | 东南大学 | Method for denoising bump-scrape acoustic emission signal based on related coefficient and EMD filtering characteristic |
CN107464226B (en) * | 2017-07-31 | 2019-10-15 | 东南大学 | A kind of image de-noising method based on improvement two-dimensional empirical mode decomposition algorithm |
CN109934050A (en) * | 2017-12-15 | 2019-06-25 | 西安艾索信息技术有限公司 | A kind of improved ocean clutter cancellation algorithm decomposed based on EMD |
CN108182415B (en) * | 2017-12-29 | 2021-05-04 | 哈尔滨工业大学 | Method for obtaining time-frequency distribution based on self-filtering frequency conversion empirical mode decomposition |
CN109520493A (en) * | 2018-10-29 | 2019-03-26 | 东南大学 | A kind of earth magnetism filtering method based on multiple dimensioned empirical mode decomposition |
CN109558857B (en) * | 2018-12-10 | 2022-05-31 | 西安邮电大学 | Chaotic signal noise reduction method |
CN109813417A (en) * | 2019-01-18 | 2019-05-28 | 国网江苏省电力有限公司检修分公司 | A kind of shunt reactor method for diagnosing faults based on improvement EMD |
CN111368778B (en) * | 2020-03-13 | 2022-05-27 | 哈尔滨理工大学 | Weak signal noise stripping method based on intelligent optimization algorithm |
CN112115764A (en) * | 2020-07-24 | 2020-12-22 | 南京航空航天大学 | Mask signal frequency selection method based on energy index |
CN115932389A (en) * | 2022-11-25 | 2023-04-07 | 国网浙江省电力有限公司嘉善县供电公司 | HHT (Hilbert-Huang transform) -based voltage quality disturbance detection method |
CN116509415B (en) * | 2023-04-21 | 2024-01-26 | 山东省人工智能研究院 | Noise reduction method based on unitized morphology features of signal components |
CN116776168B (en) * | 2023-08-22 | 2023-11-21 | 惠州帝恩科技有限公司 | Intelligent analysis method and system for production data of reagent tubes |
CN117093827B (en) * | 2023-10-16 | 2024-01-30 | 欣灵电气股份有限公司 | Intelligent fire control water supply data processing system based on Internet of things |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103984866A (en) * | 2014-05-20 | 2014-08-13 | 浙江师范大学 | Signal denoising method based on local mean value decomposition |
CN104007315A (en) * | 2014-04-29 | 2014-08-27 | 昆明理工大学 | Improved empirical mode decomposition processing method |
-
2015
- 2015-11-27 CN CN201510844356.8A patent/CN105488341B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104007315A (en) * | 2014-04-29 | 2014-08-27 | 昆明理工大学 | Improved empirical mode decomposition processing method |
CN103984866A (en) * | 2014-05-20 | 2014-08-13 | 浙江师范大学 | Signal denoising method based on local mean value decomposition |
Non-Patent Citations (5)
Title |
---|
EMD阈值滤波在光纤陀螺漂移信号去噪中的应用;崔冰波等;《光学学报》;20150210;第35卷(第2期);第0207001-1-0207001-6页 * |
Ensemble Empirical Mode Decomposition:A Noise Assisted Data Analysis Method;Zhaohua Wu et al.;《Advances in Adaptive Data Analysis》;20091231;第1-41页 * |
On empirical mode decomposition and its algorithms;Gabriel Rilling et al.;《in Proc. IEEE-EURASIP Workshop on Nonlinear Sig》;20030630;第1-5页 * |
基于现代时频分析的环境激励模态参数识别方法研究;何启源;《中国博士学位论文全文数据库 信息科技辑》;20091215;第I138-50页 * |
运用改进掩膜信号法的经验模态分解;赵玲等;《振动、测试与诊断》;20100815;第30卷(第8期);第434-439页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105488341A (en) | 2016-04-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105488341B (en) | A kind of denoising method based on mixing empirical mode decomposition | |
CN102054480B (en) | Method for separating monaural overlapping speeches based on fractional Fourier transform (FrFT) | |
CN108491355A (en) | A kind of ultrasonic signal noise-reduction method based on CEEMD and wavelet packet | |
CN106798554A (en) | A kind of denoising method of noisy IMF components and electrocardiosignal | |
CN103529436A (en) | Method for carrying out separation and time-frequency analysis on respiration and heartbeat signals in non-contact life detection on basis of HHT (Hilbert Huang Transform) | |
CN110151175A (en) | Surface electromyogram signal noise-eliminating method based on CEEMD and improvement wavelet threshold | |
CN103778921A (en) | Method for eliminating nonuniform noise in speech collected by radar | |
CN101991418A (en) | Method for improving respiratory rate detection accuracy | |
CN114469124B (en) | Method for identifying abnormal electrocardiosignals in movement process | |
CN103258543B (en) | Method for expanding artificial voice bandwidth | |
CN105845149A (en) | Predominant pitch acquisition method in acoustical signal and system thereof | |
CN104729677B (en) | A kind of time-domain digital weighted method of nonstationary noise signal | |
JP2016006536A (en) | Complex acoustic resonance speech analysis system | |
CN117609702A (en) | Pipeline leakage acoustic emission signal denoising method, system, equipment and medium | |
CN105446940B (en) | A kind of amplitude rectification method based on synchronous compression transformation reconstruct | |
CN108761202A (en) | The harmonic detecting method that pole symmetric mode decomposition and Hilbert transform are combined | |
CN108806721A (en) | signal processor | |
JP2016500847A (en) | Digital processor based complex acoustic resonance digital speech analysis system | |
CN103971012B (en) | Lamb wave signal denoising method based on John Saris model and fractional differentiation | |
CN101449321A (en) | Out-of-band signal generator and frequency band expander | |
CN107577986B (en) | Respiration and heartbeat component extraction method, electronic equipment and storage medium | |
Keizer et al. | The acceleration of sea-level rise along the coast of the Netherlands started in the 1960s | |
RU2395158C1 (en) | Digital signal filtration method | |
CN111551785A (en) | Frequency and harmonic detection method based on unscented Kalman filter | |
Hou et al. | Convolutional recurrent MetriCGAN with spectral dimension compression for full-band speech enhancement |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |