CN101261292A - Base wave and harmonic detection method based on fiver item Rife-Vincent(1)window double spectral line interpolation FFT - Google Patents

Abstract

The invention discloses a fundamental wave and harmonic wave detecting method based on five Rife-Vincent (I) windows doublet interpolation FFT (fast Fourier transform): detected signals are processed with analog-to-digital conversion and low-pass filtering, the five Rife-Vincent (I) windows which adopt a low sidelobe level and a rapid decay rate of the sidelobe carry out window-adding FFT operation on the detected signals, and the doublet interpolation calculating method is used for detecting the quantity of the fundamental wave and each harmonic wave of the detected signals. The fundamental wave and harmonic wave detecting method includes the following aspects: (1) analog-to-digital conversion and digital low-pass filtering; (2) FFT calculation based on the five Rife-Vincent (I) windows; (3) the doublet interpolation and the fitting of polynomial to obtain the parameter of the fundamental wave and the harmonic wave.

Description

Based on 5 Rife-Vincent(I) first-harmonic and the harmonic detecting method of window double spectral line interpolation FFT

Affiliated technical field

The present invention relates to a kind of method that is used for signal first-harmonic and harmonic parameters detection.Specifically be a kind of first-harmonic and harmonic parameters detection method, belong to the signal processing technology field, also can be used for the analysis of spectrum of other signal based on 5 Rife-Vincent (I) window double spectral line interpolation FFTs (fast fourier transform).

Background technology

With the electric power signal detection is example, and along with the development of Power Electronic Technique, non-linear equipment is used widely in electric system, and the harmonic problem that non-linear equipment brings is on the rise to power system security, threat stable, that economical operation constitutes.The accurate measurement of first-harmonic and harmonic wave can be electrical network electric energy metrical, harmonic trend calculating, equipment networking detection, Harmonious Waves in Power Systems compensation and inhibition etc. scientific basis is provided.

First-harmonic and harmonic component detection method generally have: (1) based on the FFT method of frequency-domain analysis, when being characterized in the mains frequency fluctuation, accuracy of detection is lower; (2) the humorous wave method of filter detection is characterized in being difficult to obtain the ideal frequency characteristic; (3) based on the method for instantaneous reactive power theory, be characterized in that calculated amount is big, handle complicated; (4) based on neural network theory and method of wavelet, be characterized in that calculated amount is big, be difficult in embedded system, realize.

Based on first-harmonic, the algorithm of harmonics analysis of FFT, be easy to go up realization at DSP (digital signal processor), be current most widely used a kind of harmonic analysis method.In electric system, line voltage, current distortion cause fundamental frequency to change, because non-synchronous sampling causes fft algorithm to have problems such as spectral leakage and fence effect, the precision that makes first-harmonic and harmonic wave detect is affected.First-harmonic and each harmonic wave accuracy of detection are the difficult problems in electric measurement signal analysis and the quality of power supply management in the spectral leakage of minimizing fft algorithm and fence effect influence, the raising electrical measurement.

Existing patent document " method of measuring factory harmonic wave and measuring instrument " (200310105446.2), " method and the gauge of metering factory harmonic wave " (200310105446.2), " Harmonious Waves in Power Systems measuring fixed amount method and gauge " (98110414.2) etc.The purpose of its invention is the dynamic perfromance that harmonic voltage, electric current in the electric system are measured in accurate quantification ground.The weak point of prior art is: prior art provides the electrical equipment method for designing of some first-harmonics and harmonic measure, but because harmonic wave has features such as diversity, randomness and polymorphism, spectral leakage and fence effect that the fundamental frequency fluctuation causes still exist, and therefore real-time, high-precision first-harmonic and harmonic measure analysis are difficult to realize.

Summary of the invention

For overcoming the deficiency of prior art, the object of the present invention is to provide double spectral line interpolation FFT method, first-harmonic and each harmonic component in quick, the accurate detection signal of this method energy based on 5 Rife-Vincent (I) window.

If comprise the time domain measured signal of multinomial integral frequency harmonizing wave be

Wherein, M is the item number of contained harmonic wave; r _mRepresent the number of times of harmonic wave; When m=0, r ₀=1, f ₀, A ₀With

Be respectively frequency, amplitude and the initial phase angle of first-harmonic; When m ≠ 0, r _m, A _mWith

Be respectively number of times, amplitude and the initial phase angle of m item harmonic wave.Signal x (t) is f through over-sampling rate _sData acquisition system (DAS) after, obtain discrete series

According to the definition of Fourier transform, the continuous Fourier transform of signal x (n) is

$X\left(f\right)=\underset{n=-\&infin;}{\overset{\&infin;}{\mathrm{\&Sigma;}}}x\left(n\right)e^{-\mathrm{<figure-callout\; id="2"\; label="j"\; filenames="A20081003106600131.png"\; state="\{\{state\}\}">j</figure-callout>}2\mathrm{\&pi;fn}}$

When realizing on embedded system or on the computing machine, data are always time-limited, promptly be equivalent to signal x (n) by a window function w (n) (n=0,1 ... N-1) be punctured into the long sequence of N point.After the windowing, signal spectrum is become with f by original spectral line ₀Be the continuous spectrum of middle mind-set both sides expansion, the phase mutual interference between the signal spectrum secondary lobe, energy is leaked to whole frequency band, is spectral leakage.

To X (f) with discrete frequency interval delta f=f _s/ N carries out discretize, can obtain X (k Δ f) (k=0,1 ... N-1 represents the sequence number of discrete spectrum line).Because there is necessarily fluctuation in mains frequency and non-constant, and may has a harmonic wave (being non-integer harmonics) in the measured signal, thereby be difficult to strictness and accomplish synchronized sampling.Non-synchronous sampling causes crest frequency to be difficult to just in time be positioned on the discrete spectral line frequency, i.e. f ₀=k ₀Δ f (k ₀Be non-integer), the observed deviation to signal peak value point that causes is fence effect.

Based on above analysis, first-harmonic and harmonic detecting method that the present invention proposes based on 5 Rife-Vincent (I) window, it is characterized in that: analog to digital conversion and low-pass filtering treatment are carried out in measured signal, adopt 5 Rife-Vincent (I) window that sidelobe level is low, the rate of decay is fast that signal is carried out windowing FFT computing, utilization double spectral line interpolation algorithm detects first-harmonic and each harmonic content in the measured signal.It comprises the following steps:

(1) analog to digital conversion and low-pass filtering

High-speed AD converter is converted to digital quantity with measured signal (analog quantity) x (t), and sample frequency is f _S, this digital quantity detects the input quantity of handling as first-harmonic and harmonic wave.According to the frequency analysis accuracy requirement, the higher hamonic wave number of times and the frequency f that determine to need detect _MAX, the design cutoff frequency is for being a bit larger tham f _MAXLowpass digital filter, input quantity obtains removing the input quantity x (n) of high frequency component signal through behind the wave digital lowpass filter.

In the formula, M is the number of times of contained higher harmonics.

(2) make up discrete 5 Rife-Vincent (I) window

5 Rife-Vincent (I) window is a kind of cosine composite window, and its discrete expression is

$w\left(n\right)=\underset{m=0}{\overset{M-1}{\mathrm{\&Sigma;}}}{(-1)}^m{a}_m\cos (2\mathrm{\&pi;n}\&CenterDot;m/N)$

In the formula, M is the item number of 5 Rife-Vincent (I) window function; N is the data length of FFT computing, also be window function length (n=1,2 ... N-1); a _mSatisfy constraint condition

$\underset{m=0}{\overset{M-1}{\mathrm{\&Sigma;}}}{(-1)}^m{a}_m=0$

The coefficient of 5 Rife-Vincent (I) window is respectively a ₀=1, a ₁=1.6, a ₂=0.8, a ₃=0.22857, a ₄=0.02857, its discrete time-domain characteristic as shown in Figure 2.Fig. 3 has provided the amplitude versus frequency characte of 5 Rife-Vincent (I) windows, 5 Rife-Vincent (I) window has good sidelobe performance, side lobe attenuation speed reaches 30 decibels/octave, sidelobe level has reached-74.5 decibels, adopts 5 Rife-Vincent (I) window that input signal is carried out windowing process.

(3) based on the FFT computing of 5 Rife-Vincent (I) window

According to the requirement of accuracy of detection, determine the data length N of FFT computing.For ease of Implementation of Embedded System, the present invention considers only to contain the situation that single harmonic component exists, and the input signal amount x (n) that removes high fdrequency component is carried out windowing FFT computing, and the expression formula that can obtain the discrete Fourier transform (DFT) of signal after the windowing is

Wherein, W () is the frequency spectrum function of w (n).

(4) double peak interpolation algorithm

Under the non-synchronous sampling situation, discrete crest frequency f ₀=k ₀Δ f is not on the discrete spectral line frequency, establishes calculative each frequency peak point k ₀Near the maximum and time maximum spectral line difference k of amplitude ₁And k ₂Obviously, k ₁≤ k ₀≤ k ₂(k ₂=k ₁+ 1), the amplitude of these two spectral lines is respectively y ₁=| X (k ₁Δ f) | and y ₂=| X (k ₂Δ f) |, establish $\mathrm{\&beta;}=\frac{y_2-y_1}{{\mathrm{<figure-callout\; id="2"\; label="y"\; filenames="A20081003106600131.png"\; state="\{\{state\}\}">y</figure-callout>}}_2+{\mathrm{<figure-callout\; id="1"\; label="y"\; filenames="A20081003106600131.png"\; state="\{\{state\}\}">y</figure-callout>}}_1},$ α=k ₀-k ₁-0.5, the α span is [0.5,0.5] as can be known, and can be able to α is independent variable, is the funtcional relationship β=h (α) of dependent variable with β

$\mathrm{\&beta;}=\frac{\left|W\left(\frac{2\mathrm{\&pi;}(-\mathrm{\&alpha;}+0.5)}{N}\right)\right|-\left|W\left(\frac{2\mathrm{\&pi;}(-\mathrm{\&alpha;}-0.5)}{N}\right)\right|}{\left|W\left(\frac{2\mathrm{\&pi;}(-\mathrm{\&alpha;}+0.5)}{N}\right)\right|+\left|W\left(\frac{2\mathrm{\&pi;}(-\mathrm{\&alpha;}-0.5)}{N}\right)\right|}$

Its inverse function is α=h ^-1(β).Can obtain parameter alpha by β, then the frequency computation part formula of measured signal is

f ₀＝k ₀Δf＝(α+k ₁+0.5)Δf

The amplitude computing formula of measured signal is

$A=\frac{A_1\left|W\left(\frac{2\mathrm{\&pi;}{(k}_1-k_0)}{N}\right)\right|{+A}_2\left|W\left(\frac{2\mathrm{\&pi;}(k_2-k_0)}{N}\right)\right|}{\left|W\left(\frac{2\mathrm{\&pi;}(k_1-k_0)}{N}\right)\right|+\left|W\left(\frac{2\mathrm{\&pi;}(k_2-k_0)}{N}\right)\right|}$

The initial phase computing formula of measured signal

Can calculate the actual frequency and the parameters of measured signal first-harmonic and harmonic wave thus.

Analog to digital conversion of the present invention is meant by high speed analog-to-digital conversion is carried out in measured signal (can be voltage or current signal etc.), measured signal is converted to digital quantity.

Low-pass filtering of the present invention is meant the filtering algorithm that measured signal medium-high frequency signal can be carried out filtering.According to analog to digital conversion speed and frequency analysis requirement, determine that measured signal carries out the high reps of frequency analysis and each parameter of lowpass digital filter, low-pass filtering treatment is carried out in measured signal, can eliminate the influence of higher hamonic wave to the frequency analysis result.

Signal data length N of the present invention is the data length that carries out a FFT conversion, and N takes all factors into consideration definite according to the travelling speed of harmonic wave accuracy of detection and computing machine or embedded system device.

5 Rife-Vincent of the present invention (I) window is a class cosine combination window function, under the non-synchronous sampling situation, the FFT computing of not windowing will exist serious spectral leakage and fence effect, accuracy of detection is low, the spectral leakage and the fence effect that add 5 Rife-Vincent (I) window fft algorithm can significantly reduce, and can improve first-harmonic and harmonic wave accuracy of detection.

Double spectral line interpolation FFT algorithm of the present invention is meant under the non-synchronous sampling situation, there are deviation in discrete spectral line and actual frequency point, by seeking near two peak value spectral lines of actual frequency point, adopt the method for fitting of a polynomial to ask for frequency values, amplitude and the initial phase at true spectral line place.

Compared with the prior art the present invention has following advantage: 1, the present invention can carry out fast the first-harmonic and the harmonic component of measured signal, accurately detect; 2, the present invention can overcome the influence of measured signal frequency jitter, and the realization frequency is accurately measured; 3, Implementation of Embedded System is convenient in the present invention, can be continuously, for a long time measured signal is detected.

Description of drawings

Fig. 1 is first-harmonic of the present invention and harmonic measure theory diagram

Fig. 2 is 5 Rife-Vincent (I) window discrete time-domain performance plot

Fig. 3 is 5 Rife-Vincent (I) window amplitude versus frequency characte figure

Fig. 4 is a FFT double peak interpolation algorithm synoptic diagram

Fig. 5 is first-harmonic of the present invention and harmonic measure FB(flow block)

Embodiment

The present invention proposes a kind of first-harmonic and harmonic detecting method based on 5 Rife-Vincent (I) window.Be described further below in conjunction with accompanying drawing, but not as qualification of the present invention.

Present embodiment is the harmonic parameters test section that can be used for realizing the three-phase multifunctional harmonic electric energy table of Harmonious Waves in Power Systems electric energy metrical.

Setting in the present embodiment harmonic wave, to detect high reps be 21 times, the 71M6513H that analog to digital converter adopts TDK company to produce, and the BF533 processor that DSP adopts ADI company to produce, major parameter is selected as follows:

(1) sampling rate: f _s=2520Hz;

(2) electrical network fundamental frequency: about 50Hz, fluctuate;

(3) FFT operational data length: N=512.

The process flow block diagram of present embodiment as shown in Figure 1, sending into DSP (BF533) after the signal process A/D conversion (analog to digital conversion) handles, BF533 finishes and makes up discrete 5 Rife-Vincent (I) window, carries out windowed interpolation FFT computing, obtains the each harmonic parameter.

The expression formula of discrete 5 Rife-Vincent (I) window is as follows in the present embodiment

$w\left(n\right)=\underset{m=0}{\overset{M-1}{\mathrm{\&Sigma;}}}{(-1)}^m{a}_m\cos (2\mathrm{\&pi;n}\&CenterDot;m/N)$

In the formula, M is the item number of 5 Rife-Vincent (I) window function; N is the data length of FFT computing, also be window function length (n=1,2 ... N-1); a _mSatisfy constraint condition

$\underset{m=0}{\overset{M-1}{\mathrm{\&Sigma;}}}{(-1)}^m{a}_m=0$

The coefficient of 5 Rife-Vincent (I) window is respectively a ₀=1, a ₁=1.6, a ₂=0.8, a ₃=0.22857, a ₄=0.02857, its discrete time-domain characteristic as shown in Figure 2.Fig. 3 has provided the amplitude versus frequency characte of 5 Rife-Vincent (I) windows, 5 Rife-Vincent (I) window has good sidelobe performance, side lobe attenuation speed reaches 30 decibels/octave, sidelobe level has reached-74.5 decibels, adopts 5 Rife-Vincent (I) window that input signal is carried out windowing process.

In the present embodiment synoptic diagram of double peak interpolation algorithm as shown in Figure 4, k ₁And k ₂Represent near two spectral lines that actual frequency point is respectively, with fundamental frequency f _s=50Hz is example, then k ₁And k ₂To be positioned at 40～60Hz frequency range, and find out two peak value spectral lines in this scope, then actual frequency is o'clock between two spectral lines.The amplitude of two peak value spectral lines is respectively y ₁=| X (k ₁Δ f) | and y ₂=| X (k ₂Δ f) |.

If $\mathrm{\&beta;}=\frac{y_2-y_1}{{\mathrm{<figure-callout\; id="2"\; label="y"\; filenames="A20081003106600131.png"\; state="\{\{state\}\}">y</figure-callout>}}_2+{\mathrm{<figure-callout\; id="1"\; label="y"\; filenames="A20081003106600131.png"\; state="\{\{state\}\}">y</figure-callout>}}_1},$ α=k ₀-k ₁-0.5, the α span is [0.5,0.5] as can be known.

The spectrum amplitude function of 5 Rife-Vincent (I) window is

$W\left(w\right)=\underset{m=0}{\overset{M-1}{\mathrm{\&Sigma;}}}{(-1)}^m\frac{a_m}{2}[W_R(w-\frac{2\mathrm{\&pi;}}{N}m)+W_R(w+\frac{2\mathrm{\&pi;}}{N}m)]$

In the formula, $W_R\left(w\right)=\frac{\mathrm{<figure-callout\; id="2"\; label="sin\; Nw"\; filenames="A20081003106600131.png"\; state="\{\{state\}\}">sin</figure-callout>}\frac{\mathrm{Nw}}{}\frac{}{2}}{\sin \frac{w}{2}}{e}^{-j\frac{N-1}{2}w}.$

Order $w=\frac{2\mathrm{\&pi;}}{N}k,$ Can get

$W\left(\frac{2\mathrm{\&pi;}}{N}k\right)=\sin \mathrm{\&pi;k}\&CenterDot;e^{-\mathrm{j\&pi;k}}\left[\underset{m=0}{\overset{M-1}{\mathrm{\&Sigma;}}}{(-1)}^m\frac{a_{\mathrm{<figure-callout\; id="2"\; label="m"\; filenames="A20081003106600131.png"\; state="\{\{state\}\}">m</figure-callout>}}}{2}\frac{\mathrm{<figure-callout\; id="2"\; label="sin"\; filenames="A20081003106600131.png"\; state="\{\{state\}\}">sin</figure-callout>}\frac{2\mathrm{\&pi;k}}{N}}{\sin \frac{\mathrm{\&pi;}(k-m)}{N}\sin \frac{\mathrm{\&pi;}(k+m)}{N}}\right]$

With k=-α ± 0.5 substitution following formula, because of |-α ± 0.5|≤1, and N is generally bigger, can get

$\left|W\left(\frac{2\mathrm{\&pi;}(-\mathrm{\&alpha;}\&PlusMinus;0.5)}{N}\right)\right|\&ap;\left|\sin \mathrm{\&pi;}(-\mathrm{\&alpha;}\&PlusMinus;0.5)\right[\underset{m=0}{\overset{M-1}{\mathrm{\&Sigma;}}}{(-1)}^m\frac{a_m}{\mathrm{\&pi;}}\frac{N(-\mathrm{\&alpha;}\&PlusMinus;0.5)}{{(-\mathrm{\&alpha;}\&PlusMinus;0.5)}^2-m^2}]|$

Carry out fitting of a polynomial and approach, obtain based on the computing formula of 5 Rife-Vincent (I) window be

α＝4.49999982β+0.00001618β ³-0.0001433β ⁵

The frequency computation part formula is

f ₀＝k ₀Δf＝(α+k ₁+0.5)Δf

The amplitude computing formula is

$A=\frac{({\mathrm{<figure-callout\; id="1"\; label="y"\; filenames="A20081003106600131.png"\; state="\{\{state\}\}">y</figure-callout>}}_1+y_2)}{N}(1.05700918+0.2106604{\mathrm{\&alpha;}}^2+0.0229521{\mathrm{\&alpha;}}^4)$

The phase calculation formula is

In the formula, i=1,2.

The first-harmonic of present embodiment and harmonic parameters detect error result shown in table 1, table 2.Wherein, aE-b represents a * 10 ^-b

Table 1 first-harmonic～21 subharmonic amplitude relative errors (%)

Overtone order	1	2	3	4	5	6	7
								Error	2.7E-5	2.8E-5	5.2E-5	4.4E-6	5.1E-5	-5.7E-5	-8.6E-6
Overtone order	8	9	10	11	12	13	14
								Error	4.8E-5	-7.4E-5	-1.9E-5	4.3E-5	-7.3E-5	-4.0E-5	-9.8E-7
Overtone order	15	16	17	18	19	20	21
								Error	-4.0E-5	-7.2E-7	4.4E-6	-1.1E-4	-5.1E-5	-1.3E-4	1.2E-5

Table 2 first-harmonic～21 subharmonic initial phase relative errors (%)

Overtone order	1	2	3	4	5	6	7
								Error	7.1E-3	-2.1E-3	5.2E-5	7.4E-5	8.4E-5	-4.7E-4	-2.6E-5
Overtone order	8	9	10	11	12	13	14
								Error	1.4E-3	-3.8E-4	8.9E-4	4.7E-3	-1.7E-2	1.2E-4	6.8E-3
Overtone order
		15	16	17	18	19	20	21
Error									-4.6E-5	-5.0E-3	-2.9E-6	-4.5E-3	5.0E-3	2.9E-2	4.0E-5

Present embodiment has also adopted other classical Rife-Vincent window function to test, the result shows, employing is during based on the first-harmonic of 5 Rife-Vincent (I) window and harmonic detecting method, humorous wave amplitude, the initial phase error of calculation will exceed 2～4 orders of magnitude of other classical Rife-Vincent window, and computational accuracy obviously improves.

Claims (5)

1. a first-harmonic and harmonic parameters detection method that is used for signal Processing based on 5 Rife-Vincent (I) window double spectral line interpolation FFT, it is characterized in that: analog to digital conversion and low-pass filtering treatment are carried out in measured signal, adopt 5 Rife-Vincent (I) window that sidelobe level is low, the rate of decay is fast that signal is carried out windowing FFT processing, utilization double spectral line interpolation algorithm detects first-harmonic and each harmonic content in the measured signal.

2. analog to digital conversion according to claim 1 is meant by high speed analog-to-digital conversion is carried out in measured signal (can be voltage or current signal etc.), measured signal is converted to digital quantity.

3. low-pass filtering according to claim 1 is meant the filtering algorithm that measured signal medium-high frequency signal can be carried out filtering.According to analog to digital conversion speed and frequency analysis requirement, determine that measured signal carries out the high reps of frequency analysis and each parameter of lowpass digital filter, low-pass filtering treatment is carried out in measured signal, can eliminate the influence of higher hamonic wave to the frequency analysis result.

4. 5 Rife-Vincent according to claim 1 (I) window is a class cosine combination window function, under the non-synchronous sampling situation, the FFT computing of not windowing will exist serious spectral leakage and fence effect, accuracy of detection is low, the spectral leakage and the fence effect that add 5 Rife-Vincent (I) window fft algorithm can significantly reduce, can effectively improve first-harmonic and harmonic wave accuracy of detection, and calculated amount is less, is beneficial to realization.

5. the double spectral line interpolation FFT algorithm of stating according to claim 1 is meant under the non-synchronous sampling situation, there are deviation in discrete spectral line and actual frequency spectral line, by seeking near two peak value spectral lines of actual frequency point, adopt the method for fitting of a polynomial to ask for frequency values, amplitude and the initial phase at true spectral line place.

Images