CN104913830A - Radar level instrument signal processing method based on curve fitting - Google Patents
Radar level instrument signal processing method based on curve fitting Download PDFInfo
- Publication number
- CN104913830A CN104913830A CN201510252502.8A CN201510252502A CN104913830A CN 104913830 A CN104913830 A CN 104913830A CN 201510252502 A CN201510252502 A CN 201510252502A CN 104913830 A CN104913830 A CN 104913830A
- Authority
- CN
- China
- Prior art keywords
- rising edge
- echo
- time
- slope
- crest
- 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.)
- Pending
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention relates to a radar level instrument signal processing method based on curve fitting. The method includes the steps of calculating the slope of each point of an echo rising edge, judging whether the duration time of the rising edge meets the minimum keeping time or not, determining a certain rising edge as a target estimation position through the maximum slope, fitting a parabola through the least square method according to data of a 60% to 90% section of a crest of the rising edge at the target estimation position, calculating the peak value and corresponding time, and finally accurately positioning a level instrument. By means of the method, the purpose of positioning the multi-diameter multi-echo crest is effectively achieved, the influences of the measurement environment and a target object itself on an echo rear edge are eliminated, and the positioning accuracy is improved.
Description
Technical field
The present invention relates to field of industrial measurement, specifically a kind of radar levelmeter signal processing method based on curve.
Background technology
Thing position is industrial important parameter.The level measurement instruments such as existing static pressure method, float-ball type, ultrasound wave, have respective features and application scope.Radar levelmeter, by advanced treatment technology, has unique advantage, demonstrates the performance of its brilliance especially, play more and more important effect in the industrial production under the harsh measurement conditions such as high temperature, high pressure, steam, corrosive medium are strong.
Industry thing position, liquid level gauge operating mode complexity, cause the problem of measuring accuracy under the measuring condition as chaff interference, guided wave bar coating, foam, steam and turbulent flow.Typical algorithm is detection peak, and the prerequisite of algorithm application is that echo waveform is enough good.In actual application, waveform multipath effect often has multiple echo, and such as, wall echo causes the peak value of more than, is exactly in addition badly to determine peak value.
Existing placement technology is by carrying out linear interpolation realization to echo, along the impact being subject to measurement environment and destination object itself after peak value and echo, linear interpolation is difficult to compound actual waveform characteristic, and this makes to be difficult to determine correct echo time of arrival, and precision improvement is limited in one's ability.
Summary of the invention
For the deficiencies in the prior art, curve is fused in radar levelmeter location algorithm by the present invention, proposes a kind of radar levelmeter signal processing method based on curve.
The technical scheme that the present invention is adopted for achieving the above object is: a kind of radar levelmeter signal processing method based on curve, comprises the following steps:
Step 1: using the reflection echo under aimless original state as Noise Background N;
Step 2: compared with Noise Background N by the reflection echo X of actual measurement and disappear mutually, disappeared echo S mutually;
Step 3: to the smoothing filtering of echo S that disappears mutually, with the small amplitude fluctuations of filtering, obtains filtered echo S
w, calculate echo S
win the slope k often put of each rising edge, determine rising edge trough S
wtnto crest S
wcninitial time t
snwith end time t
en, i.e. the duration T of rising edge
hn=t
en-t
snif, rising edge duration T
hnbe less than minimum hold time T
hmin, then this rising edge is invalid;
Step 4: in the effective situation of rising edge, using the rising edge at the maximum some place of slope as target state estimator position;
Step 5: according to the crest S of target state estimator position
wcnthe data in 60% to 90% interval, utilize least square fitting para-curve;
Step 6: obtain peak value by the Fitting Calculation, and the time t of correspondence;
Step 7: R=(t-t accurately located by level meter
c) c/2, wherein, R be the thing distance of positions from, c is the light velocity, t
cfor the clock initial correction time.
The present invention has following beneficial effect and advantage:
1. the present invention is directed to multipath effect in actual application and produce multiple echo, by the feature determination target echo of rising edge;
2. along the impact being subject to measurement environment and destination object itself after the present invention utilizes least square fitting para-curve to solve peak value and echo;
3. the present invention utilizes least square fitting para-curve to promote the positioning precision of level meter.
Accompanying drawing explanation
Fig. 1 is method flow diagram of the present invention;
Fig. 2 (a) shows radar levelmeter echo situation 1;
Fig. 2 (b) shows radar levelmeter echo situation 2;
Fig. 2 (c) shows radar levelmeter echo situation 3;
Fig. 3 is curve schematic diagram of the present invention.
Embodiment
Below in conjunction with drawings and Examples, the present invention is described in further detail.
As shown in Figure 1, be schematic flow sheet of the present invention.First calculate the slope that echo rising edge is often put, and judge whether the rising edge duration meets minimum hold time; Then determine that certain rising edge is as target state estimator position by maximum slope; Again by the data separate least square fitting para-curve in 60% to 90% interval of the rising edge crest of target state estimator position, calculate the time of peak value and correspondence, finally obtain level meter and accurately locate.Concrete steps are as follows:
Step 1: using the reflection echo under aimless original state as Noise Background N;
Step 2: compared with Noise Background N by the reflection echo X of actual measurement and disappear mutually, disappeared echo S mutually;
Step 3: to the smoothing filtering of echo S that disappears mutually, the small amplitude fluctuations of filtering, obtains filtered echo S
w, calculate echo S
win the slope k often put of the n-th rising edge, determine rising edge trough S
wtnto crest S
wcninitial time t
snwith end time t
en, i.e. the duration T of rising edge
hn=t
en-t
snif, rising edge duration T
hnbe less than minimum hold time T
hmin, by invalid for this rising edge;
Step 4: maximum slope k
maxbe the slope of the n-th rising edge point, then using the n-th rising edge as target state estimator position;
As shown in Fig. 2 (a), there are 3 rising edges in radar levelmeter echo, the duration of 3 rising edges is all greater than minimum hold time, maximum slope k
maxbe the slope of the 2nd rising edge point, then using the 2nd rising edge as target state estimator position;
As shown in Fig. 2 (b), there are 3 rising edges in radar levelmeter echo, and the duration of the 3rd rising edge is all less than minimum hold time, then the 3rd rising edge is invalid, maximum slope k
maxbe the slope of the 2nd rising edge point, then using the 2nd rising edge as target state estimator position;
As shown in Fig. 2 (c), there are 3 rising edges in radar levelmeter echo, and the duration of the 3rd rising edge is all less than minimum hold time, then the 3rd rising edge is invalid, and the slope of the 3rd rising edge point is maximum slope k
max, because the 3rd rising edge is invalid, after excluding invalid rising edge, maximum slope k
maxbe the slope of the 2nd rising edge point, then using the 2nd rising edge as target state estimator position;
Step 5: the n-th rising edge crest S
wcnthe data separate least square fitting para-curve in 60% to 90% interval;
Step 6: obtain peak value by the Fitting Calculation, and the time t of correspondence;
As shown in Figure 3, radar n-th rising edge point A is rising edge crest S with the value of some B
wcn60% and 90%, by the data separate least square fitting para-curve in this interval, obtain dotted line, dotted line fixed point C is the peak value that curve calculates, and the corresponding time is t;
Step 7: R=(t-t accurately located by level meter
c) c/2, wherein, R be the thing distance of positions from, c is the light velocity, t
cfor the clock initial correction time.
Claims (1)
1., based on a radar levelmeter signal processing method for curve, it is characterized in that: comprise the following steps:
Step 1: using the reflection echo under aimless original state as Noise Background N;
Step 2: compared with Noise Background N by the reflection echo X of actual measurement and disappear mutually, disappeared echo S mutually;
Step 3: to the smoothing filtering of echo S that disappears mutually, with the small amplitude fluctuations of filtering, obtains filtered echo S
w, calculate echo S
win the slope k often put of each rising edge, determine rising edge trough S
wtnto crest S
wcninitial time t
snwith end time t
en, i.e. the duration T of rising edge
hn=t
en-t
snif, rising edge duration T
hnbe less than minimum hold time T
hmin, then this rising edge is invalid;
Step 4: in the effective situation of rising edge, using the rising edge at the maximum some place of slope as target state estimator position;
Step 5: according to the crest S of target state estimator position
wcnthe data in 60% to 90% interval, utilize least square fitting para-curve;
Step 6: obtain peak value by the Fitting Calculation, and the time t of correspondence;
Step 7: R=(t-t accurately located by level meter
c) c/2, wherein, R be the thing distance of positions from, c is the light velocity, t
cfor the clock initial correction time.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510252502.8A CN104913830A (en) | 2015-05-15 | 2015-05-15 | Radar level instrument signal processing method based on curve fitting |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510252502.8A CN104913830A (en) | 2015-05-15 | 2015-05-15 | Radar level instrument signal processing method based on curve fitting |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104913830A true CN104913830A (en) | 2015-09-16 |
Family
ID=54083091
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510252502.8A Pending CN104913830A (en) | 2015-05-15 | 2015-05-15 | Radar level instrument signal processing method based on curve fitting |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104913830A (en) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106443646A (en) * | 2016-11-21 | 2017-02-22 | 重庆兆洲科技发展有限公司 | Ultrasonic ranging system and echo processing method and device |
CN107505027A (en) * | 2016-06-14 | 2017-12-22 | Vega格里沙贝两合公司 | Microwave reflection barrier |
CN108426520A (en) * | 2018-05-15 | 2018-08-21 | 天津大学 | The detection device in pipeline stolen hole and the detection method at stolen hole center, diameter |
CN110596662A (en) * | 2019-10-21 | 2019-12-20 | 富临精工先进传感器科技(成都)有限责任公司 | Distance deviation correction method of MIMO radar |
CN111812606A (en) * | 2020-06-03 | 2020-10-23 | 西安电子科技大学 | Guided wave radar-based level extraction method |
CN113296066A (en) * | 2021-04-27 | 2021-08-24 | 北京锐达仪表有限公司 | Echo correction method, device, readable storage medium, electronic device and system |
CN113639816A (en) * | 2021-07-06 | 2021-11-12 | 燕山大学 | Non-invasive method for monitoring the level of a liquid in a container and system for implementing said method |
CN114637021A (en) * | 2022-05-18 | 2022-06-17 | 四川吉埃智能科技有限公司 | Sub-centimeter-level full-waveform laser radar ranging method and device |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5969666A (en) * | 1997-06-06 | 1999-10-19 | Endress + Hauser GmbH +Co. | Radar-based method of measuring the level of a material in a container |
US20070201309A1 (en) * | 2006-02-21 | 2007-08-30 | Robert Brown | Method for processing an echo profile and a pulse-echo ranging system |
CN102661773A (en) * | 2012-05-15 | 2012-09-12 | 合肥工业大学 | System and method for processing echo signal of guided wave type radar level meter |
CN203177928U (en) * | 2012-09-21 | 2013-09-04 | 清华大学 | Echo sampling unit and pulse type radar material level meter with the same |
CN103542911A (en) * | 2013-11-01 | 2014-01-29 | 合肥工业大学 | Guided wave radar level meter echo signal processing system and method based on first-order derivative |
-
2015
- 2015-05-15 CN CN201510252502.8A patent/CN104913830A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5969666A (en) * | 1997-06-06 | 1999-10-19 | Endress + Hauser GmbH +Co. | Radar-based method of measuring the level of a material in a container |
US20070201309A1 (en) * | 2006-02-21 | 2007-08-30 | Robert Brown | Method for processing an echo profile and a pulse-echo ranging system |
CN102661773A (en) * | 2012-05-15 | 2012-09-12 | 合肥工业大学 | System and method for processing echo signal of guided wave type radar level meter |
CN203177928U (en) * | 2012-09-21 | 2013-09-04 | 清华大学 | Echo sampling unit and pulse type radar material level meter with the same |
CN103542911A (en) * | 2013-11-01 | 2014-01-29 | 合肥工业大学 | Guided wave radar level meter echo signal processing system and method based on first-order derivative |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107505027A (en) * | 2016-06-14 | 2017-12-22 | Vega格里沙贝两合公司 | Microwave reflection barrier |
US10704948B2 (en) | 2016-06-14 | 2020-07-07 | Vega Grieshaber Kg | Reflection microwave barrier |
CN106443646A (en) * | 2016-11-21 | 2017-02-22 | 重庆兆洲科技发展有限公司 | Ultrasonic ranging system and echo processing method and device |
CN106443646B (en) * | 2016-11-21 | 2019-03-22 | 重庆兆洲科技发展有限公司 | A kind of ultrasonic ranging system, echo processing techniques and device |
CN108426520A (en) * | 2018-05-15 | 2018-08-21 | 天津大学 | The detection device in pipeline stolen hole and the detection method at stolen hole center, diameter |
CN110596662A (en) * | 2019-10-21 | 2019-12-20 | 富临精工先进传感器科技(成都)有限责任公司 | Distance deviation correction method of MIMO radar |
CN111812606A (en) * | 2020-06-03 | 2020-10-23 | 西安电子科技大学 | Guided wave radar-based level extraction method |
CN111812606B (en) * | 2020-06-03 | 2023-12-22 | 西安电子科技大学 | Material level extraction method based on guided wave radar |
CN113296066A (en) * | 2021-04-27 | 2021-08-24 | 北京锐达仪表有限公司 | Echo correction method, device, readable storage medium, electronic device and system |
CN113296066B (en) * | 2021-04-27 | 2023-06-23 | 北京锐达仪表有限公司 | Echo correction method, device, readable storage medium, electronic equipment and system |
CN113639816A (en) * | 2021-07-06 | 2021-11-12 | 燕山大学 | Non-invasive method for monitoring the level of a liquid in a container and system for implementing said method |
CN114637021A (en) * | 2022-05-18 | 2022-06-17 | 四川吉埃智能科技有限公司 | Sub-centimeter-level full-waveform laser radar ranging method and device |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104913830A (en) | Radar level instrument signal processing method based on curve fitting | |
CN106643939B (en) | The method for calculating ultrasonic propagation time for ultrasonic flowmeter | |
CN100565241C (en) | Ultrasonic echo frontier inspection based on modulation domain measurement | |
CN107860430B (en) | Time difference measurement method of ultrasonic gas flowmeter based on time difference method | |
CN104501889B (en) | Ultrasonic flow detection method based on cross-correlation time difference method | |
CN102297712A (en) | Method for measuring propagation time of ultrasonic echo | |
CN110297250B (en) | Initial grazing angle solving method based on Taylor expansion, sound ray bending correction method and equipment | |
CN106706056B (en) | A kind of compensation method of heavy caliber ultrasonic water meter flow measurement | |
RU2012133315A (en) | VARIABILITY TRACKING | |
CN105180853A (en) | Electromagnetic ultrasonic metal material thickness measurement method | |
CN104515560A (en) | Ultrasonic flow meter and dirt detection and error compensation method thereof | |
CN104748704A (en) | Thin-walled structure ultrasonic resonance thickness measurement frequency spectrum analysis interpolation correction method | |
US11346694B2 (en) | Method for measuring the velocity of a fluid using ultrasound | |
KR101764870B1 (en) | Signal processing system for ultrasonic floemeter | |
CN109470320B (en) | Step filtering method for ultrasonic water meter | |
CN112304376A (en) | Ultrasonic flowmeter flow measuring method based on data fusion | |
CN105785123A (en) | Radar signal frequency calculation method based on apFFT phase difference | |
KR101386641B1 (en) | Distance estimation device and method using the difference of wave speed between waves | |
CN114355360B (en) | Echo frequency estimation method of Doppler log | |
CN104502627B (en) | Method is solved based on transmitting Design of Signal and the phase ambiguity of processing in ADCP | |
CN110346783B (en) | Sound line correction method based on power series equivalent ranging model | |
JP4392463B2 (en) | Film thickness measurement method | |
JP4904099B2 (en) | Pulse signal propagation time measurement device and ultrasonic flow measurement device | |
JP7174574B2 (en) | Ultrasonic flowmeter and determination method of zero-crossing time in ultrasonic flowmeter | |
CN109031260B (en) | LFM signal time delay measurement method based on fractional Fourier modulation rate analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150916 |