CN105425034A - Power utilization characteristic spectrum analysis method for specific AC components - Google Patents

Abstract

The invention relates to a power utilization characteristic spectrum analysis method for specific AC components. The method comprises that the sampling period is set as T second, the voltage and current of a power utilization device are sampled at the same time, and voltage and current sampling data at the same time is obtained; the number of the specific AC components is N, the frequency values of the AC components are f1, f2 to fN respectively, f1, f2 to fN are positive numbers that are increased progressively in order and not greater than 2pi/T, the DC components of the voltage and current as well as sine and cosine components of the AC voltage and current are directly obtained via a preestimation-correction algorithm; and amounts, including DC power as well as active power, reactive power, apparent power and the power factor of the AC component, reflecting the power utilization characteristic are calculated according to formulas. The method only includes arithmetic operation, is easy to program in computer, and has second-order precision and second-order convergence speed.

Description

A kind of electricity consumption characteristic spectrum analytical approach of specifying Alternating Component

Technical field

The invention belongs to quality of power supply assessment technique field, be specifically related to a kind of electricity consumption characteristic spectrum analytical approach of specifying Alternating Component, specifically comprise the method for solving of the active power of the power of flip-flop in alternating current and the alternating component of multiple given frequency, reactive power, applied power and power factor.

Background technology

In engineering, often need evaluating by electrical characteristics circuit and electronic electric equipment, for this reason, the active power to DC power and alternating component (given frequency), reactive power, applied power and power factor is needed to measure, be typically employed in tested voltage and the long-pending electrodymamometer metering method carrying out integral operation of electric current in fixed time length at present, such as, based on the method for Fourier transform.Under being only all the condition of positive integer in integral time with the ratio in the cycle of all alternating components of measured signal, adopt these class methods could realize the accurate measurement of electric power, and the DC power of equipment under test and circuit and the power parameter of each frequency alternating component can not be obtained simultaneously.

Summary of the invention

The object of the invention is, according to actual needs, a kind of electricity consumption characteristic spectrum analytical approach of specifying Alternating Component is provided, specifically, provide a kind of software programming simple, can directly apply in computer system and there is the electrodymamometer metering method of degree of precision and very fast speed of convergence.

Realizing technical scheme of the present invention is, a kind of electricity consumption characteristic spectrum analytical approach of specifying Alternating Component, and described method sets the sampling period as T second, samples simultaneously, obtain the sampled data of the voltage and current of synchronization to the voltage and current of consumer; Specify alternating component number to be N, the frequency values of each alternating component is increase progressively successively and be all not more than the positive number f of 2 π/T ₁, f ₂..., f _n, by estimating-correcting algorithm directly obtains the DC component of voltage and current, and the sinusoidal component of alternating voltage and electric current and cosine component; Then calculate the amount reflecting that it uses electrical characteristics, comprising: the active power of DC power and alternating component, reactive power, applied power and power factor.

Described method, for often organizing voltage sample data us and current sampling data is, performs in the steps below, and wherein the band bandwidth v of positive number is not more than 2 π/T;

[1] setup parameter T, N, f ₁, f ₂..., f _n, v value, setting voltage DC component u ₀, voltage sinusoidal component su ₁, su ₂..., su _n, voltage cosine component cu ₁, cu ₂..., cu _n; Estimated frequency, current dc component i ₀, current sinusoidal component si ₁, si ₂..., si _nwith electric current cosine component ci ₁, ci ₂..., ci _ninitial value, setting sampled data group number K, the initial value of loop control variable m;

[2] setting voltage DC component estimates initial value u ₀, current dc component estimates initial value i ₀, make subscript n respectively value be 1,2 ..., N, setting voltage sinusoidal component estimates initial value su _n, voltage cosine component estimates initial value cu _n, current sinusoidal component estimates initial value si _n, electric current cosine component estimates initial value ci _n;

[3] voltage sample data us and current sampling data is is read;

[4] discreet value is calculated: utilize formula (1) to obtain the discreet value of voltage DC component current dc component discreet value make subscript n respectively value be 1,2 ..., N, circulation execution formula (2), obtains the discreet value of voltage sinusoidal component the discreet value of voltage cosine component the discreet value of current sinusoidal component with the discreet value of electric current cosine component

$\left\{\begin{array}{c}{u}_0^p=v\·(us-u_0-T\·{\underset{\‾}{u}}_0-\underset{m=1}{\overset{N}{\Σ}}({\mathrm{su}}_m+T\·{\underset{\‾}{su}}_m\left)\right)\\ i_0^p=v\·(is-i_0-T\·{\underset{\‾}{i}}_0-\underset{m=1}{\overset{N}{\Σ}}({\mathrm{si}}_m+T\·{\underset{\‾}{si}}_m\left)\right)\end{array}\right.---\left(1\right)$

$\left\{\begin{array}{c}{\mathrm{su}}_n^p=u_0^p+f_n\·({\mathrm{cu}}_n+T\·{\underset{\‾}{cu}}_n)\\ {\mathrm{cu}}_n^p=-f_n\·({\mathrm{su}}_n+T\·{\underset{\‾}{su}}_n)\\ {\mathrm{si}}_n^p=i_0^p+f_n\·({\mathrm{ci}}_n+T\·{\underset{\‾}{ci}}_n)\\ {\mathrm{ci}}_n^p=-f_n\·({\mathrm{si}}_n+T\·{\underset{\‾}{si}}_n)\end{array}\right.---\left(2\right)$

[5] corrected value is calculated: obtain voltage DC component correction value according to formula (3) current dc component corrected value make subscript n respectively value be 1,2 ..., N, circulation execution formula (4) is respectively to voltage sinusoidal component correction value voltage cosine component corrected value current sinusoidal component correction value electric current cosine component corrected value

$\left\{\begin{array}{c}{u}_0^c=u_0+0.5T({\underset{\‾}{u}}_n+u_0^p)\\ i_0^c=i_0+0.5T({\underset{\‾}{i}}_0+i_0^p)\end{array}\right.---\left(3\right)$

$\left\{\begin{array}{c}{\mathrm{su}}_n^c={\mathrm{su}}_n+0.5T({\underset{\‾}{su}}_n+{\mathrm{su}}_n^p)\\ {\mathrm{cu}}_n^c={\mathrm{cu}}_n+0.5T({\underset{\‾}{cu}}_n+{\mathrm{cu}}_n^p)\\ {\mathrm{si}}_n^c={\mathrm{si}}_n+0.5T({\underset{\‾}{si}}_n+{\mathrm{si}}_n^p)\\ {\mathrm{ci}}_n^c={\mathrm{ci}}_n+0.5T({\underset{\‾}{ci}}_n+{\mathrm{ci}}_n^p)\end{array}\right.---\left(4\right)$

[6] corrected value is made to be estimated value: to obtain voltage DC component correction value u according to formula (5) ₀, current dc component corrected value i ₀, make subscript n respectively value be 1,2 ..., N, circulation execution formula (6) is respectively to voltage sinusoidal component correction value su _n, voltage cosine component corrected value cu _n, current sinusoidal component correction value si _n, electric current cosine component corrected value ci _n;

$\left\{\begin{array}{c}{u}_0=u_0^c\\ i_0=i_0^c\end{array}\right.---\left(5\right)$

$\left\{\begin{array}{c}{\mathrm{su}}_n={\mathrm{su}}_n^c\\ {\mathrm{cu}}_n={\mathrm{cu}}_n^c\\ {\mathrm{si}}_n={\mathrm{si}}_n^c\\ {\mathrm{ci}}_n={\mathrm{ci}}_n^c\end{array}\right.---\left(6\right)$

[7] initial value is estimated in renewal: according to the voltage DC component u obtained ₀, voltage sinusoidal component su ₁, su ₂..., su _n, voltage cosine component cu ₁, cu ₂..., cu _n, current dc component i ₀, current sinusoidal component si ₁, si ₂..., si _n, electric current cosine component ci ₁, ci ₂..., ci _nvalue, utilize formula (7) to upgrade voltage DC component and estimate initial value u ₀, current dc component estimates initial value i ₀, make subscript n respectively value be 1,2 ..., N, circulation execution formula (8) upgrades voltage sinusoidal component and estimates initial value su _n, voltage cosine component estimates initial value cu _n, current sinusoidal component estimates initial value si _nwith electric current cosine component increment ci _n;

$\left\{\begin{array}{c}{\underset{\‾}{u}}_0=v\·(us-u_0-\underset{m=1}{\overset{N}{\Σ}}{\mathrm{su}}_m)\\ {\underset{\‾}{i}}_0=v\·(is-i_0-\underset{m=1}{\overset{N}{\Σ}}{\mathrm{si}}_m)\end{array}\right.---\left(7\right)$

$\left\{\begin{array}{c}{\underset{\‾}{su}}_n={\underset{\‾}{u}}_0+f_n\·{\mathrm{cu}}_n\\ {\underset{\‾}{cu}}_n=-f_n\·{\mathrm{su}}_n\\ {\underset{\‾}{si}}_n={\underset{\‾}{i}}_0+f_n\·{\mathrm{ci}}_n\\ {\underset{\‾}{ci}}_n=-f_n\·{\mathrm{si}}_n\end{array}\right.---\left(8\right)$

[8] data are exported: according to the voltage DC component u obtained ₀, voltage sinusoidal component su ₁, su ₂..., su _n, voltage cosine component cu ₁, cu ₂..., cu _n, current dc component i ₀, current sinusoidal component si ₁, si ₂..., si _n, electric current cosine component ci ₁, ci ₂..., ci _nvalue, utilize formula (9) by DC power export be saved to P ₀[m], make subscript n respectively value be 1,2 ..., N, circulation execution formula (10), exports the active power of each alternating component, reactive power, applied power and power factor successively and is saved to P _n[m], Q _n[m], S _n[m]

P ₀[m]＝u ₀·i ₀(9)

[9] make m=m+1, if m≤K, then go to step [3], perform successively downwards to step [8], otherwise all sampled datas are processed complete, terminate.

The invention has the beneficial effects as follows,

1, the present invention is directed to the discrete time data obtained after timing sampling, adopt combined interior homotopy to calculate electric power, be convenient to utilize computing machine to realize.

2, the present invention directly obtains voltage DC component and current dc component with combined interior homotopy, and obtains DC power, has widened range of application.

3, the present invention directly obtains voltage sinusoidal component corresponding to known frequency, voltage cosine component, current sinusoidal component and electric current cosine component with combined interior homotopy, and obtain the active power of this frequency, reactive power, applied power and power factor, do not need to calculate sine function and cosine function, structure is simple, operand is little.

4, the invention belongs to second order finite impulse response (FIR) algorithm, have second order accuracy and second order convergence speed, comparatively infinite impulse response algorithm is easy to realize more, has higher precision and speed of convergence faster than single order finite impulse response (FIR) algorithm.

5, the present invention does not require that specified frequency keeps specific numerical relation, is convenient to harmonic wave and the m-Acetyl chlorophosphonazo composition of analytical voltage and electric current, and calculates harmonic power and m-Acetyl chlorophosphonazo power.

Accompanying drawing explanation

Fig. 1 is the schematic flow sheet of the electricity consumption characteristic spectrum analytical approach of the appointment Alternating Component of the embodiment of the present invention 1;

Fig. 2 is the main program flow schematic diagram of the electricity consumption characteristic spectrum analytical approach of the appointment Alternating Component of the embodiment of the present invention 2;

Fig. 3 is the interruption traffic program circuit schematic diagram of the electricity consumption characteristic spectrum analytical approach of the appointment Alternating Component of the embodiment of the present invention 2;

Fig. 4 is the DC power actual value curve over time of l-G simulation test;

Fig. 5 is the interchange active power 1 actual value curve over time of l-G simulation test;

Fig. 6 is the interchange reactive power 1 actual value curve over time of l-G simulation test;

Fig. 7 is the interchange active power 2 actual value curve over time of l-G simulation test;

Fig. 8 is the interchange reactive power 2 actual value curve over time of l-G simulation test;

Fig. 9 is the interchange active power 3 actual value curve over time of l-G simulation test;

Figure 10 is the interchange reactive power 3 actual value curve over time of l-G simulation test;

Figure 11 is the power tracking total error change curve of the electricity consumption characteristic spectrum analytical approach of appointment Alternating Component of the present invention;

Figure 12 is that the power factor of the electricity consumption characteristic spectrum analytical approach of appointment Alternating Component of the present invention follows the tracks of total error change curve.

Embodiment

Below in conjunction with accompanying drawing, specific embodiments of the invention are described in detail.

The electricity consumption characteristic spectrum analytical approach of appointment Alternating Component of the present invention comprises off-line analysis and on-line analysis two kinds of embodiments.

Embodiment 1:

Fig. 1 is the schematic flow sheet of the electricity consumption characteristic spectrum analytical approach of the appointment Alternating Component of embodiment 1, the electricity consumption characteristic spectrum analytical approach of the appointment Alternating Component of embodiment 1 is off-line analysis embodiment, and it is applicable to analyze the discrete-time series of the measured signals such as the voltage and current that the timing sampling be on request kept in memory device obtains.

The feature of off-line analysis embodiment is that sampling separates with analytical calculation and carries out.First with T second for the sampling period, timing sampling measured signal, is stored in gained sampled data in memory device, is formed with the discrete-time series of K group data, and then carries out analytical calculation to discrete-time series.

As shown in Figure 1, in the electricity consumption characteristic spectrum analytical approach of the appointment Alternating Component of embodiment 1, analytical calculation process comprises initialization step and data processing step, and wherein data processing step comprises the parts such as combined interior homotopy step, data output and cycle control step.

In initialization step (i.e. following step [1] ~ [2]), first set discrete-time series data amount check K, the location number m of data in discrete-time series that pre-treatment is worked as in setting is 1; Then setup parameter T, N, f ₁, f ₂..., f _n, v value; Setting voltage DC component u ₀, voltage sinusoidal component su ₁, su ₂..., su _n, voltage cosine component cu ₁, cu ₂..., cu _n, direct-current component i ₀, current sinusoidal component si ₁, si ₂..., si _nwith electric current cosine component ci ₁, ci ₂..., ci _ninitial value; Reset voltage DC component and estimate initial value u ₀, current dc component estimates initial value i ₀, voltage sinusoidal component estimates initial value su ₁, su ₂..., su _n, voltage cosine component estimates initial value cu ₁, cu ₂..., cu _n, current sinusoidal component estimates initial value si ₁, si ₂..., si _n, electric current cosine component estimates initial value ci ₁, ci ₂..., ci _n; Then enter data processing step to perform.

In combined interior homotopy step, first read the m group voltage data in discrete-time series and current data, respectively as the voltage sample data us and the current sampling data is that work as pre-treatment; Then combined interior homotopy step (i.e. following step [3] ~ [7]) is performed successively;

In output data (following step [8]), according to estimating-correction process after voltage DC component u ₀, voltage sinusoidal component su ₁, su ₂..., su _n, voltage cosine component cu ₁, cu ₂..., cu _n, current dc component i ₀, current sinusoidal component si ₁, si ₂..., si _n, electric current cosine component ci ₁, ci ₂..., ci _n, calculate the active power of DC power and alternating component, reactive power, applied power and power factor.For ease of Algorithm Analysis, the numerical value of DC power, active power, reactive power, applied power, power factor is saved in m the position exporting series; Then calculate and estimate initial value and enter cycle control step.

In cycle control step (following step [9]), first location number m is increased by 1, then judge whether to return execution data combined interior homotopy step according to the value of the group number K of data in location number m and discrete-time series.If m≤K, return and perform data combined interior homotopy step; If m > is K, represents that in discrete-time series, all sampled datas are all processed complete, then stop running, terminate offline analytic process.

Embodiment 2:

Fig. 2 and Fig. 3 is the schematic flow sheet of the electricity consumption characteristic spectrum analytical approach of the appointment Alternating Component of embodiment 2.Fig. 2 is the main program flow schematic diagram of the electricity consumption characteristic spectrum analytical approach of the appointment Alternating Component of embodiment 2, and Fig. 3 is the interruption traffic program circuit schematic diagram of the electricity consumption characteristic spectrum analytical approach of the appointment Alternating Component of embodiment 2.

The electricity consumption characteristic spectrum analytical approach of the appointment Alternating Component of embodiment 2 is on-line analysis embodiment, it is applicable to carry out real-time analysis process to the often group sampled data of measured signal, the feature of on-line analysis embodiment is sampled while analytical calculation, namely often to sample one group of data, just carry out an analytical calculation.

The electricity consumption characteristic spectrum analytical approach of the appointment Alternating Component of embodiment 2 comprises master routine and Interruption service routine step.

Master routine comprises initialization step and realizes other the major cycle step of function.In initialization step (i.e. following step [1] ~ [2]), first setup parameter T, N, f ₁, f ₂..., f _n, v value, setting sampled data group number K, the initial value of loop control variable m; Setting voltage DC component u ₀, voltage sinusoidal component su ₁, su ₂..., su _n, voltage cosine component cu ₁, cu ₂..., cu _n, current dc component i ₀, current sinusoidal component si ₁, si ₂..., si _n, electric current cosine component ci ₁, ci ₂..., ci _ninitial value; Setting voltage DC component estimates initial value u ₀, current dc component estimates initial value i ₀, voltage sinusoidal component estimates initial value su ₁, su ₂..., su _n, voltage cosine component estimates initial value cu ₁, cu ₂..., cu _n, current sinusoidal component estimates initial value si ₁, si ₂..., si _n, electric current cosine component estimates initial value ci ₁, ci ₂..., ci _n; Then the timing resetting timer is T second, and T requires the sampling period for meeting Shannon's sampling theorem, and the Interruption of open system.

In Interruption service step, first preserve the currency of each register of interrupt spot; Then tested voltage and current signal is sampled simultaneously, obtain voltage sample data us and current sampling data is, perform following step [3] ~ [7] successively; Recover the value of each register of interrupt spot again; Then interrupt turning back to master routine and perform following step [8] ~ [9].

On-line analysis embodiment produces Interruption event by timer circulation, causes interrupt service routine to circulate and performs.

The algorithm steps described in electricity consumption characteristic spectrum analytical approach of the appointment Alternating Component of embodiment 1 and embodiment 2 is as follows:

No matter be off-line analysis or on-line analysis, for voltage sample data us and the current sampling data is of measured signal, perform all in the steps below.

Algorithm comprises the following steps successively:

[1] setup parameter T, N, f ₁, f ₂..., f _n, v value, setting voltage DC component u ₀, voltage sinusoidal component su ₁, su ₂..., su _n, voltage cosine component cu ₁, cu ₂..., cu _n, estimated frequency, current dc component i ₀, current sinusoidal component si ₁, si ₂..., si _nwith electric current cosine component ci ₁, ci ₂..., ci _ninitial value, setting sampled data group number K, the initial value of loop control variable m;

[2] setting voltage DC component estimates initial value u ₀, current dc component estimates initial value i ₀, make subscript n respectively value be 1,2 ..., N, setting voltage sinusoidal component estimates initial value su _n, voltage cosine component estimates initial value cu _n, current sinusoidal component estimates initial value si _n, electric current cosine component estimates initial value ci _n;

[3] voltage sample data us and current sampling data is is read;

[4] discreet value is calculated: utilize formula (1) to obtain the discreet value of voltage DC component current dc component discreet value make subscript n respectively value be 1,2 ..., N, circulation execution formula (2), obtains the discreet value of voltage sinusoidal component the discreet value of voltage cosine component the discreet value of current sinusoidal component with the discreet value of electric current cosine component

$\left\{\begin{array}{c}{u}_0^p=v\·(us-u_0-T\·{\underset{\‾}{u}}_0-\underset{m=1}{\overset{N}{\Σ}}({\mathrm{su}}_m+T\·{\underset{\‾}{su}}_m\left)\right)\\ i_0^p=v\·(is-i_0-T\·{\underset{\‾}{i}}_0-\underset{m=1}{\overset{N}{\Σ}}({\mathrm{si}}_m+T\·{\underset{\‾}{si}}_m\left)\right)\end{array}\right.---\left(1\right)$

$\left\{\begin{array}{c}{\mathrm{su}}_n^p=u_0^p+f_n\·({\mathrm{cu}}_n+T\·{\underset{\‾}{cu}}_n)\\ {\mathrm{cu}}_n^p=-f_n\·({\mathrm{su}}_n+T\·{\underset{\‾}{su}}_n)\\ {\mathrm{si}}_n^p=i_0^p+f_n\·({\mathrm{ci}}_n+T\·{\underset{\‾}{ci}}_n)\\ {\mathrm{ci}}_n^p=-f_n\·({\mathrm{si}}_n+T\·{\underset{\‾}{si}}_n)\end{array}\right.---\left(2\right)$

[5] corrected value is calculated: obtain voltage DC component correction value according to formula (3) current dc component corrected value make subscript n respectively value be 1,2 ..., N, circulation execution formula (4) is respectively to voltage sinusoidal component correction value voltage cosine component corrected value current sinusoidal component correction value electric current cosine component corrected value

$\left\{\begin{array}{c}{u}_0^c=u_0+0.5T({\underset{\‾}{u}}_n+u_0^p)\\ i_0^c=i_0+0.5T({\underset{\‾}{i}}_0+i_0^p)\end{array}\right.---\left(3\right)$

$\left\{\begin{array}{c}{\mathrm{su}}_n^c={\mathrm{su}}_n+0.5T({\underset{\‾}{su}}_n+{\mathrm{su}}_n^p)\\ {\mathrm{cu}}_n^c={\mathrm{cu}}_n+0.5T({\underset{\‾}{cu}}_n+{\mathrm{cu}}_n^p)\\ {\mathrm{si}}_n^c={\mathrm{si}}_n+0.5T({\underset{\‾}{si}}_n+{\mathrm{si}}_n^p)\\ {\mathrm{ci}}_n^c={\mathrm{ci}}_n+0.5T({\underset{\‾}{ci}}_n+{\mathrm{ci}}_n^p)\end{array}\right.---\left(4\right)$

[6] corrected value is made to be estimated value: to obtain voltage DC component correction value u according to formula (5) ₀, current dc component corrected value i ₀, make subscript n respectively value be 1,2 ..., N, circulation execution formula (6) is respectively to voltage sinusoidal component correction value su _n, voltage cosine component corrected value cu _n, current sinusoidal component correction value si _n, electric current cosine component corrected value ci _n;

$\left\{\begin{array}{c}{u}_0=u_0^c\\ i_0=i_0^c\end{array}\right.---\left(5\right)$

$\left\{\begin{array}{c}{\mathrm{su}}_n={\mathrm{su}}_n^c\\ {\mathrm{cu}}_n={\mathrm{cu}}_n^c\\ {\mathrm{si}}_n={\mathrm{si}}_n^c\\ {\mathrm{ci}}_n={\mathrm{ci}}_n^c\end{array}\right.---\left(6\right)$

[7] initial value is estimated in renewal: according to the voltage DC component u obtained ₀, voltage sinusoidal component su ₁, su ₂..., su _n, voltage cosine component cu ₁, cu ₂..., cu _n, current dc component i ₀, current sinusoidal component si ₁, si ₂..., si _n, electric current cosine component ci ₁, ci ₂..., ci _nvalue, utilize formula (7) to upgrade voltage DC component and estimate initial value u ₀, current dc component estimates initial value i ₀, make subscript n respectively value be 1,2 ..., N, circulation execution formula (8) upgrades voltage sinusoidal component and estimates initial value su _n, voltage cosine component estimates initial value cu _n, current sinusoidal component estimates initial value si _nwith electric current cosine component increment ci _n;

$\left\{\begin{array}{c}{\underset{\‾}{u}}_0=v\·(us-u_0-\underset{m=1}{\overset{N}{\Σ}}{\mathrm{su}}_m)\\ {\underset{\‾}{i}}_0=v\·(is-i_0-\underset{m=1}{\overset{N}{\Σ}}{\mathrm{si}}_m)\end{array}\right.---\left(7\right)$

$\left\{\begin{array}{c}{\underset{\‾}{su}}_n={\underset{\‾}{u}}_0+f_n\·{\mathrm{cu}}_n\\ {\underset{\‾}{cu}}_n=-f_n\·{\mathrm{su}}_n\\ {\underset{\‾}{si}}_n={\underset{\‾}{i}}_0+f_n\·{\mathrm{ci}}_n\\ {\underset{\‾}{ci}}_n=-f_n\·{\mathrm{si}}_n\end{array}\right.---\left(8\right)$

[8] data are exported: according to the voltage DC component u obtained ₀, voltage sinusoidal component su ₁, su ₂..., su _n, voltage cosine component cu ₁, cu ₂..., cu _n, current dc component i ₀, current sinusoidal component si ₁, si ₂..., si _n, electric current cosine component ci ₁, ci ₂..., ci _nvalue, utilize formula (9) by DC power export be saved to P ₀[m], make subscript n respectively value be 1,2 ..., N, circulation execution formula (10), exports the active power of each alternating component, reactive power, applied power and power factor successively and is saved to P _n[m], Q _n[m], S _n[m] and

P ₀[m]＝u ₀·i ₀(9)

[9] make m=m+1, if m≤K, then go to step [3], perform successively downwards to step [8], otherwise all sampled datas are processed complete, terminate.

In above-mentioned two kinds of embodiments, sampling period T, interchange number N, frequency f ₁, f ₂..., f _n, parameter v numerical value, all can require to set according to the priori of measured signal and signal analysis.Under first sampling period T will meet Shannon's sampling theorem requirement, the sampling period is less, and analysis precision is higher.For on-line analysis mode, by the restriction of real-time, exchange number N and sampling period T and also should meet (N+1) Δ t < T, wherein Δ t represents the maximum time required for execution Interruption step.

The electricity consumption characteristic spectrum analytical approach of appointment Alternating Component of the present invention belongs to Two Order Method, error between the steady-state value of each combined interior homotopy variable and its actual value, relevant to the quadratic power of sampling period T, therefore the sampling period is less, and analysis precision is higher.

Because data in computer system are limited wordlengths, saturated for avoiding occurring in combined interior homotopy process, to DC component, each alternating component just, cosine component also can carry out amplitude limiting processing.

The physical significance of parameter v is equivalent to frequency f successively ₁, f ₂..., f _nthe bandwidth of the passband at place, according to Shannon's sampling theorem, must meet v≤2 π/T.The speed of convergence of v numerical values recited to estimated amplitude has major effect, and v value is larger, and estimated value more quickly converges to actual value, but increases the harmful effect of interference to estimated accuracy simultaneously.

If frequency f ₁, f ₂..., f _nequal the actual frequency of measured signal alternating component respectively, as the voltage DC component u of combined interior homotopy variable ₀, voltage sinusoidal component su ₁, su ₂..., su _n, voltage cosine component cu ₁, cu ₂..., cu _n, current dc component i ₀, current sinusoidal component si ₁, si ₂..., si _n, electric current cosine component ci ₁, ci ₂..., ci _ncan converge to respective actual value respectively, the initial value therefore for these variablees is not particularly limited.Preferably, all 0 is set as.

For estimating initial value, voltage DC component estimates initial value u ₀, current dc component estimates initial value i ₀, voltage sinusoidal component estimates initial value su ₁, su ₂..., su _n, voltage cosine component estimates initial value cu ₁, cu ₂..., cu _n, current sinusoidal component estimates initial value si ₁, si ₂..., si _n, electric current cosine component estimates initial value ci ₁, ci ₂..., ci _n, be not particularly limited.Preferably, all 0 is set as.

If tested voltage signal is expressed as u (t)=U ₀+ U ₁sin (w ₁t+ α ₁)+U ₂sin (w ₂t+ α ₂)+... + U _nsin (w _nt+ α _n), tested current signal representation is i (t)=I ₀+ I ₁sin (w ₁t+ β ₁)+I ₂sin (w ₂t+ β ₂)+... + I _nsin (w _nt+ β _n), and assigned frequency f ₁, f ₂..., f _nequal actual frequency w respectively successively ₁, w ₂..., w _n, then after said method process, voltage DC component u ₀, current dc component i ₀converge to U respectively ₀, I ₀, voltage sinusoidal component su _n, voltage cosine component cu _n, current sinusoidal component si _n, electric current cosine component ci _nconverge to U respectively _nsin (w _nt+ α _n), U _ncos (w _nt+ α _n), I _nsin (w _nt+ β _n), I _ncos (w _nt+ β _n), DC power P ₀converge to U ₀i ₀, active-power P _n, reactive power Q _n, applied power S _nconverge to 0.5U respectively _ni _ncos (α _n-β _n), 0.5U _ni _nsin (α _n-β _n), 0.5U _ni _n, as applied power S _nwhen being not equal to 0, power factor converge to cos (α _n-β _n), wherein subscript n respectively value be 1,2 ..., N.

The validity of the electricity consumption characteristic spectrum analytical approach of appointment Alternating Component of the present invention is described below in conjunction with example.

Such as: tested voltage signal is u=U ₀+ U ₁sin (w ₁t+ α ₁)+U ₂sin (w ₂t+ α ₂)+U ₃sin (w ₃t+ α ₃), tested current signal is i=I ₀+ I ₁sin (w ₁t+ β ₁)+I ₂sin (w ₂t+ β ₂)+I ₃sin (w ₃t+ β ₃), wherein three a-c cycle w ₁, w ₂, w ₃numerical value equals 20 π, 45 π, 71 π respectively, unit is radian per second, all the other parameters in time t change as shown in Table 1, wherein π is circular constant, and the change of the actual value of corresponding DC power, active power, reactive power is as shown in Fig. 4 to Figure 10.

For effect of the electricity consumption characteristic spectrum analytical approach of appointment Alternating Component of the present invention is described by change curve, definition power tracking total error err _pfor:

err _p＝|U ₀I ₀-P ₀|+|0.5U ₁I ₁cos(α ₁-β ₁)-P ₁|+|0.5U ₂I ₂cos(α ₂-β ₂)-P ₂|+|0.5U ₃I ₃cos(α ₃-β ₃)-P ₃|

+|0.5U ₁I ₁sin(α ₁-β ₁)-Q ₁|+|0.5U ₂I ₂sin(α ₂-β ₂)-Q ₂|+|0.5U ₃I ₃sin(α ₃-β ₃)-Q ₃|

+|0.5U ₁I ₁-S ₁|+|0.5U ₂I ₂-S ₂|+|0.5U ₃I ₃-S ₃|

Definition power factor follows the tracks of total error err _cfor:

The t change in time of table 1 measured signal parameters

If sampling period T=0.1 millisecond, first with T second for the sampling period, timing is sampled to tested voltage and current signal simultaneously, obtains sampled data and forms discrete-time series, again according to the off-line analysis embodiment shown in Fig. 1, coding is simulation run in a computer.Setting N=3, assigned frequency f ₁, f ₂, f ₃value equal w successively ₁, w ₂, w ₃, setting voltage DC component u ₀, voltage sinusoidal component su ₁, su ₂, su ₃, voltage cosine component cu ₁, cu ₂, cu ₃, current dc component i ₀, current sinusoidal component si ₁, si ₂, si ₃, electric current cosine component ci ₁, ci ₂, ci ₃initial value be 0.Setting voltage DC component estimates initial value u ₀, current dc component estimates initial value i ₀, voltage sinusoidal component estimates initial value su ₁, su ₂, su ₃, voltage cosine component estimates initial value cu ₁, cu ₂, cu ₃, current sinusoidal component estimates initial value si ₁, si ₂, si ₃, electric current cosine component estimates initial value ci ₁, ci ₂, ci ₃, be 0.

Figure 11 and Figure 12 respectively depict simulation run gained power tracking total error err under different v values _pand power factor follows the tracks of total error err _c.Illustrate when the amplitude of voltage and current and phase angle generation saltus step, the DC power that the electricity consumption characteristic spectrum analytical approach of appointment Alternating Component of the present invention is measured, the active power of each alternating component, reactive power, applied power and power factor can converge to respective actual value respectively.Meanwhile, speed of convergence is all decided by the numerical value of parameter v, and the numerical value of parameter v is larger, and speed of convergence is faster.

Above embodiment is only exemplary embodiment of the present invention, and be not used in restriction the present invention, protection scope of the present invention is defined by the claims.Those skilled in the art can in essence of the present invention and protection domain, and make various amendment or equivalent replacement to the present invention, this amendment or equivalent replacement also should be considered as dropping in protection scope of the present invention.

Claims (2)

1. specify an electricity consumption characteristic spectrum analytical approach for Alternating Component, it is characterized in that, described method sets the sampling period as T second, samples simultaneously, obtain the sampled data of the voltage and current of synchronization to the voltage and current of consumer; Specify alternating component number to be N, the frequency values of each alternating component for increase progressively successively, and is all not more than the positive number f of 2 π/T ₁, f ₂..., f _n, by estimating-correcting algorithm directly obtains the DC component of voltage and current, and the sinusoidal component of alternating voltage and electric current and cosine component; Then calculate the amount reflecting that it uses electrical characteristics, comprising: the active power of DC power and alternating component, reactive power, applied power and power factor.

2. a kind of electricity consumption characteristic spectrum analytical approach of specifying Alternating Component according to claim 1, it is characterized in that, described method, for often organizing voltage sample data us and current sampling data is, performs in the steps below, and wherein the band bandwidth v of positive number is not more than 2 π/T;

[1] setup parameter T, N, f ₁, f ₂..., f _n, v value, setting voltage DC component u ₀, voltage sinusoidal component su ₁, su ₂..., su _n, voltage cosine component cu ₁, cu ₂..., cu _n; Estimated frequency, current dc component i ₀, current sinusoidal component si ₁, si ₂..., si _nwith electric current cosine component ci ₁, ci ₂..., ci _ninitial value, setting sampled data group number K, the initial value of loop control variable m;

[2] setting voltage DC component estimates initial value u ₀, current dc component estimates initial value i ₀, make subscript n respectively value be 1,2 ..., N, setting voltage sinusoidal component estimates initial value su _n, voltage cosine component estimates initial value cu _n, current sinusoidal component estimates initial value si _n, electric current cosine component estimates initial value ci _n;

[3] voltage sample data us and current sampling data is is read;

[4] discreet value is calculated: utilize formula (1) to obtain the discreet value of voltage DC component current dc component discreet value make subscript n respectively value be 1,2 ..., N, circulation execution formula (2), obtains the discreet value of voltage sinusoidal component the discreet value of voltage cosine component the discreet value of current sinusoidal component with the discreet value of electric current cosine component

$\left\{\begin{array}{c}{u}_0^p=v\&CenterDot;(us-u_0-T\&CenterDot;{\underset{\&OverBar;}{u}}_0-\underset{m=1}{\overset{N}{\&Sigma;}}({\mathrm{su}}_m+T\&CenterDot;{\underset{\&OverBar;}{su}}_m\left)\right)\\ i_0^p=v\&CenterDot;(is-i_0-T\&CenterDot;{\underset{\&OverBar;}{i}}_0-\underset{m=1}{\overset{N}{\&Sigma;}}({\mathrm{si}}_m+T\&CenterDot;{\underset{\&OverBar;}{si}}_m\left)\right)\end{array}\right.---\left(1\right)$

$\left\{\begin{array}{c}{\mathrm{su}}_n^p=u_0^p+f_n\&CenterDot;({\mathrm{cu}}_n+T\&CenterDot;{\underset{\&OverBar;}{cu}}_n)\\ {\mathrm{cu}}_n^p=-f_n\&CenterDot;({\mathrm{su}}_n+T\&CenterDot;{\underset{\&OverBar;}{su}}_n)\\ {\mathrm{si}}_n^p=i_0^p+f_n\&CenterDot;({\mathrm{ci}}_n+T\&CenterDot;{\underset{\&OverBar;}{ci}}_n)\\ {\mathrm{ci}}_n^p=-f_n\&CenterDot;({\mathrm{si}}_n+T\&CenterDot;{\underset{\&OverBar;}{si}}_n)\end{array}\right.---\left(2\right)$

[5] corrected value is calculated: obtain voltage DC component correction value according to formula (3) current dc component corrected value make subscript n respectively value be 1,2 ..., N, circulation execution formula (4) is respectively to voltage sinusoidal component correction value voltage cosine component corrected value current sinusoidal component correction value electric current cosine component corrected value

$\left\{\begin{array}{c}{u}_0^c=u_0+0.5T({\underset{\&OverBar;}{u}}_0+u_0^p)\\ i_0^c=i_0+0.5T({\underset{\&OverBar;}{i}}_0+i_0^p)\end{array}\right.---\left(3\right)$

[6] corrected value is made to be estimated value: to obtain voltage DC component correction value u according to formula (5) ₀, current dc component corrected value i ₀, make subscript n respectively value be 1,2 ..., N, circulation execution formula (6) is respectively to voltage sinusoidal component correction value su _n, voltage cosine component corrected value cu _n, current sinusoidal component correction value si _n, electric current cosine component corrected value ci _n;

$\left\{\begin{array}{c}{u}_0=u_0^c\\ i_0=i_0^c\end{array}\right.---\left(5\right)$

[7] initial value is estimated in renewal: according to the voltage DC component u obtained ₀, voltage sinusoidal component su ₁, su ₂..., su _n, voltage cosine component cu ₁, cu ₂..., cu _n, current dc component i ₀, current sinusoidal component si ₁, si ₂..., si _n, electric current cosine component ci ₁, ci ₂..., ci _nvalue, utilize formula (7) to upgrade voltage DC component and estimate initial value u ₀, current dc component estimates initial value i ₀, make subscript n respectively value be 1,2 ..., N, circulation execution formula (8) upgrades voltage sinusoidal component and estimates initial value su _n, voltage cosine component estimates initial value cu _n, current sinusoidal component estimates initial value si _nwith electric current cosine component increment ci _n;

$\left\{\begin{array}{c}{\underset{\&OverBar;}{u}}_0=v\&CenterDot;(us-u_0-\underset{m=1}{\overset{N}{\&Sigma;}}{\mathrm{su}}_m)\\ {\underset{\&OverBar;}{i}}_0=v\&CenterDot;(is-i_0-\underset{m=1}{\overset{N}{\&Sigma;}}{\mathrm{si}}_m)\end{array}\right.---\left(7\right)$

$\left\{\begin{array}{c}{\underset{\&OverBar;}{su}}_n={\underset{\&OverBar;}{u}}_0+f_n\&CenterDot;{\mathrm{cu}}_n\\ {\underset{\&OverBar;}{cu}}_n=-f_n\&CenterDot;{\mathrm{su}}_n\\ {\underset{\&OverBar;}{si}}_n={\underset{\&OverBar;}{i}}_0+f_n\&CenterDot;{\mathrm{ci}}_n\\ {\underset{\&OverBar;}{ci}}_n=-f_n\&CenterDot;{\mathrm{si}}_n\end{array}\right.---\left(8\right)$

[8] data are exported: according to the voltage DC component u obtained ₀, voltage sinusoidal component su ₁, su ₂..., su _n, voltage cosine component cu ₁, cu ₂..., cu _n, current dc component i ₀, current sinusoidal component si ₁, si ₂..., si _n, electric current cosine component ci ₁, ci ₂..., ci _nvalue, utilize formula (9) by DC power export be saved to P ₀[m], make subscript n respectively value be 1,2 ..., N, circulation execution formula (10), exports the active power of each alternating component, reactive power, applied power and power factor successively and is saved to P _n[m], Q _n[m], S _n[m] and

P ₀[m]＝u ₀·i ₀(9)

[9] make m=m+1, if m≤K, then go to step [3], perform successively downwards to step [8], otherwise all sampled datas are processed complete, terminate.