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 PDF

Info

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
Application number
CN2011100464591A
Other languages
Chinese (zh)
Inventor
樊体峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
SHENZHEN YESTUNED INDUSTRY Co Ltd
Original Assignee
SHENZHEN YESTUNED INDUSTRY Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by SHENZHEN YESTUNED INDUSTRY Co Ltd filed Critical SHENZHEN YESTUNED INDUSTRY Co Ltd
Priority to CN2011100464591A priority Critical patent/CN102135567A/en
Publication of CN102135567A publication Critical patent/CN102135567A/en
Pending legal-status Critical Current

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

A kind of electric power system alternating current sampling real-time frequency is followed the tracks of and harmonic measuring method
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:
θ m ( k ) = tan - 1 { Im [ X m ( k ) ] Re [ X m ( k ) ] } - - - ( 1 - 1 )
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:
θ n ( k ) = tan - 1 { Im [ X n ( k ) ] Re [ X n ( k ) ] } - - - ( 1 - 2 )
Utilize these two phasing degree can obtain the tracking frequencies deviation
f Δ = θ n ( k ) - θ m ( k ) 2 π ( n - m ) - - - ( 1 - 3 )
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:
x ( t ) = sin ( 2 π f 0 t ) + 1 3 sin ( 3 × 2 π f 0 t ) + 1 5 sin ( 5 × 2 π f 0 t ) + 1 7 sin ( 7 × 2 π f 0 t )
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.
CN2011100464591A 2011-02-25 2011-02-25 Real-time frequency tracking and harmonic measuring method for AC sampling of power system Pending CN102135567A (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (3)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
李振然等: "利用富里叶变换实现电力***频率的测量", 《广西大学学报》 *
牟龙华等: "基于傅里叶变换的精确频率测量算法", 《电力***自动化》 *
王柏林等: "电网谐波分析的频率估计与跟踪", 《电测与仪表》 *

Cited By (21)

* Cited by examiner, † Cited by third party
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