CN106771546A - A kind of fixed frequency gathers measuring method with the quick amplitude of waveform signal - Google Patents
A kind of fixed frequency gathers measuring method with the quick amplitude of waveform signal Download PDFInfo
- Publication number
- CN106771546A CN106771546A CN201710062397.0A CN201710062397A CN106771546A CN 106771546 A CN106771546 A CN 106771546A CN 201710062397 A CN201710062397 A CN 201710062397A CN 106771546 A CN106771546 A CN 106771546A
- Authority
- CN
- China
- Prior art keywords
- sampling
- amplitude
- denoted
- matching
- sampled
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
- G01R19/25—Arrangements for measuring currents or voltages or for indicating presence or sign thereof using digital measurement techniques
- G01R19/2506—Arrangements for conditioning or analysing measured signals, e.g. for indicating peak values ; Details concerning sampling, digitizing or waveform capturing
- G01R19/2509—Details concerning sampling, digitizing or waveform capturing
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
A kind of quick amplitude the invention discloses fixed frequency and waveform signal gathers measuring method, comprises the following steps:Step A:If N number of 1~N of sampled point, the corresponding sampling time is denoted as t1, t2, t3..., tN, step B:Each sampling point sampling 5 times, the corresponding sampling time is denoted as ti1=ti2 Δ t, ti2=tiΔ t, ti3=ti, ti4=ti+ Δ t, ti5=t+2 Δ t, (i=1,2,3 ... N) step C:Minimum value and maximum in 5 sampled datas of removal, remaining three data are averaged and are denoted as Yi, corresponding to ti, the sampling time that really intermediate time of 5 samplings is put as this;N times are repeated, N group sampled datas (Y is obtainedi, ti);Step D:According to matching fitting algorithm principle to obtaining N group sampled datas (Yi, ti) matching the Fitting Calculation is carried out, so as to obtain real curve expression formula, and then extract amplitude.An order of magnitude is improve the invention enables sampling rate, patrol UAV is fully met to circuit or equipment tracking and the winged control requirement of obstacle.
Description
Technical field
The present invention relates to a kind of waveform signal amplitude acquisition technique field, more particularly to a kind of fixed frequency and waveform signal
Quick amplitude collection measuring method.
Background technology
In industrial control system, or instrument and equipment manufacture field, it is often necessary to known fixed frequency and the letter of waveform
Number carry out amplitude measurement, collection.Usual way is continuous quickly amplitude measurement, collection, and is continued more than a cycle
Time, according to certainty of measurement requirement, needed comprising enough sampling numbers in an entire cycle time.The one of this method
Individual major defect is, completes amplitude measurement, collection, it is necessary to be more than or equal to the time of a cycle, strongly limit signal
Amplitude measurement and the speed for gathering.For inspection and sampling that relatively large data point is required in the short time, this is particularly problematic,
For example, for 50 hertz of measurements of alternating magnetic field, using 10 × 10 inductance unit matrix inductive sensory appts, completing a pointwise
Scanography, sampling, it is necessary to be detected, sampled, every unit is needed according to the method described above one by one to 100 sensing unit devices
No less than a cycle time 20ms, 100 units add up to be needed no less than 100 × 5=5000ms times, i.e., once completely to sweep
Retouch, acquisition time is not less than 5s, if element number is more or the signal period is longer, the real-time that can have a strong impact on system refers to
Mark, or even cause that engineering demand can not be met.
The content of the invention
A kind of quick amplitude it is an object of the invention to provide fixed frequency and waveform signal gathers measuring method, can not
Measured signal cycle length is depended on, quick detection and the collection of amplitude are carried out for known fixed frequency and waveform signal.
The present invention uses following technical proposals:
A kind of fixed frequency and the quick amplitude collection measuring method of waveform signal, comprise the following steps:
Step A:If N number of 1~N of sampled point, the corresponding sampling time is denoted as t1, t2, t3..., tN,
Step B:Each sampling point sampling 5 times, the corresponding sampling time is denoted as ti1=ti- 2 Δ t, ti2=ti- Δ t, ti3=
ti, ti4=ti+ Δ t, ti5=t+2 Δ t, (i=1,2,3 ... N)
Step C:Minimum value and maximum in 5 sampled datas of removal, remaining three data are averaged and are denoted as
Yi, corresponding to ti, the sampling time that really intermediate time of 5 samplings is put as this;N times are repeated, N groups is obtained and is adopted
Sample data (Yi, ti);
Step D:According to matching fitting algorithm principle to obtaining N group sampled datas (Yi, ti) matching the Fitting Calculation is carried out, from
And real curve expression formula is obtained, and then extract amplitude.
Matching fitting algorithm described in the step D specifically includes following steps:
If ideal signal curve:
The value of N number of sampled point is expressed as:
Ideal curve is designated as with the deviation of each sampled point:
According to least square fitting principle so thatIt is minimum value, even if alsoTo AmWithSingle order local derviation
Number is equal to zero, and obtaining two equations is:
Wherein ω, tiAnd YiIt is known quantity, solves above-mentioned equation, you can obtain AmWithObtain amplitude.
The present invention carries out multi-point sampling to measured signal according to acquisition precision needs, and as N is sampled for each point
5 times, due to the frequency and waveform of measured signal be, it is known that so, using the N groups data (correspondence N number of sampled point) obtained by sampling
Curve Matching the Fitting Calculation is carried out, so as to obtain the amplitude information of measured signal.In actual applications, high-speed AD converter is adopted
Sample speed is 10MHZ, i.e. the cycle is 1/10 microsecond, for 50HZSinusoidal signal is carried out continuously 10 point samplings (according to digital filtering reason
By, sampling number is determined by certainty of measurement, takes 10~15), each point sampling 5 times (really in the surrounding sample of the sampled point,
The number of times can not be too big), then take about (× 5 times × 1/10)=5 microseconds at 10 points, it is contemplated that Curve Matching the Fitting Calculation is spent
Time, add 150 continuous sampling times, 50HZSinusoidal signal amplitude detection amounts to time-consuming about 10 microseconds with sampling, will pass
System method needs time-consuming about 20 milliseconds of a cycle, has been compressed to about ten delicate so that sampling rate improves an order of magnitude,
Patrol UAV is fully met to circuit or equipment tracking and the winged control requirement of obstacle.
Brief description of the drawings
Fig. 1 is flow chart of the invention.
Specific embodiment
Preferred embodiments of the invention will be below described in detail;It should be appreciated that preferred embodiments are only for saying
The bright present invention, rather than in order to limit the scope of the invention.As shown in figure 1, the present invention comprises the following steps:
Step A:If N number of 1~N of sampled point, the corresponding sampling time is denoted as t1, t2, t3..., tN,
Step B:Each sampling point sampling 5 times, the corresponding sampling time is denoted as ti1=ti- 2 Δ t, ti2=ti- Δ t, ti3=
ti, ti4=ti+ Δ t, ti5=t+2 Δ t, (i=1,2,3 ... N) are really using 5 intermediate times of sampling as adopting that this is put
The sample time
Step C:Minimum value and maximum in 5 sampled datas of removal, remaining three data are averaged and are denoted as
Yi, corresponding to ti, obtain N group sampled datas;
Step D:Matching the Fitting Calculation is carried out according to matching fitting algorithm principle, so as to obtain real curve expression formula, is carried
Take amplitude.
Described matching fitting algorithm specifically includes following steps:
If ideal signal curve:
The value of N number of sampled point is expressed as:
Ideal curve is designated as with the deviation of each sampled point:
According to least square fitting principle so thatIt is minimum value, even if alsoTo AmWithSingle order local derviation
Number is equal to zero, and obtaining two equations is:
Wherein ω, tiAnd YiIt is known quantity, solves above-mentioned equation, you can obtain AmWithObtain amplitude.
In view of the time that Curve Matching the Fitting Calculation is spent, 150 continuous sampling times, 50HZ sinusoidal signals are added
Amplitude detection amounts to time-consuming about 10 microseconds with sampling, and conventional method is needed into time-consuming about 20 milliseconds of a cycle, has been compressed to about
Ten is delicate so that sampling rate improves an order of magnitude.
Processed by passive type array magnetic induction antenna assembly, array scanning and acquisition control circuit, sampling and attitude data
Device is constituted, and wherein core is the high-speed A/D conversion circuit in array scanning and acquisition control circuit, and is sampled and appearance
Sampled-data processing algorithm in state data processor.
The high-speed AD converter part for using sampling rate to reach MHZ grades in array scanning and acquisition control circuit, to tested
Signal carries out multiple spot high-speed sampling, and sample information is transferred into sampling and attitude data processor below.
Sampling is with attitude data processor according to specially designed high speed data processing algorithm, including digital filtering, curve
Matching the Fitting Calculation, magnitude extraction etc., calculate the amplitude for being sampled signal.Digital filtering is mutually auxiliary with Curve Matching the Fitting Calculation
Coordinate, it is indispensable, lack digital filtering, magnitude extraction precision does not ensure that lack Curve Matching the Fitting Calculation, amplitude cannot
Extract.
The method is suitable for the quick detection of the signal amplitude to given frequency and waveform, such as sine wave, triangle with device
Ripple, square wave etc., are particularly suitable for the quick detection application of low frequency signal.
If high-speed AD converter sampling rate is 10MHZ, i.e. the cycle is 1/10 microsecond, for 50HZSinusoidal signal is continuous
(according to Digit filter theory, sampling number is determined by certainty of measurement, takes 10~15), each point sampling 5 times to carry out 10 point samplings
(really in the surrounding sample of the sampled point, the number of times can not be too big), then take about (× 5 times × 1/10)=5 microseconds at 10 points,
Due to the frequency (50H of measured signalZ) and waveform (sine) be, it is known that so, it is possible to use 10 groups of data obtained by sampling are (right
Answer 10 sampled points) Curve Matching the Fitting Calculation is carried out, so as to obtain the amplitude information of measured signal.Intend in view of Curve Matching
It is total to calculate the spent time, add 150 continuous sampling times, 50HZSinusoidal signal amplitude detection takes about altogether with sampling
10 microseconds, time-consuming about 20 milliseconds of a cycle is needed by conventional method, has been compressed to about ten delicate so that sampling rate improves
An order of magnitude, fully meets patrol UAV to circuit or equipment tracking and the winged control requirement of obstacle.
Claims (2)
1. a kind of fixed frequency gathers measuring method with the quick amplitude of waveform signal, it is characterised in that:Comprise the following steps:
Step A:If N number of 1~N of sampled point, the corresponding sampling time is denoted as t1, t2, t3..., tN,
Step B:Each sampling point sampling 5 times, the corresponding sampling time is denoted as ti1=ti- 2 Δ t, ti2=ti- Δ t, ti3=ti,
ti4=ti+ Δ t, ti5=t+2 Δ t, (i=1,2,3 ... N)
Step C:Minimum value and maximum in 5 sampled datas of removal, remaining three data are averaged and are denoted as Yi, it is right
Should be in ti, the sampling time that really intermediate time of 5 samplings is put as this;N times are repeated, N group hits is obtained
According to (Yi, ti);
Step D:According to matching fitting algorithm principle to obtaining N group sampled datas (Yi, ti) matching the Fitting Calculation is carried out, so as to obtain
Real curve expression formula is obtained, and then extracts amplitude.
2. fixed frequency according to claim 1 gathers measuring method with the quick amplitude of waveform signal, it is characterised in that:
Matching fitting algorithm described in the step D specifically includes following steps:
If ideal signal curve:
The value of N number of sampled point is expressed as:
Ideal curve is designated as with the deviation of each sampled point:
According to least square fitting principle so thatIt is minimum value, even if alsoTo AmWithFirst-order partial derivative be equal to
Zero, obtaining two equations is:
Wherein ω, tiAnd YiIt is known quantity, solves above-mentioned equation, you can obtain AmWithObtain amplitude.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710062397.0A CN106771546A (en) | 2017-01-23 | 2017-01-23 | A kind of fixed frequency gathers measuring method with the quick amplitude of waveform signal |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710062397.0A CN106771546A (en) | 2017-01-23 | 2017-01-23 | A kind of fixed frequency gathers measuring method with the quick amplitude of waveform signal |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106771546A true CN106771546A (en) | 2017-05-31 |
Family
ID=58955460
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710062397.0A Pending CN106771546A (en) | 2017-01-23 | 2017-01-23 | A kind of fixed frequency gathers measuring method with the quick amplitude of waveform signal |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106771546A (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107328987A (en) * | 2017-06-27 | 2017-11-07 | 国网福建省电力有限公司 | A kind of processing method and processing device of opto-electronic current transducer wave distortion |
CN110045180A (en) * | 2019-05-14 | 2019-07-23 | 中南民族大学 | A kind of method and system being most worth measurement for waveform |
CN111130504A (en) * | 2019-12-23 | 2020-05-08 | 成都理工大学 | Data acquisition method capable of self-adjusting sampling frequency |
CN112461376A (en) * | 2020-12-07 | 2021-03-09 | 爱沃泰科技(深圳)有限公司 | Self-adaptive infrared sensor signal sampling algorithm |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04230869A (en) * | 1990-05-03 | 1992-08-19 | Landis & Gyr Betrieps Ag | Method and apparatus for obtaining estimated value of instantaneous value of parameter of sine-shaped signal |
CN2403031Y (en) * | 1999-12-14 | 2000-10-25 | 天津德威电子有限公司 | Sine wave high-speed vector measuring device |
-
2017
- 2017-01-23 CN CN201710062397.0A patent/CN106771546A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04230869A (en) * | 1990-05-03 | 1992-08-19 | Landis & Gyr Betrieps Ag | Method and apparatus for obtaining estimated value of instantaneous value of parameter of sine-shaped signal |
CN2403031Y (en) * | 1999-12-14 | 2000-10-25 | 天津德威电子有限公司 | Sine wave high-speed vector measuring device |
Non-Patent Citations (2)
Title |
---|
佛朗哥·马洛贝蒂: "《数据转换器》", 31 July 2013, 西安交通大学出版社 * |
童敏明等: "《传感器原理与检测技术》", 31 January 2014, 机械工业出版社 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107328987A (en) * | 2017-06-27 | 2017-11-07 | 国网福建省电力有限公司 | A kind of processing method and processing device of opto-electronic current transducer wave distortion |
CN110045180A (en) * | 2019-05-14 | 2019-07-23 | 中南民族大学 | A kind of method and system being most worth measurement for waveform |
CN110045180B (en) * | 2019-05-14 | 2019-11-22 | 中南民族大学 | A kind of method and system being most worth measurement for waveform |
CN111130504A (en) * | 2019-12-23 | 2020-05-08 | 成都理工大学 | Data acquisition method capable of self-adjusting sampling frequency |
CN111130504B (en) * | 2019-12-23 | 2023-03-03 | 成都理工大学 | Data acquisition method capable of self-adjusting sampling frequency |
CN112461376A (en) * | 2020-12-07 | 2021-03-09 | 爱沃泰科技(深圳)有限公司 | Self-adaptive infrared sensor signal sampling algorithm |
CN112461376B (en) * | 2020-12-07 | 2024-05-24 | 爱沃泰科技(深圳)有限公司 | Adaptive infrared sensor signal sampling algorithm |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106771546A (en) | A kind of fixed frequency gathers measuring method with the quick amplitude of waveform signal | |
CN104545887B (en) | The recognition methods of artifact ecg wave form and device | |
CN102680860A (en) | Automatic fault-point locating method for traveling-wave based fault location of high-voltage electric power lines | |
CN102664588A (en) | Digital demodulation method and device for frequency modulation and amplitude modulation signals | |
CN203054093U (en) | Apparatus for extracting DC (direct current) signal from mixed signal by using synchronous sampling | |
CN104391177B (en) | EMUs side harmonics test system and method | |
CN105866558B (en) | A kind of zonule lightning monitoring localization method and system | |
CN106771582A (en) | The method of testing and tester of high-frequency pulse signal | |
CN112525201A (en) | Underwater target tracking method based on electromagnetic field characteristic multi-information fusion | |
CN105737793B (en) | Rolling angle measurement unit and measuring method | |
CN105403778A (en) | Measuring method and system of multi-loop sequential work | |
CN204594846U (en) | A kind of mist droplet deposition sensor | |
CN105929284A (en) | Wireless tester and testing method for live-line measuring zinc oxide arrester | |
CN110161376B (en) | Traveling wave fault time extraction algorithm | |
CN103487786A (en) | Structural impact wave-velocity-free positioning method based on two-dimensional linear arrays and spatial filters | |
CN103674177B (en) | A kind of signal processing method and device | |
CN102401856B (en) | Real-time frequency measuring method for power frequency voltage signal based on fuzzy clustering method | |
CN105004363B (en) | Multi thresholds small wave sensor performance online test device and method are based under strong jamming | |
CN106093181A (en) | Parameter of magnetic characteristic, magnetic Barkhausen, magnetosonic launch detection method and device | |
CN106610293B (en) | A kind of indoor orientation method and system based on intensity difference | |
Liu et al. | A pointer meter reading recognition method based on improved ORB algorithm for substation inspection robot | |
CN112130003A (en) | Device and method for removing electromagnetic interference signals in same frequency band | |
CN106443152A (en) | High-precision rural network low-voltage transient residual current detection method | |
CN106782711A (en) | Signal extraction and error filtering method and device | |
CN105372710B (en) | A kind of conducting wire detection 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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170531 |
|
RJ01 | Rejection of invention patent application after publication |