CN103983849A - Real-time high-accuracy power harmonic analysis method - Google Patents
Real-time high-accuracy power harmonic analysis method Download PDFInfo
- Publication number
- CN103983849A CN103983849A CN201410189931.0A CN201410189931A CN103983849A CN 103983849 A CN103983849 A CN 103983849A CN 201410189931 A CN201410189931 A CN 201410189931A CN 103983849 A CN103983849 A CN 103983849A
- Authority
- CN
- China
- Prior art keywords
- frequency
- real
- analysis method
- signal
- time high
- 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
Landscapes
- Measurement Of Resistance Or Impedance (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention provides a real-time high-accuracy power harmonic analysis method in the technical field of power harmonic detection. The real-time high-accuracy power harmonic analysis method includes the following steps: a power grid signal of an electrical power system is sampled with a fixed sampling frequency to obtain a first sampling signal; a fundamental frequency of the power grid signal is obtained according output in the mode that the first sampling signal passes through a digital filter bank; a variable sample frequency directly proportional with the fundamental frequency is obtained through the fundamental frequency, a to-be-detected power grid signal is sampled with the variable sample frequency to obtain a second sampling signal, analog-digital conversion is carried out on the second sampling signal, and the second sampling signal is input into a digital signal processor; Fourier transformation is carried out on the second sampling signal in an FFT mode, and amplitudes and phase angles of needed harmonic components are calculated. By means of the real-time high-accuracy power harmonic analysis method, the measurement accuracy on harmonious waves of the electrical power system is higher, and the requirement of the electrical power system for the instantaneity is met.
Description
Technical field
The present invention relates to the frequency analysis field in electric system, relate in particular to a kind of Electric Power Harmonic Analysis method of real-time high-precision.
Background technology
Along with the development of electronic technology and power electronic devices, the application of power electronic equipment is more and more extensive, and these non-linear power equipments make to have produced a large amount of harmonic waves in electrical network, causes the quality of power supply to decline.The accurate real-time analysis of harmonic wave is of great significance the improvement tool of the quality of power supply.
At present, for the analysis of harmonic wave, often adopt frequency domain analysis, to mainly contain FFT be Fourier transform, based on neural network theory, vector and wavelet transformation etc.Fourier transform is comparative maturity and practical harmonic analysis method, becomes the main algorithm of current frequency analysis with simple and practical feature.
Fft algorithm, the in the situation that of non-synchronous sampling, can produce spectral leakage, causes the inaccurate of measurement data; At present generally adopt the method for windowing to reduce spectrum leakage, but this method can increase the calculated amount of algorithm, affects frequency analysis speed.
So the frequency that need to provide a kind of real-time is high, calculated amount is little method to follow the tracks of input signal, to realize the synchronized sampling of signal, eliminates frequency synchronization error.
Summary of the invention
Fundamental purpose of the present invention is to provide a kind of Electric Power Harmonic Analysis method of real-time high-precision, can obtain fast the fundamental frequency of electric system by the method, improves the precision of frequency analysis.
For achieving the above object, the invention provides a kind of Electric Power Harmonic Analysis method of real-time high-precision, the method includes the steps of:
(1) with a fixed sampling frequency f
selectric system power network signal is sampled and obtained the first sampled signal;
(2) the first sampled signal is passed through to digital filter bank, according to the output of digital filter bank, obtain the fundamental frequency f of power network signal
0;
(3) by fundamental frequency f
0the variable sampling frequency f that obtains being proportional
x, utilize described variable sampling frequency f
xpower network signal to be measured is sampled and obtained the second sampled signal, by the second sampled signal supplied with digital signal processor after analog to digital conversion;
(4) take FFT mode to carry out Fourier transform to the second sampled signal, calculate amplitude and the phase angle of required harmonic component.
The fixed sampling frequency f that samples adopted to electric system power network signal in described step (1)
sfor fundamental frequency f
02 ~ 4 times.
In described step (2), need the output valve of comparative figures bank of filters, using centre frequency corresponding to the wave filter of output valve maximum as fundamental frequency f
0.
The digital filter bank of described step (2) is the logical FIR bank of filters of band.
The centre frequency of the logical FIR bank of filters of described band covers work frequency.
Described step (3) is utilized described variable sampling frequency f
xto the sampling time interval T of power network signal to be measured, be fundamental frequency f
0with variable sampling frequency f
xthe inverse of the product of one-period sampling number N, i.e. T=1/(f
0* N).
The invention has the beneficial effects as follows: can obtain rapidly and accurately the frequency of power frequency component, realize synchronized sampling, solve the problem of non-synchronous sampling generation error, meet electric system and detect the requirement to real-time.
Accompanying drawing explanation
Fig. 1 is that the Electric Power Harmonic Analysis method of real-time high-precision realizes block diagram.
Fig. 2 is FIR filter graph architecture.
Embodiment
The present invention proposes a kind of Electric Power Harmonic Analysis method of real-time high-precision, is characterized in obtaining in real time fast base wave frequency, realizes synchronized sampling.
The sampled point that frequency analysis is used need to be with power network signal with frequently, and first the present invention carries out frequency measurement by a road power network signal, and then another road power network signal is started to sampling in cycle starting point.
Measure the frequency of power network signal, first according to required measuring accuracy, set the logical FIR wave filter quantity of band that forms the logical FIR bank of filters of band, then adopt as shown in Figure 1 a fixed sampling frequency f
sanalog to digital converter ADCDui mono-road power network signal sample and obtain the first sampled signal; The first sampled signal is led to FIR bank of filters by band, the output valve of more every logical FIR wave filter of band, wherein the logical corresponding centre frequency of FIR wave filter of the band of output valve maximum is fundamental frequency f
0, for taking abundant counting, to guarantee certain frequency resolution, the fixed sampling frequency f that adopts
sfor fundamental frequency f
02 ~ 4 times; By fundamental frequency f
0by clock circuit, obtain a fundamental frequency f
0the clock signal of integral multiple M is as variable sampling frequency f
x, according to test power network signal overtone order, determine the value of M, take and test 30 subharmonic as example, M is greater than 60; The sampling time interval T of this clock signal is fundamental frequency f
0with variable sampling frequency f
xthe inverse of the product of one-period sampling number N, i.e. T=1/(f
0* N); Sample frequency using this clock signal as another variable frequency analog to digital converter ADC again, another road power network signal is sampled and obtained the second sampled signal, by the second sampled signal supplied with digital signal processor, adopt FFT mode to carry out Fourier transform to the second sampled signal, calculate amplitude and the phase angle of required harmonic component.
Claims (6)
1. an Electric Power Harmonic Analysis method for real-time high-precision, is characterized in that, the method includes the steps of:
(1) with a fixed sampling frequency f
selectric system power network signal is sampled and obtained the first sampled signal;
(2) the first sampled signal is passed through to digital filter bank, according to the output of digital filter bank, obtain the fundamental frequency f of power network signal
0;
(3) by fundamental frequency f
0the variable sampling frequency f that obtains being proportional
x, utilize described variable sampling frequency f
xpower network signal to be measured is sampled and obtained the second sampled signal, by the second sampled signal supplied with digital signal processor after analog to digital conversion;
(4) take FFT mode to carry out Fourier transform to the second sampled signal, calculate amplitude and the phase angle of required harmonic component.
2. the Electric Power Harmonic Analysis method of a kind of real-time high-precision according to claim 1, is characterized in that, the fixed sampling frequency f that samples adopted to electric system power network signal in described step (1)
sfor fundamental frequency f
02 ~ 4 times.
3. the Electric Power Harmonic Analysis method of a kind of real-time high-precision according to claim 1, is characterized in that, needs the output valve of comparative figures bank of filters in described step (2), using centre frequency corresponding to the wave filter of output valve maximum as fundamental frequency f
0.
4. according to the Electric Power Harmonic Analysis method of a kind of real-time high-precision described in claim 1 or 3, it is characterized in that, the digital filter bank of described step (2) is the logical FIR bank of filters of band.
5. the Electric Power Harmonic Analysis method of a kind of real-time high-precision according to claim 4, is characterized in that, the centre frequency of the logical FIR bank of filters of described band covers work frequency.
6. the Electric Power Harmonic Analysis method of a kind of real-time high-precision according to claim 1, is characterized in that, described step (3) is utilized described variable sampling frequency f
xto the sampling time interval T of power network signal to be measured, be fundamental frequency f
0with variable sampling frequency f
xthe inverse of the product of one-period sampling number N, i.e. T=1/(f
0* N).
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410189931.0A CN103983849B (en) | 2014-05-07 | 2014-05-07 | A kind of Electric Power Harmonic Analysis method of real-time high-precision |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410189931.0A CN103983849B (en) | 2014-05-07 | 2014-05-07 | A kind of Electric Power Harmonic Analysis method of real-time high-precision |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103983849A true CN103983849A (en) | 2014-08-13 |
| CN103983849B CN103983849B (en) | 2016-06-29 |
Family
ID=51275895
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410189931.0A Active CN103983849B (en) | 2014-05-07 | 2014-05-07 | A kind of Electric Power Harmonic Analysis method of real-time high-precision |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103983849B (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105116202A (en) * | 2015-07-14 | 2015-12-02 | 云南电网有限责任公司玉溪供电局 | Method of rapidly obtaining electric signal fundamental wave amplitude and multiple harmonic amplitudes |
| CN105973349A (en) * | 2016-04-19 | 2016-09-28 | 温州科技职业学院 | On-line prediction method for inherent frequency of Coriolis flowmeter |
| CN106372800A (en) * | 2016-08-31 | 2017-02-01 | 国网北京市电力公司 | Method and device for acquiring subway power supply quality |
| CN106918741A (en) * | 2017-03-02 | 2017-07-04 | 浙江大学 | It is applied to the adaptively sampled phase difference correction method of frequency wide swings power network |
| CN111257637A (en) * | 2020-02-21 | 2020-06-09 | 深圳市中电电力技术股份有限公司 | Ultra-high harmonic measurement method and system |
| CN117169590A (en) * | 2023-08-15 | 2023-12-05 | 嘉兴市科讯电子有限公司 | Power harmonic analysis method and device based on software variable sampling rate |
| CN117330834A (en) * | 2023-10-09 | 2024-01-02 | 国网山东省电力公司东营市东营区供电公司 | Dual-sampling-based power grid harmonic content testing method and system |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4918381A (en) * | 1989-04-06 | 1990-04-17 | Hewlett-Packard Company | Automated method for determining total harmonic distortion |
| CN101368987A (en) * | 2008-09-27 | 2009-02-18 | 武汉大学 | Power system mesomeric state harmonic wave and/or inter-harmonic wave measuring method |
| CN101915874A (en) * | 2010-07-20 | 2010-12-15 | 北海市深蓝科技发展有限责任公司 | Harmonic wave detection method based on Fourier transformation |
| CN103105532A (en) * | 2013-01-22 | 2013-05-15 | 广西电网公司电力科学研究院 | Harmonic electric energy measuring system of improved Fourier transform and control method thereof |
-
2014
- 2014-05-07 CN CN201410189931.0A patent/CN103983849B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4918381A (en) * | 1989-04-06 | 1990-04-17 | Hewlett-Packard Company | Automated method for determining total harmonic distortion |
| CN101368987A (en) * | 2008-09-27 | 2009-02-18 | 武汉大学 | Power system mesomeric state harmonic wave and/or inter-harmonic wave measuring method |
| CN101915874A (en) * | 2010-07-20 | 2010-12-15 | 北海市深蓝科技发展有限责任公司 | Harmonic wave detection method based on Fourier transformation |
| CN103105532A (en) * | 2013-01-22 | 2013-05-15 | 广西电网公司电力科学研究院 | Harmonic electric energy measuring system of improved Fourier transform and control method thereof |
Non-Patent Citations (1)
| Title |
|---|
| 杨富康: "基于FFT的电网谐波检测方法的研究", 《中国优秀硕士论文全文数据库》 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105116202A (en) * | 2015-07-14 | 2015-12-02 | 云南电网有限责任公司玉溪供电局 | Method of rapidly obtaining electric signal fundamental wave amplitude and multiple harmonic amplitudes |
| CN105973349A (en) * | 2016-04-19 | 2016-09-28 | 温州科技职业学院 | On-line prediction method for inherent frequency of Coriolis flowmeter |
| CN105973349B (en) * | 2016-04-19 | 2019-04-09 | 温州科技职业学院 | A kind of on-line prediction method of Coriolis flowmeter intrinsic frequency |
| CN106372800A (en) * | 2016-08-31 | 2017-02-01 | 国网北京市电力公司 | Method and device for acquiring subway power supply quality |
| CN106918741A (en) * | 2017-03-02 | 2017-07-04 | 浙江大学 | It is applied to the adaptively sampled phase difference correction method of frequency wide swings power network |
| CN106918741B (en) * | 2017-03-02 | 2019-04-23 | 浙江大学 | Adaptively sampled phase difference correction method applied to frequency wide swings power grid |
| CN111257637A (en) * | 2020-02-21 | 2020-06-09 | 深圳市中电电力技术股份有限公司 | Ultra-high harmonic measurement method and system |
| CN111257637B (en) * | 2020-02-21 | 2022-05-20 | 深圳市中电电力技术股份有限公司 | Ultra-high harmonic measurement method and system |
| CN117169590A (en) * | 2023-08-15 | 2023-12-05 | 嘉兴市科讯电子有限公司 | Power harmonic analysis method and device based on software variable sampling rate |
| CN117169590B (en) * | 2023-08-15 | 2024-03-29 | 嘉兴市科讯电子有限公司 | Power harmonic analysis method and device based on software variable sampling rate |
| CN117330834A (en) * | 2023-10-09 | 2024-01-02 | 国网山东省电力公司东营市东营区供电公司 | Dual-sampling-based power grid harmonic content testing method and system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103983849B (en) | 2016-06-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103983849B (en) | A kind of Electric Power Harmonic Analysis method of real-time high-precision | |
| CN102288807B (en) | Method for measuring electric network voltage flicker | |
| CN104007316B (en) | A kind of High Precision Frequency method under lack sampling speed and measuring instrument thereof | |
| CN103941087B (en) | The frequency measurement method of the high-frequency cosine signal under lack sampling speed and device thereof | |
| CN203909138U (en) | General three-phase AC signal synchronous sampling device | |
| CN102435844A (en) | Sinusoidal signal phasor calculating method being independent of frequency | |
| CN109633262A (en) | Three phase harmonic electric energy gauging method, device based on composite window multiline FFT | |
| CN101915874A (en) | Harmonic wave detection method based on Fourier transformation | |
| CN103869162A (en) | Dynamic signal phasor measurement method based on time domain quasi-synchronization | |
| CN103457603B (en) | A kind of method that ADC dynamic parameter is tested based on average frequency spectrum | |
| CN103941091A (en) | Power system HHT harmonious wave detection method based on improved EMD end point effect | |
| CN103969508A (en) | Real-time high-precision power harmonic analysis method and device | |
| CN102645585A (en) | Synchronous rapid measuring method and device for ultralow frequency sinusoidal signal phase difference | |
| CN203287435U (en) | A micro electrical network harmonic wave and inter-harmonic wave test apparatus based on an STM32F107VCT6 | |
| CN105487034A (en) | 0.05-level electronic transformer verification method and system | |
| CN113064021B (en) | Measurement and control device and method for realizing power electronic power grid higher harmonic suppression | |
| CN103543331B (en) | A kind of method calculating electric signal harmonic wave and m-Acetyl chlorophosphonazo | |
| CN110967658A (en) | Analog input merging unit calibrator tracing method based on digital differential method | |
| CN101937017A (en) | Dynamic direct-current removing method for intelligent electric meter during alternating-current sampling | |
| CN109030957B (en) | Dielectric loss measuring method | |
| RU2435168C1 (en) | Method for harmonic analysis of periodic multifrequency signal | |
| CN112180161A (en) | Harmonic inter-harmonic wave group measuring method under asynchronous high sampling rate sampling condition | |
| CN103592513B (en) | Electric power signal harmonic analysis method and device | |
| CN106872777B (en) | Harmonic and inter-harmonic separation analysis method | |
| CN103267896B (en) | The measuring method of initial phase angle of periodic signal |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |