CN106526312A - Electric energy metering method based on R-V(III) window FFT double-peak interpolation - Google Patents
Electric energy metering method based on R-V(III) window FFT double-peak interpolation Download PDFInfo
- Publication number
- CN106526312A CN106526312A CN201610882937.5A CN201610882937A CN106526312A CN 106526312 A CN106526312 A CN 106526312A CN 201610882937 A CN201610882937 A CN 201610882937A CN 106526312 A CN106526312 A CN 106526312A
- Authority
- CN
- China
- Prior art keywords
- frequency
- iii
- formula
- electric energy
- amplitude
- 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
- G01R22/00—Arrangements for measuring time integral of electric power or current, e.g. electricity meters
- G01R22/06—Arrangements for measuring time integral of electric power or current, e.g. electricity meters by electronic methods
- G01R22/10—Arrangements for measuring time integral of electric power or current, e.g. electricity meters by electronic methods using digital techniques
Abstract
The invention discloses an electric energy metering method based on R-V(III) window FFT double-peak interpolation. The method comprises the steps of based on an Rife-Vincent(III) window FFT double-peak-spectrum linear interpolation correction formula, frequency shifting property of FFT and frequency spectrum of a window function Rife-Vincent(III), obtaining a universal amplitude, phase and frequency interpolation correction formula of a signal fundamental wave, a harmonic wave and an interharmonic wave at a highest spectrum line spectrum and a secondary highest spectrum line spectrum at the left side and the right side of the fundamental wave, the harmonic wave and the interharmonic wave; calculating the amplitudes, the phases and the frequency parameters of the fundamental wave, the harmonic wave and the interharmonic wave of the metered voltage signal and the current signal according to the universal amplitude, the phase and the frequency correction formula; and furthermore calculating the electric energy of the fundamental wave, the electric energy of the harmonic wave, the electric energy of the interharmonic wave and the total electric energy. The electric energy metering method has an advantage of relatively high main lobe performance and relatively high side lobe performance.
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 R-V (III) window FFT double peak interpolations.
Background technology
With increasing nonlinear-load, particularly rolling mill, elevator, electric arc furnace, electric railway etc. it is large-scale
Application of the nonlinear-load in electric power system, generates a large amount of electric harmonics and m-Acetyl chlorophosphonazo in power supply grid.A large amount of electric power
The presence of harmonic wave and m-Acetyl chlorophosphonazo, not only has a strong impact on the quality of power supply of power supply grid, jeopardizes the safe and stable operation of electric power system,
Serious electric power accident is even resulted in, huge economic loss is caused, and is affected 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 which be measured correctly, exactly,
It is directly connected to the economic benefit of power supply and demand both sides.Traditional way of energy measuring, is directly to define by electric energy, 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
Electromotor and photovoltaic plant generate electricity by way of merging two or more grid systems, and electric arc furnace, electric welding machine, rolling mill, mine hoist, elevator etc. are non-linear negative
The presence of lotus, in power supply grid, serves not only as load and wants power consumption, and as harmonic source, can inject to power supply grid
A large amount of harmonic powers.If by traditional full way of energy measuring, although these Harmonics source customers can inject a large amount of harmonic waves to electrical network
Interference, but also pay tariff less；And vast common electricity consumer, not only by the interference and harm of harmonic wave, and more than hand over
The electricity charge.Meanwhile, this full way of energy measuring is not accounted for because of non-synchronous sampling and caused spectral leakage and fence effect yet,
Its metric results there may be relatively large deviation.
And another kind of so-called fundamental energy metering method, although the electricity consumption fairness problem of conventional power user can be solved,
Making them 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 which is electric to power supply
The electric harmonic pollution of net injection is but still not affected by any punishment.
It can be seen that, whether using full way of energy measuring or using fundamental energy metering method, how all fail solution
The measurement problem of a large amount of harmonic waves, the harmonic wave that m-Acetyl chlorophosphonazo is caused, m-Acetyl chlorophosphonazo electric energy is injected to power supply grid to Harmonics source customer, is more talked
Do not go up how economic punishment Harmonics source customer, how to administer harmonic source.
How the parameters such as amplitude, the phase and frequency of fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo is more accurately extracted, paper is " based on more
Application of the unimodal spectral line interpolation algorithm that formula is approached in m-Acetyl chlorophosphonazo analysis " (Xiao is first brave etc. electric power network technique, 2008), be given
It is based on Hanning windows, Hamming windows, Blackman-Harris windows, the unimodal spectral line interpolation of Rife-Vincent (III) window
Correction formula；Paper " carrying out the innovatory algorithm of harmonic analysis in power system using FFT " (Pang Hao etc. Proceedings of the CSEE,
2003), give based on rectangular window, Triangle windows, Hanning windows, Hamming windows, Blackman-Harris windows, Rife-
The double peak of spectral lines interpolation correction formula of Vincent (III) window；Paper " the Electric Power Harmonic Analysis algorithms based on three spectral line interpolation FFTs "
(cattle victory lock etc. Proceedings of the CSEE, 2012), " three spectral line interpolation FFTs based on five Lai Fu-Vincent window are humorous for paper
Wave analysiss " (Zhu Junwei etc. electrical technology, 2015), paper " the three peak interpolation harmonic wave Algorithm Analysis based on Nuttall windows " (Zhai
It is auspicious vast etc. protecting electrical power system and control, 2015), sets forth based on Hanning windows, Hamming windows, Blackman windows,
Blackman-Harris windows, Rife-Vincent windows, three peak spectral line interpolation correction formulas of Nuttall windows；Paper is " based on four
The power grid harmonic wave detection method of spectral line interpolation FFT " (Hao Zhu etc. protecting electrical power system and control, 2014) " it is based on paper
The Electric Power Harmonic Analysis of four spectral line interpolation FFT of Nuttall windows " (Li get Min etc. protecting electrical power system and control, 2016), respectively
Give based on Triangle windows, Hanning windows, Hamming windows, Blackman, Nuttall window four peak spectral line interpolation amendments
Formula.
But the problem that these existing algorithms are still suffered from is：1) in each class difference correction algorithm, based on different window function shapes
Into the counting accuracy of specific algorithm differ greatly；2), between inhomogeneity interpolation correction algorithm, unimodal spectral line interpolation amendment is public
Formula is more complicated, easily affected by spectrum leakage, double peak of spectral lines algorithm near the peak value frequency to leaking spectral line information using insufficient,
Three peak spectral line algorithms do not consider the comprehensive information of the symmetrical spectral line of peak value frequency, four peak spectral line precision highests but algorithm amount of calculation
Greatly increase, it is difficult to meet the requirements at the higher level such as real-time electric energy metrical.
The content of the invention
It is contemplated that at least solving one of above-mentioned technical problem.
For this purpose, it is an object of the invention to proposing a kind of based on the bimodal spectral lines of R-V (i.e. Rife-Vincent) (III) window FFT
The electric energy gauging method of interpolation amendment, can weigh the counting accuracy and computation complexity of electric energy metrical algorithm.
To achieve these goals, embodiment of the invention discloses that a kind of electricity based on R-V (III) window FFT double peak interpolations
Energy metering method, comprises the following steps：S1：Based on Rife-Vincent (III) window FFT double peak of spectral lines interpolation correction formulas, FFT
Frequency displacement and window function Rife-Vincent (III) frequency spectrum in fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo crest frequency point arranged on left and right sides
Maximum spectral line and time maximum breadth of spectral line are worth to the general amplitude of signal fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo, phase and frequency interpolation and repair
Positive formula；S2：According to the general amplitude, phase and frequency correction formula calculate tested voltage signal, the fundamental wave of current signal,
The amplitude of harmonic wave and m-Acetyl chlorophosphonazo, phase and frequency parameter；S3：According to the tested voltage signal, the fundamental wave of current signal, harmonic wave
Fundamental energy, harmonic electric energy, m-Acetyl chlorophosphonazo electric energy and total electric energy is calculated respectively with the amplitude of m-Acetyl chlorophosphonazo, phase and frequency parameter.
Electric energy gauging method based on R-V (III) window FFT double peak interpolations according to embodiments of the present invention, Rife-Vincent
(III) main lobe width of window function is 16 π/N, and side lobe peak is -73.9dB, and side lobe attenuation speed has preferable up to 12dB/otc
Master, side lobe performance.
In addition, the electric energy gauging method based on R-V (III) window FFT double peak interpolations according to the above embodiment of the present invention, also
There can be following additional technical characteristic：
Further, step S1 is further included：S101：Based on Rife-Vincent (III) window FFT double peak of spectral lines interpolation
The frequency spectrum of correction formula, the frequency displacement of FFT and window function Rife-Vincent (III) is in fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo crest frequency
Point arranged on left and right sides maximum spectral line and time maximum spectral line amplitude give expression to the formula of bimodal breadth of spectral line value ratio parameter β；S102：It is right
Ratio parameter β carries out fitting of a polynomial and approaches, and fits the universal frequency based on the bimodal spectral lines of Rife-Vincent (III) window FFT
The formula of correction factor α；S103：The public affairs of general amplitude correction coefficient g (α) are gone out according to the amplitude formula fitting of measured signal
Formula；S104：General amplitude is obtained according to the formula of the formula and universal frequency correction factor α of general amplitude correction coefficient g (α)
Correction formula, general phase only pupil filter formula and frequency interpolation correction formula.
Further, step S2 is further included：S201：Rife-Vincent is applied respectively to voltage, current signal
(III) window function, and carry out FFT to obtain the frequency spectrum of the Rife-Vincent (III) window function；S202：Using peak value
Detection technique, detects tested voltage, the fundamental wave of current signal, harmonic wave and m-Acetyl chlorophosphonazo peak point arranged on left and right sides most respectively
The amplitude of big spectral line and time maximum spectral line；S203：Using bimodal breadth of spectral line value ratio parameter ratio formula, electricity is calculated respectively
The value of the bimodal breadth of spectral line value ratio parameter β at pressure, the fundamental wave of current signal, harmonic wave and m-Acetyl chlorophosphonazo peak point；S204：According to
The value of the bimodal breadth of spectral line value ratio parameter β at voltage, the fundamental wave of current signal, harmonic wave and m-Acetyl chlorophosphonazo peak point, and according to
Rife-Vincent (III) window FFT double peak of spectral lines interpolation correction formula calculates measured signal fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo
The value of the frequency correction factor α at peak point；S205：According to correction factor α be worth to measured signal fundamental wave, harmonic wave and
The value of amplitude correction coefficient g (α) at harmonic spike point；S206：According to the value of correction factor α and being worth for correction factor g (α)
To amplitude, the phase and frequency of tested voltage, the fundamental wave of current signal, harmonic wave and m-Acetyl chlorophosphonazo.
Further, the computing formula of universal frequency correction factor α is：
+ 0.18754082 β of α=3.04072904 β3+0.08804904·β5+0.05199160·β7。
Further, the computing formula of general amplitude correction coefficient g (α) is：
G (α)=1.08537118+0.30341579 α2+0.04563736·α4+0.00514598·α6。
Further, the correction formula of general amplitude is：
A=(y1+y2)·g(α)/N；
Wherein, y1And y2Respectively the maximum spectral line of spectrum curve peak point arranged on left and right sides and time maximum spectral line are respectively the
k1Bar and kth2The amplitude of bar spectral line, N are sampled data length.
The additional aspect and advantage of the present invention will be set forth in part in the description, and partly will become from the following description
Obtain substantially, or recognized by the practice of the present invention.
Description of the drawings
The above-mentioned and/or additional aspect and advantage of the present invention will become from the description with reference to accompanying drawings below to embodiment
It is substantially and easy to understand, wherein：
Fig. 1 is the flow chart of the electric energy gauging method based on R-V (III) window FFT double peak interpolations of the present invention；
Fig. 2 is that the electrical 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 maximum spectral line of both sides and time maximum spectral line.
Specific embodiment
It is exemplary below with reference to the embodiment of Description of Drawings, is only used for explaining of the invention, and it is not intended that
Limitation of the present invention.
In describing 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 implies relative importance.
With reference to explained below and accompanying drawing, it will be clear that in terms of these and other of embodiments of the invention.In these descriptions
In accompanying drawing, some particular implementations in embodiments of the invention are specifically disclosed, represent the enforcement for implementing the present 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, the present invention
Embodiment includes all changes, modification and the equivalent fallen in the range of the spirit and intension of attached claims.
Fig. 1 is the flow process of the electric energy gauging method based on R-V (III) window FFT double peak interpolations of one embodiment of the invention
Figure.As shown in figure 1, a kind of electric energy gauging method based on R-V (III) window FFT double peak interpolations, comprises the following steps：
S1：Based on Rife-Vincent (III) window FFT double peak of spectral lines interpolation correction formulas, the frequency displacement of FFT and window function
The frequency spectrum of Rife-Vincent (III) is in fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo crest frequency point arranged on left and right sides maximum spectral line and secondary maximum
Breadth of spectral line is worth to the general amplitude of signal fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo, phase and frequency interpolation correction formula.
In one embodiment of the invention, step S1 is further included：
S101：Based on Rife-Vincent (III) window FFT double peak of spectral lines interpolation correction formulas, the frequency displacement of FFT and window letter
Number Rife-Vincent (III) frequency spectrum fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo crest frequency point arranged on left and right sides maximum spectral line and it is secondary most
Big spectral line amplitude gives expression to the formula of bimodal breadth of spectral line value ratio parameter β.
S102：Reduced value parameter beta carries out fitting of a polynomial and approaches, and fits double based on Rife-Vincent (III) windows FFT
The formula of the universal frequency correction factor α of peak spectral line.
S103：The formula of general amplitude correction coefficient g (α) is gone out according to the amplitude formula fitting of measured signal.
S104：Led to according to the formula of the formula and universal frequency correction factor α of general amplitude correction coefficient g (α)
With amplitude correction formula, general phase only pupil filter formula and frequency interpolation 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 Rife-Vincent (III) window function to voltage, current signal respectively, and carry out FFT to obtain
To the frequency spectrum of the Rife-Vincent (III) window function；
S202：Using peak detection technology, tested voltage, the fundamental wave of current signal, harmonic wave and humorous is detected respectively
The amplitude of the maximum spectral line of crest value point arranged on left and right sides and time maximum spectral line；
S203：Using bimodal breadth of spectral line value ratio parameter ratio formula, calculate respectively voltage, the fundamental wave of current signal,
The value of the bimodal breadth of spectral line value ratio parameter β at harmonic wave and m-Acetyl chlorophosphonazo peak point；
S204：According to the bimodal breadth of spectral line value ratio at voltage, the fundamental wave of current signal, harmonic wave and m-Acetyl chlorophosphonazo peak point
The value of parameter beta, and measured signal base is calculated according to the Rife-Vincent (III) window FFT double peak of spectral lines interpolation correction formula
The value of the frequency correction factor α at ripple, harmonic wave and m-Acetyl chlorophosphonazo peak point；
S205：Repaiied according to the amplitude being worth at measured signal fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo peak point of correction factor α
The value of positive coefficient g (α)；
S206：According to the value of correction factor α and correction factor g (α) be worth to tested voltage, the fundamental wave of current signal,
The amplitude of 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：
+ 0.18754082 β of α=3.04072904 β3+0.08804904·β5+0.05199160·β7。
In one embodiment of the invention, the computing formula of universal frequency correction factor α is：
The correction formula of general amplitude is：
A=(y1+y2)·g(α)/N；
Wherein, y1And y2Respectively the maximum spectral line of spectrum curve peak point arranged on left and right sides and time maximum spectral line are respectively the
k1Bar and kth2The amplitude of bar spectral line, N are sampled data length.
For making it is further understood that the present invention, will be described further by following examples.
(1) Rife-Vincent (III) window FFT double peak of spectral lines interpolation correction formulas are based on
The time-domain expression of Rife-Vincent (III) window function：
W (n)=1.0-1.43596 cos (2 π n/N)+0.49754 cos (4 π n/N) -0.06158 cos (6
π·n/N)
Universal frequency correction factors formula α,
+ 0.18754082 β of α=3.04072904 β3+0.08804904·β5+0.05199160·β7 (1)
General amplitude correction coefficient formula g (α),
G (α)=1.08537118+0.30341579 α2+0.04563736·α4+0.00514598·α6 (2)
General amplitude correction formula：
A=(y1+y2)·g(α)/N (3)
General phase only pupil filter formula：
Universal frequency correction formula：
f0=k0Δ f=(α+0.5+ki- 1) fs/N, (i=1,2) (5)
In formula, when the actual frequency of fundametal compoment are less than 50Hz, i takes 1；Otherwise i takes 2.
(2) derivation
Assume the signal containing multiple harmonic (m-Acetyl chlorophosphonazo), with sample frequency fsAfter being sampled, discrete serieses x is obtained
(n)：
In formula, m represents overtone order；AmAnd θmThe amplitude and phase place of m subharmonic is represented respectively；f0Represent fundamental frequency.
Use formula (6) instead Euler's formula to be expressed as：
The forms of time and space of Rife-Vincent (III) window function is w (n), and its discrete spectrum isThen x (n) adding windows
FFT expression formulas afterwards are:
If ignoring negative frequency-mf0The secondary lobe of place's spectral peak affects, in positive frequency mf0Neighbouring continuous frequency spectrum function is：
Formula (9) is sampled, the expression formula that its discrete Fourier transform is obtained is：
In formula, Δ f=fs/N is discrete frequency intervals；N is sampled data length.
With fundamental frequency f0=k0As a example by Δ f.As mains frequency has undulatory property, and it is generally not accurate 50Hz, such as
This, is sampled to measured signal with fixed frequency, and fence effect will occur, i.e. the Frequency point of its highest discrete spectral line is very
The frequency of difficult exactly spectrum curve peak value, that is, k0It is generally not integer.If the maximum of spectrum curve peak point arranged on left and right sides
Spectral line and time maximum spectral line are respectively kth1Bar and kth2Bar spectral line, corresponding amplitude are respectivelyWithAs shown in Figure 2.Introduce parameter alpha=k0-k1- 0.5, α ∈ [- 0.5,0.5] are obtained, bimodal breadth of spectral line value is made
Ratio parameter ratio beta：
According to fitting of a polynomial regression formulas polyfit, (β, α n), obtain inverse function α=f-1(β), that is, obtain general frequency
The formula of rate correction factor α, is shown in formula (1).
If the correction formula of measured signal fundamental voltage amplitude A is：
G (α)=2 N/ (| W [2 π (- α -0.5)/N] |+| W [2 π (- α+0.5)/N] |) is made, fitting of a polynomial is brought into and is approached
(α, g (α) 7), can obtain general amplitude correction coefficient formula g (α), see formula (2) formula polyfit.Bring formula (2) into formula
(12) general amplitude correction formula, is then obtained, formula (3) is seen.
Can be seen that by formula (3), the actual fundamental wave peak value of measured signal is the maximum and time most using peak point arranged on left and right sides
What the weighted amendment of amplitude of big spectral line was obtained, therefore the algorithm is referred to as double peak of spectral lines interpolation correction algorithm.
(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
The fundamental wave of tested actual signal, the amplitude of each harmonic and m-Acetyl chlorophosphonazo component, phase and frequency, and then sufficiently accurately can 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, WnFor n-th harmonic or the active energy of m-Acetyl chlorophosphonazo, when n takes 1, W1For fundamental energy；When n is taken more than 1
Integer when, WnFor harmonic electric energy；When the negated integers of n, WnFor 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；Primitive periods of the T for 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 (14)
Wherein, WjFor fundamental energy；WposxFor positive harmonic electric energy；WoppxFor reverse harmonic electric energy；WposjxPositive m-Acetyl chlorophosphonazo
Electric energy；WoppjxReverse m-Acetyl chlorophosphonazo electric energy；CjFor fundamental energy weight coefficient；CposxFor positive harmonic electric energy weight coefficient；CoppxFor
Reverse harmonic electric energy weight coefficient；CposjxFor positive m-Acetyl chlorophosphonazo electric energy weight coefficient；CoppjxFor the weighting of reverse m-Acetyl chlorophosphonazo electric energy it is
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 are as shown in table 1, wherein containing multiple harmonic composition, specially
In formula, f0For fundamental frequency, fluctuate generally between 49.7Hz～50.3Hz, sample frequency fs=4000Hz, according to
IEC61000-4-7 standards, the spectrum analyses time window of harmonic measure are 10 cycles by Uniform provisions, i.e. 200ms.
1 voltage of table, the fundamental wave of current signal and harmonic parameters
Frequency is respectively set to 49.7Hz, 50.3Hz and 50Hz, and the simulation result of corresponding electric energy is shown in Table 2.
The simulation result of 2 electric energy of table
By simulation result as can be seen that the derived proposition of the present invention based on Rife-Vincent (III) window FFT Bimodal Spectras
The electric energy metrical algorithm of line interpolation amendment, really with higher accuracy.By respectively with plus Rectangle windows,
Triangle windows, Hanning windows, Hamming windows, Blackman windows, Blackman-Harris windows, Nutall windows, Kaiser windows
And the arithmetic speed of the FFT double peak of spectral lines interpolation electric energy gauging methods of various combination Cosine Windows is compared, what this patent was proposed
Arithmetic speed ratio based on the electric energy metrical algorithm of Rife-Vincent (III) window is very fast；And with unimodal, three peaks, four peak spectral line bases
Be compared in Rife-Vincent (III) window FFT interpolation electric energy metrical algorithms, this patent propose algorithm arithmetic speed also compared with
It is moderate.Simulation result shows that the method has preferably weighed the counting accuracy and computation complexity of electric energy metrical algorithm, can make
One of alternative for digitalized electrical energy meter electric energy metrical algorithm.
In addition, the electric energy gauging method based on R-V (III) window FFT double peak interpolations of the embodiment of the present invention other constitute with
And effect is all known for a person skilled in the art, in order to reduce redundancy, is not repeated.
In the description of this specification, reference term " one embodiment ", " some embodiments ", " example ", " specifically show
Example ", or the description of " some examples " etc. mean specific features with reference to the embodiment or example description, structure, material or spy
Point is contained at least one embodiment or example of the present invention.In this manual, to the schematic representation of above-mentioned term not
Identical embodiment or example are referred to necessarily.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
These embodiments can be carried out with various changes, modification, replacement and modification in the case of the principle and objective that depart from the present invention, this
The scope of invention is limited by claim and its equivalent.
Claims (6)
1. one kind is based on the electric energy gauging method of R-V (III) window FFT double peak interpolations, it is characterised in that comprise the following steps：
S1：Based on Rife-Vincent (III) window FFT double peak of spectral lines interpolation correction formulas, the frequency displacement of FFT and window function Rife-
The frequency spectrum of Vincent (III) is in fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo crest frequency point arranged on left and right sides maximum spectral line and secondary maximum breadth of spectral line
It is worth to general amplitude, the phase and frequency interpolation correction formula of signal fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo；
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 is calculated respectively.
2. the electric energy gauging method based on R-V (III) window FFT double peak interpolations according to claim 1, it is characterised in that step
Rapid S1 is further included：
S101：Based on Rife-Vincent (III) window FFT double peak of spectral lines interpolation correction formulas, the frequency displacement of FFT and window function
The frequency spectrum of Rife-Vincent (III) is in fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo crest frequency point arranged on left and right sides maximum spectral line and secondary maximum
Spectral line amplitude gives expression to the formula of bimodal breadth of spectral line value ratio parameter β；
S102：Reduced value parameter beta carries out fitting of a polynomial and approaches, and fits based on Rife-Vincent (III) window FFT Bimodal Spectras
The formula of the universal frequency correction factor α of line；
S103：The formula of general amplitude correction coefficient g (α) is gone out according to the amplitude formula fitting of measured signal；
S104：General width is obtained according to the formula of the formula and universal frequency correction factor α of general amplitude correction coefficient g (α)
Value correction formula, general phase only pupil filter formula and frequency interpolation correction formula.
3. the electric energy gauging method based on R-V (III) window FFT double peak interpolations according to claim 2, it is characterised in that step
Rapid S2 is further included：
S201：Apply Rife-Vincent (III) window function to voltage, current signal respectively, and FFT is carried out to obtain
State the frequency spectrum of Rife-Vincent (III) window function；
S202：Using peak detection technology, tested voltage, the fundamental wave of current signal, harmonic wave and m-Acetyl chlorophosphonazo peak is detected respectively
The amplitude of the maximum spectral line and time maximum spectral line of value point arranged on left and right sides；
S203：Using bimodal breadth of spectral line value ratio parameter ratio formula, voltage, the fundamental wave of current signal, harmonic wave is calculated respectively
And the value of the bimodal breadth of spectral line value ratio parameter β at m-Acetyl chlorophosphonazo peak point；
S204：According to the bimodal breadth of spectral line value ratio parameter at voltage, the fundamental wave of current signal, harmonic wave and m-Acetyl chlorophosphonazo peak point
The value of β, and measured signal fundamental wave, humorous is calculated according to the Rife-Vincent (III) window FFT double peak of spectral lines interpolation correction formula
The value of the frequency correction factor α at ripple and m-Acetyl chlorophosphonazo peak 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：Tested voltage, the fundamental wave of current signal, harmonic wave be worth to according to the value of correction factor α and correction factor g (α)
Amplitude, phase and frequency with m-Acetyl chlorophosphonazo.
4. the electric energy gauging method based on R-V (III) window FFT double peak interpolations according to claim 2, it is characterised in that logical
Computing formula with frequency correction factor α is：
+ 0.18754082 β of α=3.04072904 β3+0.08804904·β5+0.05199160·β7。
5. the electric energy gauging method based on R-V (III) window FFT double peak interpolations according to claim 4, it is characterised in that logical
Computing formula with amplitude correction coefficient g (α) is：
G (α)=1.08537118+0.30341579 α2+0.04563736·α4+0.00514598·α6。
6. the electric energy gauging method based on R-V (III) window FFT double peak interpolations according to claim 5, it is characterised in that logical
With the correction formula of amplitude it is：
A=(y1+y2)·g(α)/N；
Wherein, y1And y2Respectively the maximum spectral line of spectrum curve peak point arranged on left and right sides and time maximum spectral line are respectively kth1Bar
And kth2The amplitude of bar spectral line, N are sampled data length.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610882937.5A CN106526312A (en) | 2016-10-10 | 2016-10-10 | Electric energy metering method based on R-V(III) window FFT double-peak interpolation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610882937.5A CN106526312A (en) | 2016-10-10 | 2016-10-10 | Electric energy metering method based on R-V(III) window FFT double-peak interpolation |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106526312A true CN106526312A (en) | 2017-03-22 |
Family
ID=58331258
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610882937.5A Pending CN106526312A (en) | 2016-10-10 | 2016-10-10 | Electric energy metering method based on R-V(III) window FFT double-peak interpolation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106526312A (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107179439A (en) * | 2017-06-26 | 2017-09-19 | 贵州电网有限责任公司贵阳供电局 | Electric energy gauging method based on five peak interpolations of Cosine Window FFT tri- |
| 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 |
| CN109100674A (en) * | 2018-05-28 | 2018-12-28 | 国网江西省电力有限公司电力科学研究院 | A kind of metering method based on distribution type electric energy, error detection method and device |
| CN109581045A (en) * | 2018-12-24 | 2019-04-05 | 中国船舶重工集团公司第七〇九研究所 | A kind of m-Acetyl chlorophosphonazo power measurement method meeting IEC standard frame |
| CN109683016A (en) * | 2019-01-18 | 2019-04-26 | 江苏理工学院 | A kind of harmonic amplitude analysis method |
| CN110095650A (en) * | 2019-05-05 | 2019-08-06 | 三峡大学 | The complicated harmonic detecting analysis method of four spectral line interpolation FFTs based on five Rife-Vincent (I) windows |
| CN115598416A (en) * | 2022-09-16 | 2023-01-13 | 珠海多创科技有限公司（Cn） | Method and system for processing station area sampling signal, storage medium and computer equipment |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH01143587A (en) * | 1987-11-30 | 1989-06-06 | Matsushita Electric Ind Co Ltd | Digital convergence correcting device |
| CN101261292A (en) * | 2008-04-14 | 2008-09-10 | 湖南大学 | Base wave and harmonic detection method based on fiver item Rife-Vincent(1)window double spectral line interpolation FFT |
| CN103353550A (en) * | 2013-04-24 | 2013-10-16 | 武汉大学 | Method for measuring signal frequency and harmonic parameters of electric power system |
| CN103399204A (en) * | 2013-07-31 | 2013-11-20 | 武汉大学 | Rife-Vincent (II) window interpolation FFT (Fast Fourier Transform)-based harmonic and inter-harmonic detection method |
| CN104897960A (en) * | 2015-06-15 | 2015-09-09 | 中南民族大学 | Harmonic rapid analysis method and system based on windowing four-spectral-line interpolation FFT |
| CN105486921A (en) * | 2016-01-22 | 2016-04-13 | 武汉大学 | Kaiser third-order mutual convolution window triple-spectrum-line interpolation harmonic wave and inter-harmonic wave detection method |
-
2016
- 2016-10-10 CN CN201610882937.5A patent/CN106526312A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH01143587A (en) * | 1987-11-30 | 1989-06-06 | Matsushita Electric Ind Co Ltd | Digital convergence correcting device |
| CN101261292A (en) * | 2008-04-14 | 2008-09-10 | 湖南大学 | Base wave and harmonic detection method based on fiver item Rife-Vincent(1)window double spectral line interpolation FFT |
| CN103353550A (en) * | 2013-04-24 | 2013-10-16 | 武汉大学 | Method for measuring signal frequency and harmonic parameters of electric power system |
| CN103399204A (en) * | 2013-07-31 | 2013-11-20 | 武汉大学 | Rife-Vincent (II) window interpolation FFT (Fast Fourier Transform)-based harmonic and inter-harmonic detection method |
| CN104897960A (en) * | 2015-06-15 | 2015-09-09 | 中南民族大学 | Harmonic rapid analysis method and system based on windowing four-spectral-line interpolation FFT |
| CN105486921A (en) * | 2016-01-22 | 2016-04-13 | 武汉大学 | Kaiser third-order mutual convolution window triple-spectrum-line interpolation harmonic wave and inter-harmonic wave detection method |
Non-Patent Citations (2)
| Title |
|---|
| 张立强等: "基于双CPU电能质量监测仪的谐波分析算法研究", 《电测与仪表》 * |
| 徐文佳 等: "计及间谐波的电能计量方法", 《电网技术》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107179439A (en) * | 2017-06-26 | 2017-09-19 | 贵州电网有限责任公司贵阳供电局 | Electric energy gauging method based on five peak interpolations of Cosine Window FFT tri- |
| 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 |
| CN109100674A (en) * | 2018-05-28 | 2018-12-28 | 国网江西省电力有限公司电力科学研究院 | A kind of metering method based on distribution type electric energy, error detection method and device |
| CN109100674B (en) * | 2018-05-28 | 2020-08-07 | 国网江西省电力有限公司电力科学研究院 | Metering method based on distributed electric energy, error detection method and device |
| CN109581045A (en) * | 2018-12-24 | 2019-04-05 | 中国船舶重工集团公司第七〇九研究所 | A kind of m-Acetyl chlorophosphonazo power measurement method meeting IEC standard frame |
| CN109581045B (en) * | 2018-12-24 | 2020-12-11 | 中国船舶重工集团公司第七一九研究所 | Inter-harmonic power metering method meeting IEC standard framework |
| CN109683016A (en) * | 2019-01-18 | 2019-04-26 | 江苏理工学院 | A kind of harmonic amplitude analysis method |
| CN110095650A (en) * | 2019-05-05 | 2019-08-06 | 三峡大学 | The complicated harmonic detecting analysis method of four spectral line interpolation FFTs based on five Rife-Vincent (I) windows |
| CN115598416A (en) * | 2022-09-16 | 2023-01-13 | 珠海多创科技有限公司（Cn） | Method and system for processing station area sampling signal, storage medium and computer equipment |
| CN115598416B (en) * | 2022-09-16 | 2024-01-30 | 珠海多创科技有限公司 | Processing method, system, storage medium and computer equipment for area sampling signal |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106526312A (en) | Electric energy metering method based on R-V(III) window FFT double-peak interpolation | |
| 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 | |
| Sun et al. | Harmonic contribution evaluation based on the distribution-level PMUs | |
| CN107179439A (en) | Electric energy gauging method based on five peak interpolations of Cosine Window FFT tri- | |
| CN102662106A (en) | Method for measuring electric energy of harmonic power network | |
| Li et al. | Research of harmonic distortion power for harmonic source detection | |
| CN106154040B (en) | One kind is for calculating power grid single-point equivalent time Equivalent Harmonic impedance methodologies | |
| CN105842555B (en) | A kind of Multi-harmonic Sources divisions of responsibility method based on measured data | |
| CN104897961A (en) | Three spectral line interpolation FFT harmonic wave analysis method and system based on multiplication window function | |
| CN107359612A (en) | A kind of comprehensive estimation method of quality of power supply to power distribution network energy consumption | |
| CN103105529A (en) | Harmonic wave electric energy measuring system based on parameter analysis and control method thereof | |
| CN101227092A (en) | Method for recognizing harmonic source in electric power distribution network | |
| CN110045186A (en) | System harmonic impedance estimation method and system at electric system points of common connection | |
| CN103778341B (en) | Method for estimating probabilistic load flow of harmonic pollution user | |
| Chen et al. | A modified harmonic pricing scheme for customers based on quantifying the harmonic comprehensive contribution | |
| Orts-Grau et al. | Discussion on useless active and reactive powers contained in the IEEE standard 1459 | |
| CN102185314A (en) | Power quality comprehensive evaluation method | |
| CN110333409A (en) | The online test method of distribution transformer rated capacity | |
| Cho et al. | A waveform distortion evaluation method based on a simple half-cycle RMS calculation | |
| CN105866585A (en) | Method for harmonic source identification and liability sharing based on distortion power | |
| CN105006820A (en) | Quantitative analysis method for influence of electric energy quality on comprehensive energy consumption of power transmission line | |
| CN101667730A (en) | Method for monitoring and managing electrical energy quality of electric distribution network | |
| Berrisford | The harmonic impact project—IEEE-1459 power definitions trialed in revenue meters | |
| CN106771552B (en) | A kind of distortion power metering method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | 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 |
Application publication date: 20170322 |
|
| RJ01 | Rejection of invention patent application after publication |