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
 frequency
 signal
 transformation
 quick
 broad sense
 Prior art date
Links
 239000011159 matrix materials Substances 0.000 claims abstract description 81
 238000005070 sampling Methods 0.000 claims abstract description 62
 230000000875 corresponding Effects 0.000 claims abstract description 36
 238000006243 chemical reactions Methods 0.000 claims abstract description 25
 238000000354 decomposition reactions Methods 0.000 claims abstract description 12
 238000001914 filtration Methods 0.000 claims abstract description 12
 238000000034 methods Methods 0.000 claims description 26
 238000001228 spectrum Methods 0.000 claims description 22
 230000001131 transforming Effects 0.000 claims description 15
 230000001186 cumulative Effects 0.000 claims description 14
 230000001066 destructive Effects 0.000 claims description 14
 239000000047 products Substances 0.000 claims description 12
 238000004422 calculation algorithm Methods 0.000 claims description 11
 238000004458 analytical methods Methods 0.000 claims description 7
 235000013399 edible fruits Nutrition 0.000 claims description 5
 280000584129 Justice companies 0.000 claims description 4
 238000005265 energy consumption Methods 0.000 abstract 2
 238000006073 displacement reactions Methods 0.000 description 3
 230000037010 Beta Effects 0.000 description 2
 241000196324 Embryophyta Species 0.000 description 2
 238000004364 calculation methods Methods 0.000 description 2
 238000004088 simulation Methods 0.000 description 2
 230000001052 transient Effects 0.000 description 2
 240000004678 Panax pseudoginseng Species 0.000 description 1
 235000003140 Panax quinquefolius Nutrition 0.000 description 1
 280000398338 Seismic companies 0.000 description 1
 239000006227 byproducts Substances 0.000 description 1
 230000000694 effects Effects 0.000 description 1
 238000005516 engineering processes Methods 0.000 description 1
 235000005035 ginseng Nutrition 0.000 description 1
 235000008434 ginseng Nutrition 0.000 description 1
 238000002360 preparation methods Methods 0.000 description 1
 230000001105 regulatory Effects 0.000 description 1
 230000002441 reversible Effects 0.000 description 1
 238000000926 separation method Methods 0.000 description 1
Abstract
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)
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 CN104090159B (en)  20140716  20140716  Electric energy measuring method and device 
Country Status (1)
Country  Link 

CN (1)  CN104090159B (en) 
Families Citing this family (3)
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 
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 IP Right Grant
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 

CN104635044B (en)  Power signal frequency detection method and system based on amplitude modulation  
US9069995B1 (en)  Multiply accumulate operations in the analog domain  
Krishnanand et al.  A new realtime fast discrete Stransform for crossdifferential protection of shuntcompensated power systems  
Bergland  A guided tour of the fast Fourier transform  
Vetterli et al.  Sampling signals with finite rate of innovation  
Tomic et al.  A new power system digital harmonic analyzer  
CN102435844B (en)  Sinusoidal signal phasor calculating method being independent of frequency  
Jain et al.  Harmonics estimation in emerging power system: Key issues and challenges  
Mainardi et al.  Timefractional diffusion of distributed order  
Stenbakken  A Wideband Sampling Wattmeter1  
CN104007316B (en)  A kind of High Precision Frequency method under lack sampling speed and measuring instrument thereof  
Eldar  Compressed sensing of analog signals in shiftinvariant spaces  
CN104391178B (en)  A kind of time shift phase difference stable state harmonic signal bearing calibration based on Nuttall windows  
CN103973324B (en)  A kind of wideband digital receiver and real time spectrum processing method thereof  
CN101930426B (en)  Signal processing method, data processing method and device  
Hamid  Wavelet packet transform for RMS values and power measurements  
Hamid et al.  Method for RMS and power measurements based on the wavelet packet transform  
CN101566649A (en)  Harmonic detection method in a power system  
CN103412242B (en)  A kind of harmonic source location method based on Fast Independent Component Analysis and mutual information  
CN103308766A (en)  Harmonic analysis method based on Kaiser selfconvolution window dualspectrum line interpolation FFT (Fast Fourier Transform) and device thereof  
CN104459321B (en)  Power signal base wave phase measurement method and system  
CN105115594B (en)  Gearbox vibration signal fault signature extracting method based on Wavelet Entropy and information fusion  
Lovell et al.  The relationship between instantaneous frequency and timefrequency representations  
CN100554917C (en)  Obtain the method for system features function and signal characteristic value  
CN103792427B (en)  The method and system of real time spectral analysis are carried out to nonstationary signal 
Legal Events
Date  Code  Title  Description 

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