CN106771586A - The loop signal analysis method and device of a kind of DC control protection board - Google Patents
The loop signal analysis method and device of a kind of DC control protection board Download PDFInfo
- Publication number
- CN106771586A CN106771586A CN201611058734.0A CN201611058734A CN106771586A CN 106771586 A CN106771586 A CN 106771586A CN 201611058734 A CN201611058734 A CN 201611058734A CN 106771586 A CN106771586 A CN 106771586A
- Authority
- CN
- China
- Prior art keywords
- amplitude
- frequency
- signal
- window
- harmonic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
Abstract
The invention discloses the loop signal analysis method and device of a kind of DC control protection board, be related to signal processing technology field, for solves the problems, such as cannot effectively during inhibition analysis loop signal generation spectrum leakage and invent.Main method of the invention includes：According to preset sample frequency, uniform sampling board loop signal obtains discrete signal；According to mixing convolution window interception discrete signal, filtering signal is obtained；Filtering signal is carried out the result of Fast Fourier Transform (FFT) FFT, is defined as harmonic spectrum；In the harmonic spectrum apart from preset basic frequency less than predetermined frequency deviation, the peak value frequency of and amplitude maximum nearest apart from preset basic frequency is searched；According to preset basic frequency and the corresponding amplitude of peak value frequency, amplitude departure is calculated；According to departure, each harmonic parameter in harmonic spectrum is calculated, harmonic parameters include amplitude, phase and frequency.Present invention is mainly applied in the loop signal for analyzing DC control protection board.
Description
Technical field
The present invention relates to signal processing technology field, more particularly to a kind of loop signal analysis of DC control protection board
Method and device.
Background technology
Extra-high voltage direct-current transmission refers to the direct current transportation of ± 800kV and above.Extra-high voltage direct-current transmission is with defeated
Send that capacity is big, transmission distance remote, voltage is high, the characteristics of can be used for power system asynchronous networking, range of application is more and more wider.For
Meet the growth requirement of extra-high voltage direct-current transmission, launch the research work to the safe and stable operation of DC control and protection system
It is particularly important.DC control board is control DC output circuit plate, if DC control board faults cause direct current locking,
The safe operation of power network is had a strong impact on.So for DC control panel card sets control protection board, to ensure DC control panel
Card being capable of safe operation.
Board, the analysis method of control protection board loop signal are protected in control currently for DC control protecting platform
It is FFT (Fast Fourier Transformation, Fast Fourier Transform (FFT)) algorithm, but when loop signal frequency occurs
When changing and causing non-synchronous sampling, can there is fence effect and leakage phenomenon in fft algorithm so that the loop signal being calculated
Frequency, amplitude and phase it is all inaccurate, especially phase error is larger, it is impossible to accurate measurement loop signal, to direct current transportation
System causes potential safety hazard.
In the prior art, to weaken the shadow of fence effect and spectrum leakage to accuracy of measurement in the analysis of board loop signal
Ring, reduce because of the error produced by non-synchronous sampling, result is modified, is optimized using Interpolating Window FFT Algorithm, wherein
Conventional window function mainly has rectangular window, Hanning window, Blackman window etc., and correction algorithm is frequently with double spectral line interpolation method.Although
These methods improve computational accuracy to a certain extent, but it is bad to there is also window function performance, double spectral line interpolation computational accuracy
Inadequate problem, it is impossible to effectively produce spectrum leakage during inhibition analysis loop signal.
The content of the invention
Goal of the invention of the invention is the loop signal analysis method and device for providing a kind of DC control protection board,
Can solve the problem that the problem that spectrum leakage cannot be effectively produced during inhibition analysis loop signal.
Embodiments in accordance with the present invention, there is provided a kind of DC control protects the loop signal analysis method of board, including：
According to preset sample frequency, uniform sampling board loop signal obtains discrete signal；
The discrete signal is intercepted according to mixing convolution window, filtering signal is obtained, the mixing convolution window is by cosine group
Closing window and Hanning window carries out the window function that convolution algorithm is obtained；
The filtering signal is carried out the result of Fast Fourier Transform (FFT) FFT, is defined as harmonic spectrum；
In the harmonic spectrum apart from preset basic frequency less than predetermined frequency deviation, search apart from the preset base
Plinth frequency recently and amplitude maximum peak value frequency；
According to the preset basic frequency and the corresponding amplitude of the peak value frequency, amplitude departure is calculated；
According to the amplitude departure, each harmonic parameter in the harmonic spectrum is calculated, the harmonic parameters include width
Value, phase and frequency.
Embodiments in accordance with the present invention, additionally provide the loop signal analytical equipment that a kind of DC control protects board, bag
Include：
Acquiring unit, for according to preset sample frequency, uniform sampling board loop signal to obtain discrete signal；
The acquiring unit, is additionally operable to intercept the discrete signal according to mixing convolution window, obtains filtering signal, described mixed
It is to carry out the window function that convolution algorithm is obtained by cosine combination window and Hanning window to close convolution window；
First determining unit, the result for the filtering signal to be carried out Fast Fourier Transform (FFT) FFT, is defined as harmonic wave
Frequency spectrum；
Searching unit, in the harmonic spectrum less than predetermined frequency deviation apart from preset basic frequency, looking into
Look for the peak value frequency of and amplitude maximum nearest apart from the preset basic frequency；
Computing unit, for according to the preset basic frequency and the corresponding amplitude of the peak value frequency, calculating amplitude inclined
Residual quantity；
The computing unit, is additionally operable to, according to the amplitude departure, calculate each harmonic parameter in the harmonic spectrum,
The harmonic parameters include amplitude, phase and frequency.
From above technical scheme, a kind of loop signal analysis side of direct current protecting board provided in an embodiment of the present invention
Method and device, by the way that according to preset sample frequency, uniform sampling board loop signal obtains discrete signal, then according to mixing
Convolution window intercepts discrete signal, acquisition filtering signal, then filtering signal is carried out the result of Fast Fourier Transform (FFT) FFT, it is determined that
It is harmonic frequency；Again in the harmonic spectrum apart from preset basic frequency less than predetermined frequency deviation, search apart from preset basis
Frequency recently and amplitude maximum peak value frequency, further according to preset basic frequency and the corresponding amplitude of peak value frequency, calculate amplitude
Departure, finally according to amplitude departure, calculates each harmonic parameter in harmonic spectrum.Compared with prior art, energy of the present invention
It is enough to increase the time domain scale of interception discrete signal using the weighting of convolution window function pair discrete signal is mixed, reduce discrete signal by
Frequency domain produces the energy of secondary lobe leakage, efficiently reduces the spectrum leakage produced during FFT, it is to avoid in the frequency spectrum of discrete signal
Mutual leakage interference between fundamental wave and each harmonic.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to institute in embodiment
The accompanying drawing for needing to use is briefly described, it should be apparent that, drawings in the following description are only some implementations of the invention
Example, for those of ordinary skill in the art, on the premise of not paying creative work, can also obtain according to these accompanying drawings
Obtain other accompanying drawings.
Fig. 1 is to be preferable to carry out the loop signal analysis method flow that a kind of DC control for exemplifying protects board according to one
Figure；
Fig. 2 is to be preferable to carry out the loop signal analysis method stream that another DC control that exemplifies protects board according to one
Cheng Tu；
Fig. 3 is to be preferable to carry out the amplitude-frequency response schematic diagram of the first convolution window that exemplifies according to one；
Fig. 4 is to be preferable to carry out the amplitude-frequency response schematic diagram of second convolution window that exemplifies according to one；
Fig. 5 is to be preferable to carry out the amplitude-frequency response schematic diagram of the third convolution window that exemplifies according to one；
Fig. 6 is to protect the loop signal analytical equipment of board to constitute according to an a kind of DC control for being preferable to carry out exemplifying
Block diagram；
Fig. 7 is to be preferable to carry out the loop signal analytical equipment group that another DC control that exemplifies protects board according to one
Into block diagram.
Specific embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete
Whole description, it is clear that described embodiment is only a part of embodiment of the invention, rather than whole embodiments.It is based on
Embodiment in the present invention, it is every other that those of ordinary skill in the art are obtained under the premise of creative work is not made
Embodiment, belongs to the scope of protection of the invention.
Because extra-high voltage direct-current transmission is the feeder ear of whole system, so the output stability pair of extra-high voltage direct-current transmission
There is tremendous influence in whole electric power system.So for DC control protection platform sets control protection board, control protection board energy
Enough for DC control protection platform provides a certain degree of protection, but control protection board is likely to break down.Needs pass through
Analysis protection board loop signal, judges whether control protection board breaks down, so in order to improve extra-high voltage direct-current transmission
The stability and reliability of system, the embodiment of the invention provides the loop signal analysis side that a kind of DC control protects board
Method, as shown in figure 1, including：
S101, according to preset sample frequency, uniform sampling board loop signal obtains discrete signal.
Board loop signal, is the loop signal of control protection board, can be current signal, or voltage letter
Number.No matter board loop signal is current signal or voltage signal, what the analysis method of use was just as.
Sampling refers to the process of that the continuous quantity of time-domain or spatial domain is changing into discrete magnitude.Namely acquisition time and
All continuous analog signal in amplitude, is converted on the time continuous discrete analog(ue) on discrete amplitude in the presence of sample frequency
Signal, as the sample of analog signal.Wherein sample frequency refers to the sample number of each second.Sample frequency is higher, discrete
Analog signal is closer to original analog signal.Uniform sampling is namely according to identical sample frequency equal interval sampling, uniform sampling
The frequency spectrum of sampled signal be that analog signal frequency spectrum carries out cyclic extension formation by the cycle of sample frequency.
Board loop signal is analog signal, and according to preset sample frequency, uniform sampling board loop signal obtains frequency spectrum
It is capable of the discrete signal of cyclic extension.
S102, according to mixing convolution window interception discrete signal, obtain filtering signal.
Mixing convolution window is to carry out the window function that convolution algorithm is obtained by cosine combination window and Hanning window.Hanning window is to rise remaining
One special case of porthole, Hanning window can be regarded as 3 frequency spectrum sums of rectangular time window, or perhaps 3 sinc (t) type letters
Number sum, and two in bracket respectively move π/T to left and right relative to first spectrum window, so that secondary lobe cancels each other, disappear
Remove High-frequency Interference and leakage energy.Cosine combination window can be Hamming window or Blackman window.By convolved mixtures window and discrete signal
The result of product is used as filtering signal.Filtering signal is the part signal that convolved mixtures window signal is intercepted from discrete signal.
S103, the result that filtering signal is carried out Fast Fourier Transform (FFT) FFT, are defined as harmonic spectrum.
Fft algorithm, is the fast algorithm of discrete Fourier transform, is according to the odd, even, empty, real of discrete Fourier transform
Etc. characteristic, the algorithm to DFT is improved acquisition.Fourier's analysis method is the most basic of signal analysis
Method, Fourier transformation is the core of Fourier analysis, signal is transformed to frequency domain from time-domain by it, and then study letter
Number spectrum structure and Changing Pattern.Harmonic spectrum embodies the spectrum structure and Changing Pattern of filtering signal.
S104, apart from preset basic frequency less than predetermined frequency deviation harmonic spectrum in, search apart from preset basis
Frequency recently and amplitude maximum peak value frequency.
Preset basic frequency, is theoretical frequency, it can be understood as the Harmonic Theory frequency to be extracted, such as extract 1 time it is humorous
The parameter of ripple, can be input into the theoretic frequency 50Hz of 1 subharmonic；Extract the parameter of 2 subharmonic, it is possible to be input into 2 subharmonic
Theoretic frequency 100Hz.Preset basic frequency is the corresponding multiple frequencies with multiple harmonic.According to spectrum signature, each harmonic
Maximum amplitude be different, so set predetermined frequency deviation, apart from preset basic frequency less than predetermined frequency deviation
Peak value frequency is searched in harmonic spectrum.
Search peak value frequency when, the distance with preset basic frequency be predetermined frequency deviation in the range of, find away from
From preset basic frequency recently, and amplitude maximum frequency location point.
S105, according to preset basic frequency and the corresponding amplitude of peak value frequency, calculate amplitude departure.
Amplitude departure, is exactly in harmonic spectrum, to calculate the corresponding amplitude of preset basic frequency corresponding with peak value frequency
Amplitude difference.Amplitude departure includes the departure of multiple harmonic, can form amplitude departure function.
S106, according to amplitude departure, calculate each harmonic parameter in harmonic spectrum.
Harmonic parameters include amplitude, phase and frequency.According to amplitude departure, each harmonic parameter in harmonic spectrum is counted,
Eliminate fence effect so that harmonic parameters are more accurate.
From above technical scheme, a kind of loop signal analysis side of direct current protecting board provided in an embodiment of the present invention
Method, by the way that according to preset sample frequency, uniform sampling board loop signal obtains discrete signal, then according to mixing convolution window
Interception discrete signal, acquisition filtering signal, then filtering signal is carried out the result of Fast Fourier Transform (FFT) FFT, it is defined as harmonic wave
Frequency；Again in the harmonic spectrum apart from preset basic frequency less than predetermined frequency deviation, search apart from preset basic frequency most
The peak value frequency of near and amplitude maximum, further according to preset basic frequency and the corresponding amplitude of peak value frequency, calculates amplitude departure,
Finally according to amplitude departure, each harmonic parameter in harmonic spectrum is calculated.Compared with prior art, the present invention can be using mixed
The weighting of convolution window function pair discrete signal is closed, increases the time domain scale of interception discrete signal, reduced discrete signal and produced by frequency domain
The energy of secondary lobe leakage, efficiently reduces the spectrum leakage produced during FFT, it is to avoid fundamental wave and each in the frequency spectrum of discrete signal
Mutual leakage interference between subharmonic.
The loop signal analysis method that another DC control protects board is the embodiment of the invention provides, as shown in Fig. 2
Including：
S201, according to preset sample frequency, uniform sampling board loop signal obtains discrete signal.
S101 is identical the step of method with shown in Fig. 1 for this step, repeats no more here.
S202, by the Hanning window of the first presetting length and the convolution results of the Hamming window of the second presetting length, be defined as mixing
Close convolution window.
When carrying out signal measurement and calculating, it is impossible to which the signal of endless is measured and calculated, but application window letter
Number, extracts limited time slice and is analyzed.Specially：A time slice is intercepted from signal, periodic extension is then used
Method processed, obtain virtual unlimited signal long.The time slice length for wherein being intercepted as institute's windowed function
Length.Window function with certain length, could intercept a signal for time slice from the signal of endless.
For mixing convolution window, the Hanning window that can be used or Blackman window are obtained with Hamming window convolution.Exemplary,
Cosine to being obtained after Hamming window convolution that Hanning window that discrete signal plus length are 64 and length are 64 mixes volume
Product window, Fig. 3 is the amplitude-frequency response figure of the convolution window, and its main lobe width is that side lobe attenuation speed is 24dB/oc, side lobe peak electricity
Put down as -75dB.To discrete signal plus length it is 64 Hamming window and after a length carries out convolution for 64 Blackman windows
The cosine mixing convolution window for obtaining, Fig. 4 is the amplitude-frequency response figure of the convolution window, and its main lobe width is that side lobe attenuation speed is
20dB/oc, side lobe peak level is -71dB.To discrete signal plus length it is 64 Hanning window and a length is 128
The cosine mixing convolution window that Hamming window obtained after convolution, Fig. 5 is the amplitude-frequency response figure of the convolution window, and its main lobe width is,
Side lobe attenuation speed is 24dB/oc, and side lobe peak level is -42dB.
In order to suppress the influence of spectrum leakage and fence effect to greatest extent, anti-interference is improved, selection has faster
Side lobe attenuation speed and more low sidelobe peak level, the first presetting length is equal with the second presetting length, and a length is
64 Hanning window and length be 64 Hamming window convolution after the mixing convolution window that obtains.
S203, the discrete signal is intercepted according to mixing convolution window, obtain filtering signal.
Mixing convolution window is to carry out the window function that convolution algorithm is obtained by cosine combination window and Hanning window.Wherein cosine combination
Window can be Hanning window, Hamming window or Blackman window.The determination method of mixing convolution window is detailed in step S202.From
According to mixing convolution window interception filtering signal in discrete signal.Filtering signal is some discrete signal.
S204, the result that filtering signal is carried out Fast Fourier Transform (FFT) FFT, are defined as harmonic spectrum.
S205, apart from preset basic frequency less than predetermined frequency deviation harmonic spectrum in, search apart from preset basis
Frequency recently and amplitude maximum peak value frequency.
In relatively amplitude size, can first judge centered on preset basic frequency, distance positively and negatively is
In the range of predetermined frequency deviation, harmonic spectrum is increasing function, subtraction function or having increasing to have subtracts, and is then taken according to function characteristic
Endpoint value calculate amplitude maximum peak value frequency, or according to the slope of the harmonic frequency calculate amplitude maximum peak value frequently
Point.
S206, according to preset basic frequency and the corresponding amplitude of peak value frequency, calculate amplitude departure.
Obtain the corresponding first amplitude of peak value frequency；Obtain the second amplitude apart from one unit length of peak value frequency and the
Three amplitudes；3rd amplitude and first amplitude are made the difference, again divided by the inverse function of the second amplitude, be defined as amplitude departure.According to
Polynomial approximation method, solves amplitude residual quantity.
Due to preset basic frequency peak value frequency corresponding with its, all exist in each harmonic, so preset basis frequency
Point amplitude corresponding with peak value frequency is also all exist in each harmonic, so amplitude departure is made up of multinomial.
Using polynomial approximation method, compared to general curve matching least square, amount of calculation is small, becomes to be easier.
S207, according to amplitude departure, calculate each harmonic parameter in harmonic spectrum.
Amplitude departure, calculating be the corresponding amplitude of different frequent points departure, wherein with frequency as independent variable, with width
It is dependent variable to be worth.According to amplitude departure, corresponding calculating harmonic parameters include amplitude, phase and frequency.
For the embodiment of the present invention, specific application scenarios can be with as follows, but not limited to this, including：Assuming that protection
Board loop signal is：
With fixed sampling frequency fs, carry out uniform sampling and obtain discrete signal：
Wherein n=1,23...N-1, n represent that board loop signal has nth harmonic, AiAmplitude, f for i ＆ lt harmonic wavei
For the frequency of i ＆ lt harmonic wave, N are sampling number.
The Hamming window that one length of selection is 64 Hanning window and length is 64, carries out convolution and obtains mixing convolution
Window, its time-domain expression is：
Wherein, Hn64 is the Hanning window that length is 64, and Hm64 is the Hamming window that length is 64, M be in time-domain signal most
Big time, a when m is even numberm=0.27, a when m is odd numberm=0.23, n=1,2,3...127.
Discrete signal is intercepted with mixing convolution window, filtering signal is obtained：
xw(n)=x (n) wHn64-Hm64(n)
Then filtering signal is entered into line translation, obtains harmonic spectrum：
WhereinK=1,2,3 ... 127.If ignore the corresponding amplitude of peak value frequency influences on harmonic spectrum secondary lobe,
Harmonic spectrum is：
WhereinK=1,2,3 ... 127.
In FFT calculating process, when non-synchronous sampling is carried out to signal or non-integer is blocked, due to fence effect,
The frequency h of resulting i ＆ lt harmonic signaliΔ f is difficult to be located exactly on sampling frequency point, i.e. hi=fi/ Δ f is general not
It is integer.Assuming that near the amplitude maximum frequency spectrum for obtaining of being sampled near the accurate frequency of i ＆ lt harmonic wave amplitude maximum spectral line
haAnd its left and right spectral line ha+1And ha-1, make amplitude departure δ=hi-ha, then have -0.5 ＜ δ ＜ 0.5.Accurately to try to achieve side-play amount,
So as to accurately estimate each harmonic parameter, the present invention is corrected using three spectral line interpolation methods to frequency spectrum.Due to non-synchronous sampling
When, ha、ha-1And ha+1Higher magnitude is respectively provided with, shows that they contain information largely related to actual spectrum, thus this three
The information of spectral line can be used to interpolation amendment.Remember that this three spectral line amplitudes are respectively：y1=| X (ha-1) |, y2=| X (ha) |, y1
=| X (ha-1)|.According to three amplitudes of spectral line, three relations of spectral line are set up：
Due to amplitude departure and the relation of spectral line：
If sampling number is larger, the relation of amplitude departure and frequency spectrum can be reduced to α=f (δ), its inverse function δ=
f-1(α), magnitude shift amount δ is solved by solving inverse function.As mixing convolution window function wHn64-Hm64When () is real coefficient n, its width
Frequency response answers W (2 π f) for an even function, therefore function f-1(α) and f (δ) are odd function.δ is solved using Polynomial Approximation Method
=f-1(α), obtains：
δ=f-1(α)≈c1·α+c3·α3+...+c2q+1·α2q+1
Wherein, c1, c3..., c2q+1It is the 2q+1 odd term coefficient of approximating polynomial.By amendment, cosine mixing convolution
Regression formulas corresponding to window are：
After polynomial approximation method tries to achieve side-play amount, i ＆ lt harmonic frequency signal and phase can be tried to achieve：
More accurately to obtain the amplitude of harmonic signal, by y1, y2, y3Bring into for amplitude correction, composed by this three
Wire spoke value is weighted the amplitude that average computation obtains true peak point.Consider y2To being maximum amplitude spectral line, thus when calculating
Give this spectral line weight limit.Amplitude size is：
When sampling number is larger, the big I of amplitude is reduced to：
Ai=N-1(y1+2y2+y3) v (δ),
Using polynomial approximation method to being fitted, obtain：
V (δ)=1.580952135+0.48529241 δ2+0.08251422δ4+0.01263077δ6
And then the change of amplitude size can be obtained.
From above technical scheme, a kind of loop signal analysis side of direct current protecting board provided in an embodiment of the present invention
Method, by the way that according to preset sample frequency, uniform sampling board loop signal obtains discrete signal, then according to mixing convolution window
Interception discrete signal, acquisition filtering signal, then filtering signal is carried out the result of Fast Fourier Transform (FFT) FFT, it is defined as harmonic wave
Frequency；Again in the harmonic spectrum apart from preset basic frequency less than predetermined frequency deviation, search apart from preset basic frequency most
The peak value frequency of near and amplitude maximum, further according to preset basic frequency and the corresponding amplitude of peak value frequency, calculates amplitude departure,
Finally according to amplitude departure, each harmonic parameter in harmonic spectrum is calculated.Compared with prior art, the present invention can be using mixed
The weighting of convolution window function pair discrete signal is closed, increases the time domain scale of interception discrete signal, reduced discrete signal and produced by frequency domain
The energy of secondary lobe leakage, efficiently reduces the spectrum leakage produced during FFT, it is to avoid fundamental wave and each in the frequency spectrum of discrete signal
Mutual leakage interference between subharmonic.
Used as implementing for method shown in Fig. 1 and Fig. 2, the embodiment of the present invention additionally provides a kind of DC control protection board
The loop signal analytical equipment of card, as shown in fig. 6, including：
Acquiring unit 61, for according to preset sample frequency, uniform sampling board loop signal to obtain discrete signal；
Acquiring unit 61, is additionally operable to, according to mixing convolution window interception discrete signal, obtain filtering signal, and mixing convolution window is
The window function that convolution algorithm is obtained is carried out by cosine combination window and Hanning window；
First determining unit 62, the result for filtering signal to be carried out Fast Fourier Transform (FFT) FFT, is defined as harmonic wave frequently
Spectrum；
Searching unit 63, in the harmonic spectrum less than predetermined frequency deviation apart from preset basic frequency, searching
The peak value frequency of nearest and amplitude maximum apart from preset basic frequency；
Computing unit 64, for according to preset basic frequency and the corresponding amplitude of peak value frequency, calculating amplitude departure；
Computing unit 64, is additionally operable to, according to amplitude departure, calculate each harmonic parameter, harmonic parameters bag in harmonic spectrum
Include amplitude, phase and frequency.
Further, as shown in fig. 7, the device also includes：
Second determining unit 65, for according to mixing convolution window interception discrete signal, by the before obtaining filtering signal
The convolution results of the Hamming window of the Hanning window of one presetting length and the second presetting length, are defined as mixing convolution window.
Further, the first presetting length is equal with the second presetting length.
Further, as shown in fig. 7, computing unit 64, including：
Acquisition module 641, for obtaining the corresponding first amplitude of peak value frequency；
Acquisition module 641, is additionally operable to obtain the second amplitude and the 3rd amplitude apart from one unit length of peak value frequency；
Determining module 642, for the 3rd amplitude and first amplitude to be made the difference, again divided by the inverse function of the second amplitude, it is determined that
It is amplitude departure.
Further, computing unit 64, are additionally operable to, according to polynomial approximation method, solve amplitude residual quantity.
From above technical scheme, a kind of loop signal analysis dress of direct current protecting board provided in an embodiment of the present invention
Put, by the way that according to preset sample frequency, uniform sampling board loop signal obtains discrete signal, then according to mixing convolution window
Interception discrete signal, acquisition filtering signal, then filtering signal is carried out the result of Fast Fourier Transform (FFT) FFT, it is defined as harmonic wave
Frequency；Again in the harmonic spectrum apart from preset basic frequency less than predetermined frequency deviation, search apart from preset basic frequency most
The peak value frequency of near and amplitude maximum, further according to preset basic frequency and the corresponding amplitude of peak value frequency, calculates amplitude departure,
Finally according to amplitude departure, each harmonic parameter in harmonic spectrum is calculated.Compared with prior art, the present invention can be using mixed
The weighting of convolution window function pair discrete signal is closed, increases the time domain scale of interception discrete signal, reduced discrete signal and produced by frequency domain
The energy of secondary lobe leakage, efficiently reduces the spectrum leakage produced during FFT, it is to avoid fundamental wave and each in the frequency spectrum of discrete signal
Mutual leakage interference between subharmonic.
Those skilled in the art considering specification and after putting into practice invention disclosed herein, will readily occur to it is of the invention its
Its embodiment.The application is intended to any modification of the invention, purposes or adaptations, these modifications, purposes or
Person's adaptations follow general principle of the invention and including undocumented common knowledge in the art of the invention
Or conventional techniques.Description and embodiments are considered only as exemplary, and true scope and spirit of the invention are by following
Claim is pointed out.
It should be appreciated that the invention is not limited in the precision architecture being described above and be shown in the drawings, and
And can without departing from the scope carry out various modifications and changes.The scope of the present invention is only limited by appended claim.
Claims (10)
1. a kind of DC control protects the loop signal analysis method of board, it is characterised in that methods described includes：
According to preset sample frequency, uniform sampling board loop signal obtains discrete signal；
The discrete signal is intercepted according to mixing convolution window, filtering signal is obtained, the mixing convolution window is by cosine combination window
The window function that convolution algorithm is obtained is carried out with Hanning window；
The filtering signal is carried out the result of Fast Fourier Transform (FFT) FFT, is defined as harmonic spectrum；
In the harmonic spectrum apart from preset basic frequency less than predetermined frequency deviation, search apart from the preset basis frequency
Put the peak value frequency of nearest and amplitude maximum；
According to the preset basic frequency and the corresponding amplitude of the peak value frequency, amplitude departure is calculated；
According to the amplitude departure, each harmonic parameter in the harmonic spectrum is calculated, the harmonic parameters include amplitude, phase
Position and frequency.
2. method according to claim 1, it is characterised in that the discrete letter is intercepted according to mixing convolution window described
Number, before obtaining filtering signal, methods described also includes：
By the Hanning window of the first presetting length and the convolution results of the Hamming window of the second presetting length, it is defined as the mixing convolution
Window.
3. method according to claim 2, it is characterised in that first presetting length and the second presetting length phase
Deng.
4. method according to claim 1, it is characterised in that it is described according to the preset basic frequency and the peak value frequently
The corresponding amplitude of point, calculates amplitude departure, including：
Obtain the corresponding first amplitude of the peak value frequency；
Obtain the second amplitude and the 3rd amplitude apart from one unit length of the peak value frequency；
3rd amplitude is made the difference with the first amplitude, again divided by the inverse function of second amplitude, be defined as the width
Value departure.
5. the method according to claim any one of 1-4, it is characterised in that the calculating amplitude departure, including：
According to polynomial approximation method, the amplitude residual quantity is solved.
6. a kind of DC control protects the loop signal analytical equipment of board, it is characterised in that described device includes：
Acquiring unit, for according to preset sample frequency, uniform sampling board loop signal to obtain discrete signal；
The acquiring unit, is additionally operable to intercept the discrete signal according to mixing convolution window, obtains filtering signal, the mixing volume
Product window is to carry out the window function that convolution algorithm is obtained by cosine combination window and Hanning window；
First determining unit, the result for the filtering signal to be carried out Fast Fourier Transform (FFT) FFT, is defined as harmonic wave frequently
Spectrum；
Searching unit, in the harmonic spectrum less than predetermined frequency deviation apart from preset basic frequency, search away from
The peak value frequency of nearest and amplitude maximum from the preset basic frequency；
Computing unit, for according to the preset basic frequency and the corresponding amplitude of the peak value frequency, calculating amplitude departure；
The computing unit, is additionally operable to according to the amplitude departure, calculates each harmonic parameter in the harmonic spectrum, described
Harmonic parameters include amplitude, phase and frequency.
7. device according to claim 6, it is characterised in that described device also includes：
Second determining unit, for intercepting the discrete signal according to mixing convolution window described, will before obtaining filtering signal
The convolution results of the Hamming window of the Hanning window of the first presetting length and the second presetting length, are defined as the mixing convolution window.
8. device according to claim 7, it is characterised in that first presetting length and the second presetting length phase
Deng.
9. device according to claim 6, it is characterised in that the computing unit, including：
Acquisition module, for obtaining the corresponding first amplitude of the peak value frequency；
The acquisition module, is additionally operable to obtain the second amplitude and the 3rd amplitude apart from one unit length of the peak value frequency；
Determining module, for the 3rd amplitude to be made the difference with the first amplitude, again divided by the inverse function of second amplitude,
It is defined as the amplitude departure.
10. the device according to claim any one of 6-9, it is characterised in that the computing unit, is additionally operable to
According to polynomial approximation method, the amplitude residual quantity is solved.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611058734.0A CN106771586B (en) | 2016-11-24 | 2016-11-24 | A kind of the loop signal analysis method and device of DC control protection board |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611058734.0A CN106771586B (en) | 2016-11-24 | 2016-11-24 | A kind of the loop signal analysis method and device of DC control protection board |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106771586A true CN106771586A (en) | 2017-05-31 |
| CN106771586B CN106771586B (en) | 2019-07-26 |
Family
ID=58912694
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611058734.0A Active CN106771586B (en) | 2016-11-24 | 2016-11-24 | A kind of the loop signal analysis method and device of DC control protection board |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106771586B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109490627A (en) * | 2018-12-25 | 2019-03-19 | 新疆大学 | Novel electric power system harmonics analysis method based on mixing convolution window |
| CN109598094A (en) * | 2018-12-29 | 2019-04-09 | 北京化工大学 | Earthquake vector wave field finite difference numerical simulation method, equipment and system |
| CN109900959A (en) * | 2019-04-17 | 2019-06-18 | 贵州电网有限责任公司 | The extracting method of harmonic components in a kind of dynamic sinusoidal distorted signal |
| CN110096676A (en) * | 2019-04-30 | 2019-08-06 | 哈尔滨工业大学 | A kind of metering-type optical current voltage combination transformer frequency analysis device and method |
| CN114355297A (en) * | 2021-12-06 | 2022-04-15 | 电子科技大学 | Frequency control array design method based on multi-window mixed weighted frequency offset |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN202339381U (en) * | 2011-10-19 | 2012-07-18 | 广西电网公司电力科学研究院 | Harmonic electric energy metering system based on Nuttall self-convolution window weighed FFT (Fast Fourier Transform) |
| CN103207319A (en) * | 2013-03-12 | 2013-07-17 | 广东电网公司电力科学研究院 | Harmonic wave measurement method of electricity signal of digital substation under non-synchronous sampling condition |
| CN104062528A (en) * | 2014-07-04 | 2014-09-24 | 武汉大学 | Signal harmonic analysis method and system based on Hanning product window |
| CN104062500A (en) * | 2014-07-04 | 2014-09-24 | 武汉大学 | Signal harmonic analysis method and system based on Hamming product window |
| CN105717359A (en) * | 2016-02-19 | 2016-06-29 | 云南电网有限责任公司电力科学研究院 | Harmonic analysis algorithm |
-
2016
- 2016-11-24 CN CN201611058734.0A patent/CN106771586B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN202339381U (en) * | 2011-10-19 | 2012-07-18 | 广西电网公司电力科学研究院 | Harmonic electric energy metering system based on Nuttall self-convolution window weighed FFT (Fast Fourier Transform) |
| CN103207319A (en) * | 2013-03-12 | 2013-07-17 | 广东电网公司电力科学研究院 | Harmonic wave measurement method of electricity signal of digital substation under non-synchronous sampling condition |
| CN104062528A (en) * | 2014-07-04 | 2014-09-24 | 武汉大学 | Signal harmonic analysis method and system based on Hanning product window |
| CN104062500A (en) * | 2014-07-04 | 2014-09-24 | 武汉大学 | Signal harmonic analysis method and system based on Hamming product window |
| CN105717359A (en) * | 2016-02-19 | 2016-06-29 | 云南电网有限责任公司电力科学研究院 | Harmonic analysis algorithm |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109490627A (en) * | 2018-12-25 | 2019-03-19 | 新疆大学 | Novel electric power system harmonics analysis method based on mixing convolution window |
| CN109598094A (en) * | 2018-12-29 | 2019-04-09 | 北京化工大学 | Earthquake vector wave field finite difference numerical simulation method, equipment and system |
| CN109598094B (en) * | 2018-12-29 | 2020-12-04 | 北京化工大学 | Seismic vector wave field finite difference numerical simulation method, device and system |
| CN109900959A (en) * | 2019-04-17 | 2019-06-18 | 贵州电网有限责任公司 | The extracting method of harmonic components in a kind of dynamic sinusoidal distorted signal |
| CN110096676A (en) * | 2019-04-30 | 2019-08-06 | 哈尔滨工业大学 | A kind of metering-type optical current voltage combination transformer frequency analysis device and method |
| CN114355297A (en) * | 2021-12-06 | 2022-04-15 | 电子科技大学 | Frequency control array design method based on multi-window mixed weighted frequency offset |
| CN114355297B (en) * | 2021-12-06 | 2023-11-17 | 电子科技大学 | Frequency control array design method based on multi-window mixed weighting type frequency offset |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106771586B (en) | 2019-07-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106771586A (en) | The loop signal analysis method and device of a kind of DC control protection board | |
| Krishnanand et al. | Detection, classification, and location of faults in power transmission lines | |
| Wu et al. | Ultra-high-speed directional protection of transmission lines using mathematical morphology | |
| CN103308766A (en) | Harmonic analysis method based on Kaiser self-convolution window dual-spectrum line interpolation FFT (Fast Fourier Transform) and device thereof | |
| Hasheminejad et al. | Traveling‐wave‐based protection of parallel transmission lines using Teager energy operator and fuzzy systems | |
| Yusuff et al. | Determinant-based feature extraction for fault detection and classification for power transmission lines | |
| CN109782063A (en) | A kind of dynamic m-Acetyl chlorophosphonazo analysis method based on three spectral line interpolation FFT of Nuttall self-convolution window | |
| CN105137180B (en) | High-precision harmonic analysis method based on six four spectral line interpolations of Cosine Window | |
| CN106483374A (en) | A kind of harmonic wave harmonic detection method based on Nuttall double window whole phase FFT | |
| CN103353550A (en) | Method for measuring signal frequency and harmonic parameters of electric power system | |
| Wen et al. | Hanning self-convolution window and its application to harmonic analysis | |
| CN109828163A (en) | A kind of three-phase imbalance detection method for power grid | |
| Zeng et al. | Parameter estimation of power system signals based on cosine self-convolution window with desirable side-lobe behaviors | |
| Yu | Detection and correction of saturated current transformer measurements using decaying DC components | |
| Patel et al. | Fast fault detection during power swing on a hybrid transmission line using WPT | |
| CN103743949A (en) | Detection method of harmonic and inter-harmonic based on single-channel FastICA (Fast Independent Component Analysis) | |
| Tripathy et al. | Sparse S‐transform for location of faults on transmission lines operating with unified power flow controller | |
| Patel et al. | Differential voltage‐based fault detection during power swing | |
| Wang et al. | A fast algorithm to estimate phasor in power systems | |
| CN101499651B (en) | Fast acting method for microcomputer type relay protection | |
| Naseri et al. | Power transformer differential protection based on least squares algorithm with extended kernel | |
| Mishra et al. | Intrinsic time decomposition based fault location scheme for unified power flow controller compensated transmission line | |
| Li | A total least squares enhanced smart DFT technique for frequency estimation of unbalanced three-phase power systems | |
| He et al. | A traveling wave natural frequency‐based single‐ended fault location method with unknown equivalent system impedance | |
| Li et al. | Second‐order matrix pencil‐based phasor measurement algorithm for P‐class PMUs |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |