CN107390022A - Electric energy gauging method based on discrete spectrum correction - Google Patents
Electric energy gauging method based on discrete spectrum correction Download PDFInfo
- Publication number
- CN107390022A CN107390022A CN201710729394.8A CN201710729394A CN107390022A CN 107390022 A CN107390022 A CN 107390022A CN 201710729394 A CN201710729394 A CN 201710729394A CN 107390022 A CN107390022 A CN 107390022A
- Authority
- CN
- China
- Prior art keywords
- mrow
- harmonic
- msub
- wave
- spectrum
- 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
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
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Measurement Of Current Or Voltage (AREA)
Abstract
The invention discloses a kind of electric energy gauging method based on discrete spectrum correction, including the voltage and current signal of digital quantity is converted into using a length of N of window tetra- three rank windows progress time domain truncations of Nuttall；The spectrum barycenter and fundamental frequency of fundamental wave are calculated using spectrum barycenter formula；Harmonic spike position is estimated, calculates harmonic spectrum barycenter and the maximum spectrum peak position and offset of fundamental wave and each harmonic；Principal wave harmonic wave amplitude correction coefficient is determined, and calculates tested voltage, the fundamental wave of electric current, harmonic amplitude and phase；The electrical parameter of fundamental wave harmonic is calculated, completes electric energy metrical.Therefore the inventive method need not carry out peak value searching, offset need not be calculated according to the type of institute's windowed function, and mutually eliminate zero-crossing examination, it is directly calculated using fft analysis method, eliminate the window spectrum computing of complexity, calculate and analyze obtained frequency accuracy to greatly improve, and program realizes that simply real-time is good.
Description
Technical field
Present invention relates particularly to a kind of electric energy gauging method based on discrete spectrum correction.
Background technology
With the development and the improvement of people's living standards of economic technology, electric energy has become the daily production of people and life
Essential secondary energy sources in work, production and life to people bring endless facility.
Harmonious Waves in Power Systems has significant effect to the accuracy of electric energy metrical, using metering fundamental energy harmonic respectively
The metering of electric energy is a kind of relatively rational mode, and with the development of technology, requirement of the market to principal wave harmonic wave measuring accuracy is more next
It is higher.According to the current development in electric energy metrical field, the design of harmonic wave metering mainly has two kinds：A kind of is special metering
Chip, another kind are the special AD+MCU of DSP+.Former scheme cost is relatively low, is mainly used in the design of common list, three-phase meter
In；The latter's cost is of a relatively high, and main application is mainly used in regional grid, provincial power network, light with Source of Gateway Meter platform, Source of Gateway Meter
The stoichiometric point that volt, wind-powered electricity generation, water power and electric railway even load frequently change, the field changeable in face of the big load of this harmonic content
Close, the program is more flexible compared to the former function, it is easier to meets the needs of complicated.
Source of Gateway Meter does not require nothing more than higher electric energy metrical precision, and require lower fortune as a ripe product
Row power consumption.And many algorithms have been proposed currently for the high-precision measuring of fundamental wave, the amplitude of harmonic wave, phase and frequency, this its
In have quite a few algorithm because realizing complicated difficult to accomplish requirement of real-time in embedded systems, and realize it is more complicated,
System operation power consumption can be improved accordingly, thus finds a kind of measuring accuracy height, realized simply, and time and spatial complexity
It can meet that the algorithm of existing hardware platform is particularly important.
The content of the invention
It is an object of the invention to provide a kind of precision is of a relatively high, real-time is good and algorithm is simply based on discrete spectrum
The electric energy gauging method of correction.
This electric energy gauging method based on discrete spectrum correction provided by the invention, comprises the following steps：
S1. it is digital quantity by sample conversion by voltage and current signal, using a length of N of window tetra- three ranks of Nuttall
Window carries out time domain truncation to the digital quantity after conversion, obtains blocking the spectrum signal Y (k) of rear signal；The spectrum signal Y (k)
Including voltage spectrum value YuAnd current spectrum value Y (k)i(k)；
S2. the voltage spectrum value Y obtained according to step S1u(k) the spectrum barycenter of fundamental wave, is calculated using spectrum barycenter formula
ksc1With fundamental frequency f0；
S3. the fundamental frequency f obtained using step S20Estimate harmonic spike positionAnd calculated using spectrum barycenter formula
Obtain harmonic spectrum barycenter ksch, and the maximum spectrum peak position k of fundamental wave and each harmonic is further calculatedhWith offset λh；
S4. principal wave harmonic wave amplitude correction coefficient g (λ are determined according to the offset obtained in step S3h), and be calculated by
Survey voltage, the fundamental wave of electric current, harmonic amplitude and phase；
S5. tested voltage, the fundamental wave of electric current, harmonic amplitude and the phase calculation fundamental wave harmonic obtained according to step S4
Electrical parameter, so as to complete the metering of electric energy.
The spectrum signal Y (k) for obtaining blocking rear signal described in step S1, after specially being blocked using fft analysis
The spectrum signal Y (k) of signal.
Tetra- three rank windows of Nuttall of a length of N of window described in step S1, the specially main lobe of tetra- three rank windows of Nuttall
Width isThere are 8 spectral lines in main lobe.
The a length of N of described window tetra- three rank windows of Nuttall meet following formula：
N=1,2 in formula ..., N-1；bmFor window function coefficient, meet
Fundamental wave spectrum barycenter k is calculated described in step S2sc1With fundamental frequency f0, specially calculated using following steps
Obtain：
A. fundamental wave spectrum barycenter k is calculated according to spectrum barycenter formulasc1：
In formulaAndFor downward bracket function, fmin is
The minimum value of signal fundamental frequency, fmax be signal fundamental frequency maximum, Δ f=fs/ N, fs are the sampling of analog-to-digital conversion
Frequency, N are that the analysis of fft algorithm is counted；
B. the fundamental wave obtained according to step A composes barycenter ksc1Measure fundamental frequency f0：
f0=ksc1*Δf。
Harmonic spectrum barycenter k is calculated described in step S3sch, fundamental wave and each harmonic maximum spectrum peak position khWith it is inclined
Shifting amount λh, specially it is calculated using following steps：
A. using following formula estimation harmonic spectrum peak
Wherein h value is 2,3 ..., M；M is the maximum overtone order of signal；
B. obtained according to step aHarmonic spectrum barycenter k is calculated using following formulasch：
C. the harmonic spectrum barycenter k obtained according to step bsch, fundamental wave and each harmonic are calculated most using equation below
Big spectrum peak position kh：
kh=round (ksch)
Round () is the function that rounds up in formula；
D. the maximum spectrum peak position k obtained according to step chOffset λ is calculatedh：
λh=kh-ksch。
Determination principal wave harmonic wave amplitude correction coefficient g (λ described in step S4h), specially calculated using equation below：
Wherein p1~p5It is constant.
Tested voltage, the fundamental wave of electric current, harmonic amplitude and phase is calculated described in step S4, specially using as follows
Updating formula is corrected to amplitude and phase：
A=| Y (kh)|g(λh)
A is revised amplitude in formula, and θ is revised phase, angle (Y (kh)) represent Y (kh) phase value.
The electrical parameter of calculating fundamental wave harmonic described in step S5, specially calculates fundamental wave harmonic using following formula
Electrical parameter：
Define principal wave harmonic wave voltage magnitude Uh, voltage-phase θ uh, current amplitude Ih, current phase θ ih, h=1,2 ...,
63；
Phase difference：θh=θ uh-θih
Principal wave harmonic wave active power：Ph=UhIh cos(θuh-θih)
Principal wave harmonic wave reactive power：Qh=UhIh sin(θuh-θih)
Full wave voltage：
Full-wave electric current：
All-wave is active：
All-wave is idle：
Fundamental active electric flux：EP=P1T, t are the time；
Fundamental wave reactive power electric flux：EQ=Q1T, t are the time；
Harmonic wave active energy amount：T is the time；
Harmonic wave reactive energy amount：T is the time；
All-wave active energy amount：EP=Pt, t are the time；
All-wave reactive energy amount：EQ=Qt, t are the time.
This electric energy gauging method based on discrete spectrum correction provided by the invention, is accurately looked for using spectrum centroid method first
Go out signal spectrum peak value, to eliminate caused spectral leakage and fence response under non-synchronous sampling；Further to voltage, electric current
Spectrum value be corrected using amplitude correction formula and phase only pupil filter formula, finally try to achieve all-wave, fundamental wave harmonic electric parameter
And electric flux；Therefore the inventive method need not carry out peak value searching, it is not necessary to calculated and offset according to the type of institute's windowed function
Amount, and zero-crossing examination is mutually eliminated, it is directly calculated using fft analysis method, eliminates the window spectrum computing of complexity, meter
Calculate and analyze obtained frequency accuracy to greatly improve, and program realizes that simply real-time is good.
Brief description of the drawings
Fig. 1 is the method flow diagram of the inventive method.
Fig. 2 is the spectrum centroid calculation schematic diagram of the inventive method.
Fig. 3 is the correction coefficient matched curve of the inventive method.
Fig. 4 is the harmonic frequency error precision of the inventive method.
Fig. 5 is the harmonic amplitude error precision of the inventive method.
Fig. 6 is the harmonic phase error precision of the inventive method.
Embodiment
It is the method flow diagram of the inventive method as shown in Figure 1：It is provided by the invention it is this based on discrete spectrum correction
Electric energy gauging method, comprise the following steps：
S1. it is digital quantity by sample conversion by voltage and current signal, using a length of N of window tetra- three ranks of Nuttall
Window carries out time domain truncation to the digital quantity after conversion, obtains blocking the spectrum signal Y (k) of rear signal using fft analysis；The frequency
Spectrum signal Y (k) includes voltage spectrum value YuAnd current spectrum value Y (k)i(k)；
In the specific implementation, the main lobe width of preferably tetra- three rank windows of Nuttall isThere are 8 spectral lines in main lobe；
The a length of N of described window tetra- three rank windows of Nuttall meet following formula：
N=1,2 in formula ..., N-1；bmFor window function coefficient, meetb0=
0.338946, b1=0.481973, b2=0.161054, b3=0.018027；
Described application window function is to do following processing to discrete voltage u (n), current signal i (n)：U ' (n)=u (n) * w
(n), i ' (n)=i (n) * w (n)；U ' (n) and i ' (n) obtains discrete spectrum Y after fft analysisuAnd Y (k)i(k), for ease of point
Analysis, it is used uniformly Y (k) expressions；
S2. the voltage spectrum value Y obtained according to step S1u(k), using spectrum barycenter formula to calculate, (Fig. 2 is the inventive method
Spectrum centroid calculation schematic diagram) obtain the spectrum barycenter k of fundamental wavesc1With fundamental frequency f0；Specially calculated using following steps
Arrive：
A. fundamental wave spectrum barycenter k is calculated according to spectrum barycenter formulasc1：
In formulaAndFor downward bracket function, fmin is
The minimum value of signal fundamental frequency, specially 45Hz；Fmax be signal fundamental frequency maximum, specially 65Hz；Δ f=
fs/ N, fs are the sample frequency of analog-to-digital conversion, and N is that the analysis of fft algorithm is counted；
B. the fundamental wave obtained according to step A composes barycenter ksc1Measure fundamental frequency f0：
f0=ksc1*Δf
S3. the fundamental frequency f obtained using step S20Estimate harmonic spike positionAnd calculated using spectrum barycenter formula
Obtain harmonic spectrum barycenter ksch, and the maximum spectrum peak position k of fundamental wave and each harmonic is further calculatedhWith offset λh；Tool
Body is to be calculated using following steps：
A. using following formula estimation harmonic spectrum peak
Wherein h value is 1,2,3 ..., M；M is the maximum overtone order of signal；
B. obtained according to step aHarmonic spectrum barycenter k is calculated using following formulasch：
C. the harmonic spectrum barycenter k obtained according to step bsch, fundamental wave and each harmonic are calculated most using equation below
Big spectrum peak position kh：
kh=round (ksch)
Round () is the function that rounds up in formula；
D. offset λ is calculated in the maximum spectrum peak position kh obtained according to step ch：
λh=kh-ksch
S4. principal wave harmonic wave amplitude correction coefficient g (λ are determined according to the offset obtained in step S3h), and be calculated by
Survey voltage, the fundamental wave of electric current, harmonic amplitude and phase；
Calculating g (λ are carried out using equation belowh)：
Wherein p1~p5It is constant；Preferable p1=0.2073, p2=-0.01546, p3=0.9941, p4=-
0.000251, p5=2.9503；
Amplitude and phase are corrected using following updating formula：
A=| Y (kh)|g(λh)
A is revised amplitude in formula, and θ is revised phase, angle (Y (kh)) represent Y (kh) phase value；
The derivation of above-mentioned amplitude correction coefficient described briefly below, amplitude correction formula and phase only pupil filter formula：
Consider single-frequency signals y (n) with sample frequency f firstsThe discrete-time signal that uniform sampling obtains is：
A is amplitude in formula, f0For frequency,For initial phase；
Y is obtained to y (n) plus tetra- three rank windows of Nuttallw(n)=y (n) w (n), yw(n) continuous fourier transform is
Discrete sampling is carried out to above formula, and ignores negative frequency point-f0Going out the secondary lobe at frequency peak influences, and obtains signal after adding window
The expression formula of discrete Fourier transform is：
Make f0=kscΔ f, then kscFor the spectrum barycenter of signal；As shown in Fig. 2 k*For close to spectrum barycenter kscMaximum spectral peak
Position, therefore can obtain：
Make λ=k*-ksc, λ is offset, λ ∈ [- 0.5,0.5]；
Further derive that the discrete spectrum approximation to function of tetra- three rank windows of Nuttall is as follows：
Amplitude correction formula, which can to sum up be obtained, is：
In formulaFor window function offset amplitude；Curve map it is as shown in Figure 3；Therefrom
It can be seen that the curve can use fitting of a polynomial, to reduce operand, first approached using fourth order polynomial fitting formulaCurve, fitting formula are as follows：
G (λ)=p1λ4+p2|λ|3+p3λ2+p4|λ|+p5
Wherein p1=0.2073, p2=-0.01546, p3=0.9941, p4=-0.000251, p5=2.9503；
Therefore amplitude formula can be obtained is：
A=| Y (k*)|g(λ)
Phase only pupil filter formula is：
S5. tested voltage, the fundamental wave of electric current, harmonic amplitude and the phase calculation fundamental wave harmonic obtained according to step S4
Electrical parameter, so as to complete the metering of electric energy；The electrical parameter of fundamental wave harmonic is specially calculated using following formula：
Define principal wave harmonic wave voltage magnitude Uh, voltage-phase θ uh, current amplitude Ih, current phase θ ih, h=1,2 ...,
63；
Phase difference：θh=θ uh-θih
Principal wave harmonic wave active power：Ph=UhIh cos(θuh-θih)
Principal wave harmonic wave reactive power：Qh=UhIh sin(θuh-θih)
Full wave voltage：
Full-wave electric current：
All-wave is active：
All-wave is idle：
Fundamental active electric flux：EP=P1T, t are the time；
Fundamental wave reactive power electric flux：EQ=Q1T, t are the time；
Harmonic wave active energy amount：T is the time；
Harmonic wave reactive energy amount：T is the time；
All-wave active energy amount：EP=Pt, t are the time；
All-wave reactive energy amount：EQ=Qt, t are the time.
The electric energy gauging method provided below by way of verification experimental verification invention：
The inventive method be applied to Source of Gateway Meter (ammeter specifications and models are DTSD341-MA2,3 × 57.7V, 3 × 1.5 (6) A,
20000imp/kWh, 50Hz), program is implemented in meter platform Blackfin BF533 DSP process chips.
The fundamental active pulse elementary error precision of this method is as shown in Figure 1：
The fundamental active pulse elementary error precision of table 1
Harmonic pulse Setup Experiments：
Fundamental voltage Un=57.7V；Fundamental current In is 0.5*Imax, as In=2*Ib=3A；Phasor power factor
For 1.0；Fundamental frequency is 50Hz；Single harmonic component is superimposed on the basis of fundamental wave, overtone order highest is added to 48 times；Harmonic voltage
Amplitude 0.05*Un；Harmonic current 0.4*In；Harmonic power factor is 0.5L；Harmonic pulse constant 1000000imp/kWh；
It is superimposed below for Source of Gateway Meter under single harmonic component, the error that harmonic wave goes out pulse is as shown in table 2：
The harmonic wave of table 2 goes out the error signal table of pulse
From upper table it can be found that Source of Gateway Meter measured data greatly differs from each other with theoretical value, this is due to that experiment condition is limited, is made
Obtaining experimental data can not reach best, and nonetheless, upper table error information has met GB/T17215.302-2013 pair harmonic wave
Electric energy metrical precision peak demand, and entered every strict test, fully meet standard requirement.Harmonic pulse error is main
Influenceed by caused by hardware sample circuit and transformer than difference and difference, in table harmonic error data be ammeter through amplitude and
Data are surveyed after phasing.After experiment condition permission, above-mentioned data precision will have the lifting of matter.
MATLAB simulation results as shown in Fig. 3~Fig. 6, respectively correction coefficient matched curve, harmonic frequency error precision,
Harmonic amplitude error precision harmonic phase error precision；It can be seen that the present invention, which carries algorithm, can realize that harmonic wave is believed
Number high-precision measuring.
Claims (9)
1. a kind of electric energy gauging method based on discrete spectrum correction, comprises the following steps：
S1. it is digital quantity by sample conversion by voltage and current signal, using a length of N of window tetra- three rank windows pair of Nuttall
Digital quantity after conversion carries out time domain truncation, obtains blocking the spectrum signal Y (k) of rear signal；The spectrum signal Y (k) includes
Voltage spectrum value YuAnd current spectrum value Y (k)i(k)；
S2. the voltage spectrum value Y obtained according to step S1u(k) the spectrum barycenter k of fundamental wave, is calculated using spectrum barycenter formulasc1With
Fundamental frequency f0；
S3. the fundamental frequency f obtained using step S20Estimate harmonic spike positionAnd it is calculated using spectrum barycenter formula
Harmonic spectrum barycenter ksch, and the maximum spectrum peak position k of fundamental wave and each harmonic is further calculatedhWith offset λh；
S4. principal wave harmonic wave amplitude correction coefficient g (λ are determined according to the offset obtained in step S3h), and tested electricity is calculated
Pressure, the fundamental wave of electric current, harmonic amplitude and phase；
S5. the electricity ginseng of tested voltage, the fundamental wave of electric current, harmonic amplitude and the phase calculation fundamental wave harmonic obtained according to step S4
Number, so as to complete the metering of electric energy.
2. the electric energy gauging method according to claim 1 based on discrete spectrum correction, it is characterised in that described in step S1
The spectrum signal Y (k) for obtaining blocking rear signal, specially obtain blocking the spectrum signal Y (k) of rear signal using fft analysis.
3. the electric energy gauging method according to claim 2 based on discrete spectrum correction, it is characterised in that described in step S1
The a length of N of window tetra- three rank windows of Nuttall, the main lobe width of specially tetra- three rank windows of Nuttall isHave in main lobe
8 spectral lines.
4. the electric energy gauging method according to claim 3 based on discrete spectrum correction, it is characterised in that described window length
Meet following formula for N tetra- three rank windows of Nuttall：
<mrow>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mn>3</mn>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
<msub>
<mi>b</mi>
<mi>m</mi>
</msub>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>&pi;</mi>
<mi>n</mi>
<mo>*</mo>
<mi>m</mi>
<mo>/</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
</mrow>
N=1,2 in formula ..., N-1；bmFor window function coefficient, meet
5. the electric energy gauging method according to claim 4 based on discrete spectrum correction, it is characterised in that described in step S2
Fundamental wave spectrum barycenter k is calculatedsc1With fundamental frequency f0, specially it is calculated using following steps：
A. fundamental wave spectrum barycenter k is calculated according to spectrum barycenter formulasc1：
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>min</mi>
</msub>
<mo>-</mo>
<mi>B</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<msub>
<mi>k</mi>
<mi>max</mi>
</msub>
<mo>+</mo>
<mi>B</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</munderover>
<msup>
<mi>kY</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>min</mi>
</msub>
<mo>-</mo>
<mi>B</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<msub>
<mi>k</mi>
<mi>max</mi>
</msub>
<mo>+</mo>
<mi>B</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</munderover>
<msup>
<mi>Y</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
In formulaAndFor downward bracket function, fmin is signal base
The minimum value of wave frequency rate, fmax be signal fundamental frequency maximum, Δ f=fs/ N, fs be analog-to-digital conversion sample frequency, N
Counted for the analysis of fft algorithm；
B. the fundamental wave obtained according to step A composes barycenter ksc1Measure fundamental frequency f0：
f0=ksc1*Δf。
6. the electric energy gauging method according to claim 5 based on discrete spectrum correction, it is characterised in that described in step S3
Harmonic spectrum barycenter k is calculatedsch, fundamental wave and each harmonic maximum spectrum peak position khWith offset λh, specially using such as
Lower step is calculated：
A. using following formula estimation harmonic spectrum peak
Wherein h value is 2,3 ..., M；M is the maximum overtone order of signal；
B. obtained according to step aHarmonic spectrum barycenter k is calculated using following formulasch：
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<msubsup>
<mi>k</mi>
<mi>h</mi>
<mo>*</mo>
</msubsup>
<mo>-</mo>
<mi>B</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<msubsup>
<mi>k</mi>
<mi>h</mi>
<mo>*</mo>
</msubsup>
<mo>+</mo>
<mi>B</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</munderover>
<msup>
<mi>kY</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<msubsup>
<mi>k</mi>
<mi>h</mi>
<mo>*</mo>
</msubsup>
<mo>-</mo>
<mi>B</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<msubsup>
<mi>k</mi>
<mi>h</mi>
<mo>*</mo>
</msubsup>
<mo>+</mo>
<mi>B</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</munderover>
<msup>
<mi>Y</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
C. the harmonic spectrum barycenter k obtained according to step bsch, fundamental wave is calculated using equation below and the maximum of each harmonic is composed
Peak position kh：
kh=round (ksch)
Round () is the function that rounds up in formula；
D. the maximum spectrum peak position k obtained according to step chOffset λ is calculatedh：
λh=kh-ksch。
7. the electric energy gauging method according to claim 6 based on discrete spectrum correction, it is characterised in that described in step S4
Determination principal wave harmonic wave amplitude correction coefficient g (λh), specially calculated using equation below：
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&lambda;</mi>
<mi>h</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>&lambda;</mi>
<mi>h</mi>
<mn>4</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<msubsup>
<mi>&lambda;</mi>
<mi>h</mi>
<mn>3</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>p</mi>
<mn>3</mn>
</msub>
<msubsup>
<mi>&lambda;</mi>
<mi>h</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>p</mi>
<mn>4</mn>
</msub>
<msub>
<mi>&lambda;</mi>
<mi>h</mi>
</msub>
<mo>+</mo>
<msub>
<mi>p</mi>
<mn>5</mn>
</msub>
</mrow>
Wherein p1~p5It is constant.
8. the electric energy gauging method according to claim 7 based on discrete spectrum correction, it is characterised in that described in step S4
Tested voltage, the fundamental wave of electric current, harmonic amplitude and phase is calculated, specially using following updating formula to amplitude and phase
Position is corrected：
A=| Y (kh)|g(λh)
<mrow>
<mi>&theta;</mi>
<mo>=</mo>
<mi>a</mi>
<mi>n</mi>
<mi>g</mi>
<mi>l</mi>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>Y</mi>
<mo>(</mo>
<msub>
<mi>k</mi>
<mi>h</mi>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&lambda;</mi>
<mi>h</mi>
</msub>
<mi>&pi;</mi>
<mo>+</mo>
<mfrac>
<mi>&pi;</mi>
<mn>2</mn>
</mfrac>
</mrow>
A is revised amplitude in formula, and θ is revised phase, angle (Y (kh)) represent Y (kh) phase value.
9. the electric energy gauging method according to claim 8 based on discrete spectrum correction, it is characterised in that described in step S5
Calculating fundamental wave harmonic electrical parameter, specially using following formula calculate fundamental wave harmonic electrical parameter：
Define principal wave harmonic wave voltage magnitude Uh, voltage-phase θ uh, current amplitude Ih, current phase θ ih, h=1,2 ..., 63；
Phase difference：θh=θ uh-θih
Principal wave harmonic wave active power：Ph=UhIhcos(θuh-θih)
Principal wave harmonic wave reactive power：Qh=UhIhsin(θuh-θih)
Full wave voltage：
Full-wave electric current：
All-wave is active：
All-wave is idle：
Fundamental active electric flux：EP=P1T, t are the time；
Fundamental wave reactive power electric flux：EQ=Q1T, t are the time；
Harmonic wave active energy amount：T is the time；
Harmonic wave reactive energy amount：T is the time；
All-wave active energy amount：EP=Pt, t are the time；
All-wave reactive energy amount：EQ=Qt, t are the time.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710729394.8A CN107390022B (en) | 2017-08-23 | 2017-08-23 | Electric energy metering method based on discrete spectrum correction |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710729394.8A CN107390022B (en) | 2017-08-23 | 2017-08-23 | Electric energy metering method based on discrete spectrum correction |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107390022A true CN107390022A (en) | 2017-11-24 |
| CN107390022B CN107390022B (en) | 2020-04-14 |
Family
ID=60354344
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710729394.8A Active CN107390022B (en) | 2017-08-23 | 2017-08-23 | Electric energy metering method based on discrete spectrum correction |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107390022B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110795070A (en) * | 2019-05-24 | 2020-02-14 | 杭州海兴电力科技股份有限公司 | Virtual gateway table platform and construction method |
| CN112230054A (en) * | 2019-12-20 | 2021-01-15 | 青岛鼎信通讯股份有限公司 | Harmonic content and power calculation method |
| CN113358930A (en) * | 2021-08-09 | 2021-09-07 | 南京派格测控科技有限公司 | Harmonic wave test system, device and method based on signal deviation |
| CN114167198A (en) * | 2021-10-18 | 2022-03-11 | 国网山东省电力公司平原县供电公司 | Method and platform for measuring synchronous line loss data |
| CN114184839A (en) * | 2021-12-10 | 2022-03-15 | 江苏集萃智能集成电路设计技术研究所有限公司 | ADC frequency spectrum testing method for inhibiting signal source jitter |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010144857A3 (en) * | 2004-12-21 | 2011-03-24 | Montana Tech Of The University Of Montana | Method of detecting system function by measuring frequency response |
| CN102445595A (en) * | 2011-09-23 | 2012-05-09 | 重庆大学 | Real-time measuring method for time-varying power of electrical power system |
| CN104391178A (en) * | 2014-12-05 | 2015-03-04 | 国家电网公司 | Time shift phase difference steady harmonic signal correction method based on Nuttall window |
| CN105445541A (en) * | 2015-12-24 | 2016-03-30 | 合肥工业大学 | Method for adaptively calculating power under arbitrary frequencies |
| CN106405229A (en) * | 2016-08-30 | 2017-02-15 | 威胜集团有限公司 | Fundamental wave and harmonic wave electric energy metering method |
-
2017
- 2017-08-23 CN CN201710729394.8A patent/CN107390022B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010144857A3 (en) * | 2004-12-21 | 2011-03-24 | Montana Tech Of The University Of Montana | Method of detecting system function by measuring frequency response |
| CN102445595A (en) * | 2011-09-23 | 2012-05-09 | 重庆大学 | Real-time measuring method for time-varying power of electrical power system |
| CN104391178A (en) * | 2014-12-05 | 2015-03-04 | 国家电网公司 | Time shift phase difference steady harmonic signal correction method based on Nuttall window |
| CN105445541A (en) * | 2015-12-24 | 2016-03-30 | 合肥工业大学 | Method for adaptively calculating power under arbitrary frequencies |
| CN106405229A (en) * | 2016-08-30 | 2017-02-15 | 威胜集团有限公司 | Fundamental wave and harmonic wave electric energy metering method |
Non-Patent Citations (1)
| Title |
|---|
| 曾金芳: "谱质心在电力谐波分析中的应用", 《中国电机工程学报》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110795070A (en) * | 2019-05-24 | 2020-02-14 | 杭州海兴电力科技股份有限公司 | Virtual gateway table platform and construction method |
| CN110795070B (en) * | 2019-05-24 | 2023-07-25 | 杭州海兴电力科技股份有限公司 | Virtual gateway table platform and construction method |
| CN112230054A (en) * | 2019-12-20 | 2021-01-15 | 青岛鼎信通讯股份有限公司 | Harmonic content and power calculation method |
| CN113358930A (en) * | 2021-08-09 | 2021-09-07 | 南京派格测控科技有限公司 | Harmonic wave test system, device and method based on signal deviation |
| CN114167198A (en) * | 2021-10-18 | 2022-03-11 | 国网山东省电力公司平原县供电公司 | Method and platform for measuring synchronous line loss data |
| CN114167198B (en) * | 2021-10-18 | 2024-03-01 | 国网山东省电力公司平原县供电公司 | Method and platform for measuring synchronous line loss data |
| CN114184839A (en) * | 2021-12-10 | 2022-03-15 | 江苏集萃智能集成电路设计技术研究所有限公司 | ADC frequency spectrum testing method for inhibiting signal source jitter |
| CN114184839B (en) * | 2021-12-10 | 2024-03-01 | 江苏集萃智能集成电路设计技术研究所有限公司 | ADC frequency spectrum testing method for inhibiting signal source jitter |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107390022B (en) | 2020-04-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107390022A (en) | Electric energy gauging method based on discrete spectrum correction | |
| CN102288807B (en) | Method for measuring electric network voltage flicker | |
| CN109633262A (en) | Three phase harmonic electric energy gauging method, device based on composite window multiline FFT | |
| CN108318852B (en) | Square wave influence test method for intelligent electric energy meter | |
| CN106405229B (en) | A kind of fundamental wave and harmonic wave electric energy gauging method | |
| CN104360235B (en) | Method for positioning voltage sag source of complex power distribution network including DGs | |
| CN105223434A (en) | Quality of power supply mixing detection method | |
| CN104502707A (en) | Synchronized phasor measurement method for electrical power system based on cubic spline interpolation | |
| CN103399204A (en) | Rife-Vincent (II) window interpolation FFT (Fast Fourier Transform)-based harmonic and inter-harmonic detection method | |
| CN104267258A (en) | Harmonic instantaneous power calculating method utilizing incomplete S transform | |
| CN105606900A (en) | Single-phase harmonic impedance measuring method based on square wave signals | |
| CN109061345B (en) | Effective value measuring method and system suitable for power system | |
| CN105486921A (en) | Kaiser third-order mutual convolution window triple-spectrum-line interpolation harmonic wave and inter-harmonic wave detection method | |
| CN103389412B (en) | The electrical network electric parameter measurement method and device that a kind of small generator is constituted | |
| CN102095929A (en) | Method for rapidly measuring frequency of alternating-current signals | |
| CN103675447A (en) | High-precision real-time harmonic wave analysis method of electrified railway | |
| CN104375014B (en) | Single-frequency overlapping method for power capacitor noise test | |
| CN107167658B (en) | A kind of jamproof electric system fundamental frequency of high-precision and Method for Phase Difference Measurement | |
| CN109581045A (en) | A kind of m-Acetyl chlorophosphonazo power measurement method meeting IEC standard frame | |
| CN104407197A (en) | Signal phasor measurement method based on trigonometric function iteration | |
| CN108051682B (en) | Verification method of single-phase rectifier system impedance model | |
| CN102749488A (en) | Power grid harmonic wave real-time on-line monitor and method for detecting harmonic wave using same | |
| CN104360198A (en) | High-performance power quality analyzer | |
| CN107271814B (en) | A method of obtaining load nonlinear function | |
| CN106771552A (en) | A kind of distortion power metering method |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |