CN104090159B - Electric energy measuring method and device - Google Patents

Electric energy measuring method and device Download PDF

Info

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
Application number
CN201410339464.5A
Other languages
Chinese (zh)
Other versions
CN104090159A (en
Inventor
郑可
欧习洋
侯兴哲
周全
纪静
彭鹏
雷万钧
丁忠安
Original Assignee
国家电网公司
国网重庆市电力公司电力科学研究院
西安交通大学
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 国家电网公司, 国网重庆市电力公司电力科学研究院, 西安交通大学 filed Critical 国家电网公司
Priority to CN201410339464.5A priority Critical patent/CN104090159B/en
Publication of CN104090159A publication Critical patent/CN104090159A/en
Application granted granted Critical
Publication of CN104090159B publication Critical patent/CN104090159B/en

Links

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

Electric energy gauging method and device

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.
CN201410339464.5A 2014-07-16 2014-07-16 Electric energy measuring method and device CN104090159B (en)

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 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 (3)

* Cited by examiner, † Cited by third party
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

Citations (3)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (3)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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
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 real-time fast discrete S-transform for cross-differential protection of shunt-compensated 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. Time-fractional 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 shift-invariant 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 self-convolution window dual-spectrum line interpolation FFT (Fast Fourier Transform) and device thereof
CN104459321B (en) Power signal base wave phase measurement method and system
CN105115594B (en) Gear-box vibration signal fault signature extracting method based on Wavelet Entropy and information fusion
Lovell et al. The relationship between instantaneous frequency and time-frequency 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 non-stationary 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