CN104090159B  Electric energy measuring method and device  Google Patents
Electric energy measuring method and device Download PDFInfo
 Publication number
 CN104090159B CN104090159B CN201410339464.5A CN201410339464A CN104090159B CN 104090159 B CN104090159 B CN 104090159B CN 201410339464 A CN201410339464 A CN 201410339464A CN 104090159 B CN104090159 B CN 104090159B
 Authority
 CN
 China
 Prior art keywords
 signal
 frequency
 transformation
 quick
 broad sense
 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.)
 Active
Links
Abstract
The invention provides an electric energy measuring method and device. The method includes the step of conducting lowpass filtering on electric signals so as to obtain analog signals, the step of conducting A/D sampling on the analog signals so as to obtain multiple discrete signals, the step of conducting improved generalized rapid S conversion on each discrete signal so as to obtain a first result matrix corresponding to the discrete signal, the step of conducting linear decomposition on each first result matrix through the linear characteristic of the improved generalized rapid S conversion so as to obtain a second result matrix corresponding to each frequency component of each discrete signal, the step of conducting inverse improved generalized rapid S conversion on each second result matrix according to the inverse nondestructive characteristic so as to obtain a reconstructed time domain signal corresponding to each frequency component of each discrete signal, and the step of obtaining the electric energy consumption value within the A/D sampling time interval so as to obtain the electric energy consumption value within the preset time. By means of the electric energy measuring method and device, the Gauss window shape can be adjusted according to the actual condition, and operation efficiency is improved.
Description
Technical field
The application is related to technical field of electric power, particularly to a kind of electric energy gauging method and device.
Background technology
Electric energy be modern society produce and life in indispensable important energy source commodity, it produce, sell and use according to
The system that Lai Yu is made up of power plant, power pack and user tripartite.Therefore, the accuracy of electric energy metrical and the direct shadow of reasonability
Ring to sending out, for, economic interests of tripartite and transaction fairness.
At present, existing electric energy meter can only carry out electric energy metrical to steadystate distortion signal both at home and abroad, and letter that transient state is distorted
Number electric energy metrical then seem helpless.But, in industry spot, the frequent starting of motor, input of load etc. all can produce
Huge transient power, and electric energy meter cannot be realized measuring, so that power department suffers from huge losing.Therefore, deeply
The impact to electric energy metrical and how to realize accurate, rational electric energy meter in this context when being distorted of research power network signal
Measurer has important theory and realistic meaning.
The electric energy metrical algorithm commonly used at present has Stransformation (Stockwell Transform).Stransformation is
R.G.Stockwell proposed in 1996 and is used for the algorithm of seismic signal detection, and it is a kind of reversible time frequency analysis side
Method.But, Stransformation haves the shortcomings that Gauss window shape can not be adjusted according to practical situation, and due to Stransformation algorithm
Timefrequency matrix information amount huge so that computing expend time long, operation efficiency is low.
Content of the invention
For solving abovementioned technical problem, the embodiment of the present application provides a kind of electric energy gauging method and device, to reach realization
Gauss window shape can be adjusted according to practical situation, and improves the purpose of operation efficiency, and technical scheme is as follows：
A kind of electric energy gauging method, including：
To needing to carry out the signal of telecommunication in the electrical network of electric energy metrical to carry out lowpass filtering, obtain analogue signal；
A/D sampling is carried out to described analogue signal, obtains multiple discrete signals；
Respectively the quick Stransformation of broad sense is improved to various discrete signal, obtain various discrete signal each selfcorresponding
One matrix of consequence, the quick Stransformation of described improvement broad sense is to introduce the Gaussian window width adjusting factor and direct current on the basis of Stransformation
At component the Gaussian window width adjusting factor and using fast algorithm think of conceivable；
Using the linear characteristic improving the quick Stransformation of broad sense, linear decomposition is carried out to each first matrix of consequence, obtain each
Corresponding second matrix of consequence of each frequency component of individual discrete signal；
According to inversion nondestructive, the inverse improvement quick Stransformation of broad sense is carried out respectively to each second matrix of consequence, obtains each
The reconstruct timedomain signal of each frequency component of discrete signal；
Require according to electric energy metrical, each reconstruct timedomain signal is carried out with corresponding point multiplication and result of product adds up, will
The result of product accumulated value of each reconstruct timedomain signal be added obtain cumulative and, will be described cumulative and when being multiplied by A/D sampling interval
Between, obtain the power consumption value in described A/D sampling interval duration, to obtain the power consumption value in Preset Time.
Preferably, respectively the quick Stransformation of broad sense is improved to various discrete signal, obtain various discrete signal each right
The process of the first matrix of consequence answered, including：
Using discrete improvement broad sense quick Stransformation expression formula
Various discrete signal is improved with the quick Stransformation of broad sense, obtains each selfcorresponding first matrix of consequence of various discrete signal, institute
StateThe quick S of discrete improvement broad sense for AC portion becomes
Change expression formula, describedFor the quick Stransformation expression formula of discrete improvement broad sense of direct current component, a=0,
1 ... N1, n=1 ... n_{max} 1, p=width (n) ..., 0 ... width (n) 1, described N is electricity to be analyzed
The sampling number of signal, described Δ t is the sampling time interval of the signal of telecommunication to be analyzed, and (k Δ t) is each sampling instant to described h
Sampled value.
Preferably, described discrete improvement broad sense quick Stransformation expression formula is by carrying out to improvement broad sense quick Stransformation formula
Discrete transform obtains, and described improvement broad sense quick Stransformation formula is,
DescribedH (τ) is the signal of telecommunication to be analyzed, and t is time shift variable, and f is frequency
Rate, τ is the time, and α and β is the Gaussian window width adjusting factor, and γ is the Gaussian window width adjusting factor, w at DC component_{m}(t τ,
α, beta, gamma) it is to improve Generalized Gaussian window；
Wherein, to improving broad sense quick Stransformation formula and carrying out the process of discrete transform it is：
A：Using formulaCalculate the described frequency domain radius improving Generalized Gaussian window, described
For described frequency domain radius；
B：Calculate the frequency values that frequency sampling is counted out with center frequency points, described frequency sampling is counted out and is designated as n_{max}, institute
The frequency values stating center frequency points are designated as f_{n}, described center frequency points are steppedfrequency signal；
C：The discrete radius of improvement Generalized Gaussian window calculating each center frequency points isDescribed for width_{n}N () improves the discrete radius of Generalized Gaussian window,It is to round symbol；
D：The shift length calculating the frequency spectrum of each center frequency points isDescribed centre_{n}
For shift length；
E：Calculate the coverage of described each frequency domain improving Generalized Gaussian window, coverage is by f_{start}(n) and f_{end}
N () characterizes, wherein, f_{start}(0)=0, f_{start}(n+1)=centre_{n}(n)+width_{n}(n)+width_{n}(n+1), n=
0 ... ..., n_{max} 1, f_{end}(n)=f_{start}(n)+2width_{n}(n)；
F：Calculate the fast Fourier transform FFT result of the described signal of telecommunication to be analyzed, and described FFT result shifted,
Obtain the first FFT result, be designated as H [p], the direct current result in described first FFT result intermediate frequency spectrum information is located at center, negative frequency division
Amount and positive frequency component are located at both sides respectively；
G：The time domain radius calculating the improvement Generalized Gaussian window of each center frequency points isP=width_{n}(n) ,width_{n}(n)+1 ..., 0 ..., width_{n}
(n)1；
H：Described first FFT result is added with the shift length of the frequency spectrum of each center frequency points respectively, obtains each
Each selfcorresponding second FFT result of center frequency points, is designated as H [p+centre_{n}(n)]；
I：By respective for each center frequency points the second FFT result and its respective time domain radius improving Generalized Gaussian window
It is multiplied, obtain the operation result of each center frequency points, each operation result is arranged according to distributing position in FFT spectrum
Row, obtain rank results, carry out IFFT operation to described rank results, and the quick S of discrete improvement broad sense obtaining AC portion becomes
Change expression formula.
Preferably, described step B includes：
B11：The frequency making nth steppedfrequency signal is f_{n}(n is from the beginning of 0), then the relation of two neighboring mid frequency：
B12：First frequency center point is made to be f for dc point_{0}=0, and foundationObtain the frequency values f of each center frequency points_{n}, n=0 ... ..., n_{max}1；
B13：Solve inequalityObtain frequency sampling to count out n_{max}, described f_{s}Treat for described
The sample frequency of the analysis signal of telecommunication.
Preferably, using the linear characteristic improving the quick Stransformation of broad sense, each first matrix of consequence is linearly divided
Solution, obtains the process of corresponding second matrix of consequence of each frequency component of various discrete signal, including：
Using linear characteristic S [h (t)]=S [x (t)+y (t)]=S [x (t)]+S [y improving the quick Stransformation of broad sense
(t)], calculate each the first matrix of consequence certain frequency component when, only retain in described first matrix of consequence to should frequency
The part of rate component, makes the corresponding part of frequency component in addition to this frequency component in described first matrix of consequence be entered as 0；
Described h (t), x (t), y (t) are signal to be analyzed and h (t)=x (t)+y (t).
Preferably, described inversion nondestructive is：
Described inverse improve the quick Stransformation of broad sense principle be
H (v+f)=α (v, f)/W_{m}(v),
Described a (v, f) is the Fourier transformation result to described time shift variable t and the Fourier transformation to described time τ
As a result, described v is the Fourier transformation to described time shift variable t, and described f is the Fourier transformation to described τ, and H (v+f) is right
The Fourier transformation result of the signal of telecommunication to be analyzed.
Preferably, require according to electric energy metrical, corresponding point multiplication is carried out and by product to any one reconstruct timedomain signal
The process that result adds up, including：
This reconstruct timedomain signal is carried out separate, obtain each harmonic voltage, electric current discrete signal, be designated as u respectively_{n}[k]、
i_{n}[k]；
According to formulaCalculate the electric energy consuming in this reconstruct timedomain signal unit interval.
A kind of electric power meter, including：
Filtration module, for needing to carry out the signal of telecommunication in the electrical network of electric energy metrical to carry out lowpass filtering, obtaining simulation letter
Number；
Sampling module, for carrying out A/D sampling to described analogue signal, obtains multiple discrete signals；
First conversion module, for improving the quick Stransformation of broad sense to various discrete signal respectively, obtains various discrete
Each selfcorresponding first matrix of consequence of signal, the quick Stransformation of described improvement broad sense is introducing Gauss window width on the basis of Stransformation
Degree regulatory factor and DC component at the Gaussian window width adjusting factor and using fast algorithm think of conceivable；
Decomposing module, for using the linear characteristic improving the quick Stransformation of broad sense, entering line to each first matrix of consequence
Property decompose, obtain corresponding second matrix of consequence of each frequency component of various discrete signal；
Second conversion module, fast for inverse improvement broad sense is carried out according to inversion nondestructive respectively to each second matrix of consequence
Fast Stransformation, obtains the reconstruct timedomain signal of each frequency component of various discrete signal；
Computing module, for requiring according to electric energy metrical, carries out corresponding point to each reconstruct timedomain signal and is multiplied and will take advantage of
Longpending result adds up, result of product accumulated value of each reconstruct timedomain signal is added obtain cumulative and, will be described cumulative and be multiplied by
A/D sampling interval duration, obtains the power consumption value in described A/D sampling interval duration, to obtain the electricity in Preset Time
Can consumption figures.
Preferably, described first conversion module includes：
First converter unit, for using discrete improvement broad sense quick Stransformation expression formula
Various discrete signal is improved with the quick Stransformation of broad sense, obtains each selfcorresponding first knot of various discrete signal
Fruit matrix, describedWide for the discrete improvement of AC portion
The quick Stransformation expression formula of justice, describedDiscrete improvement broad sense quick Stransformation expression for direct current component
Formula, a=0,1 ... ..N1, n=1 ... ..n_{max} 1, p=width (n) ..., 0 ... width (n) 1, described N is to treat
The sampling number of the analysis signal of telecommunication, described Δ t is the sampling time interval of the signal of telecommunication to be analyzed, and (k Δ t) adopts described h for each
The sampled value in sample moment；
Described discrete improvement broad sense quick Stransformation expression formula is by carrying out discrete change to improvement broad sense quick Stransformation formula
Get in return, described improvement broad sense quick Stransformation formula is, describedH (τ) is the signal of telecommunication to be analyzed, and t is time shift variable, and f is frequency, τ
It it is the time, α and β is the Gaussian window width adjusting factor, and γ is the Gaussian window width adjusting factor, w at DC component_{m}(t τ, α, β,
It is γ) to improve Generalized Gaussian window.
Preferably, also include：
3rd conversion module, for carrying out discrete transform to improvement broad sense quick Stransformation formula；
Described 3rd conversion module includes：
First computing unit, for utilizing formulaCalculate the described frequency improving Generalized Gaussian window
Domain radius, describedFor described frequency domain radius；
Second computing unit, for calculating the frequency values that frequency sampling is counted out with center frequency points, described frequency sampling
Count out and be designated as n_{max}, the frequency values of described center frequency points are designated as f_{n}, described center frequency points are steppedfrequency signal；
3rd computing unit, the discrete radius of improvement Generalized Gaussian window for calculating each center frequency points isDescribed for width_{n}N () improves the discrete radius of Generalized Gaussian window,It is to round symbol；
4th computing unit, the shift length of the frequency spectrum for calculating each center frequency points isDescribed centre_{n}For shift length；
5th computing unit, for calculating the coverage of described each frequency domain improving Generalized Gaussian window, coverage
By f_{start}(n) and f_{end}N () characterizes, wherein, f_{start}(0)=0, f_{start}(n+1)=centre_{n}(n)+width_{n}(n)+width_{n}
(n+1), n=0 ... ..., n_{max} 1, f_{end}(n)=f_{start}(n)+2width_{n}(n)；
6th computing unit, for calculating the FFT result of the described signal of telecommunication to be analyzed, and moves to described FFT result
Position, obtains the first FFT result, is designated as H [p], and the direct current result in described first FFT result intermediate frequency spectrum information is located at center, bears
Frequency component and positive frequency component are located at both sides respectively；
7th computing unit, the time domain radius of the improvement Generalized Gaussian window for calculating each center frequency points isP=width_{n}(n) ,width_{n}(n)+1 ..., 0 ...,
width_{n}(n)1；
8th computing unit, for by described first FFT result respectively with the displacement of the frequency spectrum of each center frequency points away from
From addition, obtain each selfcorresponding second FFT result of each center frequency points, be designated as H [p+centre_{n}(n)]；
9th computing unit, for by respective for each center frequency points the second FFT result and its respective improvement broad sense
The time domain radius of Gaussian window is multiplied, and obtains the operation result of each center frequency points, by each operation result according in FFT spectrum
Distributing position is arranged, and obtains rank results, carries out IFFT operation to described rank results, obtains the discrete of AC portion and changes
Enter broad sense quick Stransformation expression formula.
Preferably, described second computing unit includes：First computation subunit, with the frequency of nth steppedfrequency signal up to the present
Rate is f_{n}(n is from the beginning of 0), then the relation of two neighboring mid frequency：
Second computation subunit, with up to the present first frequency center point be dc point be f_{0}=0, and foundationObtain the frequency values f of each center frequency points_{n}, n=0 ... ..., n_{max}1；
3rd computation subunit, for solving inequalityObtain frequency sampling to count out n_{max},
Described f_{s}Sample frequency for the described signal of telecommunication to be analyzed.
Preferably, described decomposing module is specifically for using linear characteristic S [h (t)]=S improving the quick Stransformation of broad sense
[x (t)+y (t)]=S [x (t)]+S [y (t)], when calculating certain frequency component, it is right in described first matrix of consequence only to retain
Should frequency component part, make in described first matrix of consequence the corresponding part of frequency component in addition to this frequency component assign
It is worth for 0；
Described h (t), x (t), y (t) are signal to be analyzed and h (t)=x (t)+y (t).
Preferably, described computing module includes：
Separative element, for carrying out to this reconstruct timedomain signal separating, obtains each harmonic voltage, electric current discrete signal,
It is designated as u respectively_{n}[k]、i_{n}[k]；
Tenth computing unit, for according to formulaWhen calculating this reconstruct timedomain signal unit
The electric energy of interior consumption.
Compared with prior art, the having the beneficial effect that of the application：
In this application, due to improve the quick Stransformation of broad sense introduce on the basis of Stransformation Gaussian window width adjusting because
The Gaussian window width adjusting factor at son and DC component, therefore Gauss window shape can be adjusted according to practical situation.
It is based on fast algorithm due to improving the quick Stransformation of broad sense, therefore improve computing using when improving generalized Stransform
Speed, and be combined with improve the quick Stransformation of broad sense linear characteristic carry out linear decomposition, reduce amount of calculation, improve fortune
Calculate efficiency.
Brief description
For the technical scheme being illustrated more clearly that in the embodiment of the present application, will make to required in embodiment description below
Accompanying drawing be briefly described it should be apparent that, drawings in the following description are only some embodiments of the present application, for
For those of ordinary skill in the art, without having to pay creative labor, it can also be obtained according to these accompanying drawings
His accompanying drawing.
Fig. 1 is a kind of flow chart based on the electric energy gauging method improving the quick Stransformation of broad sense that the application provides；
Fig. 2 is a kind of structural representation of the electric power meter of the quick Stransformation of improvement broad sense that the application provides.
Specific embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present application, the technical scheme in the embodiment of the present application is carried out clear, complete
Site preparation describes it is clear that described embodiment is only some embodiments of the present application, rather than whole embodiments.It is based on
Embodiment in the application, it is every other that those of ordinary skill in the art are obtained under the premise of not making creative work
Embodiment, broadly falls into the scope of the application protection.
Embodiment one
Refer to Fig. 1, based on the electric energy gauging method improving the quick Stransformation of broad sense the one of the application offer is provided
Plant flow chart, may comprise steps of：
Step S11：To needing to carry out the signal of telecommunication in the electrical network of electric energy metrical to carry out lowpass filtering, obtain analogue signal.
In the present embodiment, the signal of telecommunication includes current signal and voltage signal.
In the present embodiment, according to the setting of actual metered requirement and A/D sample rate, it is determined to the signal of telecommunication processing
Highest frequency value, then using analogue filter circuit the signal of telecommunication is filtered process.
Step S12：A/D sampling is carried out to described analogue signal, obtains multiple discrete signals.
Step S13：Respectively the quick Stransformation of broad sense is improved to various discrete signal, obtain various discrete signal each
Corresponding first matrix of consequence, the quick Stransformation of described improvement broad sense be introduce on the basis of Stransformation Gaussian window width adjusting because
Son and DC component at the Gaussian window width adjusting factor and using fast algorithm think of conceivable.
Step S14：Using the linear characteristic improving the quick Stransformation of broad sense, each first matrix of consequence is linearly divided
Solution, obtains corresponding second matrix of consequence of each frequency component of various discrete signal.
Step S15：According to inversion nondestructive, the inverse improvement quick Stransformation of broad sense is carried out respectively to each second matrix of consequence,
Obtain the reconstruct timedomain signal of each frequency component of various discrete signal.
Step S16：Require according to electric energy metrical, corresponding point multiplication is carried out and by result of product to each reconstruct timedomain signal
Cumulative, result of product accumulated value of each reconstruct timedomain signal is added obtain cumulative and, will be described cumulative and be multiplied by A/D sampling
Interval time, obtain the power consumption value in described A/D sampling interval duration, to obtain the power consumption in Preset Time
Value.
In the present embodiment, Preset Time is divided by A/D sampling interval duration, the value obtaining with the described A/D sampling interval
Power consumption value in time is multiplied, and obtains the power consumption value in Preset Time.
In this application, due to improve the quick Stransformation of broad sense introduce on the basis of Stransformation Gaussian window width adjusting because
The Gaussian window width adjusting factor at son and DC component, therefore Gauss window shape can be adjusted according to practical situation.
It is based on fast algorithm due to improving the quick Stransformation of broad sense, therefore improve computing using when improving generalized Stransform
Speed, and be combined with improve the quick Stransformation of broad sense linear characteristic carry out linear decomposition, reduce amount of calculation, improve fortune
Calculate efficiency.
Embodiment two
Another embodiment
In the present embodiment thus it is shown that respectively the quick Stransformation of broad sense is improved to various discrete signal, obtain each
The process of each selfcorresponding first matrix of consequence of discrete signal.
Due to the quick Stransformation of broad sense is improved to various discrete signal, obtain various discrete signal each selfcorresponding
The process of one matrix of consequence is identical, only improves the quick Stransformation of broad sense to any one discrete signal therefore in the present embodiment,
The process obtaining corresponding first matrix of consequence of this discrete signal is described.
In the present embodiment, the quick Stransformation of broad sense is improved to described discrete signal, obtain this discrete signal corresponding
The detailed process of the first matrix of consequence is：
Using discrete improvement broad sense quick Stransformation expression formula
Described discrete signal is improved with the quick Stransformation of broad sense, obtains the first matrix of consequence, describedThe quick Stransformation of discrete improvement broad sense for AC portion
Expression formula, describedFor the quick Stransformation expression formula of discrete improvement broad sense of direct current component, a=0,
1 ... N1, n=1 ... n_{max} 1, p=width (n) ..., 0 ... width (n) 1, described N is electricity to be analyzed
The sampling number of signal, described Δ t is the sampling time interval of the signal of telecommunication to be analyzed, and (k Δ t) is each sampling instant to described h
Sampled value.
In the present embodiment, byTo the signal of telecommunication
AC compounent carries out discrete improvement broad sense Fast transforms, obtains the first matrix of consequence of AC compounent.
ByDiscrete improvement broad sense Fast transforms are carried out to the DC component of the signal of telecommunication, obtains straight
First matrix of consequence of flow component.
Wherein, described discrete improvement broad sense quick Stransformation expression formula pass through to improvement broad sense quick Stransformation formula carry out from
Scattered conversion obtains, and described improvement broad sense quick Stransformation formula is, describedH (τ) is the signal of telecommunication to be analyzed, and t is time shift variable, and f is frequency, and τ is
Time, α and β is the Gaussian window width adjusting factor, and γ is the Gaussian window width adjusting factor, w at DC component_{m}(t τ, α, beta, gamma)
It is to improve Generalized Gaussian window.
Improve broad sense quick Stransformation formula compared to Stransformation formula of the prior artIncreased α, β and γ.
In the present embodiment, to improving broad sense quick Stransformation formula and carrying out the process of discrete transform it is：
Step A：Using formulaCalculate the described frequency domain radius improving Generalized Gaussian window, describedFor described frequency domain radius.
Step B：Calculate the frequency values that frequency sampling is counted out with center frequency points, described frequency sampling is counted out and is designated as
n_{max}, the frequency values of described center frequency points are designated as f_{n}, described center frequency points are steppedfrequency signal.
In the present embodiment, step B includes：
B11：The frequency making nth steppedfrequency signal is f_{n}(n is from the beginning of 0), then the relation of two neighboring mid frequency：
B12：First frequency center point is made to be f for dc point_{0}=0, and foundationObtain the frequency values f of each center frequency points_{n}, n=0 ... ..., n_{max}1.
B13：Solve inequalityObtain frequency sampling to count out n_{max}, described f_{s}Treat for described
The sample frequency of the analysis signal of telecommunication.
Step C：The discrete radius of improvement Generalized Gaussian window calculating each center frequency points isDescribed for width_{n}N () improves the discrete radius of Generalized Gaussian window,It is to round symbol.
Step D：The shift length calculating the frequency spectrum of each center frequency points isDescribed
centre_{n}For shift length.
Step E：Calculate the coverage of described each frequency domain improving Generalized Gaussian window, coverage is by f_{start}(n) and
f_{end}N () characterizes, wherein, f_{start}(0)=0, f_{start}(n+1)=centre_{n}(n)+width_{n}(n)+width_{n}(n+1), n=
0 ... ..., n_{max} 1, f_{end}(n)=f_{start}(n)+2width_{n}(n).
Step F：Calculate FFT (fast Fourier transform, the Fast Fourier of the described signal of telecommunication to be analyzed
Transformation) result, and described FFT result is shifted, obtains the first FFT result, is designated as H [p], and described first
Direct current result in FFT result intermediate frequency spectrum information is located at center, and negative frequency component and positive frequency component are located at both sides respectively.
Step G：The time domain radius calculating the improvement Generalized Gaussian window of each center frequency points is
.
Step H：Described first FFT result is added with the shift length of the frequency spectrum of each center frequency points respectively, obtains
Each selfcorresponding second FFT result of each center frequency points, is designated as H [p+centre_{n}(n)].
Step I：By respective for each center frequency points the second FFT result and its respective time domain improving Generalized Gaussian window
Radius is multiplied, and obtains the operation result of each center frequency points, each operation result is carried out according to distributing position in FFT spectrum
Arrangement, obtains rank results, carries out IFFT operation to described rank results, obtain the quick S of discrete improvement broad sense of AC portion
Conversion expression formula.
Embodiment three
In the present embodiment thus it is shown that utilizing the linear characteristic improving the quick Stransformation of broad sense, to each the first result square
Battle array carries out linear decomposition, obtains the process of corresponding second matrix of consequence of each frequency component of various discrete signal.
Due to carrying out linear decomposition to each first matrix of consequence, each frequency component obtaining various discrete signal is corresponding
The process of the second matrix of consequence is identical, only carries out linear decomposition to any one first matrix of consequence, obtain therefore in the present embodiment
The process of corresponding second matrix of consequence of each frequency component of the discrete signal belonging to this first matrix of consequence is described.
In the present embodiment, linear decomposition is carried out to described first matrix of consequence, obtain belonging to this first matrix of consequence
The process of corresponding second matrix of consequence of each frequency component of discrete signal is specially：
Using linear characteristic S [h (t)]=S [x (t)+y (t)]=S [x (t)]+S [y improving the quick Stransformation of broad sense
(t)], when calculating certain frequency component, only retain in described first matrix of consequence to should frequency component part, order is described
In first matrix of consequence, the corresponding part of the frequency component in addition to this frequency component is entered as 0.
Described h (t), x (t), y (t) are signal to be analyzed and h (t)=x (t)+y (t).
Now illustrate the process of linear decomposition, such as the first matrix of consequence S [a, n] of signal h (t) to be analyzed, when a is
Between, n is frequency, when needing to calculate the n of signal to be analyzed_{0}It is only necessary to retain S [a, the n of the first matrix of consequence during frequency component_{0}]
Part, except n_{0}The corresponding part of remaining frequency component outside frequency component is entered as 0, the n obtaining_{0}Second knot of frequency component
Fruit matrix is：
Example IV
In the present embodiment thus it is shown that inverse improvement extensively is carried out according to inversion nondestructive respectively to each second matrix of consequence
The quick Stransformation of justice, obtains the process of the reconstruct timedomain signal of each frequency component of various discrete signal.
Due to carrying out the inverse improvement quick Stransformation of broad sense according to inversion nondestructive respectively to each second matrix of consequence, obtain
Identical to the process of the reconstruct timedomain signal of each frequency component of various discrete signal, therefore only to any one in the present embodiment
Second matrix of consequence according to inversion nondestructive carry out inverse improve the quick Stransformation of broad sense, obtain belonging to this second matrix of consequence from
The process of the reconstruct timedomain signal of each frequency component of scattered signal is described.
To any one second matrix of consequence according to inversion nondestructive carry out inverse improve the quick Stransformation of broad sense, obtain this
The inversion nondestructive of the reconstruct timedomain signal institute foundation of each frequency component of the discrete signal belonging to two matrixs of consequence is：
To any one second matrix of consequence according to inversion nondestructive carry out inverse improve the quick Stransformation of broad sense, obtain this
The quick Stransformation of inverse improvement broad sense of the reconstruct timedomain signal institute foundation of each frequency component of the discrete signal belonging to two matrixs of consequence
Principle be：
H (v+f)=α (v, f)/W_{m}(v)
Described a (v, f) is the Fourier transformation result to described time shift variable t and the Fourier transformation to described time τ
As a result, described v is the Fourier transformation to described time shift variable t, and described f is the Fourier transformation to described τ, and H (v+f) is right
The Fourier transformation result of the signal of telecommunication to be analyzed.
In the present embodiment, H (v+f) displacement is obtained H (v), then IFFT operation is carried out to H (v), obtain h (τ).
In the present embodiment, the detailed process reconstructing the timedomain signal of some discrete signal may include steps of：
1) setting up a length isArrayWithIts
Middle arrayDeposit the improvement broad sense quick Stransformation result (volume of array of DC component and positive frequency component
Number from the beginning of 0), then its corresponding frequency range isWherein number is that the element of i (i.e. ith element) corresponds to
FrequencyIn addition set up array H [N] that length is N to deposit the FFT result of signal；
2) For n=0,2,3 ..., n_{max}1：
To pos_spe [f_{start}(n)+1] arrive pos_spe [f_{end}(n)] common 2width_{n}N () individual element carries out FFT operation, and
It is multiplied by coefficientAnd result is still left in pos_spe [f_{start}(n)+1] arrive pos_spe [f_{end}(n)] between；
This timefrequency domain improves Generalized Gaussian window sequence
P=width_{n}(n) ..., 0 ..., width_{n}(n) 1,
Place it in gauss [f_{start}(n)+1] arrive gauss [f_{end}(n)] between；
End of For.
3) by arrayMiddle element divided byIn corresponding element, and result is placed onIn.
4) by arrayMiddle element is put into corresponding position in array H [N].
5)For：
Element H [i] is taken it to be conjugated and is placed in H [Ni]；End of For.
6) carry out IFFT to H [N] to operate and be multiplied by coefficient N, finally give reconstruct timedomain signal.
It should be noted that the method that the present invention provides is applied to the electric energy metrical field under common stable state and unstable signal
Close.In order to the specific implementation process of the method for present invention offer is described, below by with the electric energy metrical under harmonic wave, mAcetyl chlorophosphonazo signal
It is analyzed as emphasis.
Make signal sampling frequencies 6.4kHz, emulate duration 0.2s.Carry out the voltage of electric energy metrical, current signal comprises harmonic wave
With mAcetyl chlorophosphonazo.So carry out including during electric energy metrical following some：First lowpass filtering treatment is carried out to voltage, current signal,
Obtain analogue signal, be that A/D sampling is prepared；Then, analogue signal is carried out with A/D sampling and obtains discrete signal, be next step
Signal processing lay the foundation；Afterwards, using improve broad sense quick Stransformation voltage, electric current are respectively processed obtain voltage,
The matrix of consequence of electric current, afterwards with linear characteristic, matrix of consequence is carried out separate, thus represented respectively fundamental wave component and its
The matrix of consequence of his composition, then carries out the inverse quick Stransformation of broad sense of improving and respectively obtains voltage, electricity to the matrix of consequence after separation
The fundamental signal of stream and distorted signal；Last foundation formulaSolve all kinds of energy values.Metric results
As shown in table 1, table 2, table 1 is harmonic wave, the electric energy metrical analysis of simulation result of mAcetyl chlorophosphonazo signal, table 2 be traditional Stransformation with quick
The operation time of Stransformation compares.
Wherein voltage signal is：
Current signal is：
Table 1
Table 2
Table 1, table 2 indicate based on improve the quick Stransformation of broad sense electric energy gauging method be capable of accurate to electric energy
Measure and measuring accuracy is higher than wavelet package transforms, simultaneously compared to Stransformation, the inventive method highly shortened operation time.
Embodiment five
In the present embodiment, there is provided the electric power meter based on the improvement quick Stransformation of broad sense that the application provides, please
Referring to Fig. 2, it illustrates a kind of structural representation of the electric power meter of the quick Stransformation of improvement broad sense of the application offer,
Included based on the electric power meter improving the quick Stransformation of broad sense：Filtration module 21, sampling module 22, the first conversion module 23,
Decomposing module 24, the second conversion module 25 and computing module 26.
Filtration module 21, for needing to carry out the signal of telecommunication in the electrical network of electric energy metrical to carry out lowpass filtering, being simulated
Signal.
Sampling module 22, for carrying out A/D sampling to described analogue signal, obtains multiple discrete signals.
First conversion module 23, for respectively the quick Stransformation of broad sense being improved to various discrete signal, obtain each from
Each selfcorresponding first matrix of consequence of scattered signal, the quick Stransformation of described improvement broad sense is introducing Gaussian window on the basis of Stransformation
At the width adjusting factor and DC component the Gaussian window width adjusting factor and using fast algorithm think of conceivable.
In the present embodiment, the first conversion module 23 includes：First converter unit, for quick using discrete improvement broad sense
Stransformation expression formula
Various discrete signal is improved with the quick Stransformation of broad sense, obtains each selfcorresponding first knot of various discrete signal
Fruit matrix, describedWide for the discrete improvement of AC portion
The quick Stransformation expression formula of justice, a=0,1 ... N1, n=1 ... n_{max} 1, p=width (n) ..., 0,
... width (n) 1, describedDiscrete improvement broad sense quick Stransformation expression for direct current component
Formula, described N is the sampling number of the signal of telecommunication to be analyzed, and described Δ t is the sampling time interval of the signal of telecommunication to be analyzed, described h (k
Δ t) is the sampled value of each sampling instant；
Described discrete improvement broad sense quick Stransformation expression formula is by carrying out discrete change to improvement broad sense quick Stransformation formula
Get in return, described improvement broad sense quick Stransformation formula is, describedH (τ) is the signal of telecommunication to be analyzed, and t is time shift variable, and f is frequency, τ
It it is the time, α and β is the Gaussian window width adjusting factor, and γ is the Gaussian window width adjusting factor, w at DC component_{m}(t τ, α, β,
It is γ) to improve Generalized Gaussian window.
Decomposing module 24, for using the linear characteristic improving the quick Stransformation of broad sense, carrying out to each first matrix of consequence
Linear decomposition, obtains corresponding second matrix of consequence of each frequency component of various discrete signal.
Second conversion module 25, for carrying out inverse improvement broad sense according to inversion nondestructive respectively to each second matrix of consequence
Quick Stransformation, obtains the reconstruct timedomain signal of each frequency component of various discrete signal.
Computing module 26, for according to electric energy metrical require, to each reconstruct timedomain signal carry out corresponding point be multiplied and incite somebody to action
Result of product adds up, result of product accumulated value of each reconstruct timedomain signal is added obtain cumulative and, will be described cumulative and take advantage of
With A/D sampling interval duration, obtain the power consumption value in described A/D sampling interval duration, to obtain in Preset Time
Power consumption value.
In the present embodiment, the first conversion module 23, decomposing module 24, the second conversion module 25 and computing module 26 are permissible
It is integrated in same digital signal processor.
Can also be included based on the electric power meter improving the quick Stransformation of broad sense shown in Fig. 2：3rd conversion module,
For discrete transform is carried out to improvement broad sense quick Stransformation formula.
Wherein, the 3rd conversion module includes：First computing unit, the second computing unit, the 3rd computing unit, the 4th calculating
Unit, the 5th computing unit, the 6th computing unit, the 7th computing unit, the 8th computing unit and the 9th computing unit.
First computing unit, for utilizing formulaCalculate the described frequency improving Generalized Gaussian window
Domain radius, describedFor described frequency domain radius.
Second computing unit, for calculating the frequency values that frequency sampling is counted out with center frequency points, described frequency sampling
Count out and be designated as n_{max}, the frequency values of described center frequency points are designated as f_{n}, described center frequency points are steppedfrequency signal.
In the present embodiment, the second computing unit specifically includes：First computation subunit, the second computation subunit and the 3rd
Computation subunit.
First computation subunit, the frequency with nth steppedfrequency signal up to the present is f_{n}(n is from the beginning of 0), then two neighboring in
The relation of frequency of heart：
Second computation subunit, with up to the present first frequency center point be dc point be f_{0}=0, and foundationObtain the frequency values f of each center frequency points_{n}, n=0 ... ..., n_{max}1.
3rd computation subunit, for solving inequalityObtain frequency sampling to count out
n_{max}, described f_{s}Sample frequency for the described signal of telecommunication to be analyzed.
3rd computing unit, the discrete radius of improvement Generalized Gaussian window for calculating each center frequency points isDescribed for width_{n}N () improves the discrete radius of Generalized Gaussian window,It is to round symbol.
4th computing unit, the shift length of the frequency spectrum for calculating each center frequency points isDescribed centre_{n}For shift length.
5th computing unit, for calculating the coverage of described each frequency domain improving Generalized Gaussian window, coverage
By f_{start}(n) and f_{end}N () characterizes, wherein, f_{start}(0)=0, f_{start}(n+1)=centre_{n}(n)+width_{n}(n)+width_{n}
(n+1), n=0 ... ..., n_{max} 1, f_{end}(n)=f_{start}(n)+2width_{n}(n).
6th computing unit, for calculating the FFT result of the described signal of telecommunication to be analyzed, and moves to described FFT result
Position, obtains the first FFT result, is designated as H [p], and the direct current result in described first FFT result intermediate frequency spectrum information is located at center, bears
Frequency component and positive frequency component are located at both sides respectively.
7th computing unit, the time domain radius of the improvement Generalized Gaussian window for calculating each center frequency points is
.
8th computing unit, for by described first FFT result respectively with the displacement of the frequency spectrum of each center frequency points away from
From addition, obtain each selfcorresponding second FFT result of each center frequency points, be designated as H [p+centre_{n}(n)].
9th computing unit, for by respective for each center frequency points the second FFT result and its respective improvement broad sense
The time domain radius of Gaussian window is multiplied, and obtains the operation result of each center frequency points, by each operation result according in FFT spectrum
Distributing position is arranged, and obtains rank results, carries out IFFT operation to described rank results, obtains the discrete of AC portion and changes
Enter broad sense quick Stransformation expression formula.
In the present embodiment, described decomposing module 24 is specifically for using the linear characteristic S [h improving the quick Stransformation of broad sense
(t)]=S [x (t)+y (t)]=S [x (t)]+S [y (t)], when calculating certain frequency component, only retain described first result square
Battle array in should frequency component part, make the frequency component in addition to this frequency component in described first matrix of consequence corresponding
Part is entered as 0；
Described h (t), x (t), y (t) are signal to be analyzed and h (t)=x (t)+y (t)
In the present embodiment, computing module 26 specifically includes：Separative element and the tenth computing unit.
Separative element, for carrying out to this reconstruct timedomain signal separating, obtains each harmonic voltage, electric current discrete signal,
It is designated as u respectively_{n}[k]、i_{n}[k].
Tenth computing unit, for according to formulaCalculate this reconstruct timedomain signal unit
The electric energy consuming in time.
It should be noted that each embodiment in this specification is all described by the way of going forward one by one, each embodiment weight
Point explanation is all difference with other embodiment, between each embodiment identical similar partly mutually referring to.
For device class embodiment, due to itself and embodiment of the method basic simlarity, so description is fairly simple, related part ginseng
See that the part of embodiment of the method illustrates.
Last in addition it is also necessary to explanation, herein, such as first and second or the like relational terms be used merely to by
One entity or operation are made a distinction with another entity or operation, and not necessarily require or imply these entities or operation
Between there is any this actual relation or order.And, term " inclusion ", "comprising" or its any other variant meaning
Covering comprising of nonexcludability, so that including a series of process of key elements, method, article or equipment not only include that
A little key elements, but also include other key elements being not expressly set out, or also include for this process, method, article or
The intrinsic key element of equipment.In the absence of more restrictions, the key element being limited by sentence " including ... ", not
Also there is other identical element in including the process of described key element, method, article or equipment in exclusion.
Above electric energy gauging method provided herein and device are described in detail, used herein specifically
Individual example is set forth to the principle of the application and embodiment, and the explanation of above example is only intended to help and understands the application
Method and its core concept；Simultaneously for one of ordinary skill in the art, according to the thought of the application, it is being embodied as
All will change in mode and range of application, in sum, this specification content should not be construed as the restriction to the application.
Claims (13)
1. a kind of electric energy gauging method is it is characterised in that include：
To needing to carry out the signal of telecommunication in the electrical network of electric energy metrical to carry out lowpass filtering, obtain analogue signal；
A/D sampling is carried out to described analogue signal, obtains multiple discrete signals；
Respectively the quick Stransformation of broad sense is improved to various discrete signal, obtain each selfcorresponding first knot of various discrete signal
Fruit matrix, the quick Stransformation of described improvement broad sense is to introduce the Gaussian window width adjusting factor and DC component on the basis of Stransformation
Place the Gaussian window width adjusting factor and using fast algorithm think of conceivable；
Using improve the quick Stransformation of broad sense linear characteristic, linear decomposition is carried out to each first matrix of consequence, obtain each from
Corresponding second matrix of consequence of each frequency component of scattered signal；
According to inversion nondestructive, the inverse improvement quick Stransformation of broad sense is carried out respectively to each second matrix of consequence, obtains various discrete
The reconstruct timedomain signal of each frequency component of signal；
Require according to electric energy metrical, each reconstruct timedomain signal is carried out with corresponding point multiplication and result of product adds up, by each
The result of product accumulated value of reconstruct timedomain signal be added obtain cumulative and, will be described cumulative and be multiplied by A/D sampling interval duration, obtain
To the power consumption value in described A/D sampling interval duration, to obtain the power consumption value in Preset Time.
2. method according to claim 1 it is characterised in that improve the quick S of broad sense respectively to various discrete signal
Conversion, obtains the process of each selfcorresponding first matrix of consequence of various discrete signal, including：
Using discrete improvement broad sense quick Stransformation expression formula
Various discrete signal is improved with the quick Stransformation of broad sense, obtains each selfcorresponding first matrix of consequence of various discrete signal, institute
StateThe quick S of discrete improvement broad sense for AC portion becomes
Change expression formula, describedFor direct current component the quick S of discrete improvement broad sense convert expression formula, a=0,
1 ... N1, n=1 ... n_{max} 1, p=width (n) ..., 0 ... width (n) 1, described N is electricity to be analyzed
The sampling number of signal, described Δ t is the sampling time interval of the signal of telecommunication to be analyzed, and (k Δ t) is each sampling instant to described h
Sampled value, described width_{n}N () is the discrete radius of improvement Generalized Gaussian window of nth center frequency points, described centre_{n}(n)
Shift length for the frequency spectrum of nth center frequency points；
Wherein, described discrete improvement broad sense quick Stransformation expression formula is by carrying out discrete change to improvement broad sense quick Stransformation formula
Get in return, described improvement broad sense quick Stransformation formula is, describedH (τ) is the signal of telecommunication to be analyzed, and t is time shift variable, and f is frequency, τ
It it is the time, α and β is the Gaussian window width adjusting factor, and γ is the Gaussian window width adjusting factor, w at DC component_{m}(tτ,α,β,
It is γ) to improve Generalized Gaussian window.
3. method according to claim 2 is it is characterised in that described discrete improvement broad sense quick Stransformation expression formula is passed through
Carry out discrete transform to improvement broad sense quick Stransformation formula to obtain, described improvement broad sense quick Stransformation formula isDescribed
H (τ) is the signal of telecommunication to be analyzed, and t is time shift variable, and f is frequency, and τ is the time, and α and β is the Gaussian window width adjusting factor, and γ is
The Gaussian window width adjusting factor, w at DC component_{m}(t τ, α, β, γ) is to improve Generalized Gaussian window；
Wherein, to improving broad sense quick Stransformation formula and carrying out the process of discrete transform it is：
A：Using formulaCalculate the described frequency domain radius improving Generalized Gaussian window, describedFor institute
State frequency domain radius；
B：Calculate the frequency values that frequency sampling is counted out with center frequency points, described frequency sampling is counted out and is designated as n_{max}, described in
The frequency values of frequency of heart point are designated as f_{n}, described center frequency points are steppedfrequency signal；
C：The discrete radius of improvement Generalized Gaussian window calculating each center frequency points isInstitute
State as width_{n}N () improves the discrete radius of Generalized Gaussian window,It is to round symbol；
D：The shift length calculating the frequency spectrum of each center frequency points isDescribed centre_{n}For position
Move distance；
E：Calculate the coverage of described each frequency domain improving Generalized Gaussian window, coverage is by f_{start}(n) and f_{end}(n) table
Levy, wherein, f_{start}(0)=0, f_{start}(n+1)=centre_{n}(n)+width_{n}(n)+width_{n}(n+1), n=0 ..., n_{max}
1, f_{end}(n)=f_{start}(n)+2width_{n}(n)；
F：Calculate the fast Fourier transform FFT result of the described signal of telecommunication to be analyzed, and described FFT result is shifted, obtain
First FFT result, is designated as H [p], the direct current result in described first FFT result intermediate frequency spectrum information be located at center, negative frequency component and
Positive frequency component is located at both sides respectively；
G：The time domain radius calculating the improvement Generalized Gaussian window of each center frequency points is
；
H：Described first FFT result is added with the shift length of the frequency spectrum of each center frequency points respectively, obtains each center
Each selfcorresponding second FFT result of Frequency point, is designated as H [p+centre_{n}(n)]；
I：By respective for each center frequency points the second FFT result and its respective time domain radius phase improving Generalized Gaussian window
Take advantage of, obtain the operation result of each center frequency points, each operation result is arranged according to distributing position in FFT spectrum,
Obtain rank results, IFFT operation is carried out to described rank results, obtain the quick Stransformation table of discrete improvement broad sense of AC portion
Reach formula.
4. method according to claim 3 is it is characterised in that described step B includes：
B11：The frequency making nth steppedfrequency signal is f_{n}(n is from the beginning of 0), then the relation of two neighboring mid frequency：
B12：First frequency center point is made to be f for dc point_{0}=0, and foundation
Obtain the frequency values f of each center frequency points_{n}, n=0 ... ..., n_{max}1；
B13：Solve inequalityObtain frequency sampling to count out n_{max}, described f_{s}For described to be analyzed
The sample frequency of the signal of telecommunication.
5. method according to claim 3 is it is characterised in that utilize the linear characteristic improving the quick Stransformation of broad sense, to each
Individual first matrix of consequence carries out linear decomposition, obtains the mistake of corresponding second matrix of consequence of each frequency component of various discrete signal
Journey, including：
Using improve the quick Stransformation of broad sense linear characteristic S [h (t)]=S [x (t)+y (t)]=S [x (t)]+S [y (t)],
Calculate each the first matrix of consequence certain frequency component when, only retain in described first matrix of consequence to should frequency component
Part, makes the corresponding part of frequency component in addition to this frequency component in described first matrix of consequence be entered as 0；
Described h (t), x (t), y (t) are signal to be analyzed and h (t)=x (t)+y (t).
6. method according to claim 3 is it is characterised in that described inversion nondestructive is：
Described inverse improve the quick Stransformation of broad sense principle be
H (υ+f)=α (υ, f)/W_{m}(υ),
Described a (v, f) is the Fourier transformation result to described time shift variable t and the Fourier transformation result to described time τ,
Described v is the Fourier transformation to described time shift variable t, and described f is the Fourier transformation to described τ, and H (υ+f) is to treat point
The Fourier transformation result of the analysis signal of telecommunication.
7. method according to claim 6, it is characterised in that requiring according to electric energy metrical, reconstructs time domain to any one
Signal carries out corresponding point multiplication the process that result of product adds up, including：
This reconstruct timedomain signal is carried out separate, obtain each harmonic voltage, electric current discrete signal, be designated as u respectively_{n}[k]、i_{n}
[k]；
According to formulaCalculate the electric energy consuming in this reconstruct timedomain signal unit interval.
8. a kind of electric power meter is it is characterised in that include：
Filtration module, for needing to carry out the signal of telecommunication in the electrical network of electric energy metrical to carry out lowpass filtering, obtaining analogue signal；
Sampling module, for carrying out A/D sampling to described analogue signal, obtains multiple discrete signals；
First conversion module, for improving the quick Stransformation of broad sense to various discrete signal respectively, obtains various discrete signal
Each selfcorresponding first matrix of consequence, the quick Stransformation of described improvement broad sense is to introduce Gauss window width on the basis of Stransformation to adjust
Section the factor and DC component at the Gaussian window width adjusting factor and using fast algorithm think of conceivable；
Decomposing module, for using the linear characteristic improving the quick Stransformation of broad sense, linearly being divided to each first matrix of consequence
Solution, obtains corresponding second matrix of consequence of each frequency component of various discrete signal；
Second conversion module, for carrying out the inverse improvement quick S of broad sense according to inversion nondestructive respectively to each second matrix of consequence
Conversion, obtains the reconstruct timedomain signal of each frequency component of various discrete signal；
Computing module, for requiring according to electric energy metrical, carries out corresponding point multiplication and ties product to each reconstruct timedomain signal
Fruit is cumulative, result of product accumulated value of each reconstruct timedomain signal is added obtain cumulative and, will be described cumulative and be multiplied by A/D and adopt
Sample interval time, obtain the power consumption value in described A/D sampling interval duration, to obtain the power consumption in Preset Time
Value.
9. device according to claim 8 is it is characterised in that described first conversion module includes：
First converter unit, for using discrete improvement broad sense quick Stransformation expression formula
Various discrete signal is improved with the quick Stransformation of broad sense, obtains each selfcorresponding first result square of various discrete signal
Battle array, describedFast for the discrete improvement broad sense of AC portion
Fast Stransformation expression formula, describedFor the quick Stransformation expression formula of discrete improvement broad sense of direct current component, a
=0,1 ... N1, n=1 ... n_{max} 1, p=width (n) ..., 0 ... width (n) 1, described N is to treat point
The sampling number of the analysis signal of telecommunication, described Δ t is the sampling time interval of the signal of telecommunication to be analyzed, and (k Δ t's described h) samples for each
The sampled value in moment, described width_{n}N () is the discrete radius of improvement Generalized Gaussian window of nth center frequency points, described
centre_{n}N () is the shift length of the frequency spectrum of nth center frequency points；
Described discrete improvement broad sense quick Stransformation expression formula is obtained by carrying out discrete transform to improvement broad sense quick Stransformation formula
Arrive, described improvement broad sense quick Stransformation formula is, describedH (τ) is the signal of telecommunication to be analyzed, and t is time shift variable, and f is frequency, τ
It it is the time, α and β is the Gaussian window width adjusting factor, and γ is the Gaussian window width adjusting factor, w at DC component_{m}(tτ,α,β,
It is γ) to improve Generalized Gaussian window.
10. device according to claim 9 is it is characterised in that also include：
3rd conversion module, for carrying out discrete transform to improvement broad sense quick Stransformation formula；
Described 3rd conversion module includes：
First computing unit, for utilizing formulaCalculate the described frequency domain half improving Generalized Gaussian window
Footpath, describedFor described frequency domain radius；
Second computing unit, for calculating the frequency values that frequency sampling is counted out with center frequency points, described frequency sampling points
Mesh is designated as n_{max}, the frequency values of described center frequency points are designated as f_{n}, described center frequency points are steppedfrequency signal；
3rd computing unit, the discrete radius of improvement Generalized Gaussian window for calculating each center frequency points isDescribed for width_{n}N () improves the discrete radius of Generalized Gaussian window,It is to round symbol；
4th computing unit, the shift length of the frequency spectrum for calculating each center frequency points is
Described centre_{n}For shift length；
5th computing unit, for calculate described improve Generalized Gaussian window each frequency domain coverage, coverage by
f_{start}(n) and f_{end}N () characterizes, wherein, f_{start}(0)=0, f_{start}(n+1)=centre_{n}(n)+width_{n}(n)+width_{n}(n+
1), n=0 ..., n_{max} 1, f_{end}(n)=f_{start}(n)+2width_{n}(n)；
6th computing unit, for calculating the FFT result of the described signal of telecommunication to be analyzed, and shifts to described FFT result, obtains
To the first FFT result, it is designated as H [p], the direct current result in described first FFT result intermediate frequency spectrum information is located at center, negative frequency component
It is located at both sides with positive frequency component respectively；
7th computing unit, the time domain radius of the improvement Generalized Gaussian window for calculating each center frequency points is
；
8th computing unit, for by the described first FFT result shift length phase with the frequency spectrum of each center frequency points respectively
Plus, obtain each selfcorresponding second FFT result of each center frequency points, be designated as H [p+centre_{n}(n)]；
9th computing unit, for by respective for each center frequency points the second FFT result and its respective improvement Generalized Gaussian
The time domain radius of window is multiplied, and obtains the operation result of each center frequency points, and each operation result is distributed according in FFT spectrum
Position is arranged, and obtains rank results, carries out IFFT operation to described rank results, and the discrete improvement obtaining AC portion is wide
The quick Stransformation expression formula of justice.
11. devices according to claim 10 are it is characterised in that described second computing unit includes：First calculates son list
Unit, the frequency for making nth steppedfrequency signal is f_{n}(n is from the beginning of 0), then the relation of two neighboring mid frequency：
Second computation subunit, is used for making first frequency center point be f for dc point_{0}=0, and foundationObtain the frequency values f of each center frequency points_{n}, n=0 ... ..., n_{max}1；
3rd computation subunit, for solving inequalityObtain frequency sampling to count out n_{max}, described
f_{s}Sample frequency for the described signal of telecommunication to be analyzed.
12. devices according to claim 9 are it is characterised in that described decomposing module is specifically for fast using improving broad sense
Linear characteristic S [h (t)] of fast Stransformation=S [x (t)+y (t)]=S [x (t)]+S [y (t)], when calculating certain frequency component,
Only retain in described first matrix of consequence to should frequency component part, make in described first matrix of consequence except this frequency is divided
The corresponding part of frequency component outside amount is entered as 0；
Described h (t), x (t), y (t) are signal to be analyzed and h (t)=x (t)+y (t).
13. devices according to claim 12 are it is characterised in that described computing module includes：
Separative element, for carrying out to this reconstruct timedomain signal separating, obtains each harmonic voltage, electric current discrete signal, respectively
It is designated as u_{n}[k]、i_{n}[k]；
Tenth computing unit, for according to formulaCalculate in this reconstruct timedomain signal unit interval
The electric energy consuming.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201410339464.5A CN104090159B (en)  20140716  20140716  Electric energy measuring method and device 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201410339464.5A CN104090159B (en)  20140716  20140716  Electric energy measuring method and device 
Publications (2)
Publication Number  Publication Date 

CN104090159A CN104090159A (en)  20141008 
CN104090159B true CN104090159B (en)  20170215 
Family
ID=51637890
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201410339464.5A Active CN104090159B (en)  20140716  20140716  Electric energy measuring method and device 
Country Status (1)
Country  Link 

CN (1)  CN104090159B (en) 
Families Citing this family (5)
Publication number  Priority date  Publication date  Assignee  Title 

CN105092965A (en) *  20150901  20151125  银川蓝擎网络信息技术股份有限公司  Voltagecontrolled oscillatorbased power metering method and power metering circuit 
CN106645929B (en) *  20160930  20190319  国网江苏省电力公司电力科学研究院  A kind of improved electric car charging nonlinearload harmonic electric energy metering method 
CN106772204B (en) *  20161125  20190607  国网冀北电力有限公司电力科学研究院  The emulation mode and device of electric energy meter metering circuit 
CN109309513B (en) *  20180911  20210611  广东石油化工学院  Adaptive reconstruction method for power line communication signals 
CN111624399A (en) *  20200612  20200904  国网上海市电力公司  Electric energy metering method and system for nonlinear load 
Citations (3)
Publication number  Priority date  Publication date  Assignee  Title 

EP1700130B1 (en) *  20031230  20120919  Austriamicrosystems AG  Energy metering system 
CN103245832A (en) *  20130516  20130814  湖南大学  Harmonic time frequency characteristic parameter estimating method based on fast S conversion and analysis meter 
CN103308804A (en) *  20130617  20130918  湖南大学  Method for extracting timefrequency parameters of power quality disturbance signals on basis of fast KS (KaiserS) transformation 
Family Cites Families (1)
Publication number  Priority date  Publication date  Assignee  Title 

KR101759932B1 (en) *  20100427  20170720  엘지전자 주식회사  An estimating method for electric power consumption of an electric device 

2014
 20140716 CN CN201410339464.5A patent/CN104090159B/en active Active
Patent Citations (3)
Publication number  Priority date  Publication date  Assignee  Title 

EP1700130B1 (en) *  20031230  20120919  Austriamicrosystems AG  Energy metering system 
CN103245832A (en) *  20130516  20130814  湖南大学  Harmonic time frequency characteristic parameter estimating method based on fast S conversion and analysis meter 
CN103308804A (en) *  20130617  20130918  湖南大学  Method for extracting timefrequency parameters of power quality disturbance signals on basis of fast KS (KaiserS) transformation 
NonPatent Citations (2)
Title 

"基于S变换解构算法的电能计量方法研究";丁忠安等;《第三届全国电能质量学术会议暨电能质量行业发展论坛论文集》;20130815;第133141页 * 
"快速S变换在电能质量分析中的应用";张志禹等;《电网技术》;20130531;第37卷(第5期);第12851290页 * 
Also Published As
Publication number  Publication date 

CN104090159A (en)  20141008 
Similar Documents
Publication  Publication Date  Title 

CN104090159B (en)  Electric energy measuring method and device  
Vaidyanathan  Ramanujan sums in the context of signal processing—Part I: Fundamentals  
CN107064628B (en)  High Precise Frequency Measurement System and method  
CN103257271A (en)  Device and method for detecting micro grid harmonic wave and interharmonics based on STM32F107VCT6  
CN100516897C (en)  Harmonic energy metering method  
CN101900761B (en)  Highaccuracy nonintegerperiod sampled harmonic analysis and measurement method  
CN203287435U (en)  A micro electrical network harmonic wave and interharmonic wave test apparatus based on an STM32F107VCT6  
CN107478896A (en)  A kind of frequency adaptive harmonic current detection method based on cascade Generalized Integrator  
CN104535836A (en)  Electric signal fundamental frequency measuring method and system  
CN103353550A (en)  Method for measuring signal frequency and harmonic parameters of electric power system  
CN102662106A (en)  Method for measuring electric energy of harmonic power network  
US9128885B2 (en)  Computationally efficient finite impulse response comb filtering  
CN101701985B (en)  Constantfrequency variable dot power grid harmonic wave detection method and admeasuring apparatus thereof  
CN110348161B (en)  Multifrequencyband dynamic phasor electromagnetic transient simulation method and system for voltage source type converter  
CN104880592A (en)  Electric energy calculating method and system under harmonic wave condition  
CN102135552A (en)  Realtime digital detection method for active power and reactive power of electricity grid  
CN105334381A (en)  Method and device for measuring AC active power  
Hamid et al.  Wavelet packet transform for RMS and power measurements  
CN102520246A (en)  Constant frequency phasor extraction method  
CN104483545A (en)  Method and system for measuring harmonics of power system  
Wang et al.  Method for frequency measurement with ApFFT based on FPGA  
CN103105532A (en)  Harmonic electric energy measuring system of improved Fourier transform and control method thereof  
Mokeev et al.  Application of spectral representations in coordinates of complex frequency for the digital filter analysis and synthesis  
CN107707259A (en)  A kind of method of analog signal sampling and reconstruct  
Yang et al.  A smart method makes DFT more precise for power system frequency estimation 
Legal Events
Date  Code  Title  Description 

C06  Publication  
PB01  Publication  
C10  Entry into substantive examination  
SE01  Entry into force of request for substantive examination  
C14  Grant of patent or utility model  
GR01  Patent grant 