CN110082592A - A kind of ultra harmonics measurement method based on flexible atom filtering - Google Patents
A kind of ultra harmonics measurement method based on flexible atom filtering Download PDFInfo
- Publication number
- CN110082592A CN110082592A CN201910458946.5A CN201910458946A CN110082592A CN 110082592 A CN110082592 A CN 110082592A CN 201910458946 A CN201910458946 A CN 201910458946A CN 110082592 A CN110082592 A CN 110082592A
- Authority
- CN
- China
- Prior art keywords
- frequency
- ultra harmonics
- faf
- ultra
- flexible
- 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.)
- Pending
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
- G01R23/165—Spectrum analysis; Fourier analysis using filters
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The present invention relates to a kind of ultra harmonics measurement method based on flexible atom filtering, this method within the scope of superelevation subband by being arranged multiple FAF, so that filtering Measurement bandwidth covers entire superelevation subband without overlapping;After the discretization expression formula of discrete ultra harmonics signal and FAF makees inner product processing, corresponding ultra harmonics frequency and amplitude are determined according to calculated result.Compared with prior art, ultra harmonics detection method process proposed by the present invention is simple, calculation amount is few, can be realized the ultra harmonics measurement of fast accurate, provides prerequisite scientific measurement foundation for the further investigation in the fields such as ultra harmonics genesis mechanism, inhibition.
Description
Technical field
The present invention relates to the analysis of Harmonious Waves in Power Systems and fields of measurement, are filtered more particularly, to one kind based on flexible atom
Ultra harmonics measurement method.
Background technique
With the ratio of distributed power generation in the power system be continuously increased and new-energy automobile it is continuous universal, it is more next
Mostly with the equipment investment of power electronic technique into electric system, such as the new energy grid-connected powers such as photovoltaic wind use
AC/DC convertor etc. in inverter, electric automobile charging pile.This kind of equipment can bring 2~150kHz after being connected to the grid
Ultra harmonics (supraharmonics), new challenge is brought to power grid.Simultaneously as in power grid ultra harmonics hair
Persistent levels rising is penetrated, measuring instrumentss instruction distortion, abnormal work such as noise increase of electric ballast are likely to result in
Situation.In addition, situations such as ultra harmonics also results in electrical equipment abnormal shutdown, and power line carrier communication fails.Therefore, exist
Under the background that ultra harmonics problem becomes increasingly conspicuous, the detection of ultra harmonics has great importance, and can be superelevation time
The further investigation in the fields such as harmonic wave genesis mechanism, inhibition provides prerequisite scientific measurement foundation.
In recent years, in Electric Power Harmonic Analysis field, about the detection method of classical harmonic wave, current utilization is more extensive
There is fourier transform method.However, Fourier transformation has some disadvantages, Fourier transformation must take into account overall signal, therefore
The certain local information of disturbing signal cannot be characterized.It is directed to ultra harmonics simultaneously, according to Nyquist's theorem, sample frequency is needed
Reach 300kHz or more, this results in the calculation amount of signal processing to greatly increase, while processing speed decline, real-time compared with
Difference.
Summary of the invention
It is an object of the present invention to overcome the above-mentioned drawbacks of the prior art and provide one kind to be based on flexible atom
The ultra harmonics measurement method of filtering.
The purpose of the present invention can be achieved through the following technical solutions:
A kind of ultra harmonics measurement method based on flexible atom filtering, includes the following steps:
S1: building flexible atom Filtering Model, the expression formula of flexible atom filter function are as follows:
In formula,For triangle orthogonal basis, t ∈ (- ∞ ,+∞), ωcFor multiple bandpass center frequency;A is scale factor;Fb
For bandwidth adjustment parameter, for adjustingTime-frequency domain bandwidth, to obtain more compact frequency domain supported collection and shorter sight
Survey window length.For the normalization coefficient that FAF function is normalized.Wherein Fb, a and ωcIt is non-cross, it realizes
The decoupling of window function and ORTHOGONAL TRIANGULAR base, this characteristic keep the filtering control strategy of FAF more flexible.
S2: multiple FAF are arranged within the scope of superelevation subband, and (Flexible Atom Filter, flexible atom filter letter
Number) so that filtering Measurement bandwidth covers entire superelevation subband without overlapping.
The expression formula of the Fourier transformation of FAF are as follows:
In formula, when frequency is located at centre frequency ωcWhen place, ψFAFThe gain of (ω) reaches maximum value 1, when frequency is not located at
When centre frequency, need to make its gain 1 by correction factor.Window of the FAF on frequency domain can be indicated by such as following formula:
Wherein,For ψFAFThe frequency window radius of the generating function (being free of scale factor) of (ω), andFor the frequency of FAF
Window radius determines the Measurement bandwidth of filter, and then determines the model resolution of algorithm.
Each attainable frequency resolution of FAF isIt is multiple by being arranged within the scope of superelevation subband
FAF, by the centre frequency ω of FAFcWithStep-length be gradually adjusted to 150kHz from 2kHz.When a takes smaller value, frequency window half
Diameter is wider, and frequency resolution is lower, but when window it is shorter, temporal resolution is high;When a takes the larger value, frequency window radius is relatively narrow, frequency point
Resolution is higher, but when window it is longer, temporal resolution is short.Bandwidth parameter FbIt is able to ascend the adjusting energy to filter time frequency window size
Power, to obtain compact frequency domain supported collection and higher model resolution.Therefore according to the resolution requirement setup parameter a of frequency
With FbSo that filtering Measurement bandwidth covers entire ultra harmonics frequency band without overlapping, and reach higher measurement accuracy.
S3: sliding-model control is carried out to system sub-harmonic wave signal, sub-harmonic wave is adopted with certain sampling interval
Sample obtains ultra harmonics ingredient, and samples to FAF function.
The detailed process of sliding-model control is carried out to system sub-harmonic wave signal are as follows:
Assuming that system voltage, electric current continuous signal are shown below:
In formula, fhIndicate harmonic wave, and fh=hf1, f1Indicate fundamental frequency.H indicates overtone order, AhAnd θhIs characterized respectively
The amplitude and phase of h subharmonic ingredient.
Through TSSampling interval ultra harmonics is sampled after obtain corresponding ultra harmonics ingredient expression formula:
In formula, Ash、θshThe respectively amplitude and phase of sub-harmonic wave, fshFor the frequency of ultra harmonics, TsFor sampling week
Phase, n are the sampled value serial number of continuous signal, and k=0,1 ..., K-1, K is signal sequence length, and sh is the number of sub-harmonic wave.
It is equally the sampling interval to FAF function using Ts since filtering calculating of the FAF to real signal is discretization
It is sampled, can be obtained:
In formula,It isDiscretization expression formula, k is sample sequence, and window is long when value range is one
Degree, i.e. k ∈ (- aFbΔt*/Ts,aFbΔt*/Ts)。
S4: being filtered calculating to the ultra harmonics after sampling by multiple FAF, and it is multiple to obtain flexible atom filtering output
The real component and imaginary of vector.
FAF is defined as the filtering calculating of x (t)With the inner product operation of x (k), inner product result is a complex phase
Amount, wherein containing the parameter informations such as modal amplitudes, frequency.Real and imaginary parts and x (k) inner product result distinguish
For the real and imaginary parts for exporting complex phase amount.
The real component of flexible atom filtering output complex vector are as follows:
The imaginary of flexible atom filtering output complex vector are as follows:
S5: the parameter information for including in the real component and imaginary of output complex vector, benefit are filtered according to flexible atom
The amplitude and frequency of ultra harmonics are obtained with relevant calculation formula, complete ultra harmonics measurement.
Calculation formula is shown below:
The amplitude of ultra harmonics:
The frequency of ultra harmonics:
In above-mentioned two formula, A (k) is the instantaneous value of ultra harmonics amplitude envelope line in corresponding band;F (k) is instantaneous
Frequency (Hz);CfFor amplitude envelope line gain-boosted op amp coefficient, with instantaneous frequency f (k) and ΨFAF
(ω) is related.
Compared with prior art, the method for the present invention within the scope of superelevation subband by being arranged multiple FAF, according to superelevation time
Harmonic signal determines corresponding ultra harmonics frequency and amplitude to the calculated result after FAF function inner product, and process is simple, calculation amount
It is few, it can be realized the ultra harmonics measurement of fast accurate, be the further investigation in the fields such as ultra harmonics genesis mechanism, inhibition
Provide prerequisite scientific measurement foundation.
Detailed description of the invention
Fig. 1 is the flow diagram of the method for the present invention.
Specific embodiment
The present invention is described in detail with specific embodiment below in conjunction with the accompanying drawings.Obviously, described embodiment is this
A part of the embodiment of invention, rather than whole embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art exist
Every other embodiment obtained under the premise of creative work is not made, all should belong to the scope of protection of the invention.
As shown in Figure 1, being specifically included the present invention relates to a kind of ultra harmonics measurement method based on flexible atom filtering
The following steps:
Step 1: construction flexible atom filter function.
Flexible atom filter function indicates are as follows:
In formula,For triangle orthogonal basis, t ∈ (- ∞ ,+∞), ωcFor multiple bandpass center frequency;A is scale factor;Fb
For bandwidth adjustment parameter, for adjustingTime-frequency domain bandwidth, to obtain more compact frequency domain supported collection and shorter sight
Survey window length.For the normalization coefficient that FAF function is normalized.Wherein Fb, a and ωcIt is non-cross.
Step 2: multiple FAF are arranged within the scope of superelevation subband, by the centre frequency ω of FAFcWithStep
It is long to be gradually adjusted to 150kHz from 2kHz, according to resolution requirement the setup parameter a and F of frequencyb。
Step 3: sliding-model control is carried out to system sub-harmonic wave signal and FAF function, with certain sampling interval pair
Sub-harmonic wave is sampled.
Carry out the particular content of sliding-model control are as follows:
Assuming that system voltage, electric current continuous signal are shown below:
In formula, fhIndicate harmonic wave, and fh=hf1, f1Indicate fundamental frequency.H indicates overtone order, AhAnd θhIs characterized respectively
The amplitude and phase of h subharmonic ingredient.
Through TSSampling interval ultra harmonics is sampled after obtain corresponding ultra harmonics ingredient:
In formula, Ash、θshThe respectively amplitude and phase of sub-harmonic wave, fshFor the frequency of ultra harmonics, TsFor sampling week
Phase, n are the sampled value serial number of continuous signal, and k=0,1 ..., K-1, K is signal sequence length, and sh is the number of sub-harmonic wave.
Step 4: being filtered calculating to the ultra harmonics after sampling by multiple FAF, it is defeated to obtain flexible atom filtering
The real component and imaginary of complex vector out.
The calculating formula of real component are as follows:
The calculating formula of imaginary are as follows:
Step 5: passing through relevant calculation according to the parameter information for including in real and imaginary parts in the inner product result of step 4
Formula obtains the amplitude and frequency of ultra harmonics.
The calculating formula of the amplitude of ultra harmonics are as follows:
The calculating formula of the frequency of ultra harmonics are as follows:
Wherein, A (k) is the instantaneous value of ultra harmonics amplitude envelope line in corresponding band;F (k) is instantaneous frequency (Hz);CfFor amplitude envelope line gain-boosted op amp coefficient, with instantaneous frequency f (k) and ΨFAF(ω) is related.
The above description is merely a specific embodiment, but scope of protection of the present invention is not limited thereto, any
The staff for being familiar with the art in the technical scope disclosed by the present invention, can readily occur in various equivalent modifications or replace
It changes, these modifications or substitutions should be covered by the protection scope of the present invention.Therefore, protection scope of the present invention should be with right
It is required that protection scope subject to.
Claims (6)
1. a kind of ultra harmonics measurement method based on flexible atom filtering, which is characterized in that this method includes the following steps:
1) flexible atom Filtering Model is constructed;
2) it is arranged within the scope of superelevation subband and multiple covers entire superelevation subband without overlapping Measurement bandwidth will be filtered
Flexible atom filter function;
3) sliding-model control is carried out to the sub-harmonic wave signal of system, sub-harmonic wave is sampled with certain sampling interval,
Ultra harmonics ingredient is obtained, and FAF function is sampled;
4) calculating is filtered to the ultra harmonics after sampling by multiple FAF, obtains flexible atom filtering output complex vector
Real component and imaginary;
5) real component and imaginary of the flexible atom filtering output complex vector obtained according to step 4), it is humorous to obtain superelevation time
The amplitude and frequency of wave complete ultra harmonics measurement.
2. a kind of ultra harmonics measurement method based on flexible atom filtering according to claim 1, which is characterized in that
The expression formula of the transmission function of the flexible atom Filtering Model are as follows:
In formula:For triangle orthogonal basis, t ∈ (- ∞ ,+∞), ωcFor multiple bandpass center frequency, a is scale factor, FbFor for
AdjustmentTime-frequency domain bandwidth bandwidth adjustment parameter,Return for what property atom filtering function was normalized
One changes coefficient, Fb, a and ωcIt is non-cross.
3. a kind of ultra harmonics measurement method based on flexible atom filtering according to claim 2, which is characterized in that
The particular content of step 2) are as follows:
The transmission function of flexible atom Filtering Model is subjected to Fourier transformation:
In formula, when frequency is located at centre frequency ωcWhen place, ψFAFThe gain of (ω) reaches maximum value 1, when frequency is not located at center frequency
When rate, need to make its gain 1 by correction factor, window of the flexible atom Filtering Model on frequency domain is expressed from the next:
Wherein,For ψFAFThe frequency window radius of the generating function of (ω),For the frequency window radius of FAF, each FAF can reach
Frequency resolution beBy the way that multiple FAF are arranged within the scope of superelevation subband, by the centre frequency ω of FAFcWithStep-length be gradually adjusted to 150kHz from 2kHz, according to resolution requirement the setup parameter a and F of frequencyb, survey filtering
Amount bandwidth covers entire ultra harmonics frequency band without overlapping.
4. a kind of ultra harmonics measurement method based on flexible atom filtering according to claim 3, which is characterized in that
The particular content of step 3) are as follows:
Assuming that system voltage, electric current continuous signal are shown below:
In formula: fhFor harmonic wave, and fh=hf1, f1For fundamental frequency, h is overtone order, AhAnd θhRespectively h subharmonic ingredient
Amplitude and phase;
Through TSSampling interval ultra harmonics is sampled after obtain corresponding ultra harmonics ingredient expression formula it is as follows:
In formula: Ash、θshThe respectively amplitude and phase of sub-harmonic wave, fshFor the frequency of ultra harmonics, TsFor sampling period, n
For the sampled value serial number of continuous signal, k=0,1 ..., K-1, K is signal sequence length, and sh is the number of sub-harmonic wave.
5. a kind of ultra harmonics measurement method based on flexible atom filtering according to claim 4, which is characterized in that
In step 4), the calculation expression of the real component of flexible atom filtering output complex vector are as follows:
The calculation expression of the imaginary of flexible atom filtering output complex vector are as follows:
6. a kind of ultra harmonics measurement method based on flexible atom filtering according to claim 5, which is characterized in that
In step 5), the expression formula of the amplitude of ultra harmonics are as follows:
The expression formula of the frequency of ultra harmonics are as follows:
In formula: A (k) is the instantaneous value of ultra harmonics amplitude envelope line in corresponding band, and f (k) is instantaneous frequency (Hz),CfFor amplitude envelope line gain-boosted op amp coefficient.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910458946.5A CN110082592A (en) | 2019-05-29 | 2019-05-29 | A kind of ultra harmonics measurement method based on flexible atom filtering |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910458946.5A CN110082592A (en) | 2019-05-29 | 2019-05-29 | A kind of ultra harmonics measurement method based on flexible atom filtering |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110082592A true CN110082592A (en) | 2019-08-02 |
Family
ID=67422526
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910458946.5A Pending CN110082592A (en) | 2019-05-29 | 2019-05-29 | A kind of ultra harmonics measurement method based on flexible atom filtering |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110082592A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111537796A (en) * | 2020-05-12 | 2020-08-14 | 华北电力大学 | Ultrahigh harmonic measurement method based on fixed-frequency asynchronous sampling |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5794203A (en) * | 1994-03-22 | 1998-08-11 | Kehoe; Thomas David | Biofeedback system for speech disorders |
CN101368987A (en) * | 2008-09-27 | 2009-02-18 | 武汉大学 | Power system mesomeric state harmonic wave and/or inter-harmonic wave measuring method |
CN102305892A (en) * | 2011-07-13 | 2012-01-04 | 上海仪器仪表研究所 | Novel power grid harmonic detection analysis device |
CN103197143A (en) * | 2013-02-28 | 2013-07-10 | 哈尔滨工业大学 | Harmonic and inter-harmonic detection method based on Hanning-window FFT algorithm and traversal filtering |
CN106154035A (en) * | 2016-06-20 | 2016-11-23 | 哈尔滨工业大学 | A kind of quickly harmonic wave and harmonic detection method |
CN109164298A (en) * | 2018-07-25 | 2019-01-08 | 陕西科技大学 | A kind of compressed sensing based ultra harmonics detection device and detection method |
-
2019
- 2019-05-29 CN CN201910458946.5A patent/CN110082592A/en active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5794203A (en) * | 1994-03-22 | 1998-08-11 | Kehoe; Thomas David | Biofeedback system for speech disorders |
CN101368987A (en) * | 2008-09-27 | 2009-02-18 | 武汉大学 | Power system mesomeric state harmonic wave and/or inter-harmonic wave measuring method |
CN102305892A (en) * | 2011-07-13 | 2012-01-04 | 上海仪器仪表研究所 | Novel power grid harmonic detection analysis device |
CN103197143A (en) * | 2013-02-28 | 2013-07-10 | 哈尔滨工业大学 | Harmonic and inter-harmonic detection method based on Hanning-window FFT algorithm and traversal filtering |
CN106154035A (en) * | 2016-06-20 | 2016-11-23 | 哈尔滨工业大学 | A kind of quickly harmonic wave and harmonic detection method |
CN109164298A (en) * | 2018-07-25 | 2019-01-08 | 陕西科技大学 | A kind of compressed sensing based ultra harmonics detection device and detection method |
Non-Patent Citations (1)
Title |
---|
刘林等: "用于次同步振荡模态提取的柔性原子滤波", 《电工技术学报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111537796A (en) * | 2020-05-12 | 2020-08-14 | 华北电力大学 | Ultrahigh harmonic measurement method based on fixed-frequency asynchronous sampling |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6636608B2 (en) | Method and apparatus for detecting voltage in a supply network | |
CN105048995B (en) | Butterworth digital filter and realize the adaptive filtering method of frequency using it | |
CN106645919B (en) | The full spectral power of power grid based on three-phase instantaneous power vibrates phasor method for synchronously measuring | |
CN107966611A (en) | A kind of supersynchronous harmonic detecting method of electric system based on vector matching time | |
CN107844670A (en) | The computational methods of sample size needed for a kind of harmonic wave statistics | |
Zhang et al. | Shifting window average method for phasor measurement at offnominal frequencies | |
CN106155981A (en) | A kind of sub-synchronous oscillation parameter detection method | |
CN104362628A (en) | Control method and device for avoiding large-wind-field resonance | |
Shitole et al. | Comparative evaluation of synchronization techniques for grid interconnection of renewable energy sources | |
CN110082592A (en) | A kind of ultra harmonics measurement method based on flexible atom filtering | |
CN102868163A (en) | Three-phase soft phase-locked ring for adaptive grid-connected converter | |
CN109828154B (en) | Three-phase power grid impedance measurement method based on sub-band composite orthogonal pulse injection | |
CN103227640A (en) | Phase detection method of harmonic power grid and phase-locked loop | |
CN112255457B (en) | Phase angle difference measuring method suitable for automatic quasi-synchronization device | |
CN104483547B (en) | The filtering method and system of electric power signal | |
CN109342816B (en) | Detection method for spectrum leakage in power quality monitoring | |
Elvira-Ortiz et al. | Study of the harmonic and interharmonic content in electrical signals from photovoltaic generation and their relationship with environmental factors | |
CN112485522B (en) | Electric energy data perception-based flat-top window function synchronous phasor measurement method and device | |
Dubey et al. | Detection of power quality disturbances in presence of DFIG wind farm using Wavelet Transform based energy function | |
CN104483563A (en) | Method and system for synchronous sampling of power signals | |
CN107271773B (en) | Rapid detection method for harmonic waves of power system | |
CN105629064B (en) | A kind of sine wave low-pass filtering extracting method | |
CN109726490A (en) | A kind of more sinusoidal signal design methods of low-frequency range for the identification of POWER SYSTEM STATE spatial model | |
Zhang et al. | Detection of Harmonic Components using the FFT and Instantaneous Reactive Power Theory | |
CN111537796A (en) | Ultrahigh harmonic measurement method based on fixed-frequency asynchronous sampling |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190802 |
|
RJ01 | Rejection of invention patent application after publication |