CN110361608B - Multi-path asynchronous sampling data processing method based on sampling point geometric scaling - Google Patents
Multi-path asynchronous sampling data processing method based on sampling point geometric scaling Download PDFInfo
- Publication number
- CN110361608B CN110361608B CN201910474411.7A CN201910474411A CN110361608B CN 110361608 B CN110361608 B CN 110361608B CN 201910474411 A CN201910474411 A CN 201910474411A CN 110361608 B CN110361608 B CN 110361608B
- Authority
- CN
- China
- Prior art keywords
- sampling
- signal
- frequency
- path
- calculating
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R1/00—Details of instruments or arrangements of the types included in groups G01R5/00 - G01R13/00 and G01R31/00
- G01R1/02—General constructional details
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Synchronisation In Digital Transmission Systems (AREA)
- Complex Calculations (AREA)
Abstract
The invention discloses a multi-channel asynchronous sampling data processing method based on sampling point geometric scalability, and relates to the field of power distribution automation. The method comprises the following steps: measuring the frequency of each signal; setting the number of sampling points per period by taking any one of the signals as a basic signal; calculating sampling frequency and sampling interval according to the frequency of the basic signal and the number of sampling points per period; calculating the cycle difference of each path of signals except the basic signals; calculating the expansion common ratio of each path of signal according to the sampling interval and the period difference; calculating an initial coefficient and a shift coefficient of each path of signal according to the expansion common ratio; carrying out data sampling on each path of signal at the sampling frequency to obtain a sampling data sequence; and carrying out point-by-point equal ratio correction on the data of the sampling sequence. The invention reduces the calculation error and improves the synchronization processing precision.
Description
Technical Field
The invention relates to the field of distribution automation, in particular to a multi-channel asynchronous sampling data processing method based on sampling point geometric expansion.
Background
Along with the wide access of the distributed power supply, the adjustment of the operation mode of the power distribution network increases the operation of the distributed power supply on/off the network and on the network besides the original operation of closing and opening the ring, wherein if the operation of closing the ring and on the network is not processed properly, the power distribution network is impacted greatly, and the reliability and the stability of the power distribution network are influenced.
The closed loop and grid-connected operation involves two or more systems, the frequency of each system is different, if sampling is carried out at a frequency suitable for one system, other systems obtain asynchronous sampling data, further, a power signal measurement error is generated, and the operation requirement cannot be met. The accurate synchronous sampling measurement of a plurality of frequency signals is generally realized by a plurality of sets of hardware systems, and the cost is too high for power distribution network equipment, so that a synchronous processing method of asynchronous sampling data of different frequency systems is needed.
Disclosure of Invention
Aiming at the defects of the prior art, the invention aims to provide a multi-channel asynchronous sampling data processing method based on sampling point geometric scaling so as to solve the problem that the accurate synchronous sampling measurement of a plurality of frequency signals in the prior art can be realized only by a plurality of sets of hardware systems.
In order to solve the problems of the prior art, the invention adopts the technical scheme that:
a multipath asynchronous sampling data processing method based on sampling point geometric scalability comprises the following steps:
measuring the frequency of each signal;
setting the number of sampling points per period by taking any one of the signals as a basic signal;
calculating sampling frequency and sampling interval according to the frequency of the basic signal and the number of sampling points per period;
calculating the cycle difference of each path of signals except the basic signals;
calculating the expansion common ratio of each path of signal according to the sampling interval and the period difference;
calculating an initial coefficient and a shift coefficient of each path of signal according to the expansion common ratio;
carrying out data sampling on each path of signal at the sampling frequency to obtain a sampling data sequence;
and carrying out point-by-point equal ratio correction on the data of the sampling sequence.
Further, the sampling frequency and the sampling interval are calculated as follows:
fs＝f0*N；
ts＝1/fs＝1/(f0*N)；
wherein f issTo sample frequency, f0Based on the frequency of the signal, N being the number of samples per cycle, tsIs the sampling interval.
Further, the method for calculating the period difference comprises the following steps: dTn＝(fn–f0)/(f0*fn)；
Wherein, dTnIs the period difference of the nth signal, fnIs the nth signal frequency.
Further, the method for calculating the expansion ratio comprises the following steps:
qn＝dTn/(N*ts)；
wherein q isnIs the expansion ratio of the nth signal.
Further, the method for calculating the start coefficient and the shift coefficient comprises;
m_idn＝trunc[-1*(i+1)*qn]；
wherein s is_idnIs the start coefficient of the nth signal, m _ idnIs a shift coefficient of the nth signal, trunc [ 2 ]]Indicating the decimal part of the removed number, and i is the sample point number.
Further, the method for geometric correction comprises the following steps:
Yn(i)＝Xn(i+m_idn)–[Xn(i+1+s_idn+m_idn)–Xn(i+s_idn+m_idn)]*[(i+1)*qn+m_idn]；
wherein, Yn(i) For data synchronous to the i-th sample point of the n-th signal, XnIs an asynchronous sampling data sequence of the nth signal.
Compared with the prior art, the invention has the following beneficial effects:
the method converts the multi-path asynchronous sampling data into the synchronous sampling data through lower sampling frequency, realizes the accurate measurement and calculation of the power signals of different frequency systems, does not need to change hardware design, avoids increasing hardware cost, solves the synchronous processing problem of the multi-path asynchronous sampling data through a software calculation mode, can self-adaptively adjust the number of sampling points, reduces calculation errors through an initial coefficient and a shift coefficient, and improves the synchronous processing precision.
Drawings
FIG. 1 is a flow chart of the present invention.
Detailed Description
The invention is further described below with reference to the accompanying drawings. The following examples are only for illustrating the technical solutions of the present invention more clearly, and the protection scope of the present invention is not limited thereby.
As shown in fig. 1, a method for synchronizing multi-channel asynchronous sampling data based on sampling point geometric expansion includes the following specific implementation:
the method comprises the following steps: measuring the frequency f of each signal and recording as f0,f1,…,fn，fnIs the nth signal frequency;
step two: setting any one of the signals as basic signal, setting the number of sampling points N per period, and taking the product of the frequency of the basic signal and the number of the sampling points per period as sampling frequency fsCalculating the sampling interval ts；
If at f0Setting the number of sampling points per period as N and the sampling frequency f for the basic signals＝f0N, sampling interval ts＝1/fs＝1/(f0*N)。
Step three: calculating the period difference dT between each path of signal and the basic signal, and recording as dT0,dT1,…,dTn，dTn＝(fn–f0)/(f0*fn)，dTnThe cycle difference of the nth signal is obtained;
step four: calculating the expansion ratio q of each path of signal according to the sampling interval and the period difference of each path of signal, and recording as q0,q1,…,qn，qn＝dTn/(N*ts)，qnThe expansion ratio of the nth signal is obtained;
step five: calculating an initial coefficient s _ id and a shift coefficient m _ id of each path of signal according to the expansion ratio, and respectively recording the initial coefficient s _ id and the shift coefficient m _ id as s _ id0,s_id1,…,s_idnAnd m _ id0,m_id1,…,m_idnWherein s is_idnIs the start coefficient of the nth signal, m _ idnThe shift coefficient of the nth signal is;
the start coefficient s _ id is related to the signal frequency, when the signal frequency fnGreater than the base signal f0When, s _ idn-1; when the signal frequency fnLess than the base signal f0Time s _ idn＝0。
The shift coefficient m _ id is related to the sampling point serial number and the expansion common ratio, m _ idn＝trunc[-1*(i+1)*qn]Wherein trunc [ 2 ], []The decimal part of the removed number is shown, and i represents the sample point number.
Step six: carrying out data sampling on each path of signal by using the sampling frequency to obtain a sampling data sequence X0(i),X1(i),…,Xn(i)，Xn(i) When the number of sampling points in each cycle is N, the sampling point number is-1, 0,1, …, N-1, N +1 (the number-1 represents the last sampling point in the previous cycle, and the number N, N +1 represents the first two sampling points in the next cycle);
step seven: according to a preset calculation formula and the expansion common ratio q, the initial coefficient s _ id and the shift coefficient m _ id of each path of signal, performing point-by-point equal ratio correction on the data of each path of sampling sequence, and generating a synchronous data sequence Y according to the following formula0(i),Y1(i),…,Yn(i)。
Yn(i)＝Xn(i+m_idn)–[Xn(i+1+s_idn+m_idn)–Xn(i+s_idn+m_idn)]*[(i+1)*qn+m_idn]。
Wherein, Yn(i) For data synchronous to the i-th sample point of the n-th signal, XnIs an asynchronous sampling data sequence of the nth signal.
The technical solution of the present invention is clearly and completely described below by taking data samples with three frequencies of 49Hz,50Hz and 51.5Hz as research examples.
1) Measuring the frequency of each signal, denoted as f0＝49Hz,f1＝50Hz,f2＝51.5Hz。
2) Setting f1Based on 50Hz signal, the number of sampling points per period is set to be N-80 points, and the sampling frequency fs4000Hz, sample interval ts＝0.25ms。
3) Calculating the difference dT, dT between the period of the signal and the base signal0＝-4.082*10-4s，dT2＝5.825*10-4s。
4) Calculating the scaling ratio q, q of the signal according to the sampling interval and the period difference of the signal0＝-0.0204，q2＝0.0291。
5) And calculating a start coefficient s _ id and a shift coefficient m _ id of each path of signal according to the expansion common ratio.
Initial coefficient s _ id of 49Hz signal0Initial coefficient s _ id of 0, 51.5Hz signal2＝-1。
When the number of sampling points per period is 80, the 49Hz signal is shifted by a coefficient m _ id0And 51.5Hz signal shift coefficient m _ id2Respectively as follows:
6) carrying out data sampling on each path of signal by using the sampling frequency to obtain a sampling data sequence X0(i),X1(i),X2(i) Sample point numbers i ═ 1,0,1, …,79,80,81, where X1(i) For synchronizing the sampling sequences, X0(i),X2(i) Is an asynchronous sample sequence.
7) Generating X according to the following formula0(i) Corresponding synchronous data sequence Y0(i)：
Y0(0)＝X0(0)–[X0(1)–X0(0)]*q0，
Y0(48)＝X0(48)–[X0(49)–X0(48)]*49*q0，
Y0(49)＝X0(50)–[X0(51)–X0(50)]*(50*q0+1)，
…
Y0(79)＝X0(80)–[X0(81)–X0(80)]*(80*q0+1)。
Generating X according to the following formula2(i) Corresponding synchronous data sequence Y2(i)：
Y2(0)＝X2(0)–[X2(0)–X2(-1)]*q2，
Y2(33)＝X2(33)–[X2(33)–X2(32)]*34*q2，
Y2(34)＝X2(33)–[X2(33)–X2(32)]*(35*q2–1)，
Y2(67)＝X2(66)–[X2(66)–X2(65)]*(68*q2–1)，
Y2(68)＝X2(66)–[X2(66)–X2(65)]*(69*q2–2)，
…
Y2(79)＝X2(77)–[X2(77)–X2(76)]*(80*q2–2)。
The invention solves the synchronous processing problem of the multi-path asynchronous sampling data in a software processing mode, does not need to change the hardware design, avoids increasing the hardware cost and solves the problem of improper loop closing and grid connection operation caused by the traditional multi-path asynchronous sampling precision error.
The above-mentioned embodiments further describe in detail the technical problems, technical solutions and advantageous effects that are solved by the present invention, and the technical solutions disclosed in the solutions of the present invention should not be limited to the technical solutions disclosed in the above-mentioned embodiments, but also include technical solutions that are formed by any combination of the above technical features. It should be noted that those skilled in the art can make various improvements and modifications without departing from the principle of the present invention, and such improvements and modifications are also considered to be within the scope of the present invention.
Claims (1)
1. A multipath asynchronous sampling data processing method based on sampling point geometric scalability is characterized by comprising the following steps:
measuring the frequency of each signal;
setting the number of sampling points per period by taking any one of the signals as a basic signal;
calculating sampling frequency and sampling interval according to the frequency of the basic signal and the number of sampling points per period;
calculating the period difference between each path of signals except the basic signals and the basic signals;
calculating the expansion common ratio of each path of signal according to the sampling interval and the period difference;
calculating an initial coefficient and a shift coefficient of each path of signal according to the expansion common ratio;
wherein the start coefficient s _ id is related to the signal frequency fnGreater than the fundamental signal frequency f0When, s _ idn-1; when the signal frequency fnLess than the fundamental signal frequency f0Time s _ idn0; the shift coefficient m _ id is related to the sampling point serial number and the expansion common ratio, m _ idn＝trunc[-1*(i+1)*qn]Wherein trunc [ 2 ], []Representing a decimal part of a removed number, i represents a sampling point serial number; q. q.snThe expansion ratio of the nth signal is obtained;
carrying out data sampling on each path of signal at the sampling frequency to obtain a sampling data sequence;
performing point-by-point equal ratio correction on the data of the sampling data sequence; the geometric correction is that a synchronous data sequence is generated according to the following formula:
Yn(i)＝Xn(i+m_idn)–[Xn(i+1+s_idn+m_idn)–Xn(i+s_idn+m_idn)]*[(i+1)*qn+m_idn]；
wherein, Yn(i) For data synchronous to the i-th sample point of the n-th signal, XnFor asynchronously sampled data sequences of the nth signal, s_idnIs the start coefficient of the nth signal, m _ idnThe shift coefficient of the nth signal is;
the calculation method of the sampling frequency and the sampling interval is as follows:
fs＝f0*N；
ts＝1/fs＝1/(f0*N)；
wherein f issTo sample frequency, f0Based on the frequency of the signal, N being the number of samples per cycle, tsIs the sampling interval;
the calculation method of the period difference comprises the following steps: dTn＝(fn–f0)/(f0*fn)；
Wherein, dTnIs the period difference of the nth signal, fnIs the nth signal frequency;
the computing method of the expansion ratio comprises the following steps:
qn＝dTn/(N*ts)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910474411.7A CN110361608B (en) | 2019-05-31 | 2019-05-31 | Multi-path asynchronous sampling data processing method based on sampling point geometric scaling |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910474411.7A CN110361608B (en) | 2019-05-31 | 2019-05-31 | Multi-path asynchronous sampling data processing method based on sampling point geometric scaling |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110361608A CN110361608A (en) | 2019-10-22 |
| CN110361608B true CN110361608B (en) | 2022-02-15 |
Family
ID=68215478
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910474411.7A Active CN110361608B (en) | 2019-05-31 | 2019-05-31 | Multi-path asynchronous sampling data processing method based on sampling point geometric scaling |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110361608B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114609997B (en) * | 2022-03-07 | 2023-06-27 | 南京国电南自维美德自动化有限公司 | Method for realizing analog input change rate protection of process control station |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4703448A (en) * | 1984-10-25 | 1987-10-27 | Nicolet Instrument Corporation | Method for measuring skew or phase difference in electronic systems |
| CN102664630A (en) * | 2012-04-28 | 2012-09-12 | 重庆大学 | Synchronous acquisition method of multichannel signals inside mechanical systems |
| CN102946251A (en) * | 2012-11-06 | 2013-02-27 | 中国科学院上海微系统与信息技术研究所 | Method for realizing multi-channel synchronous sampling by multi-channel asynchronous sampling ADC (Analog Digital Converter) |
| CN104483563A (en) * | 2014-11-27 | 2015-04-01 | 广东电网有限责任公司电力科学研究院 | Method and system for synchronous sampling of power signals |
| CN107102220A (en) * | 2017-03-14 | 2017-08-29 | 全球能源互联网研究院 | Recorder data processing method and processing device |
-
2019
- 2019-05-31 CN CN201910474411.7A patent/CN110361608B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4703448A (en) * | 1984-10-25 | 1987-10-27 | Nicolet Instrument Corporation | Method for measuring skew or phase difference in electronic systems |
| CN102664630A (en) * | 2012-04-28 | 2012-09-12 | 重庆大学 | Synchronous acquisition method of multichannel signals inside mechanical systems |
| CN102946251A (en) * | 2012-11-06 | 2013-02-27 | 中国科学院上海微系统与信息技术研究所 | Method for realizing multi-channel synchronous sampling by multi-channel asynchronous sampling ADC (Analog Digital Converter) |
| CN104483563A (en) * | 2014-11-27 | 2015-04-01 | 广东电网有限责任公司电力科学研究院 | Method and system for synchronous sampling of power signals |
| CN107102220A (en) * | 2017-03-14 | 2017-08-29 | 全球能源互联网研究院 | Recorder data processing method and processing device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110361608A (en) | 2019-10-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107563550B (en) | PMU-based power distribution network real-time distributed state estimation and PMU optimal configuration method | |
| CN110361608B (en) | Multi-path asynchronous sampling data processing method based on sampling point geometric scaling | |
| CN103941622B (en) | Occur frequently again the method for sampling pulse based on the high accuracy pulse per second (PPS) of FPGA | |
| CN104090151A (en) | Zero-crossing-point accurate detection method for AC power-grid voltage signals | |
| EP4035270A1 (en) | Dynamic successive approximation register (sar) analog-to-digital converter (adc) (sar-adc) clock delay calibration systems and methods | |
| CN202256483U (en) | GPS (Global Positioning System) second-based real-time self-adaptive evenly-spaced sampling synchronous data acquisition device | |
| CN106099932B (en) | Day-ahead planning power flow analysis method considering uncertainty time-space correlation | |
| CN109946513B (en) | Cross-interval digital metering method, device and system based on network sampling absolute delay compensation | |
| CN110780110B (en) | Fundamental voltage zero crossing point automatic detection method and system and sampling device | |
| CN113162626A (en) | High-precision multichannel synchronous high-speed data acquisition and processing system and method | |
| CN109188987A (en) | Multi channel signals high-precise synchronization control method and device based on high speed D/A | |
| CN113533815A (en) | Multi-channel sampling synchronization method based on time stamps | |
| CN111458598A (en) | Method for aligning multiple homologous recording waveforms of asynchronous sampling | |
| CN109581018A (en) | Wide-band sampling oscillograph time base error synchroballistic method and device | |
| CN106383438B (en) | One kind taming and dociling clock method based on sliding window time extension high-precision | |
| CN103124212B (en) | Digital control system on-the-spot bus time synchronous method and device based on Bayes principle | |
| CN108011683B (en) | Large-scale synthesis sensing system multistage distributes time-frequency unified approach | |
| CN106771326B (en) | A method of based on incremental optical-electricity encoder measuring speed | |
| CN111211855B (en) | Mixed clock synchronization method for distributed processing system | |
| CN102928805A (en) | Multi-meter position instrument for calibrating digital energy meters | |
| CN108596475B (en) | PMU data recovery method based on dynamic change of interpolation interval | |
| CN107870338B (en) | A kind of satellite navigation carrier wave tracing method of low update frequency | |
| CN209496263U (en) | A kind of multi-disc ADC synchronous | |
| CN116989973A (en) | Temporary impact type high-speed wind tunnel continuous variable attack angle pressure measurement test data processing method | |
| CN109765302A (en) | Delay inequality high-precision analog device and method between a kind of supersonic array signal path |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |