CN105811921B - A kind of method and wave filter for suppressing industrial frequency harmonic interference - Google Patents

A kind of method and wave filter for suppressing industrial frequency harmonic interference Download PDF

Info

Publication number
CN105811921B
CN105811921B CN201510547011.6A CN201510547011A CN105811921B CN 105811921 B CN105811921 B CN 105811921B CN 201510547011 A CN201510547011 A CN 201510547011A CN 105811921 B CN105811921 B CN 105811921B
Authority
CN
China
Prior art keywords
signal
frequency
harmonic
reference signal
domain information
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510547011.6A
Other languages
Chinese (zh)
Other versions
CN105811921A (en
Inventor
朱锋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Vivo Mobile Communication Co Ltd
Original Assignee
Vivo Mobile Communication Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Vivo Mobile Communication Co Ltd filed Critical Vivo Mobile Communication Co Ltd
Priority to CN201510547011.6A priority Critical patent/CN105811921B/en
Publication of CN105811921A publication Critical patent/CN105811921A/en
Application granted granted Critical
Publication of CN105811921B publication Critical patent/CN105811921B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The invention provides a kind of method and wave filter for suppressing industrial frequency harmonic interference, it is related to field of signal processing, solves the problems, such as that the power frequency inhibition of existing wave filter is limited, the method includes：The time-domain information of time-domain information and default multiple harmonic reference signal according to acquired original signal, obtains the convolution nuclear information that each harmonic reference signal distinguishes corresponding industrial frequency harmonic interference signal；Frequency domain information and each harmonic reference signal according to each harmonic reference signal distinguish the convolution nuclear information of corresponding industrial frequency harmonic interference signal, obtain total frequency domain information of all industrial frequency harmonic interference signals；Total frequency domain information of frequency domain information and all industrial frequency harmonic interference signals according to acquired original signal, obtain the first estimation frequency domain information that acquired original signal removes the first signal after all industrial frequency harmonic interference signals, and estimate frequency domain information according to first, obtain the first estimation time-domain information of the first signal.The solution of the present invention effectively inhibits the interference of industrial frequency harmonic signal.

Description

A kind of method and wave filter for suppressing industrial frequency harmonic interference
Technical field
The present invention relates to field of signal processing, more particularly to a kind of method and wave filter for suppressing industrial frequency harmonic interference.
Background technology
In field of signal processing, such as ECG signal, image pick-up signal and electromagnetic probe signal, Hz noise are equal Generally existing.The presence of Hz noise would generally limit the signal bandwidth and Detection results of engineer applied, be one and more difficult overcome Problem.Especially for the industrial frequency harmonic of 50Hz, it is difficult to design good narrow band filter and go to overcome.Although some scholars propose The optimal notch filter of Line enhancement, or improve, such as IIR lattice adaptive resistance-trap filter, essence is equivalent to arrowband Wave filter, but the power frequency inhibition of these wave filters is limited.
The content of the invention
The technical problem to be solved in the present invention is to provide a kind of method and wave filter for suppressing industrial frequency harmonic interference, solves existing There is the limited problem of the power frequency inhibition of wave filter.
In order to solve the above technical problems, embodiments of the invention provide a kind of method for suppressing industrial frequency harmonic interference, including：
The time-domain information of time-domain information and default multiple harmonic reference signal according to acquired original signal, obtains each institute State the convolution nuclear information that harmonic reference signal distinguishes corresponding industrial frequency harmonic interference signal；
Frequency domain information and each described harmonic reference signal according to each harmonic reference signal distinguish corresponding work The convolution nuclear information of frequency harmonic interference signals, obtains total frequency domain information of all industrial frequency harmonic interference signals；
Total frequency domain information of frequency domain information and all industrial frequency harmonic interference signals according to the acquired original signal, Obtain the first estimation frequency domain letter that the acquired original signal removes the first signal after all industrial frequency harmonic interference signals Breath, and estimate frequency domain information according to described first, obtain the first estimation time-domain information of first signal.
In order to solve the above technical problems, embodiments of the invention also provide a kind of wave filter, including：
First acquisition module, for the time-domain information according to acquired original signal and default multiple harmonic reference signals when Domain information, obtains the convolution nuclear information that each described harmonic reference signal distinguishes corresponding industrial frequency harmonic interference signal；
Second acquisition module, for the frequency domain information according to each harmonic reference signal and each described harmonic reference Signal distinguishes the convolution nuclear information of corresponding industrial frequency harmonic interference signal, obtains total frequency of all industrial frequency harmonic interference signals Domain information；
3rd acquisition module, for the frequency domain information according to the acquired original signal and all industrial frequency harmonic interference Total frequency domain information of signal, obtains the first signal after the acquired original signal removes all industrial frequency harmonic interference signals The first estimation frequency domain information, and estimate frequency domain information according to described first, obtain the first estimation time domain of first signal Information.
Above-mentioned technical proposal of the invention has the beneficial effect that：
The embodiment of the present invention suppress industrial frequency harmonic interference method, first the time-domain information according to acquired original signal and The time-domain information of default multiple harmonic reference signal, obtains each harmonic reference signal and distinguishes corresponding industrial frequency harmonic interference signal Convolution nuclear information；Then the convolution nuclear information of corresponding industrial frequency harmonic interference signal is distinguished according to each harmonic reference signal, Obtain total frequency domain information of all industrial frequency harmonic interference signals；Frequency domain information and all power frequencies finally according to acquired original signal Total frequency domain information of harmonic interference signals, obtains the first estimation after acquired original signal removes all industrial frequency harmonic interference signals Frequency domain information, and according to first estimate frequency domain information, obtain the first signal first estimation time-domain information so that effective percentage except The interference of industrial frequency harmonic signal in acquired original signal.The method has greatly value in Hz noise suppression, reduces The cumbersome degree of trapper is set, it is simple and easy to apply, and can effectively suppress industrial frequency harmonic signal.For traditional trapper, should Method has the accuracy that frequency adaptivity and time delay are estimated, the improvement degree of signal to noise ratio is good, and signal distortion is small, is suitable for Engineering it is commonly used, and with good stability.
Brief description of the drawings
Fig. 1 is the method flow diagram that the present invention suppresses industrial frequency harmonic interference；
Fig. 2 is the flow chart that the method one that present invention suppression industrial frequency harmonic is disturbed implements embodiment；
Fig. 3 is signal time-domain information schematic diagram not by industrial frequency harmonic interference；
Fig. 4 is signal amplitude spectrum schematic diagram not by industrial frequency harmonic interference；
Fig. 5 is signal phase spectrum schematic diagram not by industrial frequency harmonic interference；
Fig. 6 is the signal time-domain information schematic diagram disturbed by industrial frequency harmonic；
Fig. 7 is that the signal amplitude disturbed by industrial frequency harmonic composes schematic diagram；
Fig. 8 is that the signal phase disturbed by industrial frequency harmonic composes schematic diagram；
Fig. 9 is to suppress the signal time-domain information schematic diagram after industrial frequency harmonic is disturbed using the method for the present invention；
Figure 10 is the method singular point rejecting schematic diagram that the present invention suppresses industrial frequency harmonic interference；
Figure 11 is the structural representation of wave filter of the present invention.
Specific embodiment
To make the technical problem to be solved in the present invention, technical scheme and advantage clearer, below in conjunction with accompanying drawing and tool Body embodiment is described in detail.
The method for suppressing industrial frequency harmonic interference of the embodiment of the present invention, using power Spectral Estimation general principle, has obtained work The convolution nuclear information of frequency harmonic interference signals, and industrial frequency harmonic interference signal is eliminated in frequency domain, effectively inhibit work The interference of frequency harmonic wave, without the computing resource of complexity high；For small-sized calculating platform strong applicability, and signal to noise ratio improves degree Good, signal distortion is small, and real-time is good, with good stability.
First embodiment：
As shown in figure 1, the method for suppressing industrial frequency harmonic interference of the embodiment of the present invention, including：
Step 101, the time-domain information of time-domain information and default multiple harmonic reference signal according to acquired original signal, obtains Take the convolution nuclear information that each harmonic reference signal distinguishes corresponding industrial frequency harmonic interference signal.
Here, harmonic reference signal corresponds to the industrial frequency harmonic interference signal in Engineering Signal, the number of harmonic reference signal Amount is the quantity of industrial frequency harmonic interference signal in Engineering Signal frequency band, and the frequency of harmonic reference signal is equal to corresponding power frequency The frequency of harmonic interference signals.
Step 102, frequency domain information and each harmonic reference signal according to each harmonic reference signal distinguish corresponding work The convolution nuclear information of frequency harmonic interference signals, obtains total frequency domain information of all industrial frequency harmonic interference signals.
Here, the frequency domain of industrial frequency harmonic interference signal can accurately be obtained by the convolution nuclear information of industrial frequency harmonic interference signal Information, and then remove interference of the industrial frequency harmonic interference signal to acquired original signal in frequency domain.
Step 103, total frequency domain information of frequency domain information and all industrial frequency harmonic interference signals according to acquired original signal, First estimation frequency domain information of the first signal after the acquisition acquired original signal all industrial frequency harmonic interference signals of removing, and according to First estimates frequency domain information, obtains the first estimation time-domain information of the first signal.
Here, total frequency domain information of all industrial frequency harmonic interference signals is subtracted by the frequency domain information of acquired original signal, Interference of the power frequency interference signals to acquired original signal can be removed in frequency domain, and it is humorous to remove all power frequencies by acquired original signal First estimation frequency domain information of the first signal after interference wave signal, can further obtain the first estimation time domain letter of the first signal Breath, and then obtain effectively filtering out the signal after industrial frequency harmonic interference.
In the present embodiment, the time domain of time-domain information and default multiple harmonic reference signal first according to acquired original signal Information, obtains the convolution nuclear information that each harmonic reference signal distinguishes corresponding industrial frequency harmonic interference signal；Then according to each Harmonic reference signal distinguishes the convolution nuclear information of corresponding industrial frequency harmonic interference signal, obtains all industrial frequency harmonic interference signals Total frequency domain information；Total frequency domain information of frequency domain information and all industrial frequency harmonic interference signals finally according to acquired original signal, The first estimation frequency domain information after acquired original signal removes all industrial frequency harmonic interference signals is obtained, and estimates frequency according to first Domain information, obtains the first estimation time-domain information of the first signal, so that effective percentage is except industrial frequency harmonic letter in acquired original signal Number interference.The method, with greatly value, reduces the cumbersome degree for setting trapper in Hz noise suppression, simple easy OK, and industrial frequency harmonic signal can be effectively suppressed.For traditional trapper, the method has frequency adaptivity and time delay The accuracy of estimation, the improvement degree of signal to noise ratio is good, and signal distortion is small, is suitable for the commonly used of engineering, and with good Stability.
Second embodiment：
The method for suppressing industrial frequency harmonic interference of the embodiment of the present invention, including：
Step 1011, the time-domain information of time-domain information and each harmonic reference signal according to acquired original signal is obtained Each harmonic reference signal respectively with the crosspower spectrum information of acquired original signal.
Step 1012, according to the time-domain information of each harmonic reference signal, obtain respectively each harmonic reference signal from Power spectral information.
Step 1013, auto-power spectrum information and corresponding crosspower spectrum information according to each harmonic reference signal are obtained Each harmonic reference signal distinguishes the convolution nuclear information of corresponding industrial frequency harmonic interference signal.
It is here, orthogonal between the auto-power spectrum information of each harmonic reference signal and corresponding crosspower spectrum information, It is also orthogonal between each harmonic reference signal.Can be mutual with corresponding according to the auto-power spectrum information of each harmonic reference signal Power spectral information, accurately obtains the convolution nuclear information that each harmonic reference signal distinguishes corresponding industrial frequency harmonic interference signal.
Step 102, frequency domain information and each harmonic reference signal according to each harmonic reference signal distinguish corresponding work The convolution nuclear information of frequency harmonic interference signals, obtains total frequency domain information of all industrial frequency harmonic interference signals.
Step 103, total frequency domain information of frequency domain information and all industrial frequency harmonic interference signals according to acquired original signal, First estimation frequency domain information of the first signal after the acquisition acquired original signal all industrial frequency harmonic interference signals of removing, and according to First estimates frequency domain information, obtains the first estimation time-domain information of the first signal.
In the present embodiment, 101 are embodied as step 1011-1013 the step of by first embodiment, every by obtaining respectively The crosspower spectrum information and the auto-power spectrum information of each harmonic reference signal of individual harmonic reference signal and acquired original signal, can The convolution nuclear information that each harmonic reference signal distinguishes corresponding industrial frequency harmonic interference signal is accurately obtained, is next to power frequency The calculating of total frequency domain information of harmonic interference signals is supported there is provided data, it is ensured that the accuracy and validity of calculating.
Generally the primary signal of collection is disturbed containing industrial frequency harmonic, and acquired original signal can typically be written as：
Wherein, x (t) is the time-domain function of acquired original signal, and s (t) is that acquired original signal removes all industrial frequency harmonics The time-domain function of the first signal after interference signal, v (t) is total time-domain function of all industrial frequency harmonic interference signals, hiT () is The convolution kernel of the corresponding industrial frequency harmonic interference signal of i-th harmonic reference signal, vi(t) be i-th harmonic reference signal when Domain function, i is 1 to the integer between p, and p is the quantity of industrial frequency harmonic interference signal in Engineering Signal frequency band, and t is time variable.
General harmonic reference signal can be written as：
vi(t)=Asin (2 π fit)；
Wherein.A freely sets constant, f for oneiBe i-th frequency of harmonic reference signal, and harmonic reference signal frequency Rate is equal to the frequency of corresponding industrial frequency harmonic interference signal.
Based on above-mentioned theory, it is preferred that in above-mentioned steps 1011, specifically each harmonic reference can be obtained by equation below Signal respectively with the cross-spectral density function of acquired original signal
In above-mentioned steps 1012, the autopower spectral density letter of each harmonic reference signal can be specifically obtained by equation below Number
Wherein, vi(t+ τ) is i-th time-domain function of harmonic reference signal delay time τ, and i is 1 to the integer between p, P is the quantity (quantity of industrial frequency harmonic interference signal i.e. in Engineering Signal frequency band) of harmonic reference signal, and t is time variable, ω It is frequency variable, j and τ is customized constant.
Due to orthogonal between each harmonic reference signal, s (t) and v (t) are also orthogonal, can obtain：
Thus, in above-mentioned steps 1013, specifically each harmonic reference signal can be obtained by equation below and distinguishes corresponding The convolution kernel function H of industrial frequency harmonic interference signali(ω)：
Wherein,It is the cross-spectral density function of i-th harmonic reference signal and acquired original signal, The Power spectral density of harmonic reference signal, H described in i-thi(ω) is hiThe Fourier transform function of (t).
Further, above-mentioned steps 102 can specifically include：
Step 1021, respectively to the time-domain function v of each harmonic reference signaliT () carries out Fourier transformation, obtain each The frequency-domain function V of the harmonic reference signali(ω)；
Step 1022, total frequency-domain function V (ω) of all industrial frequency harmonic interference signals is obtained by equation below：
Further, above-mentioned steps 103 can specifically include：
Step 1031, the time-domain function to acquired original signal carries out Fourier transformation, obtains the frequency of acquired original signal Domain function X (ω)；
Step 1032, the first estimation frequency-domain function of the first signal is obtained by equation below
Step 1033, frequency-domain function is estimated to firstInverse Fourier transform is carried out, obtain the first signal first is estimated Clocking domain function
So far, by obtain each harmonic reference signal respectively with the cross-spectral density function of acquired original signalThe Power spectral density of each harmonic reference signalAnd then obtain each harmonic reference signal difference The convolution kernel function H of corresponding industrial frequency harmonic interference signali(ω)；And the frequency-domain function V for passing through each harmonic reference signali The convolution kernel function H of (ω) and each industrial frequency harmonic interference signali(ω), obtains total frequency domain of all industrial frequency harmonic interference signals Function V (ω)；The frequency-domain function X (ω) of recycling acquired original signal subtracts total frequency domain letter of all industrial frequency harmonic interference signals Number V (ω), obtains the first estimation frequency-domain function of the first signalFourier finally is carried out to the first estimation frequency-domain function Inverse transformation, has accurately obtained the first estimation time-domain function of the first signal
3rd embodiment：
The method for suppressing industrial frequency harmonic interference of the embodiment of the present invention, including：
Step 101, the time-domain information of time-domain information and default multiple harmonic reference signal according to acquired original signal, obtains Take the convolution nuclear information that each harmonic reference signal distinguishes corresponding industrial frequency harmonic interference signal.
Step 102, frequency domain information and each harmonic reference signal according to each harmonic reference signal distinguish corresponding work The convolution nuclear information of frequency harmonic interference signals, obtains total frequency domain information of all industrial frequency harmonic interference signals.
Step 103, total frequency domain information of frequency domain information and all industrial frequency harmonic interference signals according to acquired original signal, Obtain the first estimation frequency domain information that acquired original signal removes the first signal after all industrial frequency harmonic interference signals.
Step 104, frequency domain information is estimated according to first, and the first estimation amplitude spectrum and first for obtaining the first signal estimates phase Position spectrum.
Here, frequency-domain function can be estimated by first when implementingSeparately write as the form of amplitude and phase：Wherein, ω is frequency variable, and j is customized constant.
Step 105, according to the frequency domain information of acquired original signal, obtains the amplitude spectrum of acquired original signal, and to original The amplitude spectrum for gathering signal carries out singular point rejecting treatment, obtains the second estimation amplitude spectrum.
Here, the amplitude spectrum of the frequency-domain function X (ω) of acquired original signal | X (ω) | can be carried out unusual when implementing Point rejecting treatment, the second estimation amplitude spectrum X'(ω after being rejected).
Step 106, replaces with the second estimation amplitude spectrum, and keep by the first estimation amplitude spectrum of the first estimation frequency domain information First estimates that phase spectrum is constant, obtains the second estimation frequency domain information of the first signal.
Here, frequency-domain function can be estimated by first when implementingAmplitude spectrum replace with the second estimation Amplitude spectrum X'(ω), obtain the second estimation frequency-domain function
Step 107, frequency domain information is estimated according to second, obtains the second estimation time-domain information of the first signal.
Here, frequency-domain function can be estimated to second when implementingFourier transformation is carried out, the first signal is obtained Second estimates time-domain function
In the present embodiment, after getting the first estimation frequency domain information of the first signal, to the amplitude spectrum of acquired original signal Singular point rejecting treatment is carried out, and the amplitude spectrum of the first estimation frequency domain information is replaced with into the amplitude spectrum after rejecting, obtain second Estimate amplitude spectrum, further suppress interference of the industrial frequency harmonic signal to acquired original signal.
Preferably, as shown in Figure 10, above-mentioned steps 105 are specifically included：
Step 1051, the amplitude spectrum to the acquired original signal by the way of robust iterative carries out singular point rejecting, And to rejecting point assignment again by the way of smooth compensating, obtain the second estimation amplitude spectrum.
Now, the amplitude spectrum to acquired original signal by way of robust iterative carries out singular point rejecting, and to rejecting Point assignment again by the way of smooth compensating, can effectively remove what industrial frequency harmonic interference signal was composed to acquired original signal amplitude Interference, further increases the effect for suppressing Hz noise, reduces the distortion factor of signal.
Fourth embodiment：
The method for suppressing industrial frequency harmonic interference of the embodiment of the present invention, including：
Step 1011, the time-domain information of time-domain information and each harmonic reference signal according to acquired original signal, passes through Equation below obtain each harmonic reference signal respectively with the cross-spectral density function of acquired original signal
Step 1012, according to the time-domain information of each harmonic reference signal, each harmonic reference is obtained by equation below The Power spectral density of signal
Wherein, x (t) is the time-domain function of the acquired original signal, viT () is harmonic reference signal described in i-th Time-domain function, vi(t+ τ) is the time-domain function of harmonic reference signal delay time τ described in i-th, and i is 1 to the integer between p, P is the quantity of the harmonic reference signal, and t is time variable, and ω is frequency variable, and j and τ is customized constant.
Step 1013, auto-power spectrum information and corresponding crosspower spectrum information according to each harmonic reference signal, passes through Equation below obtains the convolution kernel function H that each harmonic reference signal distinguishes corresponding industrial frequency harmonic interference signali(ω)：
Wherein,The cross-spectral density letter of harmonic reference signal and the acquired original signal described in i-th Number,The Power spectral density of harmonic reference signal described in i-th.
Step 1021, respectively to the time-domain function v of each harmonic reference signaliT () carries out Fourier transformation, obtain each The frequency-domain function V of harmonic reference signali(ω)。
Step 1022, total frequency-domain function V (ω) of all industrial frequency harmonic interference signals is obtained by equation below：
Step 1031, the time-domain function to acquired original signal carries out Fourier transformation, obtains the frequency of acquired original signal Domain function X (ω).
Step 1032, the first estimation frequency-domain function of the first signal is obtained by equation below
Step 1033, frequency-domain function is estimated to firstInverse Fourier transform is carried out, obtain the first signal first is estimated Clocking domain function
Step 104, frequency domain information is estimated according to first, and the first estimation amplitude spectrum and first for obtaining the first signal estimates phase Position spectrum.
Here, frequency-domain function can be estimated by first when implementingSeparately write as the form of amplitude and phase：Wherein, ω is frequency variable, and j is customized constant.
Step 1051, the amplitude spectrum to the acquired original signal by the way of robust iterative carries out singular point rejecting, And to rejecting point assignment again by the way of smooth compensating, obtain the second estimation amplitude spectrum.
Here, the amplitude spectrum of the frequency-domain function X (ω) of acquired original signal | X (ω) | can be carried out unusual when implementing Point rejecting treatment, the second estimation amplitude spectrum X'(ω after being rejected).
Step 106, replaces with the second estimation amplitude spectrum, and keep by the first estimation amplitude spectrum of the first estimation frequency domain information First estimates that phase spectrum is constant, obtains the second estimation frequency domain information of the first signal.
Here, frequency-domain function can be estimated by first when implementingAmplitude spectrum replace with the second estimation Amplitude spectrum X'(ω), obtain the second estimation frequency-domain function
Step 107, frequency domain information is estimated according to second, obtains the second estimation time-domain information of the first signal.
Here, frequency-domain function can be estimated to second when implementingFourier transformation is carried out, the first signal is obtained Second estimates time-domain function
The primary signal of general collection is disturbed containing industrial frequency harmonic, it is assumed that when signal is not disturbed by industrial frequency harmonic, signal Time-domain information, amplitude spectrum and phase spectrum respectively as in Figure 3-5, it is assumed that when signal is disturbed by industrial frequency harmonic, the time domain of signal Information, amplitude spectrum and phase spectrum are distinguished as shown in figs 6-8, it can be seen that being compared by the signal that industrial frequency harmonic is disturbed Relatively large deviation is occurred in that in the signal not disturbed by industrial frequency harmonic, accordingly, it would be desirable to a kind of mode effectively suppresses industrial frequency harmonic interference Influence to signal.
Using the method for the embodiment of the present invention, as shown in Fig. 2 acquired original signal x (t) is obtained first, and according to original Collection signal x (t) harmony reference signal wave 1-p, obtain each harmonic reference signal respectively with the cross-power of acquired original signal Spectral density functionAnd the Power spectral density of each harmonic reference signalThen each power frequency is calculated humorous The convolution kernel H of reference signal wavei(ω), and respectively to the time-domain function v of each harmonic reference signaliT () carries out Fourier Fourier is converted, and obtains the frequency-domain function V of each harmonic reference signali(ω)；According to the volume of each industrial frequency harmonic reference signal Product core HiThe frequency-domain function V of (ω) and each harmonic reference signali(ω), obtains total frequency domain of all industrial frequency harmonic interference signals Function V (ω)：Time-domain function x (t) to acquired original signal carries out Fourier conversion, obtains the frequency domain of acquired original signal Function X (ω), and frequency-domain function X (ω) using acquired original signal subtracts total frequency domain letter of all industrial frequency harmonic interference signals Number V (ω), obtains the first estimation frequency-domain functionAmplitude spectrum is carried out to the frequency-domain function X (ω) of acquired original signal sane Estimate, obtain the second estimation amplitude spectrum X'(ω), obtain first and estimate that the phase spectrum of frequency-domain function and the second estimation amplitude spectrum enter Row combination, obtains the second estimation frequency-domain functionFourier inverse transformations are carried out to the second estimation frequency-domain function, second is obtained Estimate time-domain functionFinally, the time-domain information for filtering the signal after Hz noise is as shown in Figure 9.
From fig. 9, it can be seen that the method for the embodiment of the present invention effectively inhibits industrial frequency harmonic signal in acquired original signal Interference, make signal distortion small, Hz noise suppression on have greatly value.
The method of the embodiment of the present invention is based on Fourier transformation, without the computing resource of complexity high, for small-sized calculating Platform strong applicability；The realization of the method for the embodiment of the present invention does not exist iterative relation, and major part is all based on Fourier transformation Carry out, real-time is good；And the method for the embodiment of the present invention can be obviously improved signal to noise ratio, can almost extract and primary signal phase Same signal, substantially reduces the distortion caused by Hz noise.
5th embodiment：
As shown in figure 11, embodiments of the invention additionally provide a kind of wave filter, including：
First acquisition module 1101, for the time-domain information according to acquired original signal and default multiple harmonic reference signal Time-domain information, obtain the convolution nuclear information that each described harmonic reference signal distinguishes corresponding industrial frequency harmonic interference signal.
Second acquisition module 1102, for the frequency domain information according to each harmonic reference signal and each described harmonic wave Reference signal distinguishes the convolution nuclear information of corresponding industrial frequency harmonic interference signal, obtains all industrial frequency harmonic interference signals Total frequency domain information.
3rd acquisition module 1103, for the frequency domain information according to the acquired original signal and all industrial frequency harmonics Total frequency domain information of interference signal, obtains first after all industrial frequency harmonic interference signals of the acquired original signal removing First estimation frequency domain information of signal, and estimate frequency domain information according to described first, obtain the first estimation of first signal Time-domain information.
The wave filter of the embodiment of the present invention, has greatly value in Hz noise suppression, reduces setting trapper Cumbersome degree, it is simple and easy to apply, and can effectively suppress industrial frequency harmonic signal.For traditional trapper, the method has frequency The accuracy that rate adaptivity and time delay are estimated, the improvement degree of signal to noise ratio is good, and signal distortion is small, is suitable for the universal of engineering Using, and with good stability.
Sixth embodiment：
Embodiments of the invention additionally provide a kind of wave filter, including：
First acquisition unit 11011, for the time-domain information according to the acquired original signal and each described harmonic wave ginseng The time-domain information of signal is examined, each described harmonic reference signal is obtained and is believed with the crosspower spectrum of the acquired original signal respectively Breath.
Second acquisition unit 11012, for the time-domain information according to each harmonic reference signal, obtains each respectively The auto-power spectrum information of the harmonic reference signal.
3rd acquiring unit 11013, for the auto-power spectrum information according to each harmonic reference signal and right The crosspower spectrum information answered, obtains the convolution that each described harmonic reference signal distinguishes corresponding industrial frequency harmonic interference signal Nuclear information.
Second acquisition module 1102, for the frequency domain information according to each harmonic reference signal and each described harmonic wave Reference signal distinguishes the convolution nuclear information of corresponding industrial frequency harmonic interference signal, obtains all industrial frequency harmonic interference signals Total frequency domain information.
3rd acquisition module 1103, for the frequency domain information according to the acquired original signal and all industrial frequency harmonics Total frequency domain information of interference signal, obtains first after all industrial frequency harmonic interference signals of the acquired original signal removing First estimation frequency domain information of signal, and estimate frequency domain information according to described first, obtain the first estimation of first signal Time-domain information.
Preferably, the first acquisition unit 11011 is specifically included：
First obtain subelement 110111, for by equation below obtain each described harmonic reference signal respectively with institute State the cross-spectral density function of acquired original signal
Wherein, x (t) is the time-domain function of the acquired original signal, vi(t+ τ) is harmonic reference signal described in i-th The time-domain function of delay time T, i is 1 to the integer between p, and p is the quantity of the harmonic reference signal, and t is time variable, ω is frequency variable, and j and τ is customized constant.
Preferably, the second acquisition unit 11012 is specifically included：
Second obtains subelement 110121, for by equation below obtain each harmonic reference signal from power Spectral density function
Wherein, viT () is the time-domain function of harmonic reference signal described in i-th, vi(t+ τ) is harmonic reference described in i-th The time-domain function of signal delay time τ, i is 1 to the integer between p, and p is the quantity of the harmonic reference signal, anaplasia when t is Amount, ω is frequency variable, and j and τ is customized constant.
Preferably, the 3rd acquiring unit 11013 is specifically included：
3rd obtains subelement 110131, is corresponded to respectively for obtaining each described harmonic reference signal by equation below Industrial frequency harmonic interference signal convolution kernel function Hi(ω)：
Wherein,The cross-spectral density letter of harmonic reference signal and the acquired original signal described in i-th Number,The Power spectral density of harmonic reference signal described in i-th, i is 1 to the integer between p, and p is described The quantity of harmonic reference signal, ω is frequency variable.
Preferably, second acquisition module 1102 is specifically included：
First converter unit 11021, Fourier is carried out for the time-domain function respectively to harmonic reference signal each described Conversion, obtains the frequency-domain function V of each harmonic reference signali(ω)；
4th acquiring unit 11022, the total frequency for obtaining all industrial frequency harmonic interference signals by equation below Domain function V (ω)：
Wherein, Hi(ω) is the convolution kernel function of the corresponding industrial frequency harmonic interference signal of harmonic reference signal described in i-th, Vi(ω) is the frequency-domain function of harmonic reference signal described in i-th, and i is 1 to the integer between p, and p is the harmonic reference signal Quantity, ω is frequency variable.
Preferably, the 3rd acquisition module 1103 is specifically included：
Second converter unit 11031, Fourier transformation is carried out for the time-domain function to the acquired original signal, is obtained The frequency-domain function X (ω) of the acquired original signal；
5th acquiring unit 11032, the first estimation frequency-domain function for obtaining first signal by equation below
Wherein, V (ω) is total frequency-domain function of all industrial frequency harmonic interference signals, and ω is frequency variable；
3rd converter unit 11033, for carrying out inverse Fourier transform to the described first estimation frequency-domain function, obtains described First estimation time-domain function of the first signal.
7th embodiment：
Embodiments of the invention additionally provide a kind of wave filter, including：
First acquisition module 1101, for the time-domain information according to acquired original signal and default multiple harmonic reference signal Time-domain information, obtain the convolution nuclear information that each described harmonic reference signal distinguishes corresponding industrial frequency harmonic interference signal.
Second acquisition module 1102, for the frequency domain information according to each harmonic reference signal and each described harmonic wave Reference signal distinguishes the convolution nuclear information of corresponding industrial frequency harmonic interference signal, obtains all industrial frequency harmonic interference signals Total frequency domain information.
3rd acquisition module 1103, for the frequency domain information according to the acquired original signal and all industrial frequency harmonics Total frequency domain information of interference signal, obtains first after all industrial frequency harmonic interference signals of the acquired original signal removing First estimation frequency domain information of signal.
4th acquisition module 1104, for estimating frequency domain information according to described first, obtains the first of first signal Estimate that amplitude spectrum and first estimates phase spectrum.
Module 1105 is rejected, for the frequency domain information according to the acquired original signal, the acquired original signal is obtained Amplitude spectrum, and amplitude spectrum to the acquired original signal carries out singular point rejecting treatment, obtains the second estimation amplitude spectrum.
Replacement module 1106, for the first estimation amplitude spectrum of the described first estimation frequency domain information to be replaced with into described second Estimate amplitude spectrum, and keep described first to estimate that phase spectrum is constant, obtain the second estimation frequency domain information of first signal.
6th acquisition module 1107, for estimating frequency domain information according to described second, obtains the second of first signal Estimate time-domain information.
Wherein, the rejecting module 1105 is specifically included：
Culling unit 11051, is carried out very for the amplitude spectrum to the acquired original signal by the way of robust iterative Dissimilarity is rejected, and to rejecting point assignment again by the way of smooth compensating, obtains the second estimation amplitude spectrum.
The wave filter of the embodiment of the present invention, has greatly value in Hz noise suppression, reduces setting trapper Cumbersome degree, it is simple and easy to apply, and can effectively suppress industrial frequency harmonic signal.For traditional trapper, the method has frequency The accuracy that rate adaptivity and time delay are estimated, the improvement degree of signal to noise ratio is good, and signal distortion is small, is suitable for the universal of engineering Using, and with good stability.
It should be noted that in accompanying drawing or specification description, similar or identical part all uses identical figure number.It is attached The implementation for not illustrated in figure or being described, is form known to a person of ordinary skill in the art in art.In addition, though this Text can provide the demonstration of the parameter comprising particular value, it is to be understood that parameter is without being definitely equal to corresponding value, but be able to can connect The error margin received is similar to corresponding value in design constraint.The direction term mentioned in embodiment, for example " on ", D score, "front", "rear", "left", "right" etc., are only the directions of refer to the attached drawing.Therefore, the direction term for using is for illustrating not to use To limit the scope of the invention.
Additionally, the above-mentioned definition to each element and method is not limited in embodiment various concrete structures, the shape mentioned Shape or mode, those of ordinary skill in the art can simply be changed or be replaced to it.For example：The phase of harmonic reference signal Setting differs and is set to zero phase, and constant A can freely be set；The parameter setting of robust iterative can be configured according to concrete condition, Too large deviation will not generally be caused.
One of ordinary skill in the art will appreciate that all or part of flow in realizing above-described embodiment method, can be The hardware of correlation is instructed to complete by computer program, described program can be stored in a computer read/write memory medium In, the program is upon execution, it may include such as the flow of the embodiment of above-mentioned each method.Wherein, described storage medium can be magnetic Dish, CD, read-only memory (Read-Only Memory, ROM) or random access memory (Random Access Memory, RAM) etc..
The above is the preferred embodiment of the present invention, it is noted that for those skilled in the art For, on the premise of principle of the present invention is not departed from, some improvements and modifications can also be made, these improvements and modifications Should be regarded as protection scope of the present invention.

Claims (18)

1. it is a kind of to suppress the method that industrial frequency harmonic is disturbed, it is characterised in that including：
The time-domain information of time-domain information and default multiple harmonic reference signal according to acquired original signal, obtains each described humorous Reference signal wave distinguishes the convolution nuclear information of corresponding industrial frequency harmonic interference signal；
It is humorous that frequency domain information and each described harmonic reference signal according to each harmonic reference signal distinguish corresponding power frequency The convolution nuclear information of interference wave signal, obtains total frequency domain information of all industrial frequency harmonic interference signals；
Total frequency domain information of frequency domain information and all industrial frequency harmonic interference signals according to the acquired original signal, obtains The acquired original signal removes the first estimation frequency domain information of the first signal after all industrial frequency harmonic interference signals, and Estimate frequency domain information according to described first, obtain the first estimation time-domain information of first signal.
2. method according to claim 1, it is characterised in that the time-domain information according to acquired original signal and default The time-domain information of multiple harmonic reference signals, obtains each described harmonic reference signal and distinguishes corresponding industrial frequency harmonic interference signal Convolution nuclear information, specifically include：
The time-domain information of time-domain information and each harmonic reference signal according to the acquired original signal, obtains each institute State harmonic reference signal respectively with the crosspower spectrum information of the acquired original signal；
According to the time-domain information of each harmonic reference signal, the auto-power spectrum of each harmonic reference signal is obtained respectively Information；
The auto-power spectrum information and the corresponding crosspower spectrum information according to each harmonic reference signal, obtain every The individual harmonic reference signal distinguishes the convolution nuclear information of corresponding industrial frequency harmonic interference signal.
3. method according to claim 2, it is characterised in that the time-domain information according to the acquired original signal and The time-domain information of each harmonic reference signal, obtain each described harmonic reference signal respectively with the acquired original signal Crosspower spectrum information, specifically include：
By equation below obtain each described harmonic reference signal respectively with the cross-spectral density of the acquired original signal Function
S xv i ( ω ) = ∫ - ∞ + ∞ [ ∫ - ∞ + ∞ x ( t ) v i ( t + τ ) d t ] - j ω τ d t ;
Wherein, x (t) is the time-domain function of the acquired original signal, vi(t+ τ) is described in i-th during harmonic reference signal delay Between τ time-domain function, i is 1 to the integer between p, and p is the quantity of the harmonic reference signal, and t is time variable, and ω is frequency Rate variable, j and τ are customized constant.
4. method according to claim 2, it is characterised in that described to be believed according to the time domain of each harmonic reference signal Breath, obtains the auto-power spectrum information of each harmonic reference signal respectively, specifically includes：
The Power spectral density of each harmonic reference signal is obtained by equation below
S v i v i ( ω ) = ∫ - ∞ + ∞ [ ∫ - ∞ + ∞ v i ( t ) v i ( t + τ ) d t ] e - j ω τ d t ;
Wherein, viT () is the time-domain function of harmonic reference signal described in i-th, vi(t+ τ) is harmonic reference signal described in i-th The time-domain function of delay time T, i is 1 to the integer between p, and p is the quantity of the harmonic reference signal, and t is time variable, ω is frequency variable, and j and τ is customized constant.
5. method according to claim 2, it is characterised in that it is described according to each described harmonic reference signal from Power spectral information and the corresponding crosspower spectrum information, obtain each described harmonic reference signal and distinguish corresponding industrial frequency harmonic The convolution nuclear information of interference signal, specifically includes：
The convolution kernel letter that each described harmonic reference signal distinguishes corresponding industrial frequency harmonic interference signal is obtained by equation below Number Hi(ω)：
H i ( ω ) = S xv i ( ω ) S v i v i ( ω ) ;
Wherein,The cross-spectral density function of harmonic reference signal and the acquired original signal described in i-th,The Power spectral density of harmonic reference signal described in i-th, i is 1 to the integer between p, and p is the harmonic wave The quantity of reference signal, ω is frequency variable.
6. method according to claim 1, it is characterised in that described to be believed according to the frequency domain of each harmonic reference signal Breath and each described harmonic reference signal distinguish the convolution nuclear information of corresponding industrial frequency harmonic interference signal, obtain all works Total frequency domain information of frequency harmonic interference signals, specifically includes：
The time-domain function to harmonic reference signal each described carries out Fourier transformation respectively, obtains each described harmonic reference letter Number frequency-domain function Vi(ω)；
Total frequency-domain function V (ω) of all industrial frequency harmonic interference signals is obtained by equation below：
V ( ω ) = Σ i = 1 P H i ( ω ) V i ( ω ) ;
Wherein, Hi(ω) is the convolution kernel function of the corresponding industrial frequency harmonic interference signal of harmonic reference signal described in i-th, Vi (ω) is the frequency-domain function of harmonic reference signal described in i-th, and i is 1 to the integer between p, and p is the harmonic reference signal Quantity, ω is frequency variable.
7. method according to claim 1, it is characterised in that the frequency domain information according to the acquired original signal and Total frequency domain information of all industrial frequency harmonic interference signals, obtains the acquired original signal and removes all industrial frequency harmonics First estimation frequency domain information of the first signal after interference signal, and estimate frequency domain information according to described first, obtains described the First estimation time-domain information of one signal, specifically includes：
Time-domain function to the acquired original signal carries out Fourier transformation, obtains the frequency-domain function of the acquired original signal X(ω)；
The first estimation frequency-domain function of first signal is obtained by equation below
Wherein, V (ω) is total frequency-domain function of all industrial frequency harmonic interference signals, and ω is frequency variable；
Inverse Fourier transform is carried out to the described first estimation frequency-domain function, the first estimation time domain letter of first signal is obtained Number.
8. method according to claim 1, it is characterised in that the frequency domain information according to the acquired original signal and Total frequency domain information of all industrial frequency harmonic interference signals, obtains the acquired original signal and removes all industrial frequency harmonics After first estimation frequency domain information of the first signal after interference signal, methods described also includes：
Estimate frequency domain information according to described first, the first estimation amplitude spectrum and first for obtaining first signal estimates phase Spectrum；
According to the frequency domain information of the acquired original signal, the amplitude spectrum of the acquired original signal is obtained, and to described original The amplitude spectrum for gathering signal carries out singular point rejecting treatment, obtains the second estimation amplitude spectrum；
First estimation amplitude spectrum of the described first estimation frequency domain information is replaced with described second and estimates amplitude spectrum, and keep described First estimates that phase spectrum is constant, obtains the second estimation frequency domain information of first signal；
Estimate frequency domain information according to described second, obtain the second estimation time-domain information of first signal.
9. method according to claim 8, it is characterised in that the amplitude spectrum to the acquired original signal is carried out very Dissimilarity rejecting is processed, and obtains the second estimation amplitude spectrum, is specifically included：
Amplitude spectrum to the acquired original signal by the way of robust iterative carries out singular point rejecting, and is used to rejecting point The mode of smooth compensating assignment again, obtains the second estimation amplitude spectrum.
10. a kind of wave filter, it is characterised in that including：
First acquisition module, for the time-domain information according to acquired original signal and the time domain letter of default multiple harmonic reference signal Breath, obtains the convolution nuclear information that each described harmonic reference signal distinguishes corresponding industrial frequency harmonic interference signal；
Second acquisition module, for the frequency domain information according to each harmonic reference signal and each described harmonic reference signal The convolution nuclear information of corresponding industrial frequency harmonic interference signal, obtains total frequency domain letter of all industrial frequency harmonic interference signals respectively Breath；
3rd acquisition module, for the frequency domain information according to the acquired original signal and all industrial frequency harmonic interference signals Total frequency domain information, obtain the acquired original signal removes the first signal after all industrial frequency harmonic interference signals the One estimates frequency domain information, and estimates frequency domain information according to described first, obtains the first estimation time-domain information of first signal.
11. wave filters according to claim 10, it is characterised in that first acquisition module is specifically included：
First acquisition unit, for the time-domain information according to the acquired original signal and each harmonic reference signal when Domain information, obtain each described harmonic reference signal respectively with the crosspower spectrum information of the acquired original signal；
Second acquisition unit, for the time-domain information according to each harmonic reference signal, obtains each described harmonic wave respectively The auto-power spectrum information of reference signal；
3rd acquiring unit, for the auto-power spectrum information according to each harmonic reference signal and it is corresponding it is described mutually Power spectral information, obtains the convolution nuclear information that each described harmonic reference signal distinguishes corresponding industrial frequency harmonic interference signal.
12. wave filters according to claim 11, it is characterised in that the first acquisition unit is specifically included：
First obtain subelement, for by equation below obtain each described harmonic reference signal respectively with the acquired original The cross-spectral density function of signal
S xv i ( ω ) = ∫ - ∞ + ∞ [ ∫ - ∞ + ∞ x ( t ) v i ( t + τ ) d t ] e - j ω τ d t ;
Wherein, x (t) is the time-domain function of the acquired original signal, vi(t+ τ) is described in i-th during harmonic reference signal delay Between τ time-domain function, i is 1 to the integer between p, and p is the quantity of the harmonic reference signal, and t is time variable, and ω is frequency Rate variable, j and τ are customized constant.
13. wave filters according to claim 11, it is characterised in that the second acquisition unit is specifically included：
Second obtains subelement, the Power spectral density for obtaining each harmonic reference signal by equation below
S v i v i ( ω ) = ∫ - ∞ + ∞ [ ∫ - ∞ + ∞ v i ( t ) v i ( t + τ ) d t ] e - j ω τ d t ;
Wherein, viT () is the time-domain function of harmonic reference signal described in i-th, vi(t+ τ) is harmonic reference signal described in i-th The time-domain function of delay time T, i is 1 to the integer between p, and p is the quantity of the harmonic reference signal, and t is time variable, ω is frequency variable, and j and τ is customized constant.
14. wave filters according to claim 11, it is characterised in that the 3rd acquiring unit is specifically included：
3rd obtains subelement, and corresponding industrial frequency harmonic is distinguished for obtaining each described harmonic reference signal by equation below The convolution kernel function H of interference signali(ω)：
H i ( ω ) = S xv i ( ω ) S v i v i ( ω ) ;
Wherein,The cross-spectral density function of harmonic reference signal and the acquired original signal described in i-th,The Power spectral density of harmonic reference signal described in i-th, i is 1 to the integer between p, and p is the harmonic wave The quantity of reference signal, ω is frequency variable.
15. wave filters according to claim 10, it is characterised in that second acquisition module is specifically included：
First converter unit, Fourier transformation is carried out for the time-domain function respectively to harmonic reference signal each described, is obtained The frequency-domain function V of each harmonic reference signali(ω)；
4th acquiring unit, the total frequency-domain function V for obtaining all industrial frequency harmonic interference signals by equation below (ω)：
V ( ω ) = Σ i = 1 P H i ( ω ) V i ( ω ) ;
Wherein, Hi(ω) is the convolution kernel function of the corresponding industrial frequency harmonic interference signal of harmonic reference signal described in i-th, Vi (ω) is the frequency-domain function of harmonic reference signal described in i-th, and i is 1 to the integer between p, and p is the harmonic reference signal Quantity, ω is frequency variable.
16. wave filters according to claim 10, it is characterised in that the 3rd acquisition module is specifically included：
Second converter unit, Fourier transformation is carried out for the time-domain function to the acquired original signal, obtains described original Gather the frequency-domain function X (ω) of signal；
5th acquiring unit, the first estimation frequency-domain function for obtaining first signal by equation below
S ^ ′ ( ω ) = X ( ω ) - V ( ω ) ;
Wherein, V (ω) is total frequency-domain function of all industrial frequency harmonic interference signals, and ω is frequency variable；
3rd converter unit, for carrying out inverse Fourier transform to the described first estimation frequency-domain function, obtains first signal First estimation time-domain function.
17. wave filters according to claim 10, it is characterised in that the wave filter also includes：
4th acquisition module, for estimating frequency domain information according to described first, obtains the first estimation amplitude of first signal Spectrum and first estimates phase spectrum；
Module is rejected, for the frequency domain information according to the acquired original signal, the amplitude spectrum of the acquired original signal is obtained, And the amplitude spectrum to the acquired original signal carries out singular point rejecting treatment, obtains the second estimation amplitude spectrum；
Replacement module, for the first estimation amplitude spectrum of the described first estimation frequency domain information to be replaced with into the second estimation amplitude Spectrum, and keep described first to estimate that phase spectrum is constant, obtain the second estimation frequency domain information of first signal；
6th acquisition module, for estimating frequency domain information according to described second, obtains the second estimation time domain of first signal Information.
18. wave filters according to claim 17, it is characterised in that the rejecting module is specifically included：
Culling unit, singular point rejecting is carried out for the amplitude spectrum to the acquired original signal by the way of robust iterative, And to rejecting point assignment again by the way of smooth compensating, obtain the second estimation amplitude spectrum.
CN201510547011.6A 2015-08-31 2015-08-31 A kind of method and wave filter for suppressing industrial frequency harmonic interference Active CN105811921B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510547011.6A CN105811921B (en) 2015-08-31 2015-08-31 A kind of method and wave filter for suppressing industrial frequency harmonic interference

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510547011.6A CN105811921B (en) 2015-08-31 2015-08-31 A kind of method and wave filter for suppressing industrial frequency harmonic interference

Publications (2)

Publication Number Publication Date
CN105811921A CN105811921A (en) 2016-07-27
CN105811921B true CN105811921B (en) 2017-06-23

Family

ID=56466130

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510547011.6A Active CN105811921B (en) 2015-08-31 2015-08-31 A kind of method and wave filter for suppressing industrial frequency harmonic interference

Country Status (1)

Country Link
CN (1) CN105811921B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108332063B (en) * 2018-01-29 2020-04-24 中国科学院声学研究所 Pipeline leakage positioning method based on cross correlation
CN112672367A (en) * 2020-12-16 2021-04-16 上海擎昆信息科技有限公司 Harmonic interference elimination method and device and receiving equipment

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE19822957C1 (en) * 1998-05-22 2000-05-25 Deutsch Zentr Luft & Raumfahrt Method for the detection and suppression of interference signals in SAR data and device for carrying out the method
CN103457638B (en) * 2013-09-11 2015-04-22 大连理工大学 Restraining device and restraining method for burst impulse noise of power line communication channel
CN104702244A (en) * 2013-12-05 2015-06-10 中国科学院深圳先进技术研究院 Adaptive filter for filtering power frequency interference in electromyography signal based on EEMD (Ensemble Empirical Mode Decomposition) algorithm

Also Published As

Publication number Publication date
CN105811921A (en) 2016-07-27

Similar Documents

Publication Publication Date Title
Ur Rehman et al. Filter bank property of multivariate empirical mode decomposition
Wang et al. Parallel-type fractional zero-phase filtering for ECG signal denoising
EP3696814A1 (en) Speech enhancement method and apparatus, device and storage medium
Chen et al. Removal of power line interference from ECG signals using adaptive notch filters of sharp resolution
CN103576060A (en) Partial discharge signal denoising method based on wavelet adaptive threshold
CN110680308B (en) Electrocardiosignal denoising method based on fusion of improved EMD and threshold method
CN106680874B (en) Harmonic noise drawing method based on wave configuration feature rarefaction modeling
CN105607125A (en) Seismic data noise suppression method based on block matching algorithm and singular value decompression
CN105811921B (en) A kind of method and wave filter for suppressing industrial frequency harmonic interference
Magsi et al. Analysis of signal noise reduction by using filters
CN106680876A (en) Joint denoising method for seismic data
CN106383107A (en) Spectral signal noise reduction method based on smoothing filtering with combination of VS-LMS
CN104635264B (en) The processing method of earthquake data before superposition and equipment
Wang et al. Fractional compound integral with application to ECG signal denoising
CN109724693A (en) A kind of fusion spectrum denoising method based on stationary wavelet
Figuera et al. Spectrally adapted Mercer kernels for support vector nonuniform interpolation
Tian et al. Early fault feature extraction for rolling bearings using adaptive variational mode decomposition with noise suppression and fast spectral correlation
CN107203972A (en) A kind of Ultrasonic Image Denoising and Enhancement Method
CN105046654B (en) A kind of adaptive non-local mean noise-reduction method of electrocardiosignal based on particle group optimizing
CN108020761B (en) A kind of Denoising of Partial Discharge
CN111631707A (en) Method, device and equipment for filtering baseline drift in electrocardiosignal and storage medium
CN111239565A (en) Oil-filled casing partial discharge pulse signal processing method and system based on layered denoising model
CN110269608A (en) A kind of method, apparatus and readable storage medium storing program for executing removing signal interference
Li et al. Robust suppression of nonstationary power-line interference in electrocardiogram signals
CN104156925A (en) Processing method and system used for carrying out speckle removing and boundary enhancing on ultrasound image

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