CN104614586A - GPS-based self-adaptive sampling method for electric quantity signals - Google Patents
GPS-based self-adaptive sampling method for electric quantity signals Download PDFInfo
- Publication number
- CN104614586A CN104614586A CN201510016798.3A CN201510016798A CN104614586A CN 104614586 A CN104614586 A CN 104614586A CN 201510016798 A CN201510016798 A CN 201510016798A CN 104614586 A CN104614586 A CN 104614586A
- Authority
- CN
- China
- Prior art keywords
- sampling
- designated
- int
- interval
- pps
- 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
Abstract
The invention discloses a GPS-based self-adaptive sampling method for electric quantity signals. When a needed sampling rate is close to a rated frequency of a chip, the dispersion of each sampling interval width and the dispersion increase degree of each sampling section width are effectively inhibited, and therefore, when improving the sampling rate under the situation that the rated frequency of the chip cannot be improved, the electric quantity gathering precision descending degree reduction became possible. The GPS-based self-adaptive sampling method for electric quantity signals is capable of quickly forecasting the clock frequency of the next second according to the previous clock frequency change trend, the anti-interference performance is strong, the great sampling frequency change due to unstable input pulse per second is avoided, and especially, the regulation difference is allowed to surpass the width of one sampling point, and the application range is expanded.
Description
Technical field
The invention belongs to electric system observation and control technology field, specifically just relate to a kind of electric parameters signal adaptive method of sampling based on GPS.
Background technology
Along with the further expansion of electric system scale, various types of power supply and load are linked in system, make the control of electric system become more complicated.In order to ensure the safe and stable operation of system, need the frequency to current/voltage, amplitude, harmonic wave, phase angle, the electric parameters such as electric flux carry out high-precision measurement.In order to obtain high-precision electrical quantities measurement value, not only need when AC sampling the number controlling interior sampled point per second to ensure the precision of sampling rate, and need to improve each sampled point interval degree of uniformity.In electric parameters collection and computation process, p.s. sampled data needs order and distribute to comparably number wheel AC sampling calculate, such as, need p.s. to carry out the calculating that 50 take turns electric voltage frequency, sampling rate is 3200Hz, then 3200 point data of p.s. need sequentially uniform distribution 50 to take turns frequency computation part to this, and each takes turns 64 total sampling widths of sampled point needed for calculating in other words, and the width being called for short sampling section here also needs identical as far as possible.Therefore in order to improve the precision of the result of calculation that each is taken turns, not only need reduction 3200 sampling interval dispersions each other, and need the dispersion each other of reduction 50 sampling section width, otherwise the measured value that the sampling section that width will be made narrower calculates the a-c cycle obtained is bigger than normal, and the measured value that the sampling section of wider width calculates the a-c cycle obtained is less than normal.
Chinese patent 201410173291.4 discloses a kind of high precision pulse per second (PPS) based on FPGA and doubly to occur frequently the method for sampling pulse, is below designated as method 1.Actual synchronization pulse per second (PPS) count value and ideal synchronisation pulse per second (PPS) count value compare by the method, deviation and sample frequency are carried out accumulation calculating, judge whether impulse meter corrects according to accumulation calculating result, thus reach the total departure of actual synchronization pulse per second (PPS) is evenly distributed in actual samples pulse, reduce the error of sampling pulse, for signals collecting provides more reliable and stable sampling pulse signal.But there is the defect of several aspect in method 1, the first, the deterministic process that the judgement of the accumulation amount in accumulation calculating process and actual samples pulse export can introduce time delay; Second, this method of sampling needs the prerequisite of a hypothesis: the frequency change of crystal oscillator is very little, but in fact the frequency of crystal oscillator can along with environment temperature and device inside state consecutive variations, and the sample effect that the effect of method 1 can produce than actual crystal oscillator frequency postpones to occur for more than one second; 3rd, jamproof robustness is strong, and the synchronous pulse per second (PPS) such as when input is unstable, such as along occur once in a while a little earlier or a little later, that will cause actual samples interval also together to broaden narrowing, cause actual still inaccurate with measurement result; 4th, this method of sampling can only be applicable to difference 4, the factor of 8 grades 4000, and the difference of adjustment can not, more than the width gauge numerical value of 1 sampled point, make range of application have limitation.
In sum, the existing electric parameters signal sampling method scope of application is narrow, along with semi-conductor chip frequency reaches physics limit, and need to improve further sampling rate to obtain high resolving power measurement, during the rated frequency of the sampling rate that therefore can cause further close to chip, the dispersion of sampling interval width and point breaking up of sampling section width become rapidly large.
Summary of the invention
The object of the invention is to there is the scope of application for existing electric parameters signal sampling method narrow, measure numerical value and still there is the defect compared with big error with actual value, a kind of electric parameters signal adaptive method of sampling based on GPS is provided, expand the scope of application of sampling, when the rated frequency of required sampling rate close to chip, the degree that the dispersion of effective suppression each sampling interval width and the dispersion of each sampling section width increase, when making to improve sampling rate when the rated frequency of chip cannot improve, also the acquisition precision decline degree reducing electric parameters becomes possibility.
Technical scheme
In order to realize above-mentioned technical purpose, the present invention designs a kind of electric parameters signal adaptive method of sampling based on GPS, and it is characterized in that, it comprises following step:
The first step, receives standard second pulse, is designated as C with the clock number between the pulse per second (PPS) of clock counter continuous counter and pulse per second (PPS)
0;
Second step, uses the pulse per second (PPS) actual count value of nearest three times, is designated as C
0, C
-1and C
-2, predict that the pulse per second (PPS) interval of next second is designated as C
1;
3rd step, sampling rate is designated as S, and basic interval timer number scale is B
1=int (C
1/ S), once share remainder equally, be designated as G
1=C
1%S, if G
1be zero, then sampling interval all equals B
1, directly jump to the 6th step, otherwise sampling rate remainder is designated as F
1=S%G
1, enter the 4th step;
4th step, if F
1be zero, the sampled point space-number of once sharing equally is designated as D
1=int (S/G
1), [0,3199] after first sampled point individual sampled point basis interval timer number is B
1, at D
1-1,2D
1-1 ..., mD
1add the delay of 1 clock in-1 sampling interval, and once to share number scale be m=[1, int (S/D
1)], jump to the 6th step, otherwise enter the 5th step;
5th step, D
1=int (S/G
1)+1, the sampled point space-number that secondary is shared equally is designated as D
2=int (S/D
2), QR is designated as R
2=G
1-int (S/D
1), [0,3199] after first sampled point individual sampled point basis interval timer number is B
1, at D
1-1,2D
1-1 ..., mD
1the delay of 1 clock is added in-1 sampling interval, and m=[1, int (S/D
1)], at D
2-1,2D
2-1 ..., pD
2add the delay of 1 clock in-1 sampling interval, secondary shares number, is designated as p=[1, R
2], enter the 6th step.
6th step, carries out the sampling of S point, jumps to second step according to the sampling interval generated;
Further, C is worked as in described second step
0, C
-1and C
-2when being initialized as 0, then need C
-2after being received in actual value, just sampling action can be carried out.
Further, in described second step 2, if outside pulse per second (PPS) is unstable, then equalization method is used to calculate C
1=int ((C
0+ C
-1+ C
-2)/3); If crystal oscillator often changes along with temperature, calculate C by slope
1=int ((C
0+ 2C
-1-C
-2)/2).
Beneficial effect
A kind of electric parameters signal adaptive method of sampling based on GPS provided by the invention, expand the scope of application of sampling, when the rated frequency of required sampling rate close to chip, the degree that the dispersion of effective suppression each sampling interval width and the dispersion of each sampling section width increase, when making to improve sampling rate when the rated frequency of chip cannot improve, the acquisition precision decline degree also reducing electric parameters becomes possibility.The present invention can predict rapidly next second clock frequency according to the variation tendency of clock frequency before, there is stronger anti-interference, sample frequency can not be produced because of the pulse per second (PPS) instability of input significantly to change, the difference particularly adjusted allows, more than the width of 1 sampled point, to expand range of application.
Accompanying drawing explanation
Accompanying drawing 1 is operational flowchart of the present invention.
Embodiment
Below in conjunction with drawings and Examples, the present invention will be further described.
Embodiment
As shown in Figure 1, with the clock count value C of actual standard pulse per second (PPS)
0=2500020, C
-1=2500014, C
-2=2500018, sampling rate S=3200 point/seconds, often wheel calculates needs 64 sampled points is example,
The first step, receives standard second pulse, with the clock number C between the pulse per second (PPS) of clock counter continuous counter and pulse per second (PPS)
0, C
-1, C
-2;
Second step, predicts the count value C at next pulse per second (PPS) interval
1=int ((C
0+ C
-1+ C
-2)/3)=2500017;
3rd step, the theoretical total linear spacing of often taking turns 64 sampled points in computation period is C
1* 64/S=50000.34clocks; Basis interval timer number B
1=int (C
1/ S)=781; Once share remainder G equally
1=C
1%S=817; Because G
1non-vanishing, then F
1=S%G
1=749, enter the 4th step,
4th step, if F
1non-vanishing, enter the 5th step;
5th step, D
1=int (S/G
1)+1=4, D
2=int (S/D
2)=188, R
2=G
1-int (S/D
1)=17; The interval, sampled point basis of 0 ~ 3199 after first sampled point is all 781 clocks; Because D
1=4, so the 3rd, 7 ..., 4m-1, adds the delay of 1 clock in 3199 sampling interval, and m=1 ~ int (S/D
1); Because D
2=188, so the 187th, 375 ..., 188p-1, adds the delay of 1 clock, p=1 ~ R in 3195 sampling interval
2.
6th step, the sampling interval according to generating carries out 3200 point samplings, redirect second step.
Often to take turns 64 sampled points in computation period, in 3200 o'clock, have 50 to take turns calculating, often take turns actual total linear spacing and theoretical value average error is 0.0009%, concrete evaluation is as following table:
Invention which is intended to be protected is not limited only to the technical scheme that above-described embodiment relates to, and any enlightenment by technical scheme of the present invention, all should fall into protection scope of the present invention.
Claims (3)
1., based on an electric parameters signal adaptive method of sampling of GPS, it is characterized in that, it comprises following step:
(I) receive standard second pulse, with the clock number between the pulse per second (PPS) of clock counter continuous counter and pulse per second (PPS), be designated as C
0;
(II) use the pulse per second (PPS) actual count value of nearest three times, be designated as C
0, C
-1and C
-2, predict the pulse per second (PPS) interval of next second, be designated as C
1;
(III) sampling rate is designated as S, and basic interval timer number scale is B
1=int (C
1/ S), once share remainder equally, be designated as G
1=C
1%S, if G
1be zero, then sampling interval all equals B
1, directly jump to step (VI), otherwise sampling rate remainder is designated as F
1=S%G
1, enter step (IV);
(IV) if F
1be zero, the sampled point space-number of once sharing equally is designated as D
1=int (S/G
1), [0,3199] after first sampled point individual sampled point basis interval timer number is B
1, at D
1-1,2D
1-1 ..., mD
1add the delay of 1 clock in-1 sampling interval, and once to share number scale be m=[1, int (S/D
1)], jump to step (VI), otherwise enter step (V);
(V) D
1=int (S/G
1)+1, the sampled point space-number that secondary is shared equally is designated as D
2=int (S/D
2), QR is designated as R
2=G
1-int (S/D
1), [0,3199] after first sampled point individual sampled point basis interval timer number is B
1, at D
1-1,2D
1-1 ..., mD
1the delay of 1 clock is added in-1 sampling interval, and m=[1, int (S/D
1)], at D
2-1,2D
2-1 ..., pD
2add the delay of 1 clock in-1 sampling interval, secondary shares number, is designated as p=[1, R
2], enter step (VI).
(VI) carry out the sampling of S point according to the sampling interval generated, jump to step (II).
2. a kind of electric parameters signal adaptive method of sampling based on GPS as claimed in claim 1, is characterized in that: work as C in described step (II)
0, C
-1and C
-2when being initialized as 0, then need C
-2after being received in actual value, just sampling action can be carried out.
3. a kind of electric parameters signal adaptive method of sampling based on GPS as claimed in claim 1, is characterized in that: in described step (II), if outside pulse per second (PPS) is unstable, then uses equalization method to calculate C
1=int ((C
0+ C
-1+ C
-2)/3); If crystal oscillator often changes along with temperature, calculate C by slope
1=int ((C
0+ 2C
-1-C
-2)/2).
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510016798.3A CN104614586B (en) | 2015-01-13 | 2015-01-13 | A kind of electrical quantity signal adaptive method of sampling based on GPS |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510016798.3A CN104614586B (en) | 2015-01-13 | 2015-01-13 | A kind of electrical quantity signal adaptive method of sampling based on GPS |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104614586A true CN104614586A (en) | 2015-05-13 |
| CN104614586B CN104614586B (en) | 2019-01-18 |
Family
ID=53149122
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510016798.3A Active CN104614586B (en) | 2015-01-13 | 2015-01-13 | A kind of electrical quantity signal adaptive method of sampling based on GPS |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104614586B (en) |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6271773B1 (en) * | 1999-07-13 | 2001-08-07 | Teratec Corporation | Coherent sampling method and apparatus |
| CN100520419C (en) * | 2003-07-11 | 2009-07-29 | Nxp股份有限公司 | Improved frequency determination |
| CN101719765A (en) * | 2009-11-25 | 2010-06-02 | 中兴通讯股份有限公司 | Method and device for generating low-jitter clock |
| CN101762742A (en) * | 2009-12-16 | 2010-06-30 | 南京弘毅电气自动化有限公司 | Sampling value correcting method based on fixed-frequency sampling sequence |
| CN101848076A (en) * | 2010-05-10 | 2010-09-29 | 新邮通信设备有限公司 | Clock signal holding method and device |
| CN102354108A (en) * | 2011-06-09 | 2012-02-15 | 国网电力科学研究院 | Method for realizing dynamic sampling rate adjustment by oversampling |
| CN202256483U (en) * | 2011-10-26 | 2012-05-30 | 江西省电力公司超高压分公司 | GPS (Global Positioning System) second-based real-time self-adaptive evenly-spaced sampling synchronous data acquisition device |
| CN102778606A (en) * | 2012-07-26 | 2012-11-14 | 南京邮电大学 | Digital measurement method for instantaneous phase of alternating current |
| CN103018555A (en) * | 2012-12-24 | 2013-04-03 | 电子科技大学 | High-precision electric power parameter software synchronous sampling method |
| CN103185837A (en) * | 2013-03-25 | 2013-07-03 | 华中科技大学 | Method for measuring frequency of power system |
| CN103713552A (en) * | 2013-12-23 | 2014-04-09 | 国电南瑞科技股份有限公司 | Self-adaptive dynamic synchronous sampling control device and method based on pulse per second |
| CN103941622A (en) * | 2014-04-28 | 2014-07-23 | 国家电网公司 | Method for adopting high-accuracy pulse per second frequency multiplication to produce sampling pulse based on FPGA |
-
2015
- 2015-01-13 CN CN201510016798.3A patent/CN104614586B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6271773B1 (en) * | 1999-07-13 | 2001-08-07 | Teratec Corporation | Coherent sampling method and apparatus |
| CN100520419C (en) * | 2003-07-11 | 2009-07-29 | Nxp股份有限公司 | Improved frequency determination |
| CN101719765A (en) * | 2009-11-25 | 2010-06-02 | 中兴通讯股份有限公司 | Method and device for generating low-jitter clock |
| CN101762742A (en) * | 2009-12-16 | 2010-06-30 | 南京弘毅电气自动化有限公司 | Sampling value correcting method based on fixed-frequency sampling sequence |
| CN101848076A (en) * | 2010-05-10 | 2010-09-29 | 新邮通信设备有限公司 | Clock signal holding method and device |
| CN102354108A (en) * | 2011-06-09 | 2012-02-15 | 国网电力科学研究院 | Method for realizing dynamic sampling rate adjustment by oversampling |
| CN202256483U (en) * | 2011-10-26 | 2012-05-30 | 江西省电力公司超高压分公司 | GPS (Global Positioning System) second-based real-time self-adaptive evenly-spaced sampling synchronous data acquisition device |
| CN102778606A (en) * | 2012-07-26 | 2012-11-14 | 南京邮电大学 | Digital measurement method for instantaneous phase of alternating current |
| CN103018555A (en) * | 2012-12-24 | 2013-04-03 | 电子科技大学 | High-precision electric power parameter software synchronous sampling method |
| CN103185837A (en) * | 2013-03-25 | 2013-07-03 | 华中科技大学 | Method for measuring frequency of power system |
| CN103713552A (en) * | 2013-12-23 | 2014-04-09 | 国电南瑞科技股份有限公司 | Self-adaptive dynamic synchronous sampling control device and method based on pulse per second |
| CN103941622A (en) * | 2014-04-28 | 2014-07-23 | 国家电网公司 | Method for adopting high-accuracy pulse per second frequency multiplication to produce sampling pulse based on FPGA |
Non-Patent Citations (1)
| Title |
|---|
| 杨启勇等: "基于GPS授时的脉冲峰值采集系统", 《测控技术》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104614586B (en) | 2019-01-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Yao et al. | A novel method for phasor measurement unit sampling time error compensation | |
| CN103941622B (en) | Occur frequently again the method for sampling pulse based on the high accuracy pulse per second (PPS) of FPGA | |
| CN109283864B (en) | Time synchronization and calibration method and system for data sampling | |
| CN102118007B (en) | Electronic mutual inductor synchronization method based on PPS | |
| CN109428598B (en) | Calibration method and system for data sampling | |
| MX2011012664A (en) | Apparatus and method for estimating synchronized phasors at predetermined times referenced to a common time standard in an electrical system. | |
| CN103616814A (en) | Synchronous sampling clock closed loop correcting method and system based on FPGA | |
| CN102590784A (en) | Single-phase intelligent energy meter distributed correction method | |
| CN104635094A (en) | Method for improving PMU (power management unit) synchronous phasor measurement precision | |
| Sinha et al. | A pre-filter based PLL for three-phase grid connected applications | |
| CN106932642B (en) | Electric Power Harmonic Analysis method | |
| CN104460311A (en) | Time calibration method and device | |
| CN101216511A (en) | Digital synchronous sampling method | |
| CN103383412A (en) | Adaptive software and hardware frequency tracking and sampling method | |
| CN102902879A (en) | Synchronous phasor calculation method based on discrete Fourier transform (DFT) recursion of field programmable gate array (FPGA) hardware | |
| CN108333426A (en) | Power system frequency measurement method based on fourier algorithm | |
| CN103983847A (en) | Self-adaptive frequency tracking measurement method based on RLS (Recursive Least Squares) in synchronized phasor measurement | |
| KR101795199B1 (en) | Signal processing device | |
| CN108964092B (en) | The method of gradual strategy adjustment low-voltage network imbalance of three-phase voltage | |
| Busarello et al. | Zero-crossing detection frequency estimator method combined with a Kalman filter for non-ideal power grid | |
| CN102170143A (en) | Microcomputer automatic quasi-synchronization paralleling realization method | |
| CN104614586A (en) | GPS-based self-adaptive sampling method for electric quantity signals | |
| CN102346219B (en) | Method for detecting phases of access point voltages of voltage source inverter by using three-phase software phase-locked loop | |
| US10613127B2 (en) | Determining the frequency of an alternating signal | |
| CN102636769A (en) | Error counter for electric energy meter and error counting 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 |