CN108344501A - Resonance identification and removing method and device in a kind of application of signal correlation - Google Patents
Resonance identification and removing method and device in a kind of application of signal correlation Download PDFInfo
- Publication number
- CN108344501A CN108344501A CN201810083684.4A CN201810083684A CN108344501A CN 108344501 A CN108344501 A CN 108344501A CN 201810083684 A CN201810083684 A CN 201810083684A CN 108344501 A CN108344501 A CN 108344501A
- Authority
- CN
- China
- Prior art keywords
- function
- resonance
- frequency
- correlation
- cross
- 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.)
- Withdrawn
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
Abstract
Resonance identification and removing method and system in a kind of application of signal correlation, the method includes：Two sensors are arranged in examining system, pick up two echo signals and calculate its cross-spectral density function, using the frequency response function of two object function of broad sense cross-correlation ROTH methods weighted calculation；The amplitude and phase of frequency response function are analyzed, if there are formants in its amplitude, and there are the offsets of the integral multiple of π for the phase of respective frequencies, then judge there is resonance, read corresponding two intrinsic frequencies of formant；Two frequency response functions of vibration pickup are calculated for two intrinsic frequencies；The cross-spectral density function of frequency domain weighting is multiplied by the ratio of the frequency response function of vibration pickup to eliminate resonance, then carries out inverse Fourier transform to it and calculates broad sense cross-correlation function；According to the main peak value of cross-correlation function, time delay estimation corresponding with main peak value is exported.The method of the present invention can improve associated time delays and estimate accuracy in engineer application, robustness, practicability.
Description
Technical field
The present invention relates to field of measuring technique, more particularly to resonance identification and removing method in a kind of signal correlation application
And device.
Background technology
Associated time delays estimation is the hot spot of field of signal processing research, and key technology has been applied to noise source and target
The engineerings such as positioning, mechanical fault diagnosis are in practice.Its principle is：It measures and is formed after different propagation paths by source signal first
Two-way echo signal, then it is filtered rear basic computing cross-correlation to signal collected, and then institute is determined according to maximum correlation peak
Acquire the time difference (time delay) between signal.Under strong jamming operating mode, such as low signal-to-noise ratio, strong reflection, multi-source positioning, substantially mutually
The precision of relevant time delay estimated value is relatively low.
Document [1] (Knapp C H, Carter G.C.The generalized correlation method for
estimation of time delay[J].IEEE Transactions on Acoustic Speech on Signal
Processing,1976:320-327.) sought using broad sense cross correlation algorithm by doing Fourier transformation to two-way echo signal
Cross-spectral density function obtains cross-correlation function to doing inverse Fourier transform after its frequency domain weighting.Common broad sense cross correlation algorithm
For ROTH, WIENER, PHAT, SCOT and ML etc..Basic cross correlation algorithm is compared, broad sense cross-correlation can more effectively sharpen main peak
(maximum correlation peak).However practical time delay estimation effect is unsatisfactory, the reason is that, in the measurement of acquisition two-way echo signal
Cheng Zhong forms new coupled system when sensor is directly contacted with examining system.Once resonating, in echo signal also
Resonant component will be mixed with；Second, because of the presence of resonant component, and then filter parameter setting is affected, practical survey will be caused
Associated time delays estimation accuracy declines in amount.
Invention content
It is an object of the invention to overcome drawbacks described above existing for presently relevant delay time estimation method, it is proposed that a kind of signal
Resonance identification and removing method and device, the method can realize the identification of resonance, utilize signal correlation in correlation application
The method that resonant component is eliminated in, the echo signal after being decoupled.
To achieve the goals above, the present invention proposes resonance identification and removing method in a kind of application of signal correlation,
The method includes：
Two sensors are arranged in examining system, pick up two echo signals and calculate its cross-spectral density function,
Using the frequency response function of two object function of broad sense cross-correlation ROTH methods weighted calculation；The amplitude and phase of frequency response function are analyzed,
If there are formants in its amplitude, and there are the offsets of the integral multiple of π for the phase of respective frequencies, then judge there is resonance, read
Corresponding two intrinsic frequencies of formant；Two frequency response functions of vibration pickup are calculated for two intrinsic frequencies；By frequency domain
The cross-spectral density function of weighting is multiplied by the ratio of the frequency response function of vibration pickup to eliminate resonance, then carries out inverse Fourier transform to it
Calculate broad sense cross-correlation function；According to the main peak value of cross-correlation function, time delay estimation corresponding with main peak value is exported.
As a kind of improvement of the above method, the described method comprises the following steps：
Step 1) arranges two sensors in examining system, and pickup echo signal is respectively xk(t), k=1,2, it calculates
Its cross-spectral density function：
Wherein, * indicates conjugation；X1(ω)、X2(ω) indicates x respectively1(t)、x2(t) Fourier transformation；| | indicate width
Value, φ (ω) indicate phase；
Step 2) carries out frequency domain weighting using broad sense cross-correlation ROTH methods to cross-spectral density, obtains the frequency of two object functions
Function is rung, is represented by：
Wherein, S11(ω) indicates echo signal x1(t) autopower spectral density；
Step 3) analyzes the amplitude and phase of frequency response function, if there are formant in its amplitude, and the phase of respective frequencies
There are the offsets of the integral multiple of π, then judge there is resonance, read the corresponding frequency of formant and are denoted as natural frequency ω1And ω2；
Step 4) is directed to each intrinsic frequency, and intrinsic frequency is arranged and damping ratio, the frequency response function for calculating vibration pickup are：
Wherein, ζ is damping ratio；
The cross-spectral density function of frequency domain weighting is multiplied by the ratio of the frequency response function of vibration pickup to eliminate resonance by step 5), then
Inverse Fourier transform is carried out to it calculates broad sense cross-correlation function；
Correlation function of the echo signal after removal is resonated is expressed as：
Step 6) exports time delay estimation corresponding with main peak value according to the main peak value of cross-correlation function.
As a kind of improvement of the above method, the described method comprises the following steps：
Step 1) arranges two sensors in examining system, and pickup echo signal is respectively xk(t), k=1,2, it calculates
Its cross-spectral density function：
Wherein, * indicates conjugation；X1(ω)、X2(ω) indicates x respectively1(t)、x2(t) Fourier transformation；| | indicate width
Value, φ (ω) indicate phase；
Step 2) carries out frequency domain weighting using broad sense cross-correlation ROTH methods to cross-spectral density, obtains the frequency of two object functions
Function is rung, is represented by：
Wherein S22(ω) indicates echo signal x2(t) autopower spectral density；
Step 3) analyzes the amplitude and phase of frequency response function, if there are formant in its amplitude, and the phase of respective frequencies
There are the offsets of the integral multiple of π, then judge there is resonance, read the corresponding frequency of formant and are denoted as natural frequency ω1And ω2；
Step 4) is directed to each intrinsic frequency, and intrinsic frequency is arranged and damping ratio, the frequency response function for calculating vibration pickup are：
Wherein, ζ is damping ratio；
The cross-spectral density function of frequency domain weighting is multiplied by the ratio of the frequency response function of vibration pickup to eliminate resonance by step 5), then
Inverse Fourier transform is carried out to it calculates broad sense cross-correlation function；
Correlation function of the echo signal after removal is resonated is expressed as：
Step 6) exports time delay estimation corresponding with main peak value according to the main peak value of cross-correlation function.
Resonance identification and system is eliminated in a kind of application of signal correlation, the system comprises：
Analog signal is converted to digital signal and put by input module for inputting the echo signal picked up by sensor
Greatly, signal x is obtainedk(t), k=1,2；
Resonate cancellation module, for judging that the echo signal of input module whether there is resonant component, if there is resonance point
Amount then carries out resonance and eliminates operation；
Filter module is connect with resonance cancellation module, for being filtered to the echo signal after elimination resonant component, into
One step filters out the interference signals such as ambient noise present in echo signal and reflection sound wave；With
Computing module is connect with the identification unit in filter module and resonance cancellation module, for mesh after resonating to removal
The frequency domain weighting cross-spectral density for marking signal carries out inverse Fourier transform, obtains broad sense cross-correlation, according to maximum correlation peak output pair
The time delay estimation answered.
As a kind of improvement of above system, the resonance cancellation module includes：First identification unit and the first decoupling are single
Member；
First identification unit, connect with input module, the echo signal for reading amplification, and mutual using broad sense
Close ROTH methods calculate two object functions frequency response function expression formula be：
Wherein, S11(ω) indicates signal x1(t) autopower spectral density；S12(ω) is cross-spectral density function：
Wherein, * indicates conjugation；X1(ω)、X2(ω) indicates x respectively1(t)、x2(t) Fourier transformation；| | indicate width
Value, φ (ω) indicate phase；
If there is the formant of mutation in the amplitude of frequency response function, and phase corresponds in intrinsic frequency that there are π integral multiples
Offset then judges that, there are resonant component, the corresponding intrinsic frequency of output resonance peak is denoted as ω1And ω2；If there is no resonance point
Amount, then be directly entered filter module；
First decoupling unit, connect with the first identification unit, for intrinsic frequency and Damping calculating vibration pickup to be arranged
Frequency response function, and eliminate the resonant component in echo signal；
The frequency response function of resonator is calculated according to following formula：
Wherein, ζ is damping ratio；
Resonance, which is carried out, according to following formula eliminates operation：
As a kind of improvement of above system, the resonance cancellation module includes：Second identification unit and the second decoupling are single
Member；
Second identification unit, connect with input module, the echo signal for reading amplification, and mutual using broad sense
Close ROTH methods calculate two object functions frequency response function expression formula be：
Wherein, S22(ω) is echo signal x2(t) autopower spectral density；S21(ω) is cross-spectral density function：
Wherein, * indicates conjugation；X1(ω)、X2(ω) indicates x respectively1(t)、x2(t) Fourier transformation；| | indicate width
Value, φ (ω) indicate phase；
If there is the formant of mutation in the amplitude of frequency response function, and phase corresponds in intrinsic frequency that there are π integral multiples
Offset then judges that, there are resonant component, the corresponding intrinsic frequency of output resonance peak is denoted as ω1And ω2；If there is no resonance point
Amount, then be directly entered filter module；
Second decoupling unit, connect with the second identification unit, for intrinsic frequency and Damping calculating vibration pickup to be arranged
Frequency response function, and eliminate the resonant component in echo signal；
The frequency response function of resonator is calculated according to following formula：
Wherein, ζ is damping ratio；
Resonance, which is carried out, according to following formula eliminates operation：
Advantage of the invention is that：
1, existing Time Delay Estimation Based does not consider to resonate, and the identification of resonance may be implemented in the method for the present invention, utilizes letter
The method that resonant component is eliminated in number correlation application, the echo signal after being decoupled；
2, method of the invention can improve associated time delays and estimate accuracy in engineer application, robustness, practicability.
Description of the drawings
Fig. 1 is the flow chart of the method for the present invention；
Fig. 2 is the coupled system schematic diagram of the present invention；
Fig. 3 is the structure chart of the device of the invention；
Fig. 4, which is resonance, to be influenced to eliminate the schematic diagram of the amplitude of forward and backward frequency response function；
Fig. 5, which is resonance, to be influenced to eliminate the schematic diagram of the phase of forward and backward frequency response function；
Fig. 6, which is resonance, influences ROTH correlated results schematic diagrames after eliminating；
Fig. 7, which is resonance, influences ROTH correlated results schematic diagrames before eliminating.
Specific implementation mode
As shown in Figure 1, resonance identification and removing method in a kind of signal correlation application, include the following steps：
Step 1：The placement sensor in examining system, pickup echo signal are xk(t), (k=1,2) calculates its mutual work(
Rate spectral density function：
Or
Wherein * indicates conjugation；X1(ω)、X2(ω) indicates x respectively1(t)、x2(t) Fourier transformation；| | indicate width
Value, φ (ω) indicate phase.
Step 2：Frequency domain weighting is carried out to cross-spectral density using broad sense cross-correlation ROTH methods, obtains the frequency of two object functions
Function is rung, is represented by：
Or
The ROTH methods are respectively 1/S using weighting function11(ω)、1/S22(ω)；S11(ω)、S22(ω) difference table
Show echo signal x1(t)、x2(t) autopower spectral density.
Step 3：The amplitude and phase of frequency response function are analyzed, if there are formant in its amplitude, and the phase of respective frequencies
There are the offsets of the integral multiple of π, then judge there is resonance, read the corresponding frequency of formant and are denoted as natural frequency ω1And ω2；
Step 4：For each intrinsic frequency, setting intrinsic frequency and damping ratio, calculate the frequency of vibration pickup in formula (3)
Ringing function is：
Wherein, ζ is damping ratio；
Step 5：The cross-spectral density function (i.e. frequency response function) of frequency domain weighting is multiplied by the ratio of the frequency response function of vibration pickup
It resonates to eliminate, then carries out inverse Fourier transform to it and calculate broad sense cross-correlation function.Phase of the echo signal after removal is resonated
Function is closed to be represented by：
Or
Step 6：According to the main peak value of cross-correlation function, time delay estimation corresponding with main peak value is exported.
When sensor and examining system contact, form a new coupled system, as shown in Figure 1, by examining system and
Vibration pickup is composed in series.L (t) is source signal in figure, by different propagation path dkIt is later echo signal x by sensor pickupk
(t), (k=1,2), the cross-spectral density function between echo signal are：
Wherein*Indicate complex conjugate；Sll(ω) is the autopower spectral density of source signal；Dk(ω)、Hk(ω,dk) it is respectively to pick up
Shake device and the frequency response function of examining system, k=1, and 2.
Each vibration pickup is made of rigidity k, damping c, quality m, and corresponding frequency response function is：
WhereinFor intrinsic frequency,For damping ratio.
If resonant component caused by being coupled there are system in echo signal, the amplitude of the frequency response function of echo signal is intrinsic
There will be apparent formants in frequency, while the offset of π integral multiples then occurring in phase spectrum.By formula (1a), (1b) as it can be seen that altogether
It shakes and will also seriously affect the amplitude and phase information of cross-spectral density function.It is delayed when being sought using correlation technique, needs first
Echo signal is filtered through row, Fourier transformation then is done to two paths of signals and acquires cross-spectral density, then to its into
Inverse Fourier transform is done after row frequency domain weighting can be obtained broad sense cross-correlation function.In above-mentioned related algorithm implementation process, when
It resonates, has not only influenced the setting of filter parameter, but also change cross-spectral density, eventually lead to associated time delays estimated result
It is inaccurate.
As shown in Fig. 2, the present invention proposes that a kind of device for eliminating resonant component, described device include：
Input module, for inputting the echo signal picked up by sensor, wherein should include A/D converter and amplification
Analog signal is converted to digital signal and amplified by device；
Resonate cancellation module, connect with input module, for judging echo signal with the presence or absence of resonant component, and, if
There are resonant components then to carry out resonating removing division operation；
The resonance cancellation module includes：
Identification unit is connect with input module, for reading amplification target signal, and uses broad sense cross-correlation ROTH methods
The frequency response function expression formula for calculating two object functions is：
Or
If there is the formant of mutation in the amplitude of frequency response function, and phase corresponds in intrinsic frequency that there are π integral multiples
Offset then judges that, there are resonant component, the corresponding intrinsic frequency of output resonance peak is denoted as ω1And ω2；If there is no resonance point
Amount, then be directly entered filter module；
The decoupling unit, connect with identification unit, the frequency response letter for intrinsic frequency and Damping calculating vibration pickup to be arranged
Number, and eliminate the resonant component in echo signal；
The frequency response function of resonator is calculated according to following formula：
Resonance, which is carried out, according to following formula eliminates operation：
Or
Filter module is connect with resonance cancellation module, for being filtered to the echo signal after elimination resonant component, into
One step filters out the interference signals such as ambient noise present in echo signal and reflection sound wave；
Computing module is connect with the identification unit in filter module and resonance cancellation module, for mesh after resonating to removal
The frequency domain weighting cross-spectral density for marking signal carries out inverse Fourier transform, obtains broad sense cross-correlation, according to maximum correlation peak output pair
The time delay estimation answered.
In order to illustrate the validity of this programme, this programme is verified with specific embodiment below：
Step 1：A sensor is respectively placed at examining system two, is picked up the time domain echo signal of 10s, is input to input
Module, sample frequency 500Hz carry out it 1024 points of Fast Fourier Transform (FFT), calculate echo signal x2(t) and x1(t)
Cross-spectral density function：
Wherein * indicates complex conjugate；
Step 2：The echo signal of input is sent to identification unit, carrying out frequency domain to cross-spectral density using ROTH methods adds
Power obtains frequency response function：
Step 3：The amplitude and phase of frequency response function are analyzed, the amplitude spectrum in the present embodiment has the formant of mutation, right
The offset of the phase presence ± π of frequency is answered, as shown in figure 4, judging that there are resonant component, the corresponding frequencies of peak value in echo signal
It is 56Hz and 83Hz, is denoted as natural frequency ω respectively1And ω2；
Step 4：Identification unit judges that, there are resonant component in echo signal, output intrinsic frequency is ω1=56Hz with
And ω2=83Hz.In the present embodiment, intrinsic frequency and dampingratioζ=0.022 are set in decoupling unit, according to following formula point
Not Ji Suan resonator frequency response function：
Step 5：By the frequency response function H in the ratio of two resonator frequency response functions and ROTH methods21(ω), which is multiplied, to carry out
Operation is eliminated in resonance：
As shown in Figure 4 and Figure 5, resonant component of the amplitude of frequency response function in intrinsic frequency is eliminated after system decoupling,
The offset in phase is also eliminated simultaneously.
Step 6：The frequency response function of echo signal after decoupling is read to filter module and is filtered, is arranged in the present embodiment
The bandwidth of bandpass filter is 5-90Hz, and inverse Fourier transform is carried out after filtering and obtains cross-correlation function.Fig. 6, which provides resonance, to be influenced
ROTH correlated results after elimination exports corresponding time delay estimadon according to maximum correlation peak, what the present embodiment finally exported
Time delay is estimated as -90ms, and relatively true time delay value -92ms, time delay error is about 2%, and comparison resonance influences the correlation before eliminating
As a result, effectively removing resonance interference, peak value obviously sharpens, as shown in Figure 7.
It should be noted last that the above examples are only used to illustrate the technical scheme of the present invention and are not limiting.Although ginseng
It is described the invention in detail according to embodiment, it will be understood by those of ordinary skill in the art that, to the technical side of the present invention
Case is modified or replaced equivalently, and without departure from the spirit and scope of technical solution of the present invention, should all be covered in the present invention
Right in.
Claims (6)
1. resonance identification and removing method in a kind of signal correlation application, the method includes：
Two sensors are arranged in examining system, pick up two echo signals and calculate its cross-spectral density function, are used
The frequency response function of two object function of broad sense cross-correlation ROTH methods weighted calculation；The amplitude and phase of frequency response function are analyzed, if its
There are formants in amplitude, and there are the offsets of the integral multiple of π for the phase of respective frequencies, then judge there is resonance, read resonance
Corresponding two intrinsic frequencies in peak；Two frequency response functions of vibration pickup are calculated for two intrinsic frequencies；By the mutual of frequency domain weighting
Spectral density function is multiplied by the ratio of the frequency response function of vibration pickup to eliminate resonance, then carries out inverse Fourier transform to it and calculate broad sense
Cross-correlation function；According to the main peak value of cross-correlation function, time delay estimation corresponding with main peak value is exported.
2. resonance identification and removing method in signal correlation application according to claim 1, the method includes following
Step：
Step 1) arranges two sensors in examining system, and pickup echo signal is respectively xk(t), k=1,2, calculate its mutual work(
Rate spectral density function：
Wherein, * indicates conjugation；X1(ω)、X2(ω) indicates the Fourier transformation of x1 (t), x2 (t) respectively；| | indicate amplitude, φ
(ω) indicates phase；
Step 2) carries out frequency domain weighting using broad sense cross-correlation ROTH methods to cross-spectral density, obtains the frequency response letter of two object functions
Number, is represented by：
Wherein, S11(ω) indicates echo signal x respectively1(t) autopower spectral density；
Step 3) analyzes the amplitude and phase of frequency response function, if there are formants in its amplitude, and there are π for the phase of respective frequencies
Integral multiple offset, then judge there is resonance, read the corresponding frequency of formant and be denoted as natural frequency ω1And ω2；
Step 4) is directed to each intrinsic frequency, and intrinsic frequency is arranged and damping ratio, the frequency response function for calculating vibration pickup are：
Wherein, ζ is damping ratio；ω is angular frequency；
The cross-spectral density function of frequency domain weighting is multiplied by the ratio of the frequency response function of vibration pickup to eliminate resonance, then to it by step 5)
It carries out inverse Fourier transform and calculates broad sense cross-correlation function；
Correlation function of the echo signal after removal is resonated is expressed as：
Step 6) exports time delay estimation corresponding with main peak value according to the main peak value of cross-correlation function.
3. resonance identification and removing method in signal correlation application according to claim 1, the method includes following
Step：
Step 1) arranges two sensors in examining system, and pickup echo signal is respectively xk(t), k=1,2, calculate its mutual work(
Rate spectral density function：
Wherein, * indicates conjugation；X1(ω)、X2(ω) indicates x1 (t), x respectively2(t) Fourier transformation；| | indicate amplitude, φ
(ω) indicates phase；
Step 2) carries out frequency domain weighting using broad sense cross-correlation ROTH methods to cross-spectral density, obtains the frequency response letter of two object functions
Number, is represented by：
Wherein, S22(ω) indicates echo signal x2(t) autopower spectral density；
Step 3) analyzes the amplitude and phase of frequency response function, if there are formants in its amplitude, and there are π for the phase of respective frequencies
Integral multiple offset, then judge there is resonance, read the corresponding frequency of formant and be denoted as natural frequency ω1And ω2；
Step 4) is directed to each intrinsic frequency, and intrinsic frequency is arranged and damping ratio, the frequency response function for calculating vibration pickup are：
Wherein, ζ is damping ratio；ω is angular frequency；
The cross-spectral density function of frequency domain weighting is multiplied by the ratio of the frequency response function of vibration pickup to eliminate resonance, then to it by step 5)
It carries out inverse Fourier transform and calculates broad sense cross-correlation function；
Correlation function of the echo signal after removal is resonated is expressed as：
Step 6) exports time delay estimation corresponding with main peak value according to the main peak value of cross-correlation function.
4. resonance identification and elimination system in a kind of signal correlation application, which is characterized in that the system comprises：
Analog signal is converted to digital signal and amplified, obtained by input module for inputting the echo signal picked up by sensor
To signal xk(t), k=1,2；
Resonate cancellation module, for judging that the echo signal of input module whether there is resonant component, if there are resonant component
It carries out resonance and eliminates operation；
Filter module is connect with resonance cancellation module, for being filtered to the echo signal after elimination resonant component, further
Filter out the interference signals such as ambient noise present in echo signal and reflection sound wave；
Computing module is connect with the identification unit in filter module and resonance cancellation module, is believed for target after resonating to removal
Number frequency domain weighting cross-spectral density carry out inverse Fourier transform, obtain broad sense cross-correlation, exported according to maximum correlation peak corresponding
Time delay is estimated.
5. resonance identification and elimination system in signal correlation application according to claim 4, which is characterized in that described total
The cancellation module that shakes includes：First identification unit and the first decoupling unit；
First identification unit, connect with input module, the echo signal for reading amplification, and uses broad sense cross-correlation
ROTH methods calculate two object functions frequency response function expression formula be：
Wherein, S11(ω) indicates signal x1(t) autopower spectral density；S12(ω) is cross-spectral density function：
Wherein, * indicates conjugation；X1(ω)、X2(ω) indicates x respectively1(t)、x2(t) Fourier transformation；| | indicate amplitude, φ
(ω) indicates phase；
If there is the formant of mutation in the amplitude of frequency response function, and phase corresponds in intrinsic frequency there are the offset of π integral multiples,
Then judge that, there are resonant component, the corresponding intrinsic frequency of output resonance peak is denoted as ω1And ω2；If resonant component is not present,
It is directly entered filter module；
First decoupling unit, connect with the first identification unit, the frequency for intrinsic frequency and Damping calculating vibration pickup to be arranged
Function is rung, and eliminates the resonant component in echo signal；
The frequency response function of resonator is calculated according to following formula：
Wherein, ζ is damping ratio；
Resonance, which is carried out, according to following formula eliminates operation：
6. resonance identification and elimination system in signal correlation application according to claim 4, which is characterized in that described total
The cancellation module that shakes includes：Second identification unit and the second decoupling unit；
Second identification unit, connect with input module, the echo signal for reading amplification, and uses broad sense cross-correlation
ROTH methods calculate two object functions frequency response function expression formula be：
Wherein, S22(ω) is echo signal x2(t) autopower spectral density；S21(ω) is cross-spectral density function：
Wherein, * indicates conjugation；X1(ω)、X2(ω) indicates x respectively1(t)、x2(t) Fourier transformation；| | indicate amplitude, φ
(ω) indicates phase；
If there is the formant of mutation in the amplitude of frequency response function, and phase corresponds in intrinsic frequency there are the offset of π integral multiples,
Then judge that, there are resonant component, the corresponding intrinsic frequency of output resonance peak is denoted as ω1And ω2；If resonant component is not present,
It is directly entered filter module；
Second decoupling unit, connect with the second identification unit, the frequency for intrinsic frequency and Damping calculating vibration pickup to be arranged
Function is rung, and eliminates the resonant component in echo signal；
The frequency response function of resonator is calculated according to following formula：
Wherein, ζ is damping ratio；
Resonance, which is carried out, according to following formula eliminates operation：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810083684.4A CN108344501A (en) | 2018-01-29 | 2018-01-29 | Resonance identification and removing method and device in a kind of application of signal correlation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810083684.4A CN108344501A (en) | 2018-01-29 | 2018-01-29 | Resonance identification and removing method and device in a kind of application of signal correlation |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108344501A true CN108344501A (en) | 2018-07-31 |
Family
ID=62961756
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810083684.4A Withdrawn CN108344501A (en) | 2018-01-29 | 2018-01-29 | Resonance identification and removing method and device in a kind of application of signal correlation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108344501A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114109798A (en) * | 2021-11-26 | 2022-03-01 | 广东美的暖通设备有限公司 | Frequency determination method, frequency determination device, compressor system and storage medium |
| CN114112006A (en) * | 2021-11-26 | 2022-03-01 | 中科传启（苏州）科技有限公司 | Noise monitoring method and device and electronic equipment |
| CN116299147A (en) * | 2023-03-13 | 2023-06-23 | 中国科学院声学研究所 | One-dimensional structure internal sound source positioning method based on acoustic coherence technology |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106502199A (en) * | 2016-10-20 | 2017-03-15 | 华中科技大学 | A kind of frequency response function recognition methodss under mechanized equipment structure running status |
-
2018
- 2018-01-29 CN CN201810083684.4A patent/CN108344501A/en not_active Withdrawn
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106502199A (en) * | 2016-10-20 | 2017-03-15 | 华中科技大学 | A kind of frequency response function recognition methodss under mechanized equipment structure running status |
Non-Patent Citations (1)
| Title |
|---|
| YAN GAO ET AL.: "Improving the shape of the cross-correlation function for leak detection in a plastic water distribution pipe using acoustic signals", 《APPLIED ACOUSTICS》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114109798A (en) * | 2021-11-26 | 2022-03-01 | 广东美的暖通设备有限公司 | Frequency determination method, frequency determination device, compressor system and storage medium |
| CN114112006A (en) * | 2021-11-26 | 2022-03-01 | 中科传启（苏州）科技有限公司 | Noise monitoring method and device and electronic equipment |
| CN114109798B (en) * | 2021-11-26 | 2023-08-29 | 广东美的暖通设备有限公司 | Frequency determining method, frequency determining device, compressor system and storage medium |
| CN116299147A (en) * | 2023-03-13 | 2023-06-23 | 中国科学院声学研究所 | One-dimensional structure internal sound source positioning method based on acoustic coherence technology |
| CN116299147B (en) * | 2023-03-13 | 2023-11-28 | 中国科学院声学研究所 | One-dimensional structure internal sound source positioning method based on acoustic coherence technology |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20020039425A1 (en) | Method and apparatus for removing noise from electronic signals | |
| CN105472191B (en) | A kind of method and apparatus tracking echo delay time | |
| CN108344501A (en) | Resonance identification and removing method and device in a kind of application of signal correlation | |
| US20100026858A1 (en) | Noise extraction device using microphone | |
| CN101426169B (en) | Time-domain tracking filter quickly detecting acoustic response parameter of sounding body and system | |
| CN107644650A (en) | A kind of improvement sound localization method based on progressive serial orthogonalization blind source separation algorithm and its realize system | |
| US10524049B2 (en) | Method for accurately calculating the direction of arrival of sound at a microphone array | |
| WO2013140399A1 (en) | System and method for robust estimation and tracking the fundamental frequency of pseudo periodic signals in the presence of noise | |
| CN114061678B (en) | Digital driving method for Coriolis flowmeter | |
| TWI363567B (en) | Apparatus and method arranged to analyse a signal comprising a series of symbols | |
| JP2016006536A (en) | Complex acoustic resonance speech analysis system | |
| CN111856401A (en) | Time delay estimation method based on cross-spectrum phase fitting | |
| US20030128848A1 (en) | Method and apparatus for removing noise from electronic signals | |
| WO2004056298A1 (en) | Method and apparatus for removing noise from electronic signals | |
| CN108718438A (en) | A kind of judgment method of the response signal start time of electro-acoustic product detecting system | |
| CN112462321A (en) | Sound beacon signal interception processing method based on pulse matching accumulation | |
| Huber et al. | Glottal source shape parameter estimation using phase minimization variants | |
| JP2016500847A (en) | Digital processor based complex acoustic resonance digital speech analysis system | |
| CN210347836U (en) | Wave velocity calculating device for cable oscillatory wave partial discharge detection system | |
| CN103207011A (en) | Ultrasonic precise time delay detection method based on circuit phase shift | |
| CN102682781B (en) | Self-adaptive suppression method of time varying sinusoidal interference in audio signals | |
| CN116312447B (en) | Directional noise elimination method and system | |
| CN108335704A (en) | Vagitus detection circuit and relevant detection method | |
| TWI684781B (en) | Magnetotelluric data processing method and system thereof | |
| CN114822579B (en) | Signal estimation method based on first-order differential microphone array |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20180731 |
|
| WW01 | Invention patent application withdrawn after publication |