CN107085144A - A kind of method of quick measurement Harmonious Waves in Power Systems - Google Patents
A kind of method of quick measurement Harmonious Waves in Power Systems Download PDFInfo
- Publication number
- CN107085144A CN107085144A CN201710297370.XA CN201710297370A CN107085144A CN 107085144 A CN107085144 A CN 107085144A CN 201710297370 A CN201710297370 A CN 201710297370A CN 107085144 A CN107085144 A CN 107085144A
- 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.)
- 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
- 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 quick measurement Harmonious Waves in Power Systems, comprises the following steps：Over-sampling ADC carries out over-sampling to signal, obtains original sampling data；Carry out down-sampled to the original sampling data that over-sampling ADC is exported and preserve；Windowing FFT processing is carried out to the data after down-sampled；Calculate fundamental frequency；Calculate quasi-synchro sampling rate and preserve；The original sampling data that the quasi-synchro sampling rate obtained using calculating is exported to over-sampling ADC carries out interpolation processing；Data Jing Guo interpolation processing are carried out down-sampled and preserved；Data are handled using rectangular window FFT；Each FFT obtained after rectangular window FFT processing spectral magnitude is multiplied by penalty coefficient, amplitude compensation, output result is carried out.The inventive method improves the signal to noise ratio of inband signaling using oversampling technique, on the premise of frequency resolution is ensured, greatly reduces computation complexity and improves the degree of accuracy of fundamental frequency calculating, and passes through the occupancy of down-sampled reduction caching.
Description
Technical field
The invention belongs to electricity field, more particularly to a kind of method for calculating analysis mains by harmonics.
Background technology
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, uninterrupted power source, magnetive cord equipment and
The household electrical appliance such as television set can all produce harmonic wave so that grid supply quality is reduced.Therefore, voltage, electric current to power supply grid
Signal carries out frequency analysis, so that carrying out power network compensation and purification becomes more and more important.In order to enter soon and exactly as far as possible
How row power network compensation and purification, the real-time and precise requirements of frequency analysis also more and more higher, rapidly and accurately calculate humorous
Ripple becomes the focus of a research.
FFT (Fast Fourier Transformation, hereinafter referred to as FFT) is the master for calculating harmonic wave
Stream instrument, but be due to truncation effect, there is spectrum leakage and fence effect in FFT, in order to accurately calculate harmonic wave, current common side
Method is as follows：
1st, using windowing FFT, the parameter such as frequency, amplitude and phase in frequency domain interpolation so as to obtain fundamental wave and each harmonic wave
Information.A kind of electric power based on windows and interpolated FFT as disclosed in Publication No. CN101701982A Chinese invention patent application
System harmonicses detection method, windowing FFT computing is carried out using hanning windows, then carries out frequency domain interpolation to calculate each harmonic wave
Accurate frequency, amplitude and phase.This method is when calculating interpolation coefficient, it is necessary to use the top of FFT spectral lines and adjacent
The ratio of the amplitude on secondary peak, therefore FFT spectral resolution sufficiently small must just can guarantee that the two spectral lines do not include other
The amplitude of frequency, during such as using 512 point FFT, sample rate for 6.4KHz, completing complete FFT computing needs sampling 512
Sampling point, the sampling time of consuming is T=N/Fs=512/6400Hz=0.08s=80ms, and source electric power is common are on the market
The response time of wave filter is typically smaller than 20ms, thus while the accuracy of this method is very high, but due to needing one relatively
Small frequency resolution distinguishes fundamental wave and each harmonic wave, causes the sampling time long, it is impossible to meet practical application request.
2nd, being changed by tracking the frequency of fundamental wave realizes plesiochronous, makes the integral multiple that quasi-synchro sampling rate is fundamental frequency,
To suppress spectrum leakage.A kind of electric energy quality harmonic as disclosed in Publication No. CN103969507A Chinese invention patent application
Analysis method, by finding out frequency shift (FS), is then entered using cosine function corresponding with principal wave harmonic wave frequency and SIN function
Row correlation computations, try to achieve amplitude and initial phase angle.This method needs each fundamental wave harmonic being multiplied by corresponding cosine function respectively
And SIN function, because each sampled point will be multiplied by cosine function and SIN function, and each fundamental wave harmonic is corresponding
Cosine function and SIN function it is all different, cause amount of calculation very big, and frequency drift parameter calculating easily by noise
Interference, make fundamental frequency calculate it is inaccurate.
3rd, quasi-synchro sampling is realized using the algorithm of time domain interpolation.Such as Publication No. CN101915874A Chinese invention
A kind of harmonic detecting method based on Fourier transform of patent application publication, the flow 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), is then adopted using wave filter from original after obtaining original sampling data
Sample extracting data fundamental signal, then using the cycle of zero crossing calculating fundamental wave, resampling is confirmed according to the cycle of fundamental wave
Interval, carries out resampling and FFT processing, the cycle of the signal after resampling and fundamental wave are tight according to the interval from original sampling data
Lattice are synchronous, thus avoid common FFT spectrum leakage problem.But this method needs the original for obtaining analog-digital converter over-sampling
Beginning data are deposited in the buffer, because over-sampling rate is very high, cause occupancy caching very big, and using wave filter to original
Also filter result is deposited in the buffer so as to the fortune for carrying out the ensuing extraction primitive period when sampled data is filtered
Calculate, for the embodiment of the patent, when analog-digital converter uses 50*128*40=256KHz sample rate, if modulus turns
The bit wide of parallel operation is 16, then the cache size that preserving the signal of a primitive period needs is 128*40*16=81920bit,
When the sample rate of wave filter is 50*128*2=12.8KHz, the cache size that preserving filtered data needs is 128*2*
16=4096bit, then need 81920+4096=86016bit caching altogether.The area that the caching takes in integrated circuits
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 is very easy to
By noise jamming, cause computation of Period inaccurate.
The content of the invention
It is few and to take caching small it is an object of the invention to provide a kind of amount of calculation, it can quickly analyze in power system each
The method of harmonic parameters.
To achieve these goals, the present invention takes following technical solution：
A kind of method of quick measurement Harmonious Waves in Power Systems, 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, step is as follows：
Carry out down-sampled to the original sampling data that over-sampling ADC is exported and preserve；
Windowing FFT processing is carried out to the data after down-sampled；
Calculate fundamental frequency；
Calculate quasi-synchro sampling rate and preserve, quasi-synchro sampling rate Fsnew=Ff× N ' × M, wherein, FfFor fundamental frequency, M
For over-sampling multiple, N ' is follow-up plus rectangular window FFT points；
Sampling rate conversion step, step is as follows：
The original sampling data that the quasi-synchro sampling rate obtained using calculating is exported to over-sampling ADC carries out interpolation processing；
Data Jing Guo interpolation processing are carried out down-sampled and preserved；
Plesiochronous FFT process steps, step is as follows：
Data after being changed to sampled rate are carried out plus rectangular window FFT processing；
Judge whether to need amplitude compensation, if the direct output result if, obtained after otherwise adding rectangular window FFT processing
To each FFT spectral magnitude be multiplied by penalty coefficient, carry out output result after amplitude compensation.
More specifically, over-sampling ADC sample rate Fs=Fnyquist× Oversample, wherein Oversample were to adopt
Sample multiple, FnyquistTo meet the sample rate of 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 plus nuttall windows FFT or Jia Hanning when calculating quasi-synchro sampling rate
Window FFT processing.
More specifically, fundamental frequency is calculated using below equation：FfAfter Δ f in=Δ f × k, formula is handled for windowing FFT
The spectral resolution of signal, k refers to fundamental frequency on kth bar FFT spectral lines.
More specifically, fundamental frequencyWherein, Y (x) is represented
The amplitude of FFT xth bar spectral line, (f/ Δs f), round (x) expression refers to be rounded L=round to x, and f is signal frequency.
More specifically, the method for carrying out interpolation to data when carrying out sample rate conversion is as follows：
The number of a, the sample rate according to over-sampling ADC and quasi-synchro sampling rate ratio calculation sampled point：Step=Fs/
Ff；
B, export using the 1st ADC sampled point as the 1st interpolation sampling point, since q=2；
C, q-th of interpolation sampling point position of calculating：Phase=Step × (q-1)+1；
D, the ADC sampled datas for calculating q-th of interpolation sampling point position needs：ADC (pre)=floor (phase) and
Floor (x) in ADC (nxt)=ceil (phase), formula is represented to round x to infinitesimal, and ceil (x) is represented to x to infinite
Round greatly, pre represents re ADC sampled point of pth, nxt represents n-th xt ADC sampled point；
E, the sampling for waiting ADC n-th xt sampling point of completion, step f is performed if existing n-th xt sampling point；
F, the amplitude for calculating q-th of interpolation sampling point simultaneously export result, according to the obtained ADC of step d sampled data meter
Calculate amplitude：Y (q)=ADC (pre) × (nxt-phase)+ADC (nxt) × (phase-pre), ADC (x) represent that x-th of ADC is adopted
Sample data；
G, make q=q+1, return to step c.
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
Repay, penalty coefficient is the inverse of the interior amplitude attenuation/yield value of band of desampling fir filter.
From above technical scheme, the present invention improves the signal to noise ratio of inband signaling using oversampling technique, reduction ADC's
Accuracy requirement, reduces system cost, while over-sampling can reduce the noise of interpolation processing；By using down-sampled technology, keep away
The problem of caching takes big caused by over-sampling is exempted from；When calculating fundamental frequency, FFT resolution is improved by down-sampled
Rate, and combine energy barycenter correction method to calculate fundamental frequency, the interference of noise and m-Acetyl chlorophosphonazo is reduced, fundamental frequency is improved
The degree of accuracy of calculating, and reduce computation complexity, it is only necessary to carry out one it is simple down-sampled, then calculate windowing FFT i.e.
Can；Additionally due to sample rate is the integral multiple of fundamental frequency, therefore the FFT of the harmonious wave amplitude of calculating fundamental wave spectral resolution can
Very wide to accomplish, maximum can reach fundamental frequency, so as to reduce FFT points, amount of calculation be reduced, when shortening sampling
Between.Compared with prior art, the present invention overcomes the shortcoming that existing time domain interpolation algorithm is computationally intensive and caching occupancy is big,
Only considerably less amount of calculation, which need to be expended, can be achieved with time domain interpolation quasi-synchro sampling, can accurately be calculated inside the very short time
Go out frequency, amplitude and the first phase of fundamental wave and each harmonic wave, and it is few to take caching.
Brief description of the drawings
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 inventive method.
The embodiment of the present invention is described in more detail below in conjunction with accompanying drawing.
Embodiment
In order to which above and other objects of the present invention, feature and advantage can be become apparent from, the embodiment of the present invention cited below particularly,
It is described below in detail.
Many details are elaborated in the following description to facilitate a thorough understanding of the present invention, still the present invention can be with
It is different from other manner described here 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 following public specific embodiment.
The inventive method is improved on the basis of existing time domain interpolation algorithm, and voltage and electricity are being obtained from power network
The signals such as stream and technology when carrying out analog-to-digital conversion using over-sampling, while being counted using new quasi-synchro sampling rate calculation procedure
Sample rate is calculated, rapidly and accurately to calculate frequency, amplitude and the first equivalent information of fundamental wave and each harmonic wave.With reference to Fig. 2,
The inventive method is described in detail, it is as follows the step of the inventive method：
First, sampled data is obtained, over-sampling ADC carries out over-sampling to signal, obtains original sampling data；
Use signal the sample rate far above signal bandwidth to carry out over-sampling processing, on the one hand can be adopted by improving
The mode of sample multiple improves ADC signal to noise ratio, in theory, often improves 1 times of over-sampling rate, 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, thus ADC required during high power over-sampling sample rate not
Height, can use relatively inexpensive low Precision A/D C, for example, when needing 50 subharmonic of measurement, signal bandwidth is 50Hz × 50=
2.5KHz, meets the sample rate F of Nyquist (nyquist) theoremnyquistFor 2 times of signal bandwidth (highest frequency), i.e.,
When 2.5KHz × 2=5KHz, 64 times of over-sampling, over-sampling ADC sample rate Fs=Fnyquist× Oversample=5KHz × 64
=320KHz, Oversample are over-sampling multiple, FnyquistTo meet the sample rate of Nyquist's theorem, now band in letter
Make 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 during follow-up interpolation processing, so as to reduce the calculation error and meter of follow-up interpolation algorithm
Calculation amount；
Then, quasi-synchro sampling rate is calculated, the calculating of quasi-synchro sampling rate comprises the following steps：
Carry out down-sampled to the original sampling data that over-sampling ADC is exported and preserve；In order to reduce amount of calculation and caching
Take, and improve the accuracy rate of fundamental frequency calculating, it is necessary 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；Because CIC calculating is not related to multiplication, so the integrated circuit area that it takes is very
Few, the down-sampled multiple of cic filter is 2048 in the present embodiment, and exponent number is 8, down-sampled sample rate=320KHz/2048
=156.25Hz, due to only needing to calculate the frequency of fundamental wave, therefore the CIC interior amplitude attenuation of band can be without processing；
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 are carried out plus nuttall window FFT processing, and the coefficient of nuttall window functions 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
It is Δ f for the spectral resolution that the data after down-sampled carried out with signal after the points of windowing FFT, plus nuttall windows FFT processing
=Fs/ N=156.25/32=4.8828125Hz, memory is preserved to (signal) points N needed for windowing FFT；Except
Windowing FFT processing can be carried out with nuttall window functions outer, it would however also be possible to employ other window functions such as Hanning window, rectangular window are carried out
Windowing FFT processing, in addition, nuttall window functions can also be using forms such as other three ranks, five ranks；
Calculate fundamental frequency；The present embodiment is using conventional fixation FFT spectrum bearing calibration --- and " power enhanced " is counted
Calculate fundamental frequency, i.e. fundamental frequencyWherein, Δ f is adding window
The spectral resolution of signal after FFT processing, the Δ f=4.8828125Hz, k of the present embodiment refer to fundamental frequency in kth bar
On FFT spectral lines, k can be decimal, and Y (x) represents the amplitude of FFT xth bar spectral line, L=round (f/ Δs f), round (x) tables
Show that finger is rounded to x, the original signal frequency f of the present embodiment is 50Hz, then L=round (50/4.8828125)=10；
Calculate quasi-synchro sampling rate and preserve；In order to eliminate spectrum leakage, quasi-synchro sampling rate is the integer of fundamental frequency
Times, due to present invention uses over-sampling ADC, therefore quasi-synchro sampling rate also needs to fundamental frequency being multiplied by over-sampling multiple M,
Final quasi-synchro sampling rate is Fsnew=Ff× N ' × M, N ' are the follow-up points that data after interpolation add with rectangular window FFT,
The N ' of the present embodiment is 128, and over-sampling multiple M is 64, then Fsnew=Ff×128×64；
Then, sample rate is changed, comprised the following steps：
The original sampling data that the quasi-synchro sampling rate obtained using calculating is exported to over-sampling ADC carries out interpolation processing；
Interpolation can be calculated using linear interpolation method or Newton interpolating method, the present embodiment calculates interpolation using linear interpolation method, including following
Step：
The number of a, the sample rate according to over-sampling ADC and quasi-synchro sampling rate ratio calculation sampled point：Step=Fs/
Ff；
B, export using the 1st ADC sampled point as the 1st interpolation sampling point, since q=2；
C, q-th of interpolation sampling point position of calculating：Phase=Step × (q-1)+1；
D, the ADC sampled datas for calculating q-th of interpolation sampling point position needs：ADC (pre)=floor (phase) and
Floor (x) in ADC (nxt)=ceil (phase), formula is represented to round x to infinitesimal, and ceil (x) is represented to x to infinite
Round greatly, such as floor (2.8)=2, ceil (2.2)=3, pre represents re ADC sampled point of pth, and nxt represents n-th xt ADC
Sampled point；
E, the sampling for waiting ADC n-th xt sampling point of completion, step f is performed if existing n-th xt sampling point；
F, the amplitude for calculating q-th of interpolation sampling point simultaneously export result to CIC, carry out down-sampled processing, are obtained according to step d
The ADC arrived sampled data calculates amplitude：Y (q)=ADC (pre) × (nxt-phase)+ADC (nxt) × (phase-pre),
ADC (x) represents x-th of ADC sampled data；
G, make q=q+1, return to step c；Because Harmonics Calculation is real-time, it is necessary to ceaselessly refresh, therefore interpolation arithmetic
It is also required to constantly update, until stopping Harmonics Calculation；
It is down-sampled to the data progress Jing Guo interpolation processing, at using cascade integral comb filter (CIC) to interpolation
Data after reason carry out down-sampled, and the CIC of the present embodiment down-sampled multiple is 64, and exponent number is 6, it is down-sampled after sample rate
FsFFT=Fsnew/ 64, the result after down-sampled processing is preserved in the buffer, and the size of caching is follow-up plus rectangular window FFT points
N ', the present embodiment rectangular window FFT points N ' it is 128 points, when data bit width is 16bit, then shared caching is 128*
16=2048bit, compared to the prior art, caching take very small；
Finally, data are carried out with plesiochronous FFT processing, and output result, is comprised 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 spectral resolution of signal is after the FFT of usable random length, plus rectangular window FFT processing：By
In the sampling time usedAs Δ foutFor 49Hz, T=1/49Hz=20.408ms, it can be seen that number
According to sampling time greatly reduced relative to prior art；
Amplitude compensation, output result；Judged whether to need to carry out amplitude compensation according to down-sampled method, if sampled
Down-sampled technology not with interior amplitude attenuation or gain then without compensation, the present embodiment using CIC carry out it is down-sampled, due to
CIC technologies are present with interior amplitude attenuation, it is therefore desirable to be multiplied by each FFT obtained after rectangular window FFT processing spectral magnitude
Penalty coefficient, carries out amplitude compensation, and penalty coefficient is the inverse of the interior amplitude attenuation/yield value of band of desampling fir filter, can be led to
Cross simulation calculation signal and penalty coefficient is determined by the attenuation/gain value of wave filter (CIC) amplitude afterwards, for example, signal exists
It is 1/10 to be changed into the interior amplitude attenuation value of band at original 1/10, i.e. 100Hz by the amplitude attenuation after CIC at 100Hz frequencies,
Then penalty coefficient is 10, and signal is changed at original 1/20, i.e. 200Hz at 200Hz frequencies by the amplitude attenuation after CIC
It is 1/20 with interior amplitude attenuation value, then the penalty coefficient at 200Hz is 20, each not phase of attenuation/gain of each frequency of wave filter
Together, and different filter constructions and down-sampled multiple penalty coefficient it is different, amplitude change is also different, as a rule
CIC exponent number is more, decays more, therefore penalty coefficient is it is determined that after desampling fir filter (CIC) structure and down-sampled multiple
Need to determine penalty coefficient by emulating.Frequency spectrum interested can be only compensated, operand is reduced.
The inventive method improves the signal to noise ratio of inband signaling using oversampling technique, reduces ADC accuracy requirement, reduction system
System cost, while the operand of sample rate conversion can be reduced, adds the FFT (the present embodiment counted on a small quantity again by first down-sampled
In be 32 points) calculate fundamental frequency, on the premise of frequency resolution is ensured, greatly reduce computation complexity, and improve
The degree of accuracy that fundamental frequency is calculated, while reduce the occupancy of caching by down-sampled technology.Because fundamental frequency change is slow
Slowly, the present invention first calculates fundamental frequency when calculating quasi-synchro sampling rate with windowing FFT, without synchronized sampling, further according to fundamental wave
Frequency calculates quasi-synchro sampling rate, then can quickly calculate harmonic wave with plesiochronous FFT method.
The foregoing description of the disclosed embodiments, enables professional and technical personnel in the field to realize or using the present invention.
A variety of modifications to these embodiments will be apparent for 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, it is of the invention
Embodiment illustrated herein is not intended to be limited to, and is to fit to consistent with principles disclosed herein and features of novelty
Widest range.
Claims (9)
1. a kind of method of quick measurement Harmonious Waves in Power Systems, it is characterised in that comprise 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, step is as follows：
Carry out down-sampled to the original sampling data that over-sampling ADC is exported and preserve；
Windowing FFT processing is carried out to the data after down-sampled；
Calculate fundamental frequency；
Calculate quasi-synchro sampling rate and preserve, quasi-synchro sampling rate Fsnew=Ff× N ' × M, wherein, FfFor fundamental frequency, M was
Sampling multiple, N ' is follow-up plus rectangular window FFT points；
Sampling rate conversion step, step is as follows：
The original sampling data that the quasi-synchro sampling rate obtained using calculating is exported to over-sampling ADC carries out interpolation processing；
Data Jing Guo interpolation processing are carried out down-sampled and preserved；
Plesiochronous FFT process steps, step is as follows：
Data after being changed to sampled rate are carried out plus rectangular window FFT processing；
Judge whether to need amplitude compensation, if the direct output result if, otherwise will add what rectangular window FFT was obtained after handling
Each FFT spectral magnitude is multiplied by penalty coefficient, carries out output result after amplitude compensation.
2. the method for quick measurement 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 are over-sampling multiple, FnyquistTo meet Nyquist's theorem
Sample rate.
3. the method for quick measurement Harmonious Waves in Power Systems according to claim 1, it is characterised in that：Combed using cascade integral
Shape wave filter carries out down-sampled processing to data.
4. the method for the quick measurement 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 during sample rate plus nuttall windows FFT or add Hanning window FFT processing.
5. the method for quick measurement Harmonious Waves in Power Systems according to claim 1, it is characterised in that：Using below equation meter
Calculate fundamental frequency：FfThe spectral resolution of signal after Δ f in=Δ f × k, formula is handled for windowing FFT, k refers to fundamental frequency
On kth bar FFT spectral lines.
6. the method for quick measurement Harmonious Waves in Power Systems according to claim 5, it is characterised in that：Fundamental frequencyWherein, Y (x) represents the amplitude of FFT xth bar spectral line, L=
(f/ Δs f), round (x) expression refers to be rounded round to x, and f is signal frequency.
7. the method for the quick measurement Harmonious Waves in Power Systems according to claim 1 or 2 or 3 or 5 or 6, it is characterised in that：Enter
The method for carrying out interpolation to data when row sample rate is changed is as follows：
The number of a, the sample rate according to over-sampling ADC and quasi-synchro sampling rate ratio calculation sampled point：Step=Fs/Ff；
B, export using the 1st ADC sampled point as the 1st interpolation sampling point, since q=2；
C, q-th of interpolation sampling point position of calculating：Phase=Step × (q-1)+1；
D, the ADC sampled datas for calculating q-th of interpolation sampling point position needs：ADC (pre)=floor (phase) and ADC
(nxt) floor (x) in=ceil (phase), formula represents to round x to infinitesimal, and ceil (x) represents to take x to infinity
Whole, pre represents re ADC sampled point of pth, and nxt represents n-th xt ADC sampled point；
E, the sampling for waiting ADC n-th xt sampling point of completion, step f is performed if existing n-th xt sampling point；
F, the amplitude for calculating q-th of interpolation sampling point simultaneously export result, and width is calculated according to the obtained ADC of step d sampled data
Value：Y (q)=ADC (pre) × (nxt-phase)+ADC (nxt) × (phase-pre), ADC (x) represent x-th of ADC hits
According to；
G, make q=q+1, return to step c.
8. the method for the quick measurement Harmonious Waves in Power Systems 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 the quick measurement Harmonious Waves in Power Systems 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 of the interior amplitude attenuation/yield value of band of wave 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 true CN107085144A (en) | 2017-08-22 |
| CN107085144B 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) |
Cited By (14)
| 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 |
| CN107707259A (en) * | 2017-11-01 | 2018-02-16 | 兰州大学 | A kind of method of analog signal sampling and reconstruct |
| CN108776263A (en) * | 2018-05-02 | 2018-11-09 | 三峡大学 | Harmonic detecting method based on the peaceful self-convolution window of the high-order Chinese and improvement interpolation algorithm |
| CN110988680A (en) * | 2019-11-28 | 2020-04-10 | 西安航天动力试验技术研究所 | Time-frequency processing-based motor rotor fault visualization method |
| CN111562438A (en) * | 2020-04-28 | 2020-08-21 | 重庆金美通信有限责任公司 | Sinusoidal signal frequency estimation method and device based on FFT and phase difference |
| CN113325215A (en) * | 2021-04-28 | 2021-08-31 | 哈尔滨工业大学 | 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 |
| CN114047392A (en) * | 2021-09-26 | 2022-02-15 | 南京国电南自电网自动化有限公司 | Method for measuring amplitude and phase angle of electric signal of low-voltage protection device |
| CN114184839A (en) * | 2021-12-10 | 2022-03-15 | 江苏集萃智能集成电路设计技术研究所有限公司 | ADC frequency spectrum testing method for inhibiting signal source jitter |
| CN114252697A (en) * | 2020-09-23 | 2022-03-29 | 麦克维尔空调制冷(武汉)有限公司 | APF oversampling method, system and storage medium |
| CN116660612A (en) * | 2023-07-31 | 2023-08-29 | 青岛鼎信通讯科技有限公司 | Characteristic current detection method |
| CN116915215A (en) * | 2023-09-12 | 2023-10-20 | 青岛艾诺仪器有限公司 | Implementation method of high sampling rate variable cut-off frequency digital filter |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101216511A (en) * | 2008-01-18 | 2008-07-09 | 杭州万工科技有限公司 | Digital synchronous sampling method |
| CN101915874A (en) * | 2010-07-20 | 2010-12-15 | 北海市深蓝科技发展有限责任公司 | Harmonic wave detection method based on Fourier transformation |
| CN102128982A (en) * | 2010-12-23 | 2011-07-20 | 东莞市创锐电子技术有限公司 | Harmonic analysis method based on windowed interpolation FFT (Fast Fourier Transform) base frequency tracking technology |
| WO2012002615A1 (en) * | 2010-06-30 | 2012-01-05 | 한국전력공사 | Power quality monitoring system and method thereof |
| JP2012112762A (en) * | 2010-11-24 | 2012-06-14 | Yokogawa Electric Corp | Higher harmonic wave component measuring device |
| CN103869162A (en) * | 2014-03-05 | 2014-06-18 | 湖南大学 | Dynamic signal phasor measurement method based on time domain quasi-synchronization |
| CN104459318A (en) * | 2014-11-13 | 2015-03-25 | 广东电网有限责任公司电力科学研究院 | Power harmonic measuring method and system |
-
2017
- 2017-04-28 CN CN201710297370.XA patent/CN107085144B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101216511A (en) * | 2008-01-18 | 2008-07-09 | 杭州万工科技有限公司 | Digital synchronous sampling method |
| WO2012002615A1 (en) * | 2010-06-30 | 2012-01-05 | 한국전력공사 | Power quality monitoring system and method thereof |
| CN101915874A (en) * | 2010-07-20 | 2010-12-15 | 北海市深蓝科技发展有限责任公司 | Harmonic wave detection method based on Fourier transformation |
| JP2012112762A (en) * | 2010-11-24 | 2012-06-14 | Yokogawa Electric Corp | Higher harmonic wave 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 |
| CN103869162A (en) * | 2014-03-05 | 2014-06-18 | 湖南大学 | Dynamic signal phasor measurement method based on time domain quasi-synchronization |
| CN104459318A (en) * | 2014-11-13 | 2015-03-25 | 广东电网有限责任公司电力科学研究院 | Power harmonic measuring method and system |
Non-Patent Citations (1)
| Title |
|---|
| 王彭 等: "基于时域准同步的谐波和间谐波检测算法", 《仪器仪表学报》 * |
Cited By (20)
| 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 |
| CN107707259A (en) * | 2017-11-01 | 2018-02-16 | 兰州大学 | A kind of method of analog signal sampling and reconstruct |
| CN107707259B (en) * | 2017-11-01 | 2020-11-03 | 兰州大学 | Method for sampling and reconstructing analog signal |
| CN108776263A (en) * | 2018-05-02 | 2018-11-09 | 三峡大学 | Harmonic detecting method based on the peaceful self-convolution window of the high-order Chinese and improvement interpolation algorithm |
| CN110988680A (en) * | 2019-11-28 | 2020-04-10 | 西安航天动力试验技术研究所 | Time-frequency processing-based motor rotor fault visualization method |
| CN111562438A (en) * | 2020-04-28 | 2020-08-21 | 重庆金美通信有限责任公司 | Sinusoidal signal frequency estimation method and device based on FFT and phase difference |
| CN111562438B (en) * | 2020-04-28 | 2022-04-08 | 重庆金美通信有限责任公司 | Sinusoidal signal frequency estimation method and device based on FFT and phase difference |
| CN114252697A (en) * | 2020-09-23 | 2022-03-29 | 麦克维尔空调制冷(武汉)有限公司 | APF oversampling method, system and storage medium |
| CN114252697B (en) * | 2020-09-23 | 2024-02-13 | 麦克维尔空调制冷(武汉)有限公司 | APF oversampling method, system and storage medium |
| CN113325215A (en) * | 2021-04-28 | 2021-08-31 | 哈尔滨工业大学 | Design method of optical fiber current transformer output controller based on multistage comb filtering |
| 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 |
| CN114047392A (en) * | 2021-09-26 | 2022-02-15 | 南京国电南自电网自动化有限公司 | Method for measuring amplitude and phase angle of electric signal of low-voltage protection device |
| CN114184839A (en) * | 2021-12-10 | 2022-03-15 | 江苏集萃智能集成电路设计技术研究所有限公司 | ADC frequency spectrum testing method for inhibiting signal source jitter |
| 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 |
| CN116915215A (en) * | 2023-09-12 | 2023-10-20 | 青岛艾诺仪器有限公司 | Implementation method of high sampling rate variable cut-off frequency digital filter |
| CN116915215B (en) * | 2023-09-12 | 2023-12-08 | 青岛艾诺仪器有限公司 | Implementation method of high sampling rate variable cut-off frequency digital filter |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107085144B (en) | 2019-08-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107085144B (en) | A kind of method of rapid survey Harmonious Waves in Power Systems | |
| Lim et al. | A new algorithm for two-dimensional maximum entropy power spectrum estimation | |
| 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 | |
| CN101915874B (en) | Harmonic wave detection method based on Fourier transformation | |
| CN103308804B (en) | Based on quick K-S converting electric power quality disturbance signal time and frequency parameter extracting method | |
| CA2929382C (en) | Method of removing decaying dc component from power system fault signal | |
| CN110837001B (en) | Method and device for analyzing harmonic waves and inter-harmonic waves in electric power system | |
| CN104391178B (en) | A kind of time shift phase difference stable state harmonic signal bearing calibration based on Nuttall windows | |
| CN103454497A (en) | Phase difference measuring method based on improved windowing discrete Fourier transform | |
| CN103207319A (en) | Harmonic wave measurement method of electricity signal of digital substation under non-synchronous sampling condition | |
| CN113077806B (en) | Audio processing method and device, model training method and device, medium and equipment | |
| CN108776263B (en) | Harmonic detection method based on high-order Hanning self-convolution window and improved interpolation algorithm | |
| CN111222088A (en) | Improved method for estimating weighted power harmonic amplitude of flat-top self-convolution window | |
| CN103969508A (en) | Real-time high-precision power harmonic analysis method and device | |
| CN103543331B (en) | A kind of method calculating electric signal harmonic wave and m-Acetyl chlorophosphonazo | |
| CN103795473B (en) | Method and system for eliminating power-frequency interference | |
| CN111505375A (en) | Frequency measurement method and device | |
| CN110269642B (en) | Doppler heart rate estimation method based on fractional Fourier transform and wavelet transform | |
| Ben et al. | MEMS gyro denoising based on second generation wavelet transform | |
| Vainio et al. | A class of predictive analog filters for sensor signal processing and control instrumentation | |
| Lu | Application of DFT filter bank to power frequency harmonic measurement | |
| CN107104674A (en) | Noise power estimation apparatus in electric parameter metering process | |
| CN110808929A (en) | Real-complex conversion type signal-to-noise ratio estimation algorithm of subtraction strategy | |
| CN103078661A (en) | Spread spectrum system interference inhibition method based on iteration threshold |
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. |