Invention content
In view of the deficienciess of the prior art, it is an object of the present invention to provide a kind of, the broadband multi-frequency based on more algorithm fusions is believed
Number measurement method, fft analysis, Prony algorithm, Grandke frequency calibrating methods, digital filter are blended, and are realized
Accurate, the quick detections of each component parameters of broadband multi-frequency signal, it is existing for monitoring low-frequency oscillation of electric power system, sub-synchronous oscillation etc.
As providing effective monitoring means.
To achieve the goals above, the present invention is to realize by the following technical solutions:
A kind of broadband multi-frequency signal measurement method based on more algorithm fusions of the present invention, including following steps:
Step (1):To being sampled containing multicomponent broadband multi-frequency signal voltage, electric current for AD acquisitions;
Step (2):Each frequency content in step (1) contained by acquisition signal, including low frequency are detected by fft algorithm
Oscillating component, harmonic wave and m-Acetyl chlorophosphonazo component and power frequency component;
Step (3):For low-frequency oscillation component, after filtering out disturbing high frequency components therein by wave digital lowpass filter;
It reuses the existing Prony algorithms of extension and calculates the amplitude of each component, frequency;
Step (4):For harmonic wave, m-Acetyl chlorophosphonazo component, bandpass filtering is carried out by digital band-pass filter, obtains removing power frequency
Other frequency components except component and low-frequency oscillation component;Reuse what the FFT methods based on Grandke Spectrum Corrections calculated
Amplitude, frequency parameter;
Step (5):For power frequency component, bandpass filtering is carried out by digital band-pass filter, it is existing fixed then to use
Interval sampling phasor correcting algorithm obtains the phasor, frequency of power frequency component, frequency change rate in broadband multi-frequency signal.
In step (1), current signal resampling frequency fsFor 4KHz.
In step (1), broadband multi-frequency signal voltage, the mathematical model of current signal are as follows:
In formula (1), XiIt is each frequency component of broadband multi-frequency signal, wiFor the frequency of each frequency component,For each frequency point
The initial phase angle of amount, AiFor the amplitude of each frequency component, i is each frequency component subscript, and i=1,2,3 ..., M, M is that frequency component is total
Number, t is the time.
For low-frequency oscillation component, wiFrequency range is between 0.1Hz-2.5Hz;For power frequency component, wiFrequency range exists
Between 45Hz-55Hz.
In step (2), the long N of data window of FFT spectrum analysis is 4096, and its FFT spectrum is denoted as X (n).
In step (3), the wave digital lowpass filter is 5Hz by frequency.
In step (4), the passband frequency range of digital band-pass filter is 1Hz.
In step (5), the passband frequency range of digital band-pass filter is 10Hz, and bandpass filter centre frequency is
50Hz。
In step (5), bandpass filtering is carried out by digital band-pass filter, obtains the power frequency component of 45Hz-55Hz.
The present invention blends fft analysis, Prony algorithm, Grandke frequency calibrating methods, digital filter, real
The various frequencies such as ground state power frequency amount and higher hamonic wave, m-Acetyl chlorophosphonazo, the low-frequency oscillation component in quick analysis electric power signal are showed
Centre frequency, the amplitude of component;To monitor the low-frequency oscillation in electric system, subsynchronous/supersynchronous oscillation and harmonic wave point in time
Analysis takes measures to ensure unit safety and power system stability operation, provides a kind of effective, reliable means.
Specific implementation mode
To make the technical means, the creative features, the aims and the efficiencies achieved by the present invention be easy to understand, with reference to
Specific implementation mode, the present invention is further explained.
Referring to Fig. 1, a kind of broadband multi-frequency signal based on more algorithm fusions of the invention unifies method for fast measuring, including
Following steps:
Step (1):To being sampled containing multicomponent broadband multi-frequency signal voltage, electric current for AD acquisitions;
Step (2):Frequency content therein is detected by using fft algorithm, realizes the automatic detection and classification of frequency band;
Step (3):For low-frequency oscillation component signal, divided using High-frequency Interference therein is filtered out after wave digital lowpass filter
After amount, Prony Algorithm Analysis its low frequency component parameter, including the amplitude of each component, frequency are reused;
Step (4):For harmonic wave, m-Acetyl chlorophosphonazo component, the low-frequency component and high-frequency noise of signal are filtered out through bandpass filter
Component analyzes the electric parameter of each component, including the amplitude of each component, frequency using Grandke spectrum correcting methods;
Step (5):For power frequency component, then after bandpass filtering, using fixed sampling interval technique phasor correcting algorithm, realize
Dynamic Phasors calculate, and income analysis result parameter includes amplitude, frequency, phase.
In step (1), the current signal resampling frequency fsFor 4KHz.
Shown in the mathematical model such as formula (1) of broadband multi-frequency signal voltage, current signal described in step (1):
X in formula 1iFor each frequency component of broadband multi-frequency signal, wiFor the frequency of each component,For the initial phase angle of each component,
AiFor the amplitude of each component.For low-frequency oscillation component, wiFrequency range is usually between 0.1Hz-2.5Hz;For power frequency point
Amount, wiFrequency range is usually between 45Hz-55Hz.
The long N of data window of FFT spectrum analysis described in step (2) is 4096, and its FFT spectrum is denoted as X (n).
Low-pass filter described in step (3) is 5Hz by frequency, and step is:
Step (3-1):Low-pass filtering is carried out to the broadband multi-frequency signal.
Step (3-2):Prony algorithm, analysis is used to calculate amplitude, the frequency of each component filtered signal.
Bandpass filter passband frequency range described in step (4) is 1Hz, and step is:
Step (4-1):Bandpass filtering is carried out to the broadband multi-frequency signal, obtains each component except power frequency component and low frequency point
Other frequency components X except amounti。
Step (4-2):X is calculated using the FFT methods based on Grandke Spectrum CorrectionsiAmplitude, frequency parameter;
Bandpass filter passband frequency range described in step (5) is 10Hz, and bandpass filter centre frequency is 50Hz, and
Step is:
Step (5-1):Bandpass filtering is carried out to the broadband multi-frequency signal, obtains the power frequency component of 45Hz-55Hz.
Step (5-2):It using fixed sampling interval technique phasor correcting algorithm, realizes that Dynamic Phasors calculate, calculates broadband multi-frequency letter
The phasor, frequency of power frequency component, frequency change rate in number.
The emulation signal applied is:
X (n)=0.5*cos (2*pi*2*n/4000)+0.4*
cos(2*pi*15*n/4000)+1*cos(2*pi*50*n/4000)+0.3*cos(2*pi*75*n/4000)+
0.4*cos(2*pi*220*n/4000);
4096 fft analysis are carried out to the signal, it is as shown in Figure 2 to obtain its spectrogram.As shown in Figure 2, by FFT spectrum
After analysis, 5 frequency components can be obtained, frequency is F1=1.953Hz, F2=14.65Hz, F3=49.8Hz, F4=respectively
75.2Hz F5=219.7Hz.
Prony algorithm is reused after 10Hz low-pass filtering for the low frequency component that frequency is 1.953Hz, is obtained
To shown in parameter such as table 1 (low frequency component result of calculation):
Table 1
Serial number |
Frequency |
Amplitude |
1 |
2.003Hz |
0.5005 |
For the subsynchronous component that frequency is 14.65Hz, by selected bandpass filter, (centre frequency is
14.6Hz, passband 1Hz) after, it reuses Grandke spectrum correcting methods and calculates its amplitude and frequency parameter, as table 2 is (secondary same
Walk component result of calculation) shown in:
Table 2
Serial number |
Frequency |
Amplitude |
1 |
14.995Hz |
0.401 |
For the supersynchronous component that frequency is 75.2Hz, by selected bandpass filter, (centre frequency is 75Hz, is led to
Band is 1Hz) after, it reuses Grandke spectrum correcting methods and calculates its amplitude and frequency parameter, as (supersynchronous component calculates table 3
As a result shown in):
Table 3
Serial number |
Frequency |
Amplitude |
1 |
75.004Hz |
0.1995 |
For the m-Acetyl chlorophosphonazo component that frequency is 219.7Hz, by selected bandpass filter (centre frequency is 220Hz,
Passband is 1Hz) after, it reuses Grandke spectrum correcting methods and calculates its amplitude and frequency parameter, such as (the m-Acetyl chlorophosphonazo component meter of table 4
Calculate result) shown in:
Table 4
Serial number |
Frequency |
Amplitude |
1 |
219.997Hz |
0.4003 |
For the power frequency component of 49.8Hz, using phasor correcting algorithm, as a result such as table 5 (power frequency component result of calculation) institute
Show:
Table 5
Serial number |
Frequency |
Amplitude |
1 |
50.000Hz |
0.9998 |
The above shows and describes the basic principles and main features of the present invention and the advantages of the present invention.The technology of the industry
Personnel are it should be appreciated that the present invention is not limited to the above embodiments, and the above embodiments and description only describe this
The principle of invention, without departing from the spirit and scope of the present invention, various changes and improvements may be made to the invention, these changes
Change and improvement all fall within the protetion scope of the claimed invention.The claimed scope of the invention by appended claims and its
Equivalent thereof.