CN110598520A - Speed measuring method for preventing strong common mode noise interference in signal - Google Patents
Speed measuring method for preventing strong common mode noise interference in signal Download PDFInfo
- Publication number
- CN110598520A CN110598520A CN201910571260.7A CN201910571260A CN110598520A CN 110598520 A CN110598520 A CN 110598520A CN 201910571260 A CN201910571260 A CN 201910571260A CN 110598520 A CN110598520 A CN 110598520A
- Authority
- CN
- China
- Prior art keywords
- signal
- calculating
- noise
- correlation
- function
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P3/00—Measuring linear or angular speed; Measuring differences of linear or angular speeds
- G01P3/64—Devices characterised by the determination of the time taken to traverse a fixed distance
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
Abstract
The invention discloses a speed measuring method for preventing strong common mode noise interference in signals, belonging to the technical field of speed measurement. Firstly, filtering high-frequency noise in a signal by adopting a digital filtering algorithm and smoothing the signal. The signal is then time-domain denoised, suppressing the noise frequency components mixed with the original signal. And calculating the displacement corresponding to the peak value of the optimal correlation function to further obtain the transit time, and calculating the correlation speed through the transit time. The invention improves the performance of the related speed measurement algorithm, enables the related speed measurement algorithm to still realize the effective measurement of the related speed in a severe environment, improves the measurement precision, widens the measurement range and improves the measurement repeatability.
Description
Technical Field
The invention belongs to the technical field of speed measurement, and particularly relates to a speed measurement method for preventing strong common-mode noise interference in signals.
Background
In the past twenty years, correlation algorithms have been widely used for speed measurement in industry. The time of flight is obtained by calculating the peak value of the cross correlation function of the output signals of the two sensors arranged in front of and behind the moving direction of the measured object, and then the moving speed is obtained according to the distance between the two sensors. However, the signal from the sensor is susceptible to various external and internal noises, and under extremely severe conditions, the signal is partially or completely submerged in strong noises, so that the signal quality is remarkably reduced, and the measurement is failed. It is therefore important to eliminate noise interference. In selecting the denoising method, the type of noise mixed with the original signal needs to be considered. The most common noises are mains frequency disturbances, vibrations of the mechanical system and intrinsic random noises. The power frequency interference mainly comes from electromagnetic interference of a power line and incorrect grounding of a signal conditioning circuit. The power line noise is centered at 50Hz or 60Hz, the bandwidth does not exceed 1Hz, and most is shielded by the grounded metal of the sensor. However, the addition of a metal layer covering the electrodes and the signal conditioning circuitry does not completely shield external electrical interference, and conventional notch filters cannot be used to remove residual noise because their frequency components (50Hz or 60Hz) may be within the spectrum of the original signal. The second type of noise is caused by vibrations of the mechanical system. The disturbing frequencies of this type of noise are generally composed of high frequency components caused by the vibrations of the connecting cable, the rotating machine main frequency and its harmonics. The third category is inherently random noise, which approximates white noise with a continuous spectral distribution. Therefore, a preferred correlation velocity measurement method is needed to remove these three common noise interferences and achieve effective measurement of correlation velocity.
Disclosure of Invention
The invention aims to provide a speed measurement method for preventing strong common-mode noise interference in signals, which comprises the following steps:
filtering high-frequency noise in the signal by adopting a digital filtering algorithm;
smoothing the filtered signal;
calculating a correlation function of the smoothed signal;
calculating the displacement corresponding to the peak value of the correlation function;
determining a transit time according to the displacement;
and calculating to obtain the related speed according to the transition time.
The invention has the beneficial effects that: the technology is suitable for any measurement based on the correlation algorithm with strong periodic noise, improves the performance of the correlation velocity measurement algorithm, enables the correlation velocity measurement algorithm to still realize effective measurement of the correlation velocity in a severe environment, improves the measurement precision, widens the measurement range, and improves the measurement repeatability.
Drawings
FIG. 1 shows a filtered signal S1nAnd S2nOf (2) (where the signal S is exemplified)1nAnd S2nWith strong common mode noise).
FIG. 2 shows the filtered signal S1nAnd S2nA cross correlation function graph of (a).
FIG. 3 shows the filtered signal S1nAnd S2nPreferably a plot of the correlation function.
Fig. 4 is a flow chart of a velocity measurement method for preventing strong common mode noise interference in a signal.
Detailed Description
The invention is described below with reference to the following figures and examples: a velocity measurement method for preventing strong common mode noise interference in signals comprises the following steps: step 1, filtering high-frequency noise in a signal by adopting a digital filtering algorithm; smoothing the filtered signal; step 2, calculating a correlation function of the smoothed signal; step 3, calculating the displacement corresponding to the peak value of the correlation function; determining a transit time according to the displacement; and 4, calculating to obtain the related speed according to the transit time.
Specifically, step 1 filters the high-frequency noise in the signal by using a digital filtering algorithm, which combines a cut-off frequency method based on fourier transform and a median filter, and smoothes the signal while removing the high-frequency noise component. In one example, a 5 th order IIR butterworth filter is used to suppress noise outside of the signal frequency, and then a median filter with a window size of 50 is used to smooth the signal.
Specifically, step 2, a first noise-free signal S is calculated1(t) autocorrelation function RS1S1A first noise-free signal S1(t) and a second noise-free signal S2(t) cross-correlation function RS1S2Can be written as
Where N is the number of samples in the correlation calculation, and m (m 0.., N) is the number of delay points; if the first noiseless signal S1(t) and a second noise-free signal S2(t) is disturbed by strong common mode periodic noise n (t) associated therewith, the first noisy signal S1n(t) autocorrelation function RS1nS1nIs composed of
While the first noisy signal S1n(t) and a second noisy signal S2n(t) cross-correlation function RS1nS2nIs composed of
Wherein S is1n(k) And S2n(k) Are respectively digitized signals S1n(t) and S2n(t)。
FIGS. 1 and 2 are respectively filtered signals S1nAnd S2nAnd two filtered signals S1nAnd S2nThe cross correlation function of (a). Although the use of a digital low-pass filter eliminates high frequency noise interference, the filtered signal is still subject to noise interference associated with the original signal, with no clear dominant peak in the correlation function of the signal.
Since noise is periodic and signal dependent, R is not sufficiently low for signal-to-noise ratioS1n,RnS1,RnS2And RnnThese correlation terms do not go to zero and the correlation function is disturbed. To address this problem, the present disclosure proposes an optimized correlation function RmodIn which S is1nAnd (S)1n-S2n) Has a correlation function of
S2nAnd (S)1n-S2n) Has a correlation function of
The subtraction of equation (11) and equation (12) yields
Wherein R ismodIs S1nAnd S2nBy a noise-free signal S only1And S2The autocorrelation and cross-correlation functions of;
in particular, step 3, a preferred correlation function R is calculatedmodObtaining the transition time tau by the displacement corresponding to the peak value;
specifically, in step 4, the correlation velocity v is calculated as L/τ, where L is the distance between the two sensors.
FIG. 3 shows two filtered signals S1nAnd S2nPreferred correlation function, RS1S1And RS2S2At the same time delay, i.e. the period T of the noise. Cross correlation function RS1S2And RS2S1Are inverted in equation (13) and are located at the transit times τ and (T- τ), respectively. Preferred correlation function RmodThe influence of periodic noise has been significantly eliminated and the correlation peak becomes more significant, and from the distance L between the two sensors, the correlation velocity v-L/τ can be calculated.
It will be evident to those skilled in the art that the invention is not limited to the details of the foregoing illustrative embodiments, and that the present invention may be embodied in other specific forms without departing from the spirit or essential attributes thereof. The present embodiments are therefore to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein. Any reference sign in a claim should not be construed as limiting the claim concerned.
Claims (4)
1. A method for measuring a velocity to prevent strong common mode noise interference in a signal, wherein the method comprises:
filtering high-frequency noise in the signal by adopting a digital filtering algorithm;
smoothing the filtered signal;
calculating a correlation function of the smoothed signal;
calculating the displacement corresponding to the peak value of the correlation function;
determining a transit time according to the displacement;
and calculating to obtain the related speed according to the transition time.
2. The method of claim 1, wherein the digital filtering algorithm is determined according to a fourier transform-based cut-off frequency method and a median filter.
3. The method of claim 1, wherein computing the correlation function of the smoothed signal comprises:
calculating a first noise-free signal S1(t) autocorrelation function RS1S1,
Calculating a first noise-free signal S1(t) and a second noise-free signal S2(t) cross-correlation function RS1S2,
Wherein N is the number of samples in the correlation calculation, m is 0.. and N is the number of delay points;
when S is1(t) and S2(t) calculating a first noisy signal S when disturbed by strong common mode periodic noise n (t) associated therewith1n(t) autocorrelation function RS1nS1n,
Calculating a first noisy signal S1n(t) and a second noisy signal S2n(t) cross-correlation function RS1nS2n,
Wherein S is1n(k) And S2n(k) Are each S1n(t) and S2n(t) the digitized signal;
determination of S1nAnd (S)1n-S2n) Is related to function RS1n(S1n-S2n),
Determination of S2nAnd (S)1n-S2n) Is related to function RS2n(S1n-S2n),
Will correlate with function RS1n(S1n-S2n)And a correlation function RS2n(S1n-S2n)Subtracting to determine a first noisy signal S1n(t) and a second noisy signal S2n(t) optimized correlation function Rmod,
4. The method of claim 1, wherein calculating a correlation velocity from the time of flight comprises:
the correlation velocity v is calculated as L/τ, where L is the spacing of the two sensors and τ is the transit time.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910571260.7A CN110598520B (en) | 2019-06-28 | 2019-06-28 | Speed measuring method for preventing strong common mode noise interference in signal |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910571260.7A CN110598520B (en) | 2019-06-28 | 2019-06-28 | Speed measuring method for preventing strong common mode noise interference in signal |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110598520A true CN110598520A (en) | 2019-12-20 |
CN110598520B CN110598520B (en) | 2023-05-05 |
Family
ID=68852762
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910571260.7A Active CN110598520B (en) | 2019-06-28 | 2019-06-28 | Speed measuring method for preventing strong common mode noise interference in signal |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110598520B (en) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6028549A (en) * | 1998-05-22 | 2000-02-22 | Deutsches Zentrum Fur Luft-Und Raumfahrt E.V. | Process for the detection and suppression of interfering signals in S.A.R. data and device for performing this process |
US6151558A (en) * | 1998-02-10 | 2000-11-21 | Conant; James R | Ultrasonic marine speedometer system |
US20120136621A1 (en) * | 2010-11-25 | 2012-05-31 | Mitsubishi Electric Corporation | Velocity measurement apparatus capable of accurately measuring velocity of moving object relative to ground surface |
CN102680958A (en) * | 2012-06-11 | 2012-09-19 | 中国工程物理研究院计算机应用研究所 | Velocity measurement signal processing method for Doppler velocimeter |
US20140136142A1 (en) * | 2012-11-09 | 2014-05-15 | Airmar Technology Corporation | Speed sensor |
CN106226739A (en) * | 2016-07-29 | 2016-12-14 | 太原理工大学 | Merge the double sound source localization method of Substrip analysis |
-
2019
- 2019-06-28 CN CN201910571260.7A patent/CN110598520B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6151558A (en) * | 1998-02-10 | 2000-11-21 | Conant; James R | Ultrasonic marine speedometer system |
US6028549A (en) * | 1998-05-22 | 2000-02-22 | Deutsches Zentrum Fur Luft-Und Raumfahrt E.V. | Process for the detection and suppression of interfering signals in S.A.R. data and device for performing this process |
US20120136621A1 (en) * | 2010-11-25 | 2012-05-31 | Mitsubishi Electric Corporation | Velocity measurement apparatus capable of accurately measuring velocity of moving object relative to ground surface |
CN102680958A (en) * | 2012-06-11 | 2012-09-19 | 中国工程物理研究院计算机应用研究所 | Velocity measurement signal processing method for Doppler velocimeter |
US20140136142A1 (en) * | 2012-11-09 | 2014-05-15 | Airmar Technology Corporation | Speed sensor |
CN106226739A (en) * | 2016-07-29 | 2016-12-14 | 太原理工大学 | Merge the double sound source localization method of Substrip analysis |
Non-Patent Citations (1)
Title |
---|
仲崇庆 等: "高精度超声回波渡越时间算法研究", 《仪表技术与传感器》 * |
Also Published As
Publication number | Publication date |
---|---|
CN110598520B (en) | 2023-05-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8386551B2 (en) | Dynamic filtering device and method thereof | |
CN110850482B (en) | Transient electromagnetic signal-noise separation method based on variational modal decomposition principle | |
CN106199532B (en) | Based on mixing Fourier-wavelet analysis Gpr Signal noise-reduction method | |
US20080243497A1 (en) | Stationary-tones interference cancellation | |
CN115023597A (en) | Noise canceller, abnormality diagnosis device, and noise cancellation method | |
CN110967735A (en) | Self-adaptive ghost wave suppression method and system | |
CN110598520B (en) | Speed measuring method for preventing strong common mode noise interference in signal | |
CN107610055B (en) | Fourier transform spectrometer interferogram noise detection and suppression method | |
CN103604404B (en) | Acceleration signal measurement displacement method based on numerical integration | |
US20090144004A1 (en) | Method and Device for Spectral Analysis in Several Frequency Bands with Different Frequency Resolution | |
CN111239814B (en) | Shallow profile data mechanical interference suppression method based on same-phase axis frequency division tracking smoothing | |
Lu et al. | A noise reduction method of rolling bearing based on empirical wavelet transform and adaptive time frequency peak filtering | |
Reda et al. | An improved method for the processing of signals contaminated with strong common-mode periodic noise in correlation velocity measurement | |
CN109724693A (en) | A kind of fusion spectrum denoising method based on stationary wavelet | |
JP2005500605A (en) | High-pass filtering method for datasets | |
CN107907542B (en) | IVMD and energy estimation combined DSPI phase filtering method | |
CN110703089B (en) | Wavelet threshold denoising method for low-frequency oscillation Prony analysis | |
CN114285393B (en) | Method and device for eliminating reverse time migration noise | |
CN114841213B (en) | Silicon micro-resonance type accelerometer noise reduction method based on improved wavelet threshold function | |
CN110764135A (en) | Irregular seismic data full-band reconstruction method | |
CN113419276B (en) | Time-varying wavelet extraction method for self-adaptive phase estimation | |
Van Moer et al. | A novel spectral subtraction technique for cognitive radios | |
JP3036985B2 (en) | Spectrum calculation device | |
CN107228980B (en) | Method and device for measuring power grid frequency | |
CN109459788B (en) | Stratum quality factor calculation method and system |
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 | ||
CB03 | Change of inventor or designer information |
Inventor after: Yan Yong Inventor after: Carmel Reid Inventor after: Zheng Ge Inventor after: Hu Yonghui Inventor before: Yan Yong Inventor before: Carmel Reid Inventor before: Zheng Ge Inventor before: Hu Yonghui |
|
CB03 | Change of inventor or designer information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |