CN107085144B - A kind of method of rapid survey Harmonious Waves in Power Systems - Google Patents
A kind of method of rapid survey Harmonious Waves in Power Systems Download PDFInfo
- Publication number
- CN107085144B CN107085144B CN201710297370.XA CN201710297370A CN107085144B CN 107085144 B CN107085144 B CN 107085144B CN 201710297370 A CN201710297370 A CN 201710297370A CN 107085144 B CN107085144 B CN 107085144B
- Authority
- CN
- China
- Prior art keywords
- sampling
- sampled
- adc
- data
- fft
- 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.)
- Active
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
- G01R23/167—Spectrum analysis; Fourier analysis using filters with digital filters
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
A kind of method of rapid survey Harmonious Waves in Power Systems, comprising the following steps: over-sampling ADC carries out over-sampling to signal, obtains original sampling data；It carries out down-sampled to the original sampling data of over-sampling ADC output and saves；Windowing FFT processing is carried out to the data after down-sampled；Calculate fundamental frequency；It calculates quasi-synchro sampling rate and saves；Calculated quasi-synchro sampling rate carries out interpolation processing to the original sampling data that over-sampling ADC is exported；Data Jing Guo interpolation processing are carried out down-sampled and are saved；Data are handled using rectangular window FFT；The spectral magnitude of each FFT obtained after rectangular window FFT processing is multiplied by penalty coefficient, amplitude compensation is carried out, exports result.The method of the present invention improves the signal-to-noise ratio of inband signaling using oversampling technique, under the premise of guaranteeing frequency resolution, greatly reduces computation complexity and improves the accuracy of fundamental frequency calculating, and pass through the down-sampled occupancy for reducing caching.
Description
Technical field
The invention belongs to electricity field more particularly to a kind of methods for calculating analysis mains by harmonics.
Background technique
In recent years, with the rapid growth of Electricity Demand, electrical load is more and more, and all non-linear negative
Carry, for example, switched-mode power supply, electronic ballast for fluorescent lamp, speed-regulating actuator, uninterruptible power supply, magnetive cord equipment and
The household electrical appliance such as television set can all generate harmonic wave, so that grid supply quality reduces.Therefore, to the voltage of power supply grid, electric current
Signal carries out frequency analysis, to carry out power network compensation and purification becomes more and more important.In order to as fast as possible accurately into
Row power network compensation and purification, the real-time and precise requirements of frequency analysis are also higher and higher, how rapidly and accurately to calculate humorous
Wave becomes the hot spot of a research.
Fast Fourier Transform (Fast Fourier Transformation, hereinafter referred to as FFT) is to calculate the master of harmonic wave
Stream tool, but due to truncation effect, there are spectrum leakages and fence effect by FFT, in order to accurately calculate harmonic wave, current common side
Method is as follows:
1, using windowing FFT, in parameters such as frequency, amplitude and phase of the frequency domain interpolation to obtain fundamental wave and each harmonic wave
Information.A kind of electric power based on windows and interpolated FFT as disclosed in the Chinese invention patent application of Publication No. CN101701982A
System harmonics detection method, carries out windowing FFT operation using hanning window, carries out frequency domain interpolation then to calculate each harmonic wave
Accurate frequency, amplitude and phase.This method needs to use the top of FFT spectral line and adjacent when calculating interpolation coefficient
The ratio of the amplitude on secondary peak, therefore the spectral resolution of FFT sufficiently small must just can guarantee that the two spectral lines do not include other
The amplitude of frequency, when the use of 512 point FFT, sample rate being such as 6.4KHz, the operation for completing a complete FFT needs to sample 512
Sampling point, the sampling time of consuming are T=N/Fs=512/6400Hz=0.08s=80ms, and common are source electric power on the market
The response time of filter is usually less than 20ms, thus while the accuracy of this method is very high, but due to need one it is opposite
Small frequency resolution distinguishes fundamental wave and each harmonic wave, causes the sampling time too long, is not able to satisfy practical application request.
2, plesiochronous by tracking the frequency variation realization of fundamental wave, making quasi-synchro sampling rate is the integral multiple of fundamental frequency,
To inhibit spectrum leakage.A kind of electric energy quality harmonic as disclosed in the Chinese invention patent application of Publication No. CN103969507A
Analysis method, by finding out frequency shift (FS), then using cosine function corresponding with principal wave harmonic wave frequency and SIN function come into
Row relevant calculation, acquires amplitude and initial phase angle.This method is needed each fundamental wave harmony wave respectively multiplied by corresponding cosine function
And SIN function, since each sampled point will be multiplied by cosine function and SIN function, and each fundamental wave and harmonic wave are corresponding
Cosine function and SIN function it is all different, cause calculation amount very big, and frequency drift parameter calculating be easy by noise
Interference, make fundamental frequency calculate inaccuracy.
3, quasi-synchro sampling is realized using the algorithm of time domain interpolation.Such as the Chinese invention of Publication No. CN101915874A
A kind of harmonic detecting method based on Fourier transform of patent application publication, the process of this method is as shown in Figure 1, by over-sampling
Analog-digital converter (Analog to Digital Converter, hereinafter referred to as ADC) is to coherent signal (including electric current, voltage etc.
Signal) over-sampling is carried out, it is stored in caching (memory) after obtaining original sampling data, is then adopted using filter from original
Fundamental signal is extracted in sample data, then calculates the period of fundamental wave using zero crossing, and resampling is confirmed according to the period of fundamental wave
Interval carries out resampling from original sampling data according to the interval and FFT is handled, and the period of the signal after resampling and fundamental wave are tight
Lattice are synchronous, thus avoid the spectrum leakage problem of common FFT.But this method needs the original for obtaining analog-digital converter over-sampling
Beginning data are stored in the buffer, since over-sampling rate is very high, cause occupancy caching very big, and using filter to original
Filter result is also stored into the fortune to carry out the next extraction primitive period in the buffer when sampled data is filtered
It calculates, for the embodiment of the patent, when analog-digital converter uses the sample rate of 50*128*40=256KHz, if modulus turns
The bit wide of parallel operation is 16, then the cache size for saving the signal needs of a primitive period is 128*40*16=81920bit,
When the sample rate of filter is 50*128*2=12.8KHz, saving the cache size that filtered data need is 128*2*
16=4096bit then needs the caching of 81920+4096=86016bit in total.The area that the caching occupies in integrated circuits
It can be very big；On the other hand, this method is to calculate the primitive period by calculating the time difference of two zero crossings, therefore be very easy to
By noise jamming, the period is caused to calculate inaccuracy.
Summary of the invention
And a kind of occupancy caching few the purpose of the present invention is to provide calculation amount is small, can quickly analyze each in electric system
The method of harmonic parameters.
To achieve the goals above, the present invention takes following technical solution:
A kind of method of rapid survey Harmonious Waves in Power Systems, comprising the following steps:
Sampled data step is obtained, over-sampling ADC carries out over-sampling to signal, obtains original sampling data；
Quasi-synchro sampling rate step is calculated, steps are as follows:
It carries out down-sampled to the original sampling data of over-sampling ADC output and saves；
Windowing FFT processing is carried out to the data after down-sampled；
Calculate fundamental frequency；
It calculates quasi-synchro sampling rate and saves, quasi-synchro sampling rate Fsnew=Ff× N ' × M, wherein FfFor fundamental frequency, M
For over-sampling multiple, N ' is subsequent plus rectangular window FFT points；
Sampling rate conversion step, steps are as follows:
Calculated quasi-synchro sampling rate carries out interpolation processing to the original sampling data that over-sampling ADC is exported；
Data Jing Guo interpolation processing are carried out down-sampled and are saved；
Plesiochronous FFT processing step, steps are as follows:
Data after sampled rate conversion are carried out plus rectangular window FFT is handled；
Judge whether to need amplitude compensation, be obtained after directly exporting if as a result, rectangular window FFT processing otherwise will be added
To the spectral magnitude of each FFT be multiplied by penalty coefficient, export result after carrying out amplitude compensation.
More specifically, the sample rate F of over-sampling ADCs=Fnyquist× Oversample, wherein Oversample was to adopt
Sample multiple, FnyquistFor the sample rate for meeting Nyquist's theorem.
More specifically, down-sampled processing is carried out to data using cascade integral comb filter.
More specifically, the data after down-sampled are carried out when calculating quasi-synchro sampling rate plus nuttall window FFT or Jia Hanning
Window FFT processing.
More specifically, it is calculated using the following equation fundamental frequency: FfAfter=Δ f × k, the Δ f in formula are windowing FFT processing
The spectral resolution of signal, k refer to fundamental frequency on kth FFT spectral line.
More specifically, fundamental frequencyWherein, Y (x) is indicated
The amplitude of the xth spectral line of FFT, (f/ Δ f), round (x) expression refers to is rounded x to L=round, and f is signal frequency.
More specifically, it carries out as follows to the method for data progress interpolation when sample rate conversion:
A, according to the number of the ratio calculation sampled point of the sample rate of over-sampling ADC and quasi-synchro sampling rate: Step=Fs/
Ff；
B, it is exported using the 1st ADC sampled point as the 1st interpolation sampling point, since q=2；
C, q-th of interpolation sampling point position: phase=Step × (q-1)+1 is calculated；
D, calculate q-th of interpolation sampling point position need ADC sampled data: ADC (pre)=floor (phase) and
ADC (nxt)=ceil (phase), the floor (x) in formula indicate to be rounded x to infinitesimal, and ceil (x) is indicated to x to infinite
Big to be rounded, pre indicates that re ADC sampled point of pth, nxt indicate n-th xt ADC sampled point；
E, ADC is waited to complete the sampling of n-th xt sampling point, if existing n-th xt sampling point thens follow the steps f；
F, it calculates the amplitude of q-th of interpolation sampling point and exports result, according to the sampled data meter of the obtained ADC of step d
Calculate amplitude: Y (q)=ADC (pre) × (nxt-phase)+ADC (nxt) × (phase-pre), ADC (x) indicate that x-th of ADC is adopted
Sample data；
G, q=q+1, return step c are enabled.
More specifically, interpolation processing is carried out to data using Newton interpolating method.
More specifically, whether occurred to carry out amplitude benefit with interior amplitude attenuation or gain after down-sampled filtering according to signal
It repays, penalty coefficient is the inverse with interior amplitude attenuation/yield value of desampling fir filter.
From the above technical scheme, the present invention improves the signal-to-noise ratio of inband signaling using oversampling technique, reduces ADC's
Accuracy requirement reduces system cost, while over-sampling can reduce the noise of interpolation processing；By using down-sampled technology, keep away
Exempt from caching caused by over-sampling and occupies big problem；When calculating fundamental frequency, the resolution of FFT is improved by down-sampled
Rate, and fundamental frequency is calculated in conjunction with energy barycenter correction method, reduce the interference of noise and m-Acetyl chlorophosphonazo, improves fundamental frequency
The accuracy of calculating, and reduce computation complexity, it is only necessary to it is simple down-sampled to carry out one, then calculates windowing FFT i.e.
It can；Additionally due to sample rate is the integral multiple of fundamental frequency, therefore the spectral resolution for calculating the FFT of fundamental wave and harmonic amplitude can
Very wide to accomplish, maximum can reach fundamental frequency, to reduce the points of FFT, calculation amount be reduced, when shortening sampling
Between.Compared with prior art, computationally intensive the present invention overcomes existing time domain interpolation algorithm and caching occupies big disadvantage,
Only considerably less calculation amount, which need to be expended, can be achieved with time domain interpolation quasi-synchro sampling, can accurately calculate inside the very short time
The frequency of fundamental wave and each harmonic wave, amplitude and first phase out, and it is few to occupy caching.
Detailed description of the invention
Fig. 1 is the flow chart of quasi-synchro sampling time domain interpolation algorithm in the prior art；
Fig. 2 is the flow chart of the method for the present invention.
A specific embodiment of the invention is described in more detail below in conjunction with attached drawing.
Specific embodiment
In order to above and other objects of the present invention, feature and advantage can be become apparent from, the embodiment of the present invention is cited below particularly,
It is described below in detail.
In the following description, numerous specific details are set forth in order to facilitate a full understanding of the present invention, but the present invention can be with
It is different from other way described herein using other and implements, those skilled in the art can be without prejudice to intension of the present invention
In the case of do similar popularization, therefore the present invention is not limited by the specific embodiments disclosed below.
The method of the present invention is improved on the basis of existing time domain interpolation algorithm, is obtaining voltage and electricity from power grid
The signals such as stream simultaneously calculate step using the technology of over-sampling, while using new quasi-synchro sampling rate to count when carrying out analog-to-digital conversion
Sample rate is calculated, rapidly and accurately to calculate the frequency, amplitude and first equivalent information of fundamental wave and each harmonic wave.Below with reference to Fig. 2,
The method of the present invention is described in detail, the step of the method for the present invention is as follows:
Firstly, obtaining sampled data, over-sampling ADC carries out over-sampling to signal, obtains original sampling data；
It uses the sample rate much higher than signal bandwidth to carry out over-sampling processing signal, on the one hand can be adopted by improving
The mode of sample multiple improves the signal-to-noise ratio of ADC, and theoretically, 1 times of over-sampling rate of every raising, signal-to-noise ratio improves about 3dB, due to electricity
The fundamental wave of Force system and its frequency of multiple harmonic are relatively low, therefore the sample rate of ADC required when high power over-sampling is not
Relatively inexpensive low Precision A/D C can be used in height, for example, signal bandwidth is 50Hz × 50=when needing to measure 50 subharmonic
2.5KHz meets the sample rate F of Nyquist (nyquist) theoremnyquistIt is 2 times of signal bandwidth (highest frequency), i.e.,
2.5KHz × 2=5KHz, at 64 times of over-sampling, the sample rate F of over-sampling ADCs=Fnyquist× Oversample=5KHz × 64
=320KHz, Oversample are over-sampling multiple, FnyquistFor the sample rate for meeting Nyquist's theorem, the letter in band at this time
It makes an uproar than SNR can be improvedoptimize≈3×log2(Oversample)=3 × log264=18dB；On the other hand, due to passing through
Over-sampling can save liter sampling filter when subsequent interpolation processing, to reduce the calculating error and meter of subsequent interpolation algorithm
Calculation amount；
Then, calculate quasi-synchro sampling rate, the calculating of quasi-synchro sampling rate the following steps are included:
It carries out down-sampled to the original sampling data of over-sampling ADC output and saves；In order to reduce calculation amount and caching
It occupies, and improves the accuracy rate of fundamental frequency calculating, need to carry out down-sampled, the present embodiment to the data that over-sampling ADC is exported
Over-sampling ADC is exported using cascade integral comb filter (Cascaded Integrator-Comb, hereinafter referred to as CIC)
Original sampling data carries out down-sampled processing；Since the calculating of CIC is not related to multiplication, so the integrated circuit area that it is occupied is very
Few, the down-sampled multiple of cic filter is 2048 in the present embodiment, order 8, down-sampled sample rate=320KHz/2048
=156.25Hz, due to only needing to calculate the frequency of fundamental wave, CIC's can not have to processing with interior amplitude attenuation；
Then, windowing FFT processing is carried out to the data after down-sampled, to reduce the influence of spectral leakage；The present embodiment pair
Data after down-sampled carry out plus nuttall window FFT processing, and the coefficient of nuttall window function is w (n)=a0-a1*cos (2 π
(n/N))+a2*cos (4 π (n/N))-a3*cos (6 π (n/N)), wherein n=0,1,2 ..., N-1, a0, a1, a2, a3 are
Constant, the a0=0.3635819 of the present embodiment, a1=0.4891775, a2=0.1365995, a3=0.0106411, N=32
For the points that the data after down-sampled are carried out with windowing FFT, the spectral resolution of signal is Δ f after adding nuttall window FFT to handle
=Fs/ N=156.25/32=4.8828125Hz, memory save (signal) points N needed for windowing FFT；In addition to
It is outer that windowing FFT processing can be carried out with nuttall window function, can also be carried out using other window functions such as Hanning window, rectangular window
Windowing FFT processing, in addition, nuttall window function can also be using forms such as other three ranks, five ranks；
Calculate fundamental frequency；The present embodiment is using common fixed FFT spectrum bearing calibration --- " power enhanced " meter
Calculate fundamental frequency, i.e. fundamental frequencyWherein, Δ f is adding window
The spectral resolution of FFT treated signal, the Δ f=4.8828125Hz, k of the present embodiment refer to fundamental frequency in kth item
On FFT spectral line, k can be decimal, and Y (x) indicates the amplitude of the xth spectral line of FFT, L=round (f/ Δ f), round (x) table
Show that finger is rounded x, the original signal frequency f of the present embodiment is 50Hz, then (50/4.8828125)=10 L=round；
It calculates quasi-synchro sampling rate and saves；In order to eliminate spectrum leakage, quasi-synchro sampling rate is the integer of fundamental frequency
Times, since present invention uses over-sampling ADC, quasi-synchro sampling rate also needs for fundamental frequency to be multiplied by over-sampling multiple M,
Final quasi-synchro sampling rate is Fsnew=Ff× N ' × M, N ' are the subsequent points for carrying out adding rectangular window FFT to data after interpolation,
The N ' of the present embodiment is 128, and over-sampling multiple M is 64, then Fsnew=Ff×128×64；
Then, sample rate is converted, comprising the following steps:
Calculated quasi-synchro sampling rate carries out interpolation processing to the original sampling data that over-sampling ADC is exported；
Linear interpolation method can be used or Newton interpolating method calculates interpolation, the present embodiment calculates interpolation using linear interpolation method, including following
Step:
A, according to the number of the ratio calculation sampled point of the sample rate of over-sampling ADC and quasi-synchro sampling rate: Step=Fs/
Ff；
B, it is exported using the 1st ADC sampled point as the 1st interpolation sampling point, since q=2；
C, q-th of interpolation sampling point position: phase=Step × (q-1)+1 is calculated；
D, calculate q-th of interpolation sampling point position need ADC sampled data: ADC (pre)=floor (phase) and
ADC (nxt)=ceil (phase), the floor (x) in formula indicate to be rounded x to infinitesimal, and ceil (x) is indicated to x to infinite
It is big to be rounded, as floor (2.8)=2, ceil (2.2)=3, pre indicates that re ADC sampled point of pth, nxt indicate n-th xt ADC
Sampled point；
E, ADC is waited to complete the sampling of n-th xt sampling point, if existing n-th xt sampling point thens follow the steps f；
F, it calculates the amplitude of q-th of interpolation sampling point and exports result to CIC, carry out down-sampled processing, obtained according to step d
The sampled data of the ADC arrived calculates amplitude: Y (q)=ADC (pre) × (nxt-phase)+ADC (nxt) × (phase-pre),
ADC (x) indicates x-th of ADC sampled data；
G, q=q+1, return step c are enabled；Since Harmonics Calculation is in real time, to need ceaselessly to refresh, therefore interpolation arithmetic
It is also required to constantly update, until stopping Harmonics Calculation；
To Jing Guo interpolation processing data carry out it is down-sampled, as using cascade integral comb filter (CIC) to interpolation at
Data after reason carry out down-sampled, and the down-sampled multiple of the CIC of the present embodiment is 64, order 6, it is down-sampled after sample rate
FsFFT=Fsnew/ 64, save that down-sampled treated as a result, the size of caching is the points of subsequent plus rectangular window FFT in the buffer
N ', the points N of the present embodiment rectangular window FFT ' it is 128 points, when data bit width is 16bit, then occupied caching is 128*
16=2048bit, compared to the prior art, caching occupy very small；
Finally, carrying out plesiochronous FFT processing to data, and export result, comprising the following steps:
To interpolation processing and it is down-sampled after data, i.e., sampled rate conversion after data carry out plus rectangular window FFT processing,
The FFT that random length can be used, the spectral resolution of signal after adding rectangular window FFT to handle are as follows:By
In the sampling time usedAs Δ foutFor 49Hz, T=1/49Hz=20.408ms, it can be seen that number
According to sampling time greatly reduce compared with the existing technology；
Amplitude compensation exports result；Judge whether to need to carry out amplitude compensation according to down-sampled method, if sampled
Down-sampled technology do not have to then compensate with interior amplitude attenuation or gain, the present embodiment using CIC carry out it is down-sampled, due to
CIC technology exists with interior amplitude attenuation, it is therefore desirable to be multiplied by the spectral magnitude of each FFT obtained after rectangular window FFT processing
Penalty coefficient carries out amplitude compensation, and penalty coefficient is the inverse with interior amplitude attenuation/yield value of desampling fir filter, can lead to
It crosses simulation calculation signal and penalty coefficient is determined by the attenuation/gain value of filter (CIC) amplitude afterwards, for example, signal exists
Becoming original 1/10 by the amplitude attenuation after CIC at 100Hz frequency, i.e. at 100Hz is 1/10 with interior amplitude attenuation value,
Then penalty coefficient is 10, and signal becomes original 1/20 by the amplitude attenuation after CIC at 200Hz frequency, i.e. at 200Hz
It is 1/20 with interior amplitude attenuation value, then penalty coefficient 200Hz at is 20, each not phase of the attenuation/gain of each frequency of filter
Together, and the penalty coefficient of different filter constructions and down-sampled multiple is different, and amplitude variation is also different, usually
The order of CIC is more, decays more, therefore penalty coefficient is after determining desampling fir filter (CIC) structure and down-sampled multiple
It needs to determine penalty coefficient by emulating.Interested frequency spectrum can be only compensated, operand is reduced.
The method of the present invention improves the signal-to-noise ratio of inband signaling using oversampling technique, reduces the accuracy requirement of ADC, reduces system
System cost, while the operand of sample rate conversion can be reduced, add the FFT (the present embodiment counted on a small quantity again by first down-sampled
In be 32 points) calculate fundamental frequency, under the premise of guaranteeing frequency resolution, greatly reduce computation complexity, and improve
The accuracy that fundamental frequency calculates, while the occupancy of caching is reduced by down-sampled technology.Since fundamental frequency variation is slow
Slowly, the present invention is calculating quasi-synchro sampling rate Shi Xianyong windowing FFT calculating fundamental frequency, does not have to synchronized sampling, further according to fundamental wave
Frequency calculates quasi-synchro sampling rate, then can quickly calculate harmonic wave with the method for plesiochronous FFT.
The foregoing description of the disclosed embodiments enables those skilled in the art to implement or use the present invention.
Various modifications to these embodiments will be readily apparent to those skilled in the art, as defined herein
General Principle can be realized in other embodiments without departing from the spirit or scope of the present invention.Therefore, of the invention
It is not intended to be limited to embodiment illustrated herein, and is to fit to consistent with the principles and novel features disclosed in this article
Widest range.
Claims (9)
1. a kind of method of rapid survey Harmonious Waves in Power Systems, which comprises the following steps:
Sampled data step is obtained, over-sampling ADC carries out over-sampling to signal, obtains original sampling data；
Quasi-synchro sampling rate step is calculated, steps are as follows:
It carries out down-sampled to the original sampling data of over-sampling ADC output and saves；
Windowing FFT processing is carried out to the data after down-sampled；
Calculate fundamental frequency；
It calculates quasi-synchro sampling rate and saves, quasi-synchro sampling rate Fsnew=Ff× N ' × M, wherein FfFor fundamental frequency, M was
Multiple is sampled, N ' is subsequent plus rectangular window FFT points；
Sampling rate conversion step, steps are as follows:
Calculated quasi-synchro sampling rate carries out interpolation processing to the original sampling data that over-sampling ADC is exported；
Data Jing Guo interpolation processing are carried out down-sampled and are saved；
Plesiochronous FFT processing step, steps are as follows:
Data after sampled rate conversion are carried out plus rectangular window FFT is handled；
Judge whether to need amplitude compensation, if not having to otherwise directly output is as a result, will add rectangular window FFT to obtain after handling
The spectral magnitude of each FFT is multiplied by penalty coefficient, exports result after carrying out amplitude compensation.
2. the method for rapid survey Harmonious Waves in Power Systems according to claim 1, it is characterised in that: over-sampling ADC's adopts
Sample rate Fs=Fnyquist× Oversample, wherein Oversample is the over-sampling multiple of over-sampling ADC, FnyquistTo meet
The sample rate of Nyquist's theorem.
3. the method for rapid survey Harmonious Waves in Power Systems according to claim 1, it is characterised in that: combed using cascade integral
Shape filter carries out down-sampled processing to data.
4. the method for rapid survey Harmonious Waves in Power Systems according to claim 1 or 3, it is characterised in that: calculate plesiochronous
The data after down-sampled are carried out when sample rate plus nuttall window FFT or add Hanning window FFT processing.
5. the method for rapid survey Harmonious Waves in Power Systems according to claim 1, it is characterised in that: use following formula meter
Calculate fundamental frequency: Ff=Δ f × k, the Δ f in formula are the spectral resolution of signal after windowing FFT processing, and k refers to fundamental frequency
On kth FFT spectral line.
6. the method for rapid survey Harmonious Waves in Power Systems according to claim 5, it is characterised in that: fundamental frequencyWherein, Y (x) indicates the amplitude of the xth spectral line of FFT, L=
(f/ Δ f), round (x) expression refers to is rounded x to round, and f is signal frequency.
7. the method for rapid survey Harmonious Waves in Power Systems described according to claim 1 or 2 or 3 or 5 or 6, it is characterised in that: into
It is as follows to the method for data progress interpolation when row sample rate is converted:
A, according to the number of the ratio calculation sampled point of the sample rate of over-sampling ADC and quasi-synchro sampling rate:
Step=Fs/Fsnew；
B, it is exported using the 1st ADC sampled point as the 1st interpolation sampling point, since q=2；
C, q-th of interpolation sampling point position: phase=Step × (q-1)+1 is calculated；
D, the ADC sampled data that q-th of interpolation sampling point position needs: ADC (pre)=floor (phase) and ADC is calculated
(nxt)=ceil (phase), the floor (x) in formula indicate to be rounded x to infinitesimal, and ceil (x) expression takes x to infinity
Whole, pre indicates that re ADC sampled point of pth, nxt indicate n-th xt ADC sampled point；
E, ADC is waited to complete the sampling of n-th xt sampling point, if existing n-th xt sampling point thens follow the steps f；
F, it calculates the amplitude of q-th of interpolation sampling point and exports result, width is calculated according to the sampled data of the obtained ADC of step d
Value: Y (q)=ADC (pre) × (nxt-phase)+ADC (nxt) × (phase-pre), ADC (x) indicate x-th of ADC hits
According to；
G, q=q+1, return step c are enabled.
8. the method for rapid survey Harmonious Waves in Power Systems described according to claim 1 or 2 or 3 or 5 or 6, it is characterised in that: adopt
Interpolation processing is carried out to data with Newton interpolating method.
9. the method for rapid survey Harmonious Waves in Power Systems described according to claim 1 or 2 or 3 or 5 or 6, it is characterised in that: root
Whether it is believed that number occurring to carry out amplitude compensation with interior amplitude attenuation or gain after down-sampled filtering, penalty coefficient is down-sampled
The inverse with interior amplitude attenuation or yield value of filter.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710297370.XA CN107085144B (en) | 2017-04-28 | 2017-04-28 | A kind of method of rapid survey Harmonious Waves in Power Systems |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710297370.XA CN107085144B (en) | 2017-04-28 | 2017-04-28 | A kind of method of rapid survey Harmonious Waves in Power Systems |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107085144A CN107085144A (en) | 2017-08-22 |
| CN107085144B true CN107085144B (en) | 2019-08-20 |
Family
ID=59612273
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710297370.XA Active CN107085144B (en) | 2017-04-28 | 2017-04-28 | A kind of method of rapid survey Harmonious Waves in Power Systems |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107085144B (en) |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107644003A (en) * | 2017-09-29 | 2018-01-30 | 北京智芯微电子科技有限公司 | A kind of computing Co-processor Module of main control chip |
| CN107707259B (en) * | 2017-11-01 | 2020-11-03 | 兰州大学 | Method for sampling and reconstructing analog signal |
| CN108776263B (en) * | 2018-05-02 | 2020-07-28 | 三峡大学 | Harmonic detection method based on high-order Hanning self-convolution window and improved interpolation algorithm |
| CN110988680A (en) * | 2019-11-28 | 2020-04-10 | 西安航天动力试验技术研究所 | Time-frequency processing-based motor rotor fault visualization method |
| CN111562438B (en) * | 2020-04-28 | 2022-04-08 | 重庆金美通信有限责任公司 | Sinusoidal signal frequency estimation method and device based on FFT and phase difference |
| CN114252697B (en) * | 2020-09-23 | 2024-02-13 | 麦克维尔空调制冷(武汉)有限公司 | APF oversampling method, system and storage medium |
| CN113325215B (en) * | 2021-04-28 | 2022-04-26 | 哈尔滨工业大学 | Design method of optical fiber current transformer output controller based on multistage comb filtering |
| CN113533846A (en) * | 2021-07-13 | 2021-10-22 | 国网江苏省电力有限公司营销服务中心 | High-speed synchronous acquisition and calculation method and device for electrical load data |
| CN113608023A (en) * | 2021-08-18 | 2021-11-05 | 国网湖南省电力有限公司 | Sampling frequency self-adaptive harmonic detection method and system |
| CN113866705A (en) * | 2021-08-27 | 2021-12-31 | 威胜集团有限公司 | Harmonic correction method and device, electric energy metering equipment and readable storage medium |
| CN114184839B (en) * | 2021-12-10 | 2024-03-01 | 江苏集萃智能集成电路设计技术研究所有限公司 | ADC frequency spectrum testing method for inhibiting signal source jitter |
| CN116660612A (en) * | 2023-07-31 | 2023-08-29 | 青岛鼎信通讯科技有限公司 | Characteristic current detection method |
| CN116915215B (en) * | 2023-09-12 | 2023-12-08 | 青岛艾诺仪器有限公司 | Implementation method of high sampling rate variable cut-off frequency digital filter |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100587504C (en) * | 2008-01-18 | 2010-02-03 | 杭州万工科技有限公司 | Digital synchronous sampling method |
| KR101144276B1 (en) * | 2010-06-30 | 2012-05-11 | 한국전력공사 | Power quality monitoring system and method thereof |
| CN101915874B (en) * | 2010-07-20 | 2013-01-09 | 北海市深蓝科技发展有限责任公司 | Harmonic wave detection method based on Fourier transformation |
| JP5770454B2 (en) * | 2010-11-24 | 2015-08-26 | 横河電機株式会社 | Harmonic component measuring device |
| CN102128982A (en) * | 2010-12-23 | 2011-07-20 | 东莞市创锐电子技术有限公司 | Harmonic analysis method based on windowed interpolation FFT (Fast Fourier Transform) base frequency tracking technology |
| CN103869162B (en) * | 2014-03-05 | 2017-02-08 | 湖南大学 | Dynamic signal phasor measurement method based on time domain quasi-synchronization |
| CN104459318B (en) * | 2014-11-13 | 2017-06-16 | 广东电网有限责任公司电力科学研究院 | The measuring method and system of electric harmonic |
-
2017
- 2017-04-28 CN CN201710297370.XA patent/CN107085144B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN107085144A (en) | 2017-08-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107085144B (en) | A kind of method of rapid survey Harmonious Waves in Power Systems | |
| CN101915874B (en) | Harmonic wave detection method based on Fourier transformation | |
| CN106483374B (en) | A kind of harmonic wave harmonic detection method based on Nuttall double window whole phase FFT | |
| CN103869162B (en) | Dynamic signal phasor measurement method based on time domain quasi-synchronization | |
| CN103257271B (en) | A kind of micro-capacitance sensor harmonic wave based on STM32F107VCT6 and m-Acetyl chlorophosphonazo pick-up unit and detection method | |
| CN110837001B (en) | Method and device for analyzing harmonic waves and inter-harmonic waves in electric power system | |
| CN106405230B (en) | Frequency measurement method and device | |
| CN106771591B (en) | A kind of method for parameter estimation of Complex Power harmonic wave | |
| CN108918964A (en) | A kind of harmonic analysis method of sparsity enhancing | |
| CN203287435U (en) | A micro electrical network harmonic wave and inter-harmonic wave test apparatus based on an STM32F107VCT6 | |
| CN106018956B (en) | A kind of power system frequency computational methods of adding window spectral line interpolation | |
| CN103969508A (en) | Real-time high-precision power harmonic analysis method and device | |
| CN107942139B (en) | Synchronous sampling method for power harmonic parameter software | |
| Kaiser et al. | Estimation of power systems amplitudes, frequencies, and phase characteristics using energy operators | |
| CN113608023A (en) | Sampling frequency self-adaptive harmonic detection method and system | |
| CN103543331A (en) | Method for calculating harmonics and inter-harmonics of electric signal | |
| CN111505375A (en) | Frequency measurement method and device | |
| CN209102804U (en) | A kind of synchronous phasor measuring device applied to power distribution network | |
| CN109254202A (en) | A kind of synchronous phasor measuring device applied to power distribution network | |
| Carbone et al. | Effective frequency-domain ADC testing | |
| CN114184838A (en) | Power system harmonic detection method, system and medium based on SN mutual convolution window | |
| Abser et al. | CleanEMG: Comparing interpolation strategies for power line interference quantification in surface EMG signals | |
| CN107104674A (en) | Noise power estimation apparatus in electric parameter metering process | |
| Sousa et al. | Non-iterative frequency estimation in the DFT magnitude domain | |
| CN109444957B (en) | Hydrate reservoir earthquake prediction method and processing terminal based on decision tree evaluation |
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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20200205 Address after: 519000 Room 403, 4th floor, building 2, science and Technology Innovation Park, No.1 Gangwan, Jintang Road, hi tech Zone, Zhuhai City, Guangdong Province Patentee after: Zhuhai Taiwei Electronic Co., Ltd Address before: 519000 unit B, room 605, floor 6, main building, No. 10, Keji 1st Road, Gangwan Avenue, Tangjiawan Town, Xiangzhou District, Zhuhai City, Guangdong Province Patentee before: Zhuhai core semiconductor Co., Ltd. |