CN106814230A - Electric energy gauging method based on the peak interpolations of Kaiser windows FFT tetra- - Google Patents
Electric energy gauging method based on the peak interpolations of Kaiser windows FFT tetra- Download PDFInfo
- Publication number
- CN106814230A CN106814230A CN201611157851.2A CN201611157851A CN106814230A CN 106814230 A CN106814230 A CN 106814230A CN 201611157851 A CN201611157851 A CN 201611157851A CN 106814230 A CN106814230 A CN 106814230A
- Authority
- CN
- China
- Prior art keywords
- peak
- amplitude
- frequency
- formula
- electric energy
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R11/00—Electromechanical arrangements for measuring time integral of electric power or current, e.g. of consumption
- G01R11/48—Meters specially adapted for measuring real or reactive components; Meters specially adapted for measuring apparent energy
- G01R11/54—Meters specially adapted for measuring real or reactive components; Meters specially adapted for measuring apparent energy for measuring simultaneously at least two of the following three variables: real component, reactive component, apparent energy
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
Abstract
The invention discloses a kind of electric energy gauging method based on the peak interpolations of Kaiser windows FFT tetra-, including：Amplitude of the frequency spectrum based on the peak interpolation correction formulas of Kaiser windows FFT tetra-, the frequency displacement of FFT and window function Kaiser in fundamental wave, harmonic wave and the symmetrical spectral line of left and right four of m-Acetyl chlorophosphonazo crest frequency point, using fitting polynomial formulas, general amplitude, the phase and frequency interpolation correction formula of signal fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo are obtained；Amplitude, the phase and frequency parameter of tested voltage signal, the fundamental wave of current signal, harmonic wave and m-Acetyl chlorophosphonazo are calculated according to general amplitude, phase and frequency interpolation correction formula；And then calculate fundamental energy, harmonic electric energy, m-Acetyl chlorophosphonazo electric energy and total electric energy.The invention has the advantages that：Four peak interpolation correction algorithms, take full advantage of the symmetry of four spectral lines of amplitude Frequency point or so and excellent master, the side lobe performance of Kaiser window functions, with accuracy in computation higher.
Description
Technical field
The present invention relates to electric power system electric energy metrical field and digitalized electric energy metering field, and in particular to one kind is based on
The electric energy gauging method of the peak interpolations of Kaiser windows FFT tetra-.
Background technology
With increasing electric railway, rolling mill, electric arc furnaces, electric automobile charging pile, elevator, intelligent family
Coming into operation for the Nonlinear Dynamic electric load such as electricity, elevator, generates a large amount of electric harmonics and m-Acetyl chlorophosphonazo in power supply grid.
The interference components such as a large amount of electric harmonic, m-Acetyl chlorophosphonazo, dc-decay, Voltage Drops, cause power supply grid electric current, voltage waveform
Distortion.The presence of these interference, not only has a strong impact on power quality, jeopardizes the safe and stable operation of electric power system, or even
Cause serious electric power accident, cause huge economic loss, and influence the accuracy of electric energy metrical.
Electric energy is the foundation that grid company, sale of electricity company etc. carry out business accounting, and can it be measured correctly, exactly,
It is directly connected to the economic benefit of power supply and demand both sides.Traditional way of energy measuring, be directly by electric energy define, by voltage,
The instantaneous sampling value of current signal is multiplied, carries out time integral again to calculate the so-called full way of energy measuring of electric energy.And wind-force
Generator and photovoltaic plant generate electricity by way of merging two or more grid systems, and the Nonlinear Dynamic electric power such as electric arc furnaces, rolling mill, electric locomotive, charging pile
The a large amount of of load put into operation, and in power supply grid, serving not only as load will consume power, and as harmonic source, pollution sources, cause
While supply current, voltage waveform distortion, a large amount of harmonic waves and m-Acetyl chlorophosphonazo power etc. can be injected to power supply grid.If by traditional
Full way of energy measuring, although Harmonics source customer can inject a large amount of harmonic waves to power network disturbing, also pays tariff less；And it is vast general
Logical electricity consumer not only by harmonic wave interference and harm, and more than pay tariff；Meanwhile, this full way of energy measuring does not have yet
There is consideration because of spectral leakage and fence effect caused by non-synchronous sampling, therefore use this way of energy measuring, undoubtedly may
Bring larger error.
And another so-called fundamental energy metering method, although the electricity consumption fairness problem of conventional power user can be solved,
Making conventional power user be not necessarily harmonic electric energy to pay tariff more, but Harmonics source customer still only needs to pay the fundamental energy electricity charge, and its
The electric harmonic pollution injected to power supply grid is but not affected by any punishment.
It can be seen that, whether using full way of energy measuring or using fundamental energy metering method, how fail solution
To Harmonics source customer to the inaccurate problem of the caused electric energy metrical such as power supply grid harmonic, m-Acetyl chlorophosphonazo, let alone such as
What is disturbed by punitive measures reduced harmonic source user as far as possible to electric power system harmonic etc..
Therefore, electric power system electric current, voltage waveform distortion are caused for specific, exemplary dynamic nonlinear-load, how more
The parameters such as amplitude, the phase and frequency of fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo are extracted exactly, are before one of high accuracy electric energy metrical
Carry.Conventional method is to apply the computings such as window function, Fourier transformation, interpolation amendment to measured signal.Specifically to tested
Before signal carries out FFT, apply suitable window function, let out with suppressing the frequency spectrum caused by non-synchronous sampling and non-integer-period are blocked
Leakage；Again by interpolation amendment, error caused by compensation fence effect.Paper " calculate by the unimodal spectral line interpolation based on approximation by polynomi-als
Application of the method in m-Acetyl chlorophosphonazo analysis " (Xiao is first brave to wait electric power network techniques, 2008), gives based on Hanning windows, Hamming
Window, Blackman-Harris windows, the unimodal spectral line interpolation correction formula of Rife-Vincent windows；Paper " carries out electricity using FFT
The innovatory algorithm of Force system frequency analysis " (the Proceedings of the CSEEs such as Pang Hao, 2003), give based on rectangular window,
Triangle windows, Hanning windows, Hamming windows, the double peak interpolation correction formula of Blackman windows；Paper " is inserted based on three spectral lines
The Electric Power Harmonic Analysis algorithm of value FFT " (the Proceedings of the CSEEs, 2012 such as ox victory lock), paper " are based on five Lai Fu-texts
Three spectral line interpolation FFT frequency analyses of cent window " (the electrical technologies such as Zhu Junwei, 2015), paper " based on Nuttall windows three
Peak interpolation harmonic wave Algorithm Analysis " (the electric power system protection and controls such as Zhai Ruimiao, 2015), sets forth based on Hanning
Window, Hamming windows, Blackman windows, Blackman-Harris windows, Rife-Vincent windows, three peak spectral lines of Nuttall windows
Interpolation correction formula.
But the problem that these existing algorithms are still suffered from is：1) because the spectral performance of window function is changeless, therefore
When carrying out windowed interpolation correction algorithm to Dynamic Signal, computational accuracy by window function fix side lobe performance limited, cause by
The harmonic components surveyed in signal cannot be extracted accurately due to what is influenceed by spectral leakage；2) unimodal spectral line interpolation correction formula
It is more complicated, easily influenceed by spectrum leakage, double peak of spectral lines algorithm is to containing spectral line information using insufficient, three peaks near peak value frequency
Spectral line algorithm does not consider the comprehensive information of the symmetrical spectral line of peak value frequency, is difficult to meet the requirement of high accuracy electric energy metrical；
3) in existing interpolation correction algorithm, the poor problem of generally existing even-order harmonic electric energy accuracy in computation.
In view of this, it is necessary to a kind of electric energy gauging method based on the peak interpolations of Kaiser windows FFT tetra- is provided, to meet reality
Border application needs.
The content of the invention
It is contemplated that at least solving one of above-mentioned technical problem.
Therefore, it is an object of the invention to propose a kind of electric energy gauging method based on the peak interpolations of Kaiser windows FFT tetra-, filling
Divide excellent master, the side lobe performance of the symmetry and Kaiser window functions using four spectral lines of crest frequency point or so, to reach
The requirement of high accuracy electric energy metrical.
To achieve these goals, the technical solution adopted in the present invention is：One kind is based on the peak interpolations of Kaiser windows FFT tetra-
Electric energy gauging method, it is characterised in that comprise the following steps：
S1：Frequency spectrum based on the peak interpolation correction formulas of Kaiser windows FFT tetra-, the frequency displacement of FFT and window function Kaiser exists
The amplitude of fundamental wave, harmonic wave and the symmetrical spectral line of m-Acetyl chlorophosphonazo crest frequency point arranged on left and right sides four, using fitting of a polynomial, obtains letter
General amplitude, the phase and frequency interpolation correction formula of number fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo；
S2：The base of tested voltage signal, current signal is calculated according to the general amplitude, phase and frequency correction formula
The amplitude of ripple, harmonic wave and m-Acetyl chlorophosphonazo, phase and frequency parameter；
S3：According to amplitude, the phase and frequency of the tested voltage signal, the fundamental wave of current signal, harmonic wave and m-Acetyl chlorophosphonazo
Parameter calculates fundamental energy, harmonic electric energy, m-Acetyl chlorophosphonazo electric energy and total electric energy respectively.
Electric energy gauging method based on the peak interpolation amendments of Kaiser windows FFT tetra- according to embodiments of the present invention, Kaiser windows
Function is one group of adjustable window function, and the proportion of main lobe width and side lobe height, β are adjusted by unrestricted choice form parameter β
Value is bigger, and the main lobe width of window function frequency spectrum is wider, and side lobe peak is smaller.This patent uses the Kaiser of form parameter β=22
Window function, its main lobe width is 28 π/N, and side lobe peak is -172dB, side lobe attenuation speed 14dB/otc, with it is preferably main,
Side lobe performance.
In addition, the electric energy gauging method based on the peak interpolation amendments of Kaiser windows FFT tetra- according to the above embodiment of the present invention,
There can also be following additional technical characteristic：
Further, step S1 is further included：
S101：Frequency spectrum based on the peak interpolation correction formulas of Kaiser windows FFT tetra-, the frequency displacement of FFT and window function Kaiser
Four peak spectral line Amplitude Rations are given expression in the amplitude of fundamental wave, harmonic wave and the symmetrical spectral line of m-Acetyl chlorophosphonazo crest frequency point arranged on left and right sides four
The formula of value parameter γ；
S102：Contrast value parameter γ carries out fitting of a polynomial and approaches, and fits logical based on the peak spectral lines of Kaiser windows FFT tetra-
With the formula of frequency correction factor δ；
S103：Amplitude formula fitting according to measured signal goes out the formula of general amplitude correction coefficient g (δ)；
S104：The formula of formula and universal frequency correction factor δ according to general amplitude correction coefficient g (δ) is led to
With amplitude correction formula, general phase only pupil filter formula and universal frequency correction formula.
Further, step S2 is further included：
S201：Apply Kaiser window functions respectively to voltage, current signal, and it is described to obtain to carry out FFT
The frequency spectrum of Kaiser window functions；
S202：Using peak detection technology, tested voltage, the fundamental wave of current signal, harmonic wave and humorous are detected respectively
Four amplitudes of symmetrical spectral line of crest value point arranged on left and right sides；
S203：Using four peak spectral line amplitude ratio parameter ratio formula, calculate respectively voltage, the fundamental wave of current signal,
The value of four peak spectral line amplitude ratio parameter γ at harmonic wave and m-Acetyl chlorophosphonazo peak point；
S204：According to four peak spectral line amplitude ratios at voltage, the fundamental wave of current signal, harmonic wave and m-Acetyl chlorophosphonazo peak point
The value of parameter γ, and measured signal fundamental wave, harmonic wave and humorous are calculated according to the peak interpolation correction formulas of Kaiser windows FFT tetra-
The value of the frequency correction factor δ at crest value point；
S205：The amplitude being worth at measured signal fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo peak point according to correction factor δ is repaiied
The value of positive coefficient g (δ)；
S206：The value interpolation amendment of value and correction factor g (δ) according to correction factor δ goes out tested voltage, current signal
The amplitude of fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo, phase and frequency parameter.
Further, the computing formula of universal frequency correction factor δ is：
The γ of δ=3.1041273 γ+0.70465993+0.35296224·γ5+0.22441945·γ7+
0.16913521×γ9。
Further, the computing formula of general amplitude correction coefficient g (δ) is：
G (δ)=1.50096168+0.22288389 δ2+0.01752579·δ4+0.00098765·δ6。
Further, the correction formula of general amplitude is：
A=2 (y1+2y2+2y3+y4)·g(δ)/N；
Wherein, y1、y2、y3And y4Respectively the symmetrical spectral line of spectrum curve peak point arranged on left and right sides four be respectively kth1
Bar, kth2Bar, kth3Bar and kth4The amplitude of bar spectral line, N is sampled data length.
Compared with prior art, the present invention has following beneficial effect：It is of the invention to be inserted based on the peaks of Kaiser windows FFT tetra-
The electric energy gauging method of value, take full advantage of four spectral lines of amplitude Frequency point or so symmetry and Kaiser window functions it is excellent
Master, side lobe performance, with electric energy metrical accuracy in computation higher.
Brief description of the drawings
Of the invention above-mentioned and/or additional aspect and advantage will become from description of the accompanying drawings below to embodiment is combined
Substantially and be readily appreciated that, wherein：
Fig. 1 is the flow chart of the electric energy gauging method based on the peak interpolations of Kaiser windows FFT tetra- of the invention；
Fig. 2 is that the power network fundamental frequency of one embodiment of the invention is left and right more than or less than peak point in the case of 50Hz
The schematic diagram of the symmetrical spectral line of four of both sides.
Fig. 3 is spectral characteristic of the Kaiser window functions when form parameter β is respectively 0,4,8,22.It is other with the increase of β
Valve peak value is smaller, but main lobe width broadens.When β=22, with relatively low side lobe peak and faster side lobe attenuation speed, together
When there is moderate main lobe width.
Specific embodiment
Embodiment below with reference to Description of Drawings is exemplary, is only used for explaining of the invention, and it is not intended that
Limitation of the present invention.
In the description of the invention, it is to be understood that term " first ", " second " are only used for describing purpose, and can not
It is interpreted as indicating or implying relative importance.
With reference to following description and accompanying drawing, it will be clear that these and other aspects of embodiments of the invention.In these descriptions
In accompanying drawing, specifically disclose some particular implementations in embodiments of the invention to represent implementation implementation of the invention
Some modes of the principle of example, but it is to be understood that the scope of embodiments of the invention is not limited.Conversely, of the invention
Embodiment includes all changes, modification and the equivalent that fall into the range of the spiritual and intension of attached claims.
Fig. 1 is the stream of the electric energy gauging method based on the peak interpolation amendments of Kaiser windows FFT tetra- of one embodiment of the invention
Cheng Tu.As shown in figure 1, a kind of electric energy gauging method based on the peak interpolation amendments of Kaiser windows FFT tetra-, comprises the following steps：
S1：Frequency spectrum based on the peak interpolation correction formulas of Kaiser windows FFT tetra-, the frequency displacement of FFT and window function Kaiser exists
The amplitude of fundamental wave, harmonic wave and the symmetrical spectral line of m-Acetyl chlorophosphonazo crest frequency point arranged on left and right sides four obtains signal fundamental wave, harmonic wave and humorous
The general amplitude of ripple, phase and frequency interpolation correction formula.
In one embodiment of the invention, step S1 is further included：
S101：Frequency spectrum based on the peak interpolation correction formulas of Kaiser windows FFT tetra-, the frequency displacement of FFT and window function Kaiser
Four peak spectral line Amplitude Rations are given expression in the amplitude of fundamental wave, harmonic wave and the symmetrical spectral line of m-Acetyl chlorophosphonazo crest frequency point arranged on left and right sides four
The formula of value parameter γ.
S102：Contrast value parameter γ carries out fitting of a polynomial and approaches, and fits logical based on the peak spectral lines of Kaiser windows FFT tetra-
With the formula of frequency correction factor δ.
S103：Amplitude formula fitting according to measured signal goes out the formula of general amplitude correction coefficient g (δ).
S104：The formula of formula and universal frequency correction factor δ according to general amplitude correction coefficient g (δ) is led to
With amplitude correction formula, general phase only pupil filter formula and universal frequency correction formula.
S2：The base of tested voltage signal, current signal is calculated according to the general amplitude, phase and frequency correction formula
The amplitude of ripple, harmonic wave and m-Acetyl chlorophosphonazo, phase and frequency parameter.
In one embodiment of the invention, step S2 is further included：
S201：Apply Kaiser window functions respectively to voltage, current signal, and it is described to obtain to carry out FFT
The frequency spectrum of Kaiser window functions；
S202：Using peak detection technology, tested voltage, the fundamental wave of current signal, harmonic wave and humorous are detected respectively
The amplitude of the symmetrical spectral line of crest value point arranged on left and right sides four；
S203：Using four peak spectral line amplitude ratio parameter ratio formula, calculate respectively voltage, the fundamental wave of current signal,
The value of four peak spectral line amplitude ratio parameter γ at harmonic wave and m-Acetyl chlorophosphonazo peak point；
S204：According to four peak spectral line amplitude ratios at voltage, the fundamental wave of current signal, harmonic wave and m-Acetyl chlorophosphonazo peak point
The value of parameter γ, and measured signal fundamental wave, harmonic wave and humorous are calculated according to the peak interpolation correction formulas of Kaiser windows FFT tetra-
The value of the frequency correction factor δ at crest value point；
S205：The amplitude being worth at measured signal fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo peak point according to correction factor δ is repaiied
The value of positive coefficient g (δ)；
S206：The value interpolation amendment of value and correction factor g (δ) according to correction factor δ goes out tested voltage, current signal
The amplitude of fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo, phase and frequency.
S3：According to amplitude, the phase and frequency of the tested voltage signal, the fundamental wave of current signal, harmonic wave and m-Acetyl chlorophosphonazo
Parameter calculates fundamental energy, harmonic electric energy, m-Acetyl chlorophosphonazo electric energy and total electric energy respectively.
In one embodiment of the invention, the computing formula of universal frequency correction factor δ is：
The γ of δ=3.1041273 γ+0.70465993+0.35296224·γ5+0.22441945·γ7+
0.16913521·γ9。
In one embodiment of the invention, the computing formula of general amplitude correction coefficient g (δ) is：
G (δ)=1.50096168+0.22288389 δ2+0.01752579·δ4+0.00098765·δ6。
The correction formula of general amplitude is：
A=2 (y1+2y2+2y3+y4)·g(δ)/N；
Wherein, y1、y2、y3And y4The respectively symmetrical spectral line of spectrum curve peak point arranged on left and right sides four is respectively kth1Bar,
Kth2Bar, kth3Bar and kth4The amplitude of bar spectral line, N is sampled data length.
To make it is further understood that the present invention, will be described further by following examples.
(1) based on the peak interpolation correction formulas of Kaiser windows FFT tetra-
The time-domain expression of Kaiser window functions：
In formula, I0(β) is the 1st class deformation zero Bessel function；β is the form parameter of Kaiser functions；
In formula, α is the difference between the main lobe value and side lobe levels of Kaiser windows, by setting different β value, thus it is possible to vary
The main lobe and side lobe performance of Kaiser window functions.
Universal frequency correction factors formula δ,
General amplitude correction coefficient formula g (δ),
G (δ)=1.50096168+0.22288389 δ2+0.01752579·δ4+0.00098765·δ6 (2)
General amplitude correction formula：
A=2 (y1+2y2+2y3+y4)·g(δ)/N (3)
General phase only pupil filter formula：
Universal frequency correction formula：
f0=k0Δ f=(δ+k2-0.5)fs/N (5)
(2) derivation
Assuming that the signal containing multiple harmonic (m-Acetyl chlorophosphonazo), with sample frequency fsAfter being sampled, discrete series x is obtained
(n)：
In formula, m represents overtone order；AmAnd θmThe amplitude and phase of m subharmonic are represented respectively；f0Represent fundamental frequency.
Formula (6) is used instead Euler's formula to be expressed as：
The forms of time and space of Kaiser window functions is w (n), and its discrete spectrum isThen FFT tables after x (n) adding windows
It is up to formula:
If ignoring negative frequency-mf0Locate the secondary lobe influence of spectral peak, in positive frequency mf0Neighbouring continuous frequency spectrum function is：
Formula (9) is sampled, the expression formula that can obtain its discrete Fourier transform is：
In formula, Δ f=fs/N is discrete frequency intervals；N is sampled data length.
With fundamental frequency f0=k0As a example by Δ f.Because mains frequency has fluctuation, and it is generally not accurate 50Hz, such as
This, is sampled with fixed frequency to measured signal, and fence effect will occur, i.e. the Frequency point of its highest discrete spectral line is very
The frequency of hardly possible exactly spectrum curve peak value, that is, k0It is generally not integer.If spectrum curve peak point arranged on left and right sides four is right
Spectral line is claimed to be respectively kth1Bar, kth2Bar, kth3Bar and kth4Bar spectral line (k1≤k2≤k0≤k3≤k4), correspondence amplitude is respectivelyAs shown in Figure 2.Introduce ginseng
Number δ=k0-k1- 0.5, δ ∈ [- 0.5,0.5] are obtained, make four peak spectral line amplitude ratio parameter ratio γ=((2y3+y4)-(2y2+
y1))/(y1+2y2+2y3+y4), order：
R=2 | W [2 π (- δ+0.5)/N] |+| W [2 π (- δ+1.5)/N] | (11)
S=2 | W [2 π (- δ -0.5)/N] |+| W [2 π (- δ -1.5)/N] | (12)
According to fitting of a polynomial regression formulas polyfit (γ, δ, 9), inverse function δ=f is obtained-1(γ), that is, obtain general
The formula of frequency correction factor δ, is shown in formula (1).
Due to the kth near peak point2And kth3The amplitude information amount of bar spectral line is maximum, when then calculating, gives this two
With bigger weight, specifically, in turn, the weight of four spectral lines is respectively 1,2,2,1 to spectral line, and the correction formula of amplitude is：
G (δ)=2N/ (R+S) is made, when n is large, approximation by polynomi-als is carried out by polyfit (δ, g (δ), 7) program,
The formula of amplitude correction coefficient g (δ) is obtained, formula (2) is seen.Bring formula (2) into formula (14), then obtain general amplitude correction public
Formula, is shown in formula (3).
Be can be seen that by formula (3), the actual fundamental wave peak value of measured signal is using the symmetrical spectral line of peak point arranged on left and right sides four
Amplitude weighted amendment obtain, therefore the algorithm is referred to as four peak interpolation correction algorithms.
(3) electric energy computing formula
There are the actual fundametal compoment amplitude of measured signal, the amendment type of phase and frequency, so that it may sufficiently accurately calculate
Fundamental wave, the amplitude of each harmonic and m-Acetyl chlorophosphonazo component, the phase and frequency of tested actual signal, and then can sufficiently accurately count
Calculate fundamental wave, each harmonic and m-Acetyl chlorophosphonazo electric energy.From orthogonality of trigonometric function, different integers time voltage x current harmonic wave point
Amount does not produce active energy, therefore the electric energy calculating formula of fundamental wave, harmonic wave, m-Acetyl chlorophosphonazo is：
In formula, WnIt is n-th harmonic or the active energy of m-Acetyl chlorophosphonazo, when n takes 1, W1It is fundamental energy；When n is taken more than 1
Integer when, WnIt is harmonic electric energy；When the negated integers of n, WnIt is m-Acetyl chlorophosphonazo electric energy；Un、InRespectively voltage signal, electric current are believed
Number n-th harmonic or m-Acetyl chlorophosphonazo amplitude, when n takes 1 be fundamental voltage amplitude；αn、βnRespectively voltage signal, current signal
The initial phase of n-th harmonic or m-Acetyl chlorophosphonazo, is fundamental phase when n takes 1；T is the primitive period of line voltage electric current；When k is
Between window number, according to IEC61000-4-7 standards, typically take 10 primitive periods.
Meter and fundamental wave, harmonic wave, total electric energy computing formula of m-Acetyl chlorophosphonazo are：
W=CjWj-CposxWposx+CoppxWoppx-CposjxWposjx+CoppjxWoppjx (16)
Wherein, WjIt is fundamental energy；WposxIt is positive harmonic electric energy；WoppxIt is reverse harmonic electric energy；WposjxPositive m-Acetyl chlorophosphonazo
Electric energy；WoppjxReverse m-Acetyl chlorophosphonazo electric energy；CjIt is fundamental energy weight coefficient；CposxIt is positive harmonic electric energy weight coefficient；CoppxFor
Reverse harmonic electric energy weight coefficient；CposjxIt is positive m-Acetyl chlorophosphonazo electric energy weight coefficient；CoppjxIt is for reverse m-Acetyl chlorophosphonazo electric energy is weighted
Number.Suggestion weight coefficient value：Cj=1,0<Cposx、Cposjx<1, Coppx、Coppjx>1。
(4) electric energy emulation
Power supply grid voltage, current simulations model as shown in table 1, wherein containing multiple harmonic composition, specially
In formula, f0It is fundamental frequency, is fluctuated generally between 49.7Hz~50.3Hz, sample frequency fs=4000Hz, according to
IEC61000-4-7 standards, the spectrum analysis time window of harmonic measure is 200ms by Uniform provisions.
The voltage of table 1, the fundamental wave of current signal and harmonic parameters
Hanning windows, Hamming windows, Blackman windows, 45 ranks of Nuttall windows and Kaiser windows (β is respectively adopted
=electric energy metrical simulation calculation 22) is carried out based on four peak interpolation algorithms, the length of window function elects as at 800 points, and fundamental frequency is
50.1Hz, simulation result is given in table 2.
The simulation result of the electric energy of table 2
By simulation result as can be seen that the calculating of the electric energy metrical algorithm based on the peak interpolation amendments of Kaiser windows FFT tetra- is accurate
Exactness compares one high of other four kinds of algorithms even two orders of magnitude, while the degree of accuracy of even-order harmonic electric energy is also carried significantly
Height, illustrates that the electric energy metrical algorithm based on the peak interpolation amendments of Kaiser windows FFT tetra- proposed by the present invention has meter higher really
Calculate exactness accurately.
In addition, other structures of the electric energy gauging method based on the peak interpolation amendments of Kaiser windows FFT tetra- of the embodiment of the present invention
Into and effect be all for a person skilled in the art known, in order to reduce redundancy, do not repeat.
In the description of this specification, reference term " one embodiment ", " some embodiments ", " example ", " specifically show
The description of example " or " some examples " etc. means to combine specific features, structure, material or spy that the embodiment or example are described
Point is contained at least one embodiment of the invention or example.In this manual, to the schematic representation of above-mentioned term not
Necessarily refer to identical embodiment or example.And, the specific features of description, structure, material or feature can be any
One or more embodiments or example in combine in an appropriate manner.
Although an embodiment of the present invention has been shown and described, it will be understood by those skilled in the art that：Not
Can these embodiments be carried out with various changes, modification, replacement and modification in the case of departing from principle of the invention and objective, this
The scope of invention is by claim and its equivalent limits.
Claims (6)
1. a kind of electric energy gauging method based on the peak interpolations of Kaiser windows FFT tetra-, it is characterised in that comprise the following steps：
S1：Frequency spectrum based on the peak interpolation correction formulas of Kaiser windows FFT tetra-, the frequency displacement of FFT and window function Kaiser fundamental wave,
The amplitude of harmonic wave and the symmetrical spectral line of m-Acetyl chlorophosphonazo crest frequency point arranged on left and right sides four, using fitting of a polynomial, obtains signal base
The general amplitude of ripple, harmonic wave and m-Acetyl chlorophosphonazo, phase and frequency interpolation correction formula；
S2：The tested voltage signal, fundamental wave of current signal, humorous is calculated according to the general amplitude, phase and frequency correction formula
The amplitude of ripple and m-Acetyl chlorophosphonazo, phase and frequency parameter；
S3：Amplitude, phase and frequency parameter according to the tested voltage signal, the fundamental wave of current signal, harmonic wave and m-Acetyl chlorophosphonazo
Fundamental energy, harmonic electric energy, m-Acetyl chlorophosphonazo electric energy and total electric energy are calculated respectively.
2. the electric energy gauging method based on the peak interpolations of Kaiser windows FFT tetra- according to claim 1, it is characterised in that step
Rapid S1 is further included：
S101：Frequency spectrum based on the peak interpolation correction formulas of Kaiser windows FFT tetra-, the frequency displacement of FFT and window function Kaiser is in base
The amplitude of ripple, harmonic wave and the symmetrical spectral line of m-Acetyl chlorophosphonazo crest frequency point arranged on left and right sides four gives expression to four peak spectral line amplitude ratio ginsengs
The formula of number γ；
S102：Contrast value parameter γ carries out fitting of a polynomial and approaches, and fits the general frequency based on the peak spectral lines of Kaiser windows FFT tetra-
The formula of rate correction factor δ；
S103：Amplitude formula fitting according to measured signal goes out the formula of general amplitude correction coefficient g (δ)；
S104：The formula of formula and universal frequency correction factor δ according to general amplitude correction coefficient g (δ) obtains general width
Value correction formula, general phase only pupil filter formula and universal frequency correction formula.
3. the electric energy gauging method based on the peak interpolations of Kaiser windows FFT tetra- according to claim 2, it is characterised in that step
Rapid S2 is further included：
S201：Apply Kaiser window functions respectively to voltage, current signal, and carry out FFT to obtain the Kaiser windows
The frequency spectrum of function；
S202：Using peak detection technology, tested voltage, the fundamental wave of current signal, harmonic wave and m-Acetyl chlorophosphonazo peak are detected respectively
Four amplitudes of symmetrical spectral line of value point arranged on left and right sides；
S203：Using four peak spectral line amplitude ratio parameter ratio formula, voltage, the fundamental wave of current signal, harmonic wave are calculated respectively
And the value of four peak spectral line amplitude ratio parameter γ at m-Acetyl chlorophosphonazo peak point；
S204：According to four peak spectral line amplitude ratio parameters at voltage, the fundamental wave of current signal, harmonic wave and m-Acetyl chlorophosphonazo peak point
The value of γ, and measured signal fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo peak are calculated according to the peak interpolation correction formulas of Kaiser windows FFT tetra-
The value of the frequency correction factor δ at value point；
S205：According to the amplitude correction system being worth at measured signal fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo peak point of correction factor δ
The value of number g (δ)；
S206：The value interpolation amendment of value and correction factor g (δ) according to correction factor δ goes out the base of tested voltage, current signal
The amplitude of ripple, harmonic wave and m-Acetyl chlorophosphonazo, phase and frequency parameter.
4. the electric energy gauging method based on the peak interpolations of Kaiser windows FFT tetra- according to claim 2, it is characterised in that logical
Computing formula with frequency correction factor δ is：
The γ of δ=3.1041273 γ+0.70465993+0.35296224·γ5+0.22441945·γ7+
0.16913521·γ9。
5. the electric energy gauging method based on the peak interpolations of Kaiser windows FFT tetra- according to claim 4, it is characterised in that logical
Computing formula with amplitude correction coefficient g (δ) is：
G (δ)=1.50096168+0.22288389 δ2+0.01752579·δ4+0.00098765·δ6。
6. the electric energy gauging method based on the peak interpolations of Kaiser windows FFT tetra- according to claim 5, it is characterised in that logical
It is with the correction formula of amplitude：
A=2 (y1+2y2+2y3+y4)·g(δ)/N；
Wherein, y1、y2、y3And y4Respectively the symmetrical spectral line of spectrum curve peak point arranged on left and right sides four be respectively kth1Article,
k2Bar, kth3Bar and kth4The amplitude of bar spectral line, N is sampled data length.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611157851.2A CN106814230A (en) | 2016-12-15 | 2016-12-15 | Electric energy gauging method based on the peak interpolations of Kaiser windows FFT tetra- |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611157851.2A CN106814230A (en) | 2016-12-15 | 2016-12-15 | Electric energy gauging method based on the peak interpolations of Kaiser windows FFT tetra- |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106814230A true CN106814230A (en) | 2017-06-09 |
Family
ID=59109563
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611157851.2A Pending CN106814230A (en) | 2016-12-15 | 2016-12-15 | Electric energy gauging method based on the peak interpolations of Kaiser windows FFT tetra- |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106814230A (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107945147A (en) * | 2017-12-04 | 2018-04-20 | 云南电网有限责任公司普洱供电局 | One kind is based on the tetra- modified image interfusion methods of peak interpolation of Kaiser windows FFT |
| CN110954746A (en) * | 2019-11-27 | 2020-04-03 | 云南电网有限责任公司电力科学研究院 | Six-interpolation FFT algorithm based on four-term Nuttall cosine window |
| CN111122941A (en) * | 2019-12-04 | 2020-05-08 | 国网湖南综合能源服务有限公司 | Kaiser window function-based voltage sag characteristic quantity detection method, system and medium for improving S conversion |
| CN111221812A (en) * | 2020-04-13 | 2020-06-02 | 广东电网有限责任公司佛山供电局 | Data interpolation processing method for generator parameter identification |
| CN112073339A (en) * | 2019-06-11 | 2020-12-11 | 大唐移动通信设备有限公司 | Method and device for determining calibration information |
| CN113341223A (en) * | 2021-06-02 | 2021-09-03 | 贵州电网有限责任公司 | Method suitable for power grid harmonic analysis and harmonic source positioning |
| CN113640579A (en) * | 2021-10-13 | 2021-11-12 | 四川大学 | Harmonic measurement method based on double spectral line transformation, electronic device and storage medium |
-
2016
- 2016-12-15 CN CN201611157851.2A patent/CN106814230A/en active Pending
Non-Patent Citations (3)
| Title |
|---|
| 徐文佳 等: ""计及间谐波的电能计量方法"", 《电网技术》 * |
| 李得民 等: ""基于 Nuttall 窗四谱线插值 FFT 的电力谐波分析"", 《电力系统保护与控制》 * |
| 陈立伟 等: ""Kaiser 窗三谱线插值电力谐波分析"", 《应用科技》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107945147A (en) * | 2017-12-04 | 2018-04-20 | 云南电网有限责任公司普洱供电局 | One kind is based on the tetra- modified image interfusion methods of peak interpolation of Kaiser windows FFT |
| CN112073339A (en) * | 2019-06-11 | 2020-12-11 | 大唐移动通信设备有限公司 | Method and device for determining calibration information |
| CN112073339B (en) * | 2019-06-11 | 2021-12-28 | 大唐移动通信设备有限公司 | Method and device for determining calibration information |
| CN110954746A (en) * | 2019-11-27 | 2020-04-03 | 云南电网有限责任公司电力科学研究院 | Six-interpolation FFT algorithm based on four-term Nuttall cosine window |
| CN111122941A (en) * | 2019-12-04 | 2020-05-08 | 国网湖南综合能源服务有限公司 | Kaiser window function-based voltage sag characteristic quantity detection method, system and medium for improving S conversion |
| CN111221812A (en) * | 2020-04-13 | 2020-06-02 | 广东电网有限责任公司佛山供电局 | Data interpolation processing method for generator parameter identification |
| CN113341223A (en) * | 2021-06-02 | 2021-09-03 | 贵州电网有限责任公司 | Method suitable for power grid harmonic analysis and harmonic source positioning |
| CN113640579A (en) * | 2021-10-13 | 2021-11-12 | 四川大学 | Harmonic measurement method based on double spectral line transformation, electronic device and storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106814230A (en) | Electric energy gauging method based on the peak interpolations of Kaiser windows FFT tetra- | |
| CN106814246A (en) | Electric energy gauging method based on the Kaiser windows unimodal interpolation of FFT | |
| CN106526312A (en) | Electric energy metering method based on R-V(III) window FFT double-peak interpolation | |
| Barbaro et al. | A novel approach based on nonactive power for the identification of disturbing loads in power systems | |
| CN201352236Y (en) | Harmonic power meter based on Kaiser window double-spectrum-line interpolation FFT | |
| CN101567559B (en) | Tabular method of non-intrusive electrical load decomposition | |
| CN109030941A (en) | Tri- spectral line interpolation harmonic analysis method of Hanning involution convolution window FFT | |
| CN107179439A (en) | Electric energy gauging method based on five peak interpolations of Cosine Window FFT tri- | |
| CN102662106B (en) | Method for measuring electric energy of harmonic power network | |
| CN103116064A (en) | Method and device for detecting voltage fluctuation and flicker based on energy operator and spectrum correction | |
| Dirik et al. | A computation and metering method for harmonic emissions of individual consumers | |
| CN104897961A (en) | Three spectral line interpolation FFT harmonic wave analysis method and system based on multiplication window function | |
| Li et al. | Research of harmonic distortion power for harmonic source detection | |
| CN103675447B (en) | A kind of high-precision real time-harmonic wave analysis method of electric railway | |
| CN102938130A (en) | Comprehensive assessment method for performance indexes of power quality based on PQView | |
| CN101950012A (en) | Field tester for alternating current (AC) energy meter | |
| CN105842555A (en) | Multi-harmonic source responsibility division method based on measured data | |
| CN106154040A (en) | One is used for calculating electrical network single-point equivalent time Equivalent Harmonic impedance methodologies | |
| CN110031674A (en) | Motor train unit train energy consumption metering system and method | |
| CN109494735A (en) | A kind of low-voltage network line loss calculation model and method | |
| CN102185314A (en) | Power quality comprehensive evaluation method | |
| CN106771552B (en) | A kind of distortion power metering method | |
| CN106712021A (en) | Multi-harmonic-limit multi-period distribution method for public grid | |
| CN106709623A (en) | Power grid marketing inspection risk control method based on risk calculation model | |
| Tan et al. | High‐accurate estimation method of sub‐synchronous and super‐synchronous harmonic parameters on PMU |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170609 |