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 PDF

Info

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
Application number
CN201910571260.7A
Other languages
Chinese (zh)
Other versions
CN110598520B (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.)
Watson Energy Technology (langfang) Co Ltd
Original Assignee
Watson Energy Technology (langfang) 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 Watson Energy Technology (langfang) Co Ltd filed Critical Watson Energy Technology (langfang) Co Ltd
Priority to CN201910571260.7A priority Critical patent/CN110598520B/en
Publication of CN110598520A publication Critical patent/CN110598520A/en
Application granted granted Critical
Publication of CN110598520B publication Critical patent/CN110598520B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P3/00Measuring linear or angular speed; Measuring differences of linear or angular speeds
    • G01P3/64Devices characterised by the determination of the time taken to traverse a fixed distance
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/02Preprocessing
    • G06F2218/04Denoising

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

Speed measuring method for preventing strong common mode noise interference in signal
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.
CN201910571260.7A 2019-06-28 2019-06-28 Speed measuring method for preventing strong common mode noise interference in signal Active CN110598520B (en)

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)

* Cited by examiner, † Cited by third party
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

Patent Citations (6)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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