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
- 238000000034 method Methods 0.000 title claims abstract description 55
- 239000011159 matrix material Substances 0.000 claims abstract description 81
- 238000005070 sampling Methods 0.000 claims abstract description 62
- 238000006243 chemical reaction Methods 0.000 claims abstract description 25
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 12
- 238000001914 filtration Methods 0.000 claims abstract description 12
- 230000006872 improvement Effects 0.000 claims description 82
- 230000008569 process Effects 0.000 claims description 26
- 238000001228 spectrum Methods 0.000 claims description 22
- 230000009466 transformation Effects 0.000 claims description 15
- 230000001186 cumulative effect Effects 0.000 claims description 14
- 230000001066 destructive effect Effects 0.000 claims description 14
- 238000004422 calculation algorithm Methods 0.000 claims description 11
- 238000004458 analytical method Methods 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 7
- 235000013399 edible fruits Nutrition 0.000 claims description 5
- 230000005611 electricity Effects 0.000 claims description 5
- 230000008901 benefit Effects 0.000 claims description 3
- 238000005265 energy consumption Methods 0.000 abstract 2
- 239000000047 product Substances 0.000 description 7
- 238000006073 displacement reaction Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 241000196324 Embryophyta Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000006227 byproduct Substances 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007717 exclusion Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Landscapes
- Complex Calculations (AREA)
- Compression, Expansion, Code Conversion, And Decoders (AREA)
Abstract
The invention provides an electric energy measuring method and device. The method includes the step of conducting low-pass 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 steady-state 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 S-transformation (Stockwell Transform).S-transformation 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, S-transformation haves the shortcomings that Gauss window shape can not be adjusted according to practical situation, and due to S-transformation algorithm
Time-frequency matrix information amount huge so that computing expend time long, operation efficiency is low.
Content of the invention
For solving above-mentioned 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 low-pass filtering, obtain analogue signal；
A/D sampling is carried out to described analogue signal, obtains multiple discrete signals；
Respectively the quick S-transformation of broad sense is improved to various discrete signal, obtain various discrete signal each self-corresponding
One matrix of consequence, the quick S-transformation of described improvement broad sense is to introduce the Gaussian window width adjusting factor and direct current on the basis of S-transformation
At component the Gaussian window width adjusting factor and using fast algorithm think of conceivable；
Using the linear characteristic improving the quick S-transformation 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 non-destructive, the inverse improvement quick S-transformation of broad sense is carried out respectively to each second matrix of consequence, obtains each
The reconstruct time-domain signal of each frequency component of discrete signal；
Require according to electric energy metrical, each reconstruct time-domain signal is carried out with corresponding point multiplication and result of product adds up, will
The result of product accumulated value of each reconstruct time-domain 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 S-transformation 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 S-transformation expression formula
Various discrete signal is improved with the quick S-transformation of broad sense, obtains each self-corresponding 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 S-transformation expression formula of discrete improvement broad sense of direct current component, a=0,
1 ... N-1, n=1 ... nmax- 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 S-transformation expression formula is by carrying out to improvement broad sense quick S-transformation formula
Discrete transform obtains, and described improvement broad sense quick S-transformation 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 componentm(t- τ,
α, beta, gamma) it is to improve Generalized Gaussian window；
Wherein, to improving broad sense quick S-transformation 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 nmax, institute
The frequency values stating center frequency points are designated as fn, described center frequency points are stepped-frequency signal；
C：The discrete radius of improvement Generalized Gaussian window calculating each center frequency points isDescribed for widthnN () 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 centren
For shift length；
E：Calculate the coverage of described each frequency domain improving Generalized Gaussian window, coverage is by fstart(n) and fend
N () characterizes, wherein, fstart(0)=0, fstart(n+1)=centren(n)+widthn(n)+widthn(n+1), n=
0 ... ..., nmax- 1, fend(n)=fstart(n)+2widthn(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=-widthn(n) ,-widthn(n)+1 ..., 0 ..., widthn
(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 self-corresponding second FFT result of center frequency points, is designated as H [p+centren(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 n-th stepped-frequency signal is fn(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 point0=0, and foundationObtain the frequency values f of each center frequency pointsn, n=0 ... ..., nmax-1；
B13：Solve inequalityObtain frequency sampling to count out nmax, described fsTreat for described
The sample frequency of the analysis signal of telecommunication.
Preferably, using the linear characteristic improving the quick S-transformation 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 S-transformation 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 non-destructive is：
Described inverse improve the quick S-transformation of broad sense principle be
H (v+f)=α (v, f)/Wm(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 time-domain signal
The process that result adds up, including：
This reconstruct time-domain signal is carried out separate, obtain each harmonic voltage, electric current discrete signal, be designated as u respectivelyn[k]、
in[k]；
According to formulaCalculate the electric energy consuming in this reconstruct time-domain 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 low-pass 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 S-transformation of broad sense to various discrete signal respectively, obtains various discrete
Each self-corresponding first matrix of consequence of signal, the quick S-transformation of described improvement broad sense is introducing Gauss window width on the basis of S-transformation
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 S-transformation 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 non-destructive respectively to each second matrix of consequence
Fast S-transformation, obtains the reconstruct time-domain 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 time-domain signal and is multiplied and will take advantage of
Long-pending result adds up, result of product accumulated value of each reconstruct time-domain 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 S-transformation expression formula
Various discrete signal is improved with the quick S-transformation of broad sense, obtains each self-corresponding first knot of various discrete signal
Fruit matrix, describedWide for the discrete improvement of AC portion
The quick S-transformation expression formula of justice, describedDiscrete improvement broad sense quick S-transformation expression for direct current component
Formula, a=0,1 ... ..N-1, n=1 ... ..nmax- 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 S-transformation expression formula is by carrying out discrete change to improvement broad sense quick S-transformation formula
Get in return, described improvement broad sense quick S-transformation 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 componentm(t- τ, α, β,
It is γ) to improve Generalized Gaussian window.
Preferably, also include：
3rd conversion module, for carrying out discrete transform to improvement broad sense quick S-transformation 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 nmax, the frequency values of described center frequency points are designated as fn, described center frequency points are stepped-frequency signal；
3rd computing unit, the discrete radius of improvement Generalized Gaussian window for calculating each center frequency points isDescribed for widthnN () 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 centrenFor shift length；
5th computing unit, for calculating the coverage of described each frequency domain improving Generalized Gaussian window, coverage
By fstart(n) and fendN () characterizes, wherein, fstart(0)=0, fstart(n+1)=centren(n)+widthn(n)+widthn
(n+1), n=0 ... ..., nmax- 1, fend(n)=fstart(n)+2widthn(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=-widthn(n) ,-widthn(n)+1 ..., 0 ...,
widthn(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 self-corresponding second FFT result of each center frequency points, be designated as H [p+centren(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 S-transformation expression formula.
Preferably, described second computing unit includes：First computation subunit, with the frequency of n-th stepped-frequency signal up to the present
Rate is fn(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 f0=0, and foundationObtain the frequency values f of each center frequency pointsn, n=0 ... ..., nmax-1；
3rd computation subunit, for solving inequalityObtain frequency sampling to count out nmax,
Described fsSample 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 S-transformation 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 time-domain signal separating, obtains each harmonic voltage, electric current discrete signal,
It is designated as u respectivelyn[k]、in[k]；
Tenth computing unit, for according to formulaWhen calculating this reconstruct time-domain 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 S-transformation of broad sense introduce on the basis of S-transformation 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 S-transformation of broad sense, therefore improve computing using when improving generalized S-transform
Speed, and be combined with improve the quick S-transformation 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 S-transformation of broad sense that the application provides；
Fig. 2 is a kind of structural representation of the electric power meter of the quick S-transformation 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 S-transformation 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 low-pass 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 S-transformation of broad sense is improved to various discrete signal, obtain various discrete signal each
Corresponding first matrix of consequence, the quick S-transformation of described improvement broad sense be introduce on the basis of S-transformation 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 S-transformation 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 non-destructive, the inverse improvement quick S-transformation of broad sense is carried out respectively to each second matrix of consequence,
Obtain the reconstruct time-domain 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 time-domain signal
Cumulative, result of product accumulated value of each reconstruct time-domain 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 S-transformation of broad sense introduce on the basis of S-transformation 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 S-transformation of broad sense, therefore improve computing using when improving generalized S-transform
Speed, and be combined with improve the quick S-transformation 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 S-transformation of broad sense is improved to various discrete signal, obtain each
The process of each self-corresponding first matrix of consequence of discrete signal.
Due to the quick S-transformation of broad sense is improved to various discrete signal, obtain various discrete signal each self-corresponding
The process of one matrix of consequence is identical, only improves the quick S-transformation 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 S-transformation 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 S-transformation expression formula
Described discrete signal is improved with the quick S-transformation of broad sense, obtains the first matrix of consequence, describedThe quick S-transformation of discrete improvement broad sense for AC portion
Expression formula, describedFor the quick S-transformation expression formula of discrete improvement broad sense of direct current component, a=0,
1 ... N-1, n=1 ... nmax- 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 S-transformation expression formula pass through to improvement broad sense quick S-transformation formula carry out from
Scattered conversion obtains, and described improvement broad sense quick S-transformation 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 componentm(t- τ, α, beta, gamma)
It is to improve Generalized Gaussian window.
Improve broad sense quick S-transformation formula compared to S-transformation formula of the prior artIncreased α, β and γ.
In the present embodiment, to improving broad sense quick S-transformation 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
nmax, the frequency values of described center frequency points are designated as fn, described center frequency points are stepped-frequency signal.
In the present embodiment, step B includes：
B11：The frequency making n-th stepped-frequency signal is fn(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 point0=0, and foundationObtain the frequency values f of each center frequency pointsn, n=0 ... ..., nmax-1.
B13：Solve inequalityObtain frequency sampling to count out nmax, described fsTreat 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 widthnN () 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
centrenFor shift length.
Step E：Calculate the coverage of described each frequency domain improving Generalized Gaussian window, coverage is by fstart(n) and
fendN () characterizes, wherein, fstart(0)=0, fstart(n+1)=centren(n)+widthn(n)+widthn(n+1), n=
0 ... ..., nmax- 1, fend(n)=fstart(n)+2widthn(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 self-corresponding second FFT result of each center frequency points, is designated as H [p+centren(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 S-transformation 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 S-transformation 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 analyzed0It is only necessary to retain S [a, the n of the first matrix of consequence during frequency component0]
Part, except n0The corresponding part of remaining frequency component outside frequency component is entered as 0, the n obtaining0Second 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 non-destructive respectively to each second matrix of consequence
The quick S-transformation of justice, obtains the process of the reconstruct time-domain signal of each frequency component of various discrete signal.
Due to carrying out the inverse improvement quick S-transformation of broad sense according to inversion non-destructive respectively to each second matrix of consequence, obtain
Identical to the process of the reconstruct time-domain 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 non-destructive carry out inverse improve the quick S-transformation of broad sense, obtain belonging to this second matrix of consequence from
The process of the reconstruct time-domain signal of each frequency component of scattered signal is described.
To any one second matrix of consequence according to inversion non-destructive carry out inverse improve the quick S-transformation of broad sense, obtain this
The inversion non-destructive of the reconstruct time-domain 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 non-destructive carry out inverse improve the quick S-transformation of broad sense, obtain this
The quick S-transformation of inverse improvement broad sense of the reconstruct time-domain signal institute foundation of each frequency component of the discrete signal belonging to two matrixs of consequence
Principle be：
H (v+f)=α (v, f)/Wm(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 time-domain signal of some discrete signal may include steps of：
1) setting up a length isArrayWithIts
Middle arrayDeposit the improvement broad sense quick S-transformation 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. i-th 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 ..., nmax-1：
To pos_spe [fstart(n)+1] arrive pos_spe [fend(n)] common 2widthnN () individual element carries out FFT operation, and
It is multiplied by coefficientAnd result is still left in pos_spe [fstart(n)+1] arrive pos_spe [fend(n)] between；
This time-frequency domain improves Generalized Gaussian window sequence
P=-widthn(n) ..., 0 ..., widthn(n) -1,
Place it in gauss [fstart(n)+1] arrive gauss [fend(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 [N-i]；End of For.
6) carry out IFFT to H [N] to operate and be multiplied by coefficient N, finally give reconstruct time-domain 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, m-Acetyl 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 m-Acetyl chlorophosphonazo.So carry out including during electric energy metrical following some：First low-pass 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 S-transformation 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 S-transformation 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 m-Acetyl chlorophosphonazo signal, table 2 be traditional S-transformation with quick
The operation time of S-transformation compares.
Wherein voltage signal is：
Current signal is：
Table 1
Table 2
Table 1, table 2 indicate based on improve the quick S-transformation 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 S-transformation, 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 S-transformation 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 S-transformation of improvement broad sense of the application offer,
Included based on the electric power meter improving the quick S-transformation 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 low-pass 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 S-transformation of broad sense being improved to various discrete signal, obtain each from
Each self-corresponding first matrix of consequence of scattered signal, the quick S-transformation of described improvement broad sense is introducing Gaussian window on the basis of S-transformation
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
S-transformation expression formula
Various discrete signal is improved with the quick S-transformation of broad sense, obtains each self-corresponding first knot of various discrete signal
Fruit matrix, describedWide for the discrete improvement of AC portion
The quick S-transformation expression formula of justice, a=0,1 ... N-1, n=1 ... nmax- 1, p=-width (n) ..., 0,
... width (n) -1, describedDiscrete improvement broad sense quick S-transformation 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 S-transformation expression formula is by carrying out discrete change to improvement broad sense quick S-transformation formula
Get in return, described improvement broad sense quick S-transformation 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 componentm(t- τ, α, β,
It is γ) to improve Generalized Gaussian window.
Decomposing module 24, for using the linear characteristic improving the quick S-transformation 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 non-destructive respectively to each second matrix of consequence
Quick S-transformation, obtains the reconstruct time-domain signal of each frequency component of various discrete signal.
Computing module 26, for according to electric energy metrical require, to each reconstruct time-domain 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 time-domain 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 S-transformation of broad sense shown in Fig. 2：3rd conversion module,
For discrete transform is carried out to improvement broad sense quick S-transformation 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 nmax, the frequency values of described center frequency points are designated as fn, described center frequency points are stepped-frequency 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 n-th stepped-frequency signal up to the present is fn(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 f0=0, and foundationObtain the frequency values f of each center frequency pointsn, n=0 ... ..., nmax-1.
3rd computation subunit, for solving inequalityObtain frequency sampling to count out
nmax, described fsSample 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 widthnN () 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 centrenFor shift length.
5th computing unit, for calculating the coverage of described each frequency domain improving Generalized Gaussian window, coverage
By fstart(n) and fendN () characterizes, wherein, fstart(0)=0, fstart(n+1)=centren(n)+widthn(n)+widthn
(n+1), n=0 ... ..., nmax- 1, fend(n)=fstart(n)+2widthn(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 self-corresponding second FFT result of each center frequency points, be designated as H [p+centren(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 S-transformation expression formula.
In the present embodiment, described decomposing module 24 is specifically for using the linear characteristic S [h improving the quick S-transformation 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 time-domain signal separating, obtains each harmonic voltage, electric current discrete signal,
It is designated as u respectivelyn[k]、in[k].
Tenth computing unit, for according to formulaCalculate this reconstruct time-domain 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 low-pass filtering, obtain analogue signal；
A/D sampling is carried out to described analogue signal, obtains multiple discrete signals；
Respectively the quick S-transformation of broad sense is improved to various discrete signal, obtain each self-corresponding first knot of various discrete signal
Fruit matrix, the quick S-transformation of described improvement broad sense is to introduce the Gaussian window width adjusting factor and DC component on the basis of S-transformation
Place the Gaussian window width adjusting factor and using fast algorithm think of conceivable；
Using improve the quick S-transformation 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 non-destructive, the inverse improvement quick S-transformation of broad sense is carried out respectively to each second matrix of consequence, obtains various discrete
The reconstruct time-domain signal of each frequency component of signal；
Require according to electric energy metrical, each reconstruct time-domain signal is carried out with corresponding point multiplication and result of product adds up, by each
The result of product accumulated value of reconstruct time-domain 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 self-corresponding first matrix of consequence of various discrete signal, including：
Using discrete improvement broad sense quick S-transformation expression formula
Various discrete signal is improved with the quick S-transformation of broad sense, obtains each self-corresponding 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 ... N-1, n=1 ... nmax- 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 widthnN () is the discrete radius of improvement Generalized Gaussian window of n-th center frequency points, described centren(n)
Shift length for the frequency spectrum of n-th center frequency points；
Wherein, described discrete improvement broad sense quick S-transformation expression formula is by carrying out discrete change to improvement broad sense quick S-transformation formula
Get in return, described improvement broad sense quick S-transformation 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 componentm(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 S-transformation expression formula is passed through
Carry out discrete transform to improvement broad sense quick S-transformation formula to obtain, described improvement broad sense quick S-transformation 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 componentm(t- τ, α, β, γ) is to improve Generalized Gaussian window；
Wherein, to improving broad sense quick S-transformation 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 nmax, described in
The frequency values of frequency of heart point are designated as fn, described center frequency points are stepped-frequency signal；
C：The discrete radius of improvement Generalized Gaussian window calculating each center frequency points isInstitute
State as widthnN () 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 centrenFor position
Move distance；
E：Calculate the coverage of described each frequency domain improving Generalized Gaussian window, coverage is by fstart(n) and fend(n) table
Levy, wherein, fstart(0)=0, fstart(n+1)=centren(n)+widthn(n)+widthn(n+1), n=0 ..., nmax-
1, fend(n)=fstart(n)+2widthn(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 self-corresponding second FFT result of Frequency point, is designated as H [p+centren(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 S-transformation 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 n-th stepped-frequency signal is fn(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 point0=0, and foundation
Obtain the frequency values f of each center frequency pointsn, n=0 ... ..., nmax-1；
B13：Solve inequalityObtain frequency sampling to count out nmax, described fsFor 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 S-transformation 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 S-transformation 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 non-destructive is：
Described inverse improve the quick S-transformation of broad sense principle be
H (υ+f)=α (υ, f)/Wm(υ),
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 time-domain signal is carried out separate, obtain each harmonic voltage, electric current discrete signal, be designated as u respectivelyn[k]、in
[k]；
According to formulaCalculate the electric energy consuming in this reconstruct time-domain 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 low-pass 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 S-transformation of broad sense to various discrete signal respectively, obtains various discrete signal
Each self-corresponding first matrix of consequence, the quick S-transformation of described improvement broad sense is to introduce Gauss window width on the basis of S-transformation 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 S-transformation 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 non-destructive respectively to each second matrix of consequence
Conversion, obtains the reconstruct time-domain 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 time-domain signal
Fruit is cumulative, result of product accumulated value of each reconstruct time-domain 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 S-transformation expression formula
Various discrete signal is improved with the quick S-transformation of broad sense, obtains each self-corresponding first result square of various discrete signal
Battle array, describedFast for the discrete improvement broad sense of AC portion
Fast S-transformation expression formula, describedFor the quick S-transformation expression formula of discrete improvement broad sense of direct current component, a
=0,1 ... N-1, n=1 ... nmax- 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 widthnN () is the discrete radius of improvement Generalized Gaussian window of n-th center frequency points, described
centrenN () is the shift length of the frequency spectrum of n-th center frequency points；
Described discrete improvement broad sense quick S-transformation expression formula is obtained by carrying out discrete transform to improvement broad sense quick S-transformation formula
Arrive, described improvement broad sense quick S-transformation 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 componentm(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 S-transformation 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 nmax, the frequency values of described center frequency points are designated as fn, described center frequency points are stepped-frequency signal；
3rd computing unit, the discrete radius of improvement Generalized Gaussian window for calculating each center frequency points isDescribed for widthnN () 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 centrenFor shift length；
5th computing unit, for calculate described improve Generalized Gaussian window each frequency domain coverage, coverage by
fstart(n) and fendN () characterizes, wherein, fstart(0)=0, fstart(n+1)=centren(n)+widthn(n)+widthn(n+
1), n=0 ..., nmax- 1, fend(n)=fstart(n)+2widthn(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 self-corresponding second FFT result of each center frequency points, be designated as H [p+centren(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 S-transformation 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 n-th stepped-frequency signal is fn(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 point0=0, and foundationObtain the frequency values f of each center frequency pointsn, n=0 ... ..., nmax-1；
3rd computation subunit, for solving inequalityObtain frequency sampling to count out nmax, described
fsSample 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 S-transformation=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 time-domain signal separating, obtains each harmonic voltage, electric current discrete signal, respectively
It is designated as un[k]、in[k]；
Tenth computing unit, for according to formulaCalculate in this reconstruct time-domain signal unit interval
The electric energy consuming.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410339464.5A CN104090159B (en) | 2014-07-16 | 2014-07-16 | Electric energy measuring method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410339464.5A CN104090159B (en) | 2014-07-16 | 2014-07-16 | Electric energy measuring method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104090159A CN104090159A (en) | 2014-10-08 |
| CN104090159B true CN104090159B (en) | 2017-02-15 |
Family
ID=51637890
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410339464.5A Active CN104090159B (en) | 2014-07-16 | 2014-07-16 | Electric energy measuring method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104090159B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105092965A (en) * | 2015-09-01 | 2015-11-25 | 银川蓝擎网络信息技术股份有限公司 | Voltage-controlled oscillator-based power metering method and power metering circuit |
| CN106645929B (en) * | 2016-09-30 | 2019-03-19 | 国网江苏省电力公司电力科学研究院 | A kind of improved electric car charging nonlinear-load harmonic electric energy metering method |
| CN106772204B (en) * | 2016-11-25 | 2019-06-07 | 国网冀北电力有限公司电力科学研究院 | The emulation mode and device of electric energy meter metering circuit |
| CN109309513B (en) * | 2018-09-11 | 2021-06-11 | 广东石油化工学院 | Adaptive reconstruction method for power line communication signals |
| CN111624399A (en) * | 2020-06-12 | 2020-09-04 | 国网上海市电力公司 | Electric energy metering method and system for nonlinear load |
| CN117388569B (en) * | 2023-12-11 | 2024-03-01 | 浙江宏仁电气有限公司 | Electric energy metering method, electric energy metering box and medium under waveform distortion of power grid |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1700130B1 (en) * | 2003-12-30 | 2012-09-19 | Austriamicrosystems AG | Energy metering system |
| CN103245832A (en) * | 2013-05-16 | 2013-08-14 | 湖南大学 | Harmonic time frequency characteristic parameter estimating method based on fast S conversion and analysis meter |
| CN103308804A (en) * | 2013-06-17 | 2013-09-18 | 湖南大学 | Method for extracting time-frequency parameters of power quality disturbance signals on basis of fast K-S (Kaiser-S) transformation |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101759932B1 (en) * | 2010-04-27 | 2017-07-20 | 엘지전자 주식회사 | An estimating method for electric power consumption of an electric device |
-
2014
- 2014-07-16 CN CN201410339464.5A patent/CN104090159B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1700130B1 (en) * | 2003-12-30 | 2012-09-19 | Austriamicrosystems AG | Energy metering system |
| CN103245832A (en) * | 2013-05-16 | 2013-08-14 | 湖南大学 | Harmonic time frequency characteristic parameter estimating method based on fast S conversion and analysis meter |
| CN103308804A (en) * | 2013-06-17 | 2013-09-18 | 湖南大学 | Method for extracting time-frequency parameters of power quality disturbance signals on basis of fast K-S (Kaiser-S) transformation |
Non-Patent Citations (2)
| Title |
|---|
| "基于S变换解构算法的电能计量方法研究";丁忠安等;《第三届全国电能质量学术会议暨电能质量行业发展论坛论文集》;20130815;第133-141页 * |
| "快速S变换在电能质量分析中的应用";张志禹等;《电网技术》;20130531;第37卷(第5期);第1285-1290页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104090159A (en) | 2014-10-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104090159B (en) | Electric energy measuring method and device | |
| CN103869162B (en) | Dynamic signal phasor measurement method based on time domain quasi-synchronization | |
| CN107064628B (en) | High Precise Frequency Measurement System and method | |
| CN103257271A (en) | Device and method for detecting micro grid harmonic wave and inter-harmonics based on STM32F107VCT6 | |
| CN100516897C (en) | Harmonic energy metering method | |
| CN203287435U (en) | A micro electrical network harmonic wave and inter-harmonic wave test apparatus based on an STM32F107VCT6 | |
| CN107478896A (en) | A kind of frequency adaptive harmonic current detection method based on cascade Generalized Integrator | |
| CN103454497A (en) | Phase difference measuring method based on improved windowing discrete Fourier transform | |
| CN104535836A (en) | Electric signal fundamental frequency measuring method and system | |
| CN102662106A (en) | Method for measuring electric energy of harmonic power network | |
| CN110348161B (en) | Multi-frequency-band dynamic phasor electromagnetic transient simulation method and system for voltage source type converter | |
| CN104502675B (en) | Fundamental wave amplitude method and system of power signal | |
| CN106771592A (en) | A kind of trend of harmonic detection method of power | |
| CN102135552A (en) | Real-time digital detection method for active power and reactive power of electricity grid | |
| US20140108479A1 (en) | Computationally Efficient Finite Impulse Response Comb Filtering | |
| CN104880592A (en) | Electric energy calculating method and system under harmonic wave condition | |
| CN105334381A (en) | Method and device for measuring AC active power | |
| CN104111373A (en) | Metering method used for digital electric energy of intelligent substation | |
| Shui et al. | M-band compactly supported orthogonal symmetric interpolating scaling functions | |
| CN102520246A (en) | Constant frequency phasor extraction method | |
| CN105372471A (en) | Amplitude detection method and amplitude detection system for sine signal | |
| Mokeev et al. | Application of spectral representations in coordinates of complex frequency for the digital filter analysis and synthesis | |
| CN103105532A (en) | Harmonic electric energy measuring system of improved Fourier transform and control method thereof | |
| Wang et al. | Method for frequency measurement with ApFFT based on FPGA | |
| CN103308759A (en) | Active power measurement method and system based on digital filtering |
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 |