CN102135567A - Real-time frequency tracking and harmonic measuring method for AC sampling of power system - Google Patents
Real-time frequency tracking and harmonic measuring method for AC sampling of power system Download PDFInfo
- Publication number
- CN102135567A CN102135567A CN2011100464591A CN201110046459A CN102135567A CN 102135567 A CN102135567 A CN 102135567A CN 2011100464591 A CN2011100464591 A CN 2011100464591A CN 201110046459 A CN201110046459 A CN 201110046459A CN 102135567 A CN102135567 A CN 102135567A
- Authority
- CN
- China
- Prior art keywords
- digital
- sampling
- fast fourier
- frequency
- signal
- 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
Images
Landscapes
- Measurement Of Current Or Voltage (AREA)
Abstract
The invention discloses a real-time frequency tracking and harmonic measuring method for the AC sampling of a power system, which is characterized by comprising the steps: an AC sampling input is made to enter an analog-to-digital converter at a fixed sampling rate; the analog-to-digital converter at the fixed sampling rate respectively outputs a first segment of digital signals with the length of N points and a second segment of digital signals with the length of N points to enter a first fast Fourier module; the first fast Fourier module outputs phase angles of the highest frequency components of the first segment of digital signals with the length of N points and the second segment of digital signals with the length of N points, and sends the phase angles into a phase discriminator one after another; the phase discriminator outputs a correction frequency, and sends the correction frequency into a digital signal reconstruction filter; the digital signal reconstruction filter performs digital reconstruction sampling on the first segment of digital output signals with the length of N points, of the analog-to-digital converter at the fixed sampling rate, according to the correction frequency, and sends an obtained N-point reconstruction sampling signal into a second fast Fourier module; the phase angle output of the highest frequency component of the second fast Fourier module is sent into the phase discriminator; the third segment of digital signals with the length of N points, of the analog-to-digital converter, enter the first fast Fourier module; and the phase discriminator judges the phase angle based on accuracy requirements, and outputs a correction frequency to a digital signal reconstruction filter. After the iterations of the process reach M, the output of the second fast Fourier module is the values of respective subharmonic components of the AC sampling signal, and the output frequency of the phase discriminator is the frequency of the AC sampling signal.
Description
Technical field
The invention belongs to AC sampling and measure applied technical field, relate more specifically to electric system industrial frequency AC sampled measurements real-time frequency and follow the tracks of and the harmonic measure technology.
Background technology
Fast development along with electric system, the power network capacity constantly increases, structure is increasingly sophisticated, it is very important that the robotization of monitoring in real time in the electric system, dispatching just seems, and digital collection is the important step that realizes robotization, especially how the electric weight of the simulation of each in the acquisition system accurately and rapidly is the focus that work about electric power person pays close attention to always.According to the difference of sampled signal, can be divided into direct current sampling and AC sampling 2 big classes.The direct current sampling is the DC voltage that alternating voltage, current signal is converted into 0~5V, and the major advantage of this method is that algorithm is simple, is convenient to filtering, but invest bigger, safeguard complexity, can't realize the live signal collection, thereby be restricted in application on power system.AC sampling be of ac is converted into ± 5V (or gather, and major advantage is that real-time is good, and phase distortion is little, small investment, is convenient to safeguard by 0~5V) alternating voltage; Its shortcoming is the algorithm complexity, and precision is difficult to improve, and the A/D slewing rate is had relatively high expectations.Along with the development of microcomputer technology, AC sampling has the trend that progressively replaces the direct current sampling with its excellent ratio of performance to price.
Summary of the invention
The objective of the invention is the A/D slewing rate to be had relatively high expectations for the data acquisition that overcomes electric power power frequency component ac electric, precision is difficult to the shortcoming of conventional AC Sampling techniques such as improving, and a kind of real-time frequency-tracking and high-precision AC quantity measuring method fast is provided.
A kind of electric power system alternating current sampling real-time frequency of the present invention is followed the tracks of and harmonic measuring method is achieved by the following technical solution, and the described real-time synchronous AC method of sampling comprises that step is as follows:
1, gets AC sampling input 1 and enter fixed sample rate analog to digital converter 2.
2, fixed sample rate analog to digital converter 2 is exported the digital signal of first section N point length respectively, second section N point length digital signal enters the first fast Fourier module (FFT) 3, the digital signal of first section N point length of the first fast Fourier module (FFT), 3 outputs, the highest frequency component phase angle 8 of second section N point length digital signal are successively sent into phase detector 4, and phase detector 4 output calibration frequencies 10 are sent into digital signal reconfigurable filter 5.
3, digital signal reconfigurable filter 5 carries out the digital reconstruction sampling according to the digital output signal 12 of first section N point length of proofreading and correct 10 pairs of fixed sample rate analog to digital converters 2 of back frequency, and the N point reconstructed sample signal that obtains is sent into the second fast Fourier module 6.
4, phase detector 4 is sent in the highest frequency component phase angle of the second fast Fourier module 6 output 9, the digital signal of the 3rd section N point of analog to digital converter 2 length enters the first fast Fourier module 3, and output calibration frequency 10 was sent into digital signal reconfigurable filter 5 after phase detector 4 was judged according to accuracy requirement.
5, repeating step 2,3 and 4, reach predetermined iterations M after, 6 outputs 7 of the second fast Fourier module are the each harmonic component value of ac sampling signal, phase detector 4 output frequencies 10 are the frequency of ac sampling signal.
Description of drawings
Fig. 1 is a structural representation of the present invention.
Fig. 2 is a practical measuring examples that contains harmonic noise AC sampling input signal of the present invention.
Embodiment
Fig. 1 is the inventive method structural representation.The AC sampling input 1 of stable state electrical network is expressed as x (t), and its work frequency is expressed as f
0Because harmonic interference and other factors disturb the frequency departure that is caused to make the frequency of power network signal fluctuation occur, analog to digital converter 1 is with fixed sample rate f
0N times of speed the one-period signal of signal x (t) is sampled, digital length is the N point.Sampled signal is expressed as x
i(nT
s), n=0 wherein, 1,2, L N-1 is an integer, T
sBe sampling time interval.
Simulation power network signal x (t) is sampled, and getting two segment length on time shaft respectively all is that the numeral that N is ordered obtains x
m(n) and x
n(n), n=0,1,2, L N-1.Then with x
m(n) and x
n(n) send into first Fast Fourier Transform (FFT) (FFT) module 3, be output as X respectively
m(k) and X
n(k), k=0,1,2, L, N-1
M highest frequency component constantly is expressed as X
m(k), k is expressed as the frequency component of input signal.The initial phase that obtains the high frequency components of m moment N point input signal from fast Fourier transform module (FFT) output is:
In the formula, Im and Re represent real and imaginary part respectively.At the n of m after the moment constantly, fast Fourier transform module (FFT component) highest frequency component of N point numeral is expressed as X
n(k), its initial phase angle is:
Utilize these two phasing degree can obtain the tracking frequencies deviation
The minimum frequency deviation that obtains by (1-3) formula be frequency resolution far above the frequency resolution precision that conventional fast Fourier transform FFT obtains, have speed of convergence fast simultaneously.Fixed sampling frequency f
0With fine frequency f
ΔAddition obtains the new frequency of power network signal:
f
new=f
0+f
Δ (1-4)
So just can obtain an estimation frequency value of approaching electrical network power frequency component frequency relatively, power network signal x (t) is pressed new sample rate f with latter two periodic signal
sSampling, wherein
f
s=f
new·N (1-5)
From the output signal 13 of reconfigurable filter 5, get two sections N continuous and put digital signal, send into module 6 and be Fast Fourier Transform (FFT) FFT again, new sampling rate after (1-1) to (1-5) process of utilization obtains proofreading and correct, the iteration said process just can infinitely be approached the frequency of true power network signal, thereby the high precision real-time frequency that realizes AC sampling input signal 1 is followed the tracks of and measured.
Phase detector 2 estimations obtain new sample rate f
NewAfter, to entering analog to digital converter 2 AC sampling input 1 before by new sample rate f
NewCarry out resampling, this process realizes by reconfigurable filter 5.Reconfigurable filter 5 is a kind of digital signal reconfigurable filters, the output sampled signal 12 of analog to digital converter 2 is carried out numeral recover reconstruct.
Fig. 2 is a practical measuring examples that contains harmonic noise AC sampling input signal that adopts the inventive method.The AC sampling input 1 of stable state electrical network is:
F wherein
0Get 52, N gets 128, implements following steps according to the inventive method:
1, gets AC sampling input 1 and enter fixed sample rate analog to digital converter 2.
2, fixed sample rate analog to digital converter 2 is exported the digital signal of first section N point length respectively, second section N point length digital signal enters the first fast Fourier module (FFT) 3, the digital signal of first section N point length of the first fast Fourier module (FFT), 3 outputs, the highest frequency component phase angle 8 of second section N point length digital signal are successively sent into phase detector 4, and phase detector 4 output calibration frequencies 10 are sent into digital signal reconfigurable filter 5.
3, digital signal reconfigurable filter 5 carries out the digital reconstruction sampling according to the digital output signal 12 of first section N point length of proofreading and correct 10 pairs of fixed sample rate analog to digital converters 2 of back frequency, and the N point reconstructed sample signal that obtains is sent into the second fast Fourier module 6.
4, phase detector 4 is sent in the highest frequency component phase angle of the second fast Fourier module 6 output 9, the digital signal of the 3rd section N point of analog to digital converter 2 length enters the first fast Fourier module 3, and output calibration frequency 10 was sent into digital signal reconfigurable filter 5 after phase detector 4 was judged according to accuracy requirement.
5, repeating step 2,3 and 4, reach 2.7e-12 through frequency-tracking after 3 iteration and measuring error, and each harmonic is measured maximum error 1e-15.
Claims (4)
1. an electric power system alternating current sampling real-time frequency is followed the tracks of and harmonic measuring method, and it is characterized in that: it is as follows that the method comprising the steps of:
(1) gets AC sampling input 1 and enter fixed sample rate analog to digital converter 2.
(2) fixed sample rate analog to digital converter 2 is exported the digital signal of first section N point length respectively, second section N point length digital signal enters the first fast Fourier module (FFT) 3, the digital signal of first section N point length of the first fast Fourier module (FFT), 3 outputs, the highest frequency component phase angle 8 of second section N point length digital signal are successively sent into phase detector 4, and phase detector 4 output calibration frequencies 10 are sent into digital signal reconfigurable filter 5.
(3) digital signal reconfigurable filter 5 carries out the digital reconstruction sampling according to the digital output signal 12 of first section N point length of proofreading and correct 10 pairs of fixed sample rate analog to digital converters 2 of back frequency, and the N point reconstructed sample signal that obtains is sent into the second fast Fourier module 6.
Phase detector 4 is sent in the highest frequency component phase angle output 9 of (4) second fast Fourier modules 6, the digital signal of the 3rd section N point of analog to digital converter 2 length enters the first fast Fourier module 3, and output calibration frequency 10 was sent into digital signal reconfigurable filter 5 after phase detector 4 was judged according to accuracy requirement.
(5) repeating step 2,3 and 4, reach predetermined iterations M after, 6 outputs 7 of the second fast Fourier module are the each harmonic component value of ac sampling signal, phase detector 4 output frequencies 10 are the frequency of ac sampling signal.
2. method according to claim 1 is characterized in that: predetermined iterations M is a positive integer, generally gets the M value and be 3 to 6 proper.
3. method according to claim 1 is characterized in that: the highest frequency component phase angle of the first fast Fourier module (FFT) 3 and the second fast Fourier module 6 is sent into phase detector 4, utilizes the difference of the phase angle of the two to obtain the tracking frequencies deviation.
4. method according to claim 1 is characterized in that: digital signal reconfigurable filter 5 is a sampling rate with emending frequency 10, and the digital signal of fixed sample rate analog to digital converter 2 outputs is carried out the variable Rate sampling.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100464591A CN102135567A (en) | 2011-02-25 | 2011-02-25 | Real-time frequency tracking and harmonic measuring method for AC sampling of power system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100464591A CN102135567A (en) | 2011-02-25 | 2011-02-25 | Real-time frequency tracking and harmonic measuring method for AC sampling of power system |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102135567A true CN102135567A (en) | 2011-07-27 |
Family
ID=44295409
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011100464591A Pending CN102135567A (en) | 2011-02-25 | 2011-02-25 | Real-time frequency tracking and harmonic measuring method for AC sampling of power system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102135567A (en) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103063913A (en) * | 2012-12-07 | 2013-04-24 | 深圳市金宏威技术股份有限公司 | Frequency tracking method for Fourier transform |
CN103728523A (en) * | 2014-01-21 | 2014-04-16 | 国家电网公司 | Urban power grid yield measuring method |
CN103941086A (en) * | 2014-03-25 | 2014-07-23 | 中国科学院长春光学精密机械与物理研究所 | Ultrahigh precision frequency measurement instrument and measuring method thereof |
CN105004920A (en) * | 2015-07-10 | 2015-10-28 | 国网天津市电力公司 | Fourier correction coefficient frequency measuring method |
CN105004913A (en) * | 2015-07-21 | 2015-10-28 | 黎文安 | Real-time tracking method for alternating current sampling frequency of electric power system |
CN105044459A (en) * | 2015-07-21 | 2015-11-11 | 青岛艾诺智能仪器有限公司 | Harmonic analysis method |
CN106597095A (en) * | 2016-12-22 | 2017-04-26 | 南京因泰莱电器股份有限公司 | Double-CPU uniformly-spaced sampling frequency-tracking implementation method |
EP3287795A1 (en) * | 2016-08-24 | 2018-02-28 | Schneider Electric Industries SAS | Method for determining the frequency of an ac signal |
CN108037357A (en) * | 2017-11-29 | 2018-05-15 | 西电通用电气自动化有限公司 | A kind of frequency tracking method of Automation of Electric Systems device |
CN109030957A (en) * | 2015-05-19 | 2018-12-18 | 江苏理工学院 | Dielectric loss measurement method |
CN109030942A (en) * | 2015-05-19 | 2018-12-18 | 江苏理工学院 | Humorous phase angle analysis method |
CN110187176A (en) * | 2019-05-31 | 2019-08-30 | 中国电力科学研究院有限公司 | A kind of mains by harmonics measuring device and method |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101221201A (en) * | 2008-01-25 | 2008-07-16 | 杭州万工科技有限公司 | Method for measuring fundamental frequency |
US20100069030A1 (en) * | 2006-02-10 | 2010-03-18 | Thierry Briand | Frequency measuring broadband digital receiver |
CN101852826A (en) * | 2009-03-30 | 2010-10-06 | 西门子公司 | Harmonic analysis method for power system and device thereof |
-
2011
- 2011-02-25 CN CN2011100464591A patent/CN102135567A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100069030A1 (en) * | 2006-02-10 | 2010-03-18 | Thierry Briand | Frequency measuring broadband digital receiver |
CN101221201A (en) * | 2008-01-25 | 2008-07-16 | 杭州万工科技有限公司 | Method for measuring fundamental frequency |
CN101852826A (en) * | 2009-03-30 | 2010-10-06 | 西门子公司 | Harmonic analysis method for power system and device thereof |
Non-Patent Citations (3)
Title |
---|
李振然等: "利用富里叶变换实现电力***频率的测量", 《广西大学学报》 * |
牟龙华等: "基于傅里叶变换的精确频率测量算法", 《电力***自动化》 * |
王柏林等: "电网谐波分析的频率估计与跟踪", 《电测与仪表》 * |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103063913A (en) * | 2012-12-07 | 2013-04-24 | 深圳市金宏威技术股份有限公司 | Frequency tracking method for Fourier transform |
CN103063913B (en) * | 2012-12-07 | 2016-01-20 | 深圳市金宏威技术有限责任公司 | For the frequency tracking method of Fourier transform |
CN103728523A (en) * | 2014-01-21 | 2014-04-16 | 国家电网公司 | Urban power grid yield measuring method |
CN103941086A (en) * | 2014-03-25 | 2014-07-23 | 中国科学院长春光学精密机械与物理研究所 | Ultrahigh precision frequency measurement instrument and measuring method thereof |
CN109030957A (en) * | 2015-05-19 | 2018-12-18 | 江苏理工学院 | Dielectric loss measurement method |
CN109030957B (en) * | 2015-05-19 | 2020-12-25 | 江苏理工学院 | Dielectric loss measuring method |
CN109030942B (en) * | 2015-05-19 | 2020-12-04 | 江苏理工学院 | Harmonic phase angle analysis method |
CN109030942A (en) * | 2015-05-19 | 2018-12-18 | 江苏理工学院 | Humorous phase angle analysis method |
CN105004920A (en) * | 2015-07-10 | 2015-10-28 | 国网天津市电力公司 | Fourier correction coefficient frequency measuring method |
CN105004920B (en) * | 2015-07-10 | 2017-11-17 | 国网天津市电力公司 | Fourier's correction factor frequency measurement method |
CN105044459A (en) * | 2015-07-21 | 2015-11-11 | 青岛艾诺智能仪器有限公司 | Harmonic analysis method |
CN105044459B (en) * | 2015-07-21 | 2017-09-29 | 青岛艾诺智能仪器有限公司 | A kind of harmonic analysis method |
CN105004913A (en) * | 2015-07-21 | 2015-10-28 | 黎文安 | Real-time tracking method for alternating current sampling frequency of electric power system |
FR3055417A1 (en) * | 2016-08-24 | 2018-03-02 | Schneider Electric Industries Sas | DETERMINING THE FREQUENCY OF AN ALTERNATIVE SIGNAL |
CN107782966A (en) * | 2016-08-24 | 2018-03-09 | 施耐德电器工业公司 | Determine the frequency of AC signal |
EP3287795A1 (en) * | 2016-08-24 | 2018-02-28 | Schneider Electric Industries SAS | Method for determining the frequency of an ac signal |
US10613127B2 (en) | 2016-08-24 | 2020-04-07 | Schneider Electric Industries Sas | Determining the frequency of an alternating signal |
CN106597095B (en) * | 2016-12-22 | 2019-09-10 | 南京因泰莱电器股份有限公司 | A kind of frequency-tracking implementation method of dual processors equal interval sampling |
CN106597095A (en) * | 2016-12-22 | 2017-04-26 | 南京因泰莱电器股份有限公司 | Double-CPU uniformly-spaced sampling frequency-tracking implementation method |
CN108037357A (en) * | 2017-11-29 | 2018-05-15 | 西电通用电气自动化有限公司 | A kind of frequency tracking method of Automation of Electric Systems device |
CN110187176A (en) * | 2019-05-31 | 2019-08-30 | 中国电力科学研究院有限公司 | A kind of mains by harmonics measuring device and method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102135567A (en) | Real-time frequency tracking and harmonic measuring method for AC sampling of power system | |
de la O Serna et al. | Instantaneous Oscillating Phasor Estimates With Taylor $^ K $-Kalman Filters | |
CN103869162B (en) | Dynamic signal phasor measurement method based on time domain quasi-synchronization | |
CN103257271A (en) | Device and method for detecting micro grid harmonic wave and inter-harmonics based on STM32F107VCT6 | |
CN203287435U (en) | A micro electrical network harmonic wave and inter-harmonic wave test apparatus based on an STM32F107VCT6 | |
CN103543335B (en) | A kind of synchronous phasor measuring method | |
CN101915874A (en) | Harmonic wave detection method based on Fourier transformation | |
CN103941088A (en) | Method for quickly measuring frequency of electric power system based on three-phase signals | |
CN203909138U (en) | General three-phase AC signal synchronous sampling device | |
CN107478896A (en) | A kind of frequency adaptive harmonic current detection method based on cascade Generalized Integrator | |
CN109142863B (en) | Power system frequency measurement method and system | |
CN102236048A (en) | Method for measuring phasor frequency of electric system | |
CN102035554A (en) | System and method for metering and analyzing electric energy as well as analog-to-digital conversion circuit | |
CN103018555A (en) | High-precision electric power parameter software synchronous sampling method | |
CN102243272A (en) | High-precision harmonic analysis method of sampling data synchronization | |
CN105629060A (en) | Grid frequency measurement method and device based on optimal baseband filtering | |
CN104330623A (en) | Sine wave parameter measuring method and system in electric power system | |
CN104155517A (en) | Error compensation method and system for non-full period sampling of digital electric energy meter | |
CN102495250A (en) | Quasi-synchronous wideband reactive energy meter based on Hilbert algorithm and sampling method thereof | |
CN102809687A (en) | Digital measurement method for alternating-current frequency | |
CN103969508A (en) | Real-time high-precision power harmonic analysis method and device | |
CN114966194A (en) | ADC-based low-frequency metering three-phase electric energy meter | |
CN105044459A (en) | Harmonic analysis method | |
CN105004913A (en) | Real-time tracking method for alternating current sampling frequency of electric power system | |
CN105182072A (en) | Harmonic detection method and system based on virtual filter coefficients |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20110727 |