CN105318940B - A kind of multi-way stream gauge calibrating installation pulse counting signal reconstructing method - Google Patents
A kind of multi-way stream gauge calibrating installation pulse counting signal reconstructing method Download PDFInfo
- Publication number
- CN105318940B CN105318940B CN201510649332.7A CN201510649332A CN105318940B CN 105318940 B CN105318940 B CN 105318940B CN 201510649332 A CN201510649332 A CN 201510649332A CN 105318940 B CN105318940 B CN 105318940B
- Authority
- CN
- China
- Prior art keywords
- signal
- amplitude
- rising edge
- failing edge
- calibrating installation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F25/00—Testing or calibration of apparatus for measuring volume, volume flow or liquid level or for metering by volume
Abstract
The present invention proposes a kind of multi-way stream gauge calibrating installation pulse counting signal reconstructing method, the method includes：Obtain flowmeter and proving flowmeter output signal u to be checkedo(t), u is analyzedo(t) amplitude versus frequency characte, first trough frequency is filter cutoff frequency before selecting amplitude maximum signal, using N rank high pass digital Chebyshev filters, obtains filtered signal v (t)；Calculate v (t) T average periodv, amplitude root mean square Av, setting rising edge, failing edge threshold coefficient ku、kd, then rising edge, failing edge activation threshold value are respectively kuAv、kdAv, it is one-dimension array P to record all rising edge trigger positionsu, all failing edge trigger positions be one-dimension array Pd；According to Pu、Pd, with AvFor amplitude, the time Δ t between adjacent rising edges (failing edge) should meet condition (1 3 σ) Tv＜ Δ t ＜ (1+3 σ) Tv, wherein σ is that pulse exports periodic deviation, Reconstruction Impulse signal u (t).
Description
Technical field
The present invention relates to flowmeter pulses signal denoising field more particularly to a kind of multi-way stream gauge calibrating installation pulsimeters
Number signal reconfiguring method.
Background technology
With development of modern industry, flowmeter demand is continuously increased, and especially pulse-output type flowmeter (is stablized in flow field
In the case of, pulse-output type flowmeter instantaneous flow is proportional to pulse frequency) it is unprecedentedly applied, calibrating quick to its, accurate
It is significant.General calibrating installation, industrial personal computer drive commutator according to the lockin signal of flowmeter pulses signal, can be to flow
The pulse signal of meter realizes interception complete cycle.However, in more position devices, proving flowmeter, flowmeter to be checked, converter,
The non-ideal device of inverse transformer generates signal there are nonlinear distortion, and interrupted current, lightning etc. cause impulse disturbances, power supply, nothing
Influence of the thermal noise and Johnson noise of single-frequency noise and intra-channel caused by line electricity etc. for system is more very important.
These factors may cause the case where step-by-step counting error count.
Invention content
To solve the problems, such as that above-mentioned and defect, the present invention propose a kind of multi-way stream gauge calibrating installation step-by-step counting letter
Number reconstructing method realizes the method combined with reconstruct using filtering, reconstructs flowmeter output pulse signal, removes output signal
In single-frequency noise, impulse disturbances, white noise, improve the signal-to-noise ratio of signal, prevent from causing because of noise jamming in flow measurement
The case where step-by-step counting error count.
The purpose of the present invention is realized by technical solution below：
A kind of multi-way stream gauge calibrating installation pulse counting signal reconstructing method, which is characterized in that the method includes：
A obtains flowmeter and proving flowmeter output signal u to be checkedo(t)；
B analyzes uo(t) amplitude versus frequency characte, first trough frequency is filter cutoff frequency before selecting amplitude maximum signal,
Using N rank high pass digital Chebyshev filters, filtered signal v (t) is obtained；
C calculates v (t) T average periodv, amplitude root mean square Av, setting rising edge, failing edge threshold coefficient ku、kd, then rise
Edge, failing edge activation threshold value are respectively kuAv、kdAv；
It is one-dimension array P that D, which records all rising edge trigger positions,u, all failing edge trigger positions be one-dimension array Pd；Root
According to Pu、Pd, with AvFor amplitude, the time Δ t between adjacent rising edges (failing edge) should meet condition (1-3 σ) Tv＜ Δ t ＜ (1+3
σ)Tv, wherein σ is that pulse exports periodic deviation Reconstruction Impulse signal u (t).
Present invention has the advantages that：
Realize the reconstruct of flowmeter output pulse signal, single-frequency noise, impulse disturbances in removal output signal, white noise
Sound improves the signal-to-noise ratio of signal, the case where preventing from flow measurement leading to step-by-step counting error count because of noise jamming.
Description of the drawings
Fig. 1 is multi-way stream gauge calibrating installation impulse compensation signal reconfiguring method flow diagram of the present invention；
Fig. 2 is original signal figure；
Fig. 3 is Chebyshev filter curve figures；
Fig. 4 is reconstruct afterpulse signal graph.
Specific implementation mode
Present invention will now be described in further detail with reference to the embodiments and the accompanying drawings.
The present invention is to remove single-frequency noise based on N rank high pass digital Chebyshev filters, real based on signal reconstruction technology
Existing multi-way stream gauge calibrating installation impulse compensation signal reconstruction, as shown in Figure 1, this method comprises the following steps：
Step 10 obtains flowmeter and proving flowmeter output signal u to be checkedo(t), as shown in Figure 2；
Step 20, analysis uo(t) amplitude versus frequency characte, first trough frequency is filter cutoff before selecting amplitude maximum signal
Frequency obtains filtered signal v (t) using N rank high pass digital Chebyshev filters；
The wherein Amplitude-squared function A (Ω of N ranks high pass digital Chebyshev filters2) be：
In formula, ΩcFor effective cut-off frequecy of passband of Chebyshev filters, ε is and passband ripple associated arguments, ε ∈
(0,1)；VN(x) it is N rank Chebyshev multinomials, is：
Ω is setcFor uo(t) first trough frequency before amplitude maximum signal in amplitude versus frequency characte；In the present embodiment,
Ωc=18.03Hz, ε=0.5, N=1, uo(t) as shown in Figure 3 through Chebyshev filtered signal v (t).
Chebyshev filters passband peak error it is minimum, in stopband amplitude response monotone decreasing, than
Butterworth filter attenuation slopes bigger, the rate of decay are faster.
Step 30 calculates v (t) T average periodv=33.4ms, amplitude root mean square Av=0.9754V, setting rising edge, under
Drop is along threshold coefficient ku=0.8, kd=-0.2, then rising edge, failing edge activation threshold value are respectively kuAv=0.7803V, kdAv=-
0.1951V；
Step 40 records all rising edge trigger positions as one-dimension array Pu, all failing edge trigger positions be a dimension
Group Pd；According to Pu、Pd, with AvFor amplitude, the time Δ t between adjacent rising edges (failing edge) should meet condition (1-3 σ) Tv＜ Δs t
＜ (1+3 σ) Tv, Reconstruction Impulse signal u (t).
1, all rising edge trigger positions are one-dimension array Pu=(pu1,pu2…pui), wherein puiIt is touched for i-th of rising edge
Send out position；All failing edge trigger positions are one-dimension array Pd=(pd1,pd2…pdn), wherein pdnFor n-th of failing edge trigger bit
It sets.
2, according to Pu、Pd, with AvFor amplitude, the time Δ t between adjacent rising edges (failing edge) should meet condition (1-3 σ) Tv
＜ Δ t ＜ (1+3 σ) Tv, wherein σ is that the method for pulse output periodic deviation Reconstruction Impulse signal u (t) is：
1. determining that first edge of pulse is rising edge or failing edge, and determine the position c on first edge1；
2. for k-th along ckFor rising edge, and ck=pui, then+1 edge of kth is failing edge, and corresponding failing edge meets item
Part：
if pdn＞ ck,pdn∈Pdk
ck+1=min (Pdk)
K ∈ N in formula+、n∈N+。
For j-th along cjFor failing edge, and cj=pdi, then+1 edge of jth is rising edge, and corresponding rising edge meets item
Part：
if pdi＞ cj,pdi∈Puj
cj+1=min (Puj)
3. obtaining all pulse transition position one-dimension array C=(c1,c2……cx), Reconstruction Impulse signal u (t), expression
Formula is：
It is as shown in Figure 4 to reconstruct afterpulse signal u (t).
Although disclosed herein embodiment it is as above.But the content is only to facilitate understanding the present invention and adopting
Embodiment is not limited to the present invention.Any those skilled in the art to which this invention pertains are not departing from this
Under the premise of the disclosed spirit and scope of invention, any modification and change can be made in the implementing form and in details,
But the scope of patent protection of the present invention, still should be subject to the scope of the claims as defined in the appended claims.
Claims (3)
1. a kind of multi-way stream gauge calibrating installation pulse counting signal reconstructing method, which is characterized in that the method includes：
A obtains flowmeter and proving flowmeter output signal u to be checkedo(t)；
B analyzes uo(t) amplitude versus frequency characte, first trough frequency is filter cutoff frequency before selecting amplitude maximum signal, using N
Rank high pass digital Chebyshev filters, obtain filtered signal v (t)；
C calculates v (t) T average periodv, amplitude root mean square Av, setting rising edge, failing edge threshold coefficient ku、kd, then rising edge, under
Drop is respectively k along activation threshold valueuAv、kdAv；
It is one-dimension array P that D, which records all rising edge trigger positions,u, all failing edge trigger positions be one-dimension array Pd；According to Pu、
Pd, with AvFor amplitude, the time Δ t between adjacent rising edges or failing edge should meet condition (1-3 σ) Tv＜ Δ t ＜ (1+3 σ) Tv,
Wherein σ is that pulse exports periodic deviation, Reconstruction Impulse signal u (t).
2. multi-way stream gauge calibrating installation pulse counting signal reconstructing method as described in claim 1, which is characterized in that described
The Amplitude-squared function A (Ω of N ranks high pass digital Chebyshev filters in step B2) be：
Wherein, ΩcFor effective cut-off frequecy of passband of Chebyshev filters, ε be with passband ripple associated arguments, ε ∈ (0,
1)；VN(x) it is N rank Chebyshev multinomials, is：
3. multi-way stream gauge calibrating installation pulse counting signal reconstructing method as described in claim 1, which is characterized in that described
In step D, all rising edge trigger positions are one-dimension array Pu=(pu1,pu2…pui), wherein puiIt is triggered for i-th of rising edge
Position；All failing edge trigger positions are one-dimension array Pd=(pd1,pd2…pdn), wherein pdnFor n-th of failing edge trigger bit
It sets.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510649332.7A CN105318940B (en) | 2015-10-08 | 2015-10-08 | A kind of multi-way stream gauge calibrating installation pulse counting signal reconstructing method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510649332.7A CN105318940B (en) | 2015-10-08 | 2015-10-08 | A kind of multi-way stream gauge calibrating installation pulse counting signal reconstructing method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105318940A CN105318940A (en) | 2016-02-10 |
| CN105318940B true CN105318940B (en) | 2018-07-13 |
Family
ID=55246827
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510649332.7A Active CN105318940B (en) | 2015-10-08 | 2015-10-08 | A kind of multi-way stream gauge calibrating installation pulse counting signal reconstructing method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105318940B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115792507B (en) * | 2023-02-09 | 2023-05-23 | 昆明理工大学 | Multiple lightning stroke distinguishing method and system based on monotonicity of short-time window slope |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100403653C (en) * | 2002-09-24 | 2008-07-16 | 联发科技股份有限公司 | Method and apparatus for calibrating acceptable deviation of maximum jitter |
| EP2861955A4 (en) * | 2012-06-15 | 2016-06-08 | Handyem Inc | Method and flow cell for characterizing particles by means of non-gaussian temporal signals |
| CN102901853A (en) * | 2012-09-19 | 2013-01-30 | 四川大学 | Nuclear signal frequency conversion digitization sampling method |
| CN104266732B (en) * | 2014-10-24 | 2017-06-30 | 广州能源检测研究院 | Dual stage bit stream gauge calibrating installation step-by-step counting precision compensation method |
| CN104639109B (en) * | 2015-02-12 | 2017-03-01 | 上海凌世电子有限公司 | Spike generator |
| CN104897250B (en) * | 2015-06-25 | 2018-02-02 | 华南理工大学 | A kind of more bit stream gauge step-by-step counting compensation methodes for resisting strong harmonic wave interference |
-
2015
- 2015-10-08 CN CN201510649332.7A patent/CN105318940B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105318940A (en) | 2016-02-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bossavy et al. | Forecasting ramps of wind power production with numerical weather prediction ensembles | |
| CN103185837B (en) | Method for measuring frequency of power system | |
| CN101539596A (en) | Method for monitoring electric network frequency | |
| DE3490412T (en) | Apparatus and method for determining the size of a capacitance | |
| CN109557573B (en) | Digital alpha/beta ray discrimination method | |
| DE102014211872A1 (en) | Control unit for capacitive touch sensor with automatic gain control | |
| CN103837884A (en) | Digital nucleus pulse signal trapezoidal shaping algorithm based on time-domain analysis | |
| CN105318940B (en) | A kind of multi-way stream gauge calibrating installation pulse counting signal reconstructing method | |
| CN205374588U (en) | Leakage current collection system | |
| CN110673532B (en) | Pulse signal processing and collecting device | |
| US20210111730A1 (en) | Systems and methods for removing low frequency offset components from a digital data stream | |
| JP2011202956A5 (en) | ||
| US20220103165A1 (en) | Filters for removing disturbances from signals | |
| CN111505375A (en) | Frequency measurement method and device | |
| CN107356799A (en) | A kind of electric current delay method of sampling | |
| US11817875B2 (en) | Systems and methods for removing low frequency offset components from a digital data stream | |
| US11758308B2 (en) | Systems and methods for improving frequency response of a high-speed data acquisition device | |
| CN103308765B (en) | Wide region, high resolution frequency monitor | |
| CN108872692B (en) | PWM wave duty ratio measuring method under industrial noise condition | |
| CN110212892B (en) | High-precision electric energy meter variable threshold value integral differential pulse generation method | |
| WO2013021175A3 (en) | Device and method for photon detection | |
| WO2011137255A1 (en) | Digital fractional integrator | |
| CN103580646B (en) | A kind of method for estimating analog filter frequency response characteristic | |
| CN102420664A (en) | Simulation system of noise signal for carrier test of low-voltage power line | |
| Kumar et al. | Design of Bandpass Finite Impulse Response Filter Using Various Window Method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |