CN104062500A  Signal harmonic analysis method and system based on Hamming product window  Google Patents
Signal harmonic analysis method and system based on Hamming product window Download PDFInfo
 Publication number
 CN104062500A CN104062500A CN201410321280.6A CN201410321280A CN104062500A CN 104062500 A CN104062500 A CN 104062500A CN 201410321280 A CN201410321280 A CN 201410321280A CN 104062500 A CN104062500 A CN 104062500A
 Authority
 CN
 China
 Prior art keywords
 signal
 hamming
 harmonic
 window
 product window
 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
 238000004458 analytical method Methods 0.000 title abstract description 15
 238000006243 chemical reaction Methods 0.000 claims description 42
 230000003595 spectral Effects 0.000 claims description 37
 238000001228 spectrum Methods 0.000 claims description 23
 230000000875 corresponding Effects 0.000 claims description 11
 238000005070 sampling Methods 0.000 abstract description 10
 238000000034 method Methods 0.000 description 5
 230000000694 effects Effects 0.000 description 4
 238000010183 spectrum analysis Methods 0.000 description 3
 230000001360 synchronised Effects 0.000 description 3
 238000001514 detection method Methods 0.000 description 2
 238000010586 diagram Methods 0.000 description 2
 230000002401 inhibitory effect Effects 0.000 description 2
 239000000203 mixture Substances 0.000 description 2
 238000006467 substitution reaction Methods 0.000 description 2
 235000004418 Durio kutejensis Nutrition 0.000 description 1
 238000005259 measurement Methods 0.000 description 1
Abstract
The invention discloses a signal harmonic analysis method and system based on a Hamming product window. The method includes the steps that discrete signals are windowed through the Hamming product window, and FFT is conducted on the windowed signals and the Hamming product window; actual fundamental frequency and the amplitude and the phases of subharmonics are acquired according to a signal spectrogram obtained after FFT and the Hamming product window obtained after FFT. The Hamming product window is combined with FFT, and harmonic analysis with extremely high accuracy can be achieved in an asynchronous sampling state; meanwhile, the signal harmonic analysis method and system further have the advantages of being small in calculated amount, short in consumed time, high in practicability and the like.
Description
Technical field
The present invention relates to a kind of Harmonious Waves in Power Systems frequency spectrum analysis method, be specially a kind of signal harmonic analytical approach and system based on Hamming product window.
Background technology
Nonlinearload in electric system can cause the deterioration of the quality of power supply.In actual signal processing procedure, signal is processed degree of accuracy can be subject to the impact that the quality of power supply changes, the nonsynchronous sampling that the fluctuation of fundamental frequency causes will cause spectral leakage and hurdle effect, make electric system detect harmonic amplitude, phase angle and frequency and all have larger error, can not meet accuracy requirement.Hurdle effect can adopt method of interpolation to solve, but the measuring error causing because of spectrum leakage phenomenon still exists.To carrying out the error that direct computation of DFT leaf analysis can reduce spectral leakage after signal windowing, but process is complicated, calculated amount is large again.
For addressing the above problem, Chinese scholars has been done corresponding research, has proposed high precision harmonic detecting method based on FFT, based on Nuttall (Nuttall) window and 5 Lai FuVincents (RifeVincent) (I) the Electric Power Harmonic Analysis method, rectangle convolution window etc. of window double spectral line interpolation FFT.Yet because the frequency spectrum sidelobe performance of above window function is not ideal enough, limited to spectral leakage inhibiting effect, make frequency analysis error larger, be difficult to realize highprecision frequency measurement and frequency analysis, and computation process is complicated, length consuming time, is difficult to meet the requirement of harmonic analysis in power system realtime and accuracy.
Summary of the invention
The problem existing for prior art, the invention provides a kind of signal harmonic frequency spectrum analysis method and system based on Hamming product window that can effectively suppress spectral leakage, the method calculated amount is little, consuming time short, degree of accuracy is high, can realize signal in real time and Accurate Analysis.
The present invention is based on Hamming product window (Hamming product window) and realize, take Hamming window function as basic window function, basic window function is multiplied each other and obtains Hamming product window w
_{m}(n):
In formula (1), w
_{i}(n) represent i basic window function, i=1,2 ..., p.
The signal harmonic frequency spectrum analysis method based on Hamming product window that the present invention proposes, comprises step:
Step 1, carries out analog to digital conversion to power system signal x (t) and obtains discretize signal x (n), and described power system signal x (t) is voltage signal or current signal;
Step 2, adopts Hamming product window to discretize signal x (n) windowing, and signal after windowing and Hamming product window are carried out to FFT conversion, obtains signal spectrum figure and Hamming product window after FFT conversion; Described Hamming product window is the acquisition of multiplying each other of p Hamming window, and p value is set according to actual conditions;
Step 3, obtains actual fundamental frequency according to signal spectrum figure and Hamming product window after FFT conversion;
Step 4, according to amplitude, the phase place of the signal spectrum figure after FFT conversion and actual fundamental frequency acquisition firstharmonic and each harmonic.
Hamming window described in step 2 is cosine Hamming window.
Step 3 is specially:
By put the discrete spectral line sequence number in the nearest left and right sides from theoretical peak, be designated as respectively k
_{1}and k
_{2}, discrete spectral line k
_{1}and k
_{2}corresponding amplitude is respectively y
_{1}and y
_{2}, according to
obtain the discrete spectral line sequence number k of actual firstharmonic
_{0}thereby, obtain fundamental frequency; Wherein, W () is the Hamming product window after FFT conversion.
Step 4 is specially:
[(the i*f of signal spectrum figure after FFT conversion
_{0}a), (i*f
_{0}+ a)] look for the discrete spectral line k that i subharmonic is corresponding in frequency range
_{i1}and k
_{i2}, i is overtone order; The harmonic wave maximum deviation that a allows according to electric system arranges; Based on discrete spectral line k
_{i1}and k
_{i2}corresponding amplitude y
_{i1}and y
_{i2}, adopt amplitude correction formula and phase place correction formula to obtain amplitude and the phase angle of i subharmonic.
The signal harmonic analytic system based on Hamming product window that the present invention proposes, comprising:
Analogtodigital conversion module, is used for that power system signal x (t) is carried out to analog to digital conversion and obtains discretize signal x (n), and described power system signal x (t) is voltage signal or current signal;
Windowing and FFT change module, are used for adopting Hamming product window to discretize signal x (n) windowing, and signal after windowing and Hamming product window are carried out to FFT conversion, obtain signal spectrum figure and Hamming product window after FFT conversion; Described Hamming product window is the acquisition of multiplying each other of p Hamming window, and p value is set according to actual conditions;
Actual fundamental frequency obtains module, is used for obtaining actual fundamental frequency according to signal spectrum figure and Hamming product window after FFT conversion;
Harmonic amplitude and harmonic phase obtain module, are used for according to amplitude, the phase place of the signal spectrum figure after FFT conversion and actual fundamental frequency acquisition firstharmonic and each harmonic.
Abovementioned harmonic amplitude and harmonic phase obtain module and further comprise:
The discrete spectral line that harmonic wave is corresponding is looked for module, is used for the [(i*f of the signal spectrum figure after FFT conversion
_{0}a), (i*f
_{0}+ a)] look for the discrete spectral line k that i subharmonic is corresponding in frequency range
_{i1}and k
_{i2}, i is overtone order; The harmonic wave maximum deviation that a allows according to electric system arranges;
Amplitude and phase solution module, be used for based on discrete spectral line k
_{i1}and k
_{i2}corresponding amplitude y
_{i1}and y
_{i2}, adopt amplitude correction formula and phase place correction formula to obtain amplitude and the phase angle of i subharmonic.
Compared with prior art, the present invention has the following advantages and beneficial effect:
(1) based on Hamming product window, realize, under nonsynchronous sampling, the inhibition of spectral leakage is much better than to classical window function.
(2) by Hammng product window and FFT conversion combination, under nonsynchronous sampling state, can realize the frequency analysis of very high degree of precision.
(3) calculated amount is little, and consuming time short, practicality is high, can guarantee realtime and the degree of accuracy of frequency analysis simultaneously.
Accompanying drawing explanation
Fig. 1 is the time domain specification of different Hamming product windows;
Fig. 2 is the frequency domain characteristic of different Hamming product windows;
Fig. 3 is the concrete implementing procedure figure of the inventive method;
Fig. 4 is the spectrogram of windowing signal after FFT conversion;
Fig. 5 is the amplitude relative error comparison diagram of the inventive method and classic method;
Fig. 6 is the phase place relative error comparison diagram of the inventive method and classic method.
Embodiment
The present invention is based on Hamming product window and realize the harmonic wave detection of power system signal, Hamming product window w
_{hammingm}(n) by some Hamming window w
_{hamming}(n) acquisition of multiplying each other:
In formula (1), p is the Hamming window number multiplying each other.In concrete enforcement, p value is selected according to actual conditions and experience, and generally speaking, p value is larger, and the side lobe peak height of Hamming product window is lower, but main lobe is wider.And suitable window function requires, main lobe is narrow, secondary lobe is low.
In this concrete enforcement, adopt cosine Hamming window w
_{hamming}(n), it has very little side lobe peak:
In formula (2), N is window function length, 0≤n≤N1.
When basic window function type and number change in Hamming product window, can obtain different Hamming product windows.Fig. 1 and Fig. 2 have shown respectively time domain specification and the frequency domain characteristic of different Hamming product windows, and the characteristic parameter of each Hamming product window is in Table 1.
The spectral characteristic parameter of table 1Hamming product window
In table 1, Hamming represents cosine Hamming window, Hamming
^{2}, Hamming
^{3}, Hamming
^{4}, Hamming
^{5}represent respectively the Hamming product window that 2,3,4,5 cosine Hamming windows multiply each other and obtain.
The main lobe width (MB) of Hamming product window, side lobe attenuation speed (SRR) have following approximation relation with parameter p:
In formula (3): w
_{hammingm}(n)
_{mB}the main lobe width that represents Hamming product window, w
_{hamming}(n)
_{mB}the main lobe width that represents Hamming window, w
_{hammingm}(n)
_{sSR}the side lobe attenuation speed that represents Hamming product window.
In Table 1, when p increases, the side lobe peak height (PSL) of Hamming product window diminishes gradually, and along with basic window function amount increases, Hamming product window suppresses spectral leakage ability can be strengthened.
Fig. 3 is the particular flow sheet of the inventive method, comprises step:
Step 1, sampling power system signal x (t) also obtains discretize signal x (n).
With fixed sampling frequency 4kHz sampling power system signal x (t), signal x (t) can be voltage signal or current signal; Signal x (t) is obtained to discretize signal x (n) through analog to digital conversion, send discretize signal x (n) to merge cells after merging combination, transfer to digital information network.
Discretize signal x (n) is as follows:
x(n)＝Asin(2πnf
_{0}/f
_{s}+θ) (4)
In formula (4), n is sampled point numbering; A and θ are respectively amplitude and the phase place of harmonic wave; f
_{0}for signal x (t) crest frequency, fundamental frequency to be asked; f
_{s}for sample frequency, be 4kHz herein.
Step 2, adopts Hamming product window to carry out windowing to discretize signal x (n), and signal after windowing is carried out to the signal spectrum figure after FFT conversion acquisition conversion.
In electric energy quality monitor, discretize signal x (n) is carried out to windowing: adopt Hamming product window w
_{hammingm}(n) block discretize signal x (n), signal x after acquisition windowing
_{m}(n)=x (n) * w
_{hammingm}(n), n=0,1 ... N1, N is window length, carries out according to actual needs value, N=1000 in this concrete enforcement.
Ignore the peak secondary lobe impact frequently of negative frequency peak point place, signal after windowing is carried out to FFT conversion, after conversion, signal X (k Δ f) is:
In formula (5), j represents imaginary part; A and θ are respectively amplitude and the phase place of nth harmonic; W () is the Hamming product window after FFT conversion; Δ f is discrete frequency interval, Δ f=f
_{s}/ N, N is window length; K is independent variable, represents discrete spectral line sequence number; Crest frequency f
_{0}=k
_{0}Δ f generally there will not be on the Frequency point of discrete spectral line, k
_{0}discrete spectral line sequence number for firstharmonic.
Step 3, obtains actual fundamental frequency according to the signal spectrum figure after Fast Fourier Transform (FFT) (FFT) and Hamming product window value.
By put the discrete spectral line sequence number in the nearest left and right sides from theoretical peak, be designated as respectively k
_{1}and k
_{2}, spectral line k
_{1}and k
_{2}corresponding amplitude is respectively y
_{1}and y
_{2}, from formula (5):
Because window function W () is known, according to formula (6), can obtain unknown quantity k
_{0}thereby, obtain fundamental frequency.
Step 4, obtains amplitude and the phase angle of each harmonic based on fundamental frequency.
Because of 0≤k
_{0}k
_{1}≤ 1, for simplifying, calculate defined parameters
α=k
_{0}k
_{1}0.5, wherein, α ∈ [0.5,0.5], when window length N is larger, formula (6) can be reduced to:
α＝g
^{1}(β) (7)
Adopt polynomial approach algorithm to solve formula (7), by control, approach number of times, can control flexibly approximation accuracy.The weight of supposing spectral line is directly proportional to its amplitude, approaches in the approximate treatment formula of obtaining with polynomial expression,, containing odd term, does not obtain amplitude correction formula (8) and phase place correction formula (9):
A＝N
^{1}(y
_{1}+y
_{2})(b
_{0}+b
_{2}β
^{2}+…+b
_{2i}β
^{2i}) (8)
θ＝arg[X(k
_{i}·Δf)]+π/2arg[W(2π(k
_{i}k
_{0})/N)] (9)
In formula (8)～(9), i represents overtone order.
Therefore, adopt amplitude correction formula (8) and phase place correction formula (9) can directly obtain amplitude and the phase place of each harmonic.
The amplitude of each harmonic and the acquisition methods of phase angle are specific as follows:
[(the i*f of signal spectrum figure after FFT conversion
_{0}a), (i*f
_{0}+ a)] look for the discrete spectral line k that i subharmonic is corresponding in frequency range
_{i1}and k
_{i2}, i is overtone order, i=1,2 ..., n, during i=1, be firstharmonic; The harmonic wave maximum deviation that a allows according to electric system arranges, and this embodiment makes a=5; By discrete spectral line k
_{1}and k
_{2}corresponding amplitude y
_{1}and y
_{2}substitution amplitude correction formula (8) and phase place correction formula (9) obtain amplitude and the phase angle of i subharmonic.
To further illustrate technical solution of the present invention and technique effect by embodiment below.
Power system signal x (t) with fixed sampling frequency 4kHz sampling containing 11 subharmonic, fundamental frequency f
_{1}=50.1Hz, each harmonic amplitude and phase place are in Table 2.
The amplitude of table 2 each harmonic and phase place
Abovementioned power system signal x (t) is carried out to harmonic wave detection, and step is as follows:
Step 1, carries out analog to digital conversion to original signal x (t) and obtains discretize signal x (n), adopts Hamming product window to carry out windowing to discretize signal x (n) and obtains truncated signal x
_{m}(n)=x (n) * w
_{hammingm}(n), wherein, n=0,1 ..., N1, Hamming product window length N=1000.
Truncated signal is carried out to FFT conversion, and after conversion, the spectral characteristic of signal is shown in Fig. 4.
Step 2 is looked for the discrete spectral line k nearest from theoretic frequency point both sides in [45Hz, 55Hz] frequency range
_{1}and k
_{2}, spectral line k
_{1}and k
_{2}corresponding amplitude is designated as y
_{1}and y
_{2}, wherein:
Obtain factor alpha and β:
α＝0.1314β+0.2865β
^{3}+0.4671β
^{5}(12)
Step 3, obtains amplitude correction formula (13) and the phase angle correction formula (14) of each harmonic according to factor alpha and β value:
A＝N
^{1}·(y
_{1}+y
_{2})·(0.3065+0.6624*β
^{2}+0.9504*β
^{4}) (13)
θ＝arg[X(k
_{i}·Δf)]π·(α(1)
^{i}·0.5)...(i＝1,2) (14)
Step 4, can obtain fundamental frequency f according to formula (10)
_{1}=f
_{0}; Amplitude and the phase place of each harmonic obtain based on fundamental frequency, are specially: at [(i*f
_{0}5), (i*f
_{0}+ 5)] in frequency range, find the corresponding discrete spectral line k of i subharmonic
_{i1}and k
_{i2}, by corresponding amplitude y
_{i1}and y
_{i2}substitution formula (13)～(14), amplitude and the phase place of acquisition i subharmonic.Repeatedly repeat this step until obtain amplitude and the phase place of all harmonic waves.
Frequency analysis result to the present embodiment is carried out error analysis, and fundamental frequency relative error is in Table 3, and the amplitude of harmonic wave and the relative error of phase place are in Table 4～5.
The fundamental frequency relative error of the different product windows of table 3
The each harmonic amplitude relative error of the different product windows of table 4
The each harmonic phase place relative error of the different product windows of table 5
The present embodiment also adopts respectively traditional Hamming (Hamming) window, peaceful (Hanning) window of the Chinese and this (BlackmanHarris) window of haab to carry out frequency analysis to abovementioned power system signal x (t), and obtain amplitude relative error and the phase angle relative error that various window functions are corresponding, see respectively Fig. 5 and Fig. 6.Result shows, compares with classical window function, and Hamming multiplication window is with the obvious advantage: side lobe peak is lower.When N basic window multiplies each other, the sidelobe performance of Hamming multiplication window is obviously better than basic window, can better suppress spectral leakage, obviously improves harmonic parameters computational accuracy.
Claims (6)
1. the signal harmonic analytical approach based on Hamming product window, is characterized in that, comprises step:
Step 1, carries out analog to digital conversion to power system signal x (t) and obtains discretize signal x (n), and described power system signal x (t) is voltage signal or current signal;
Step 2, adopts Hamming product window to discretize signal x (n) windowing, and signal after windowing and Hamming product window are carried out to FFT conversion, obtains signal spectrum figure and Hamming product window after FFT conversion; Described Hamming product window is the acquisition of multiplying each other of p Hamming window, and p value is set according to actual conditions;
Step 3, obtains actual fundamental frequency according to signal spectrum figure and Hamming product window after FFT conversion;
Step 4, according to amplitude, the phase place of the signal spectrum figure after FFT conversion and actual fundamental frequency acquisition firstharmonic and each harmonic.
2. the signal harmonic analytical approach based on Hamming product window as claimed in claim 1, is characterized in that:
Described Hamming window is cosine Hamming window.
3. the signal harmonic analytical approach based on Hamming product window as claimed in claim 1, is characterized in that:
Described signal spectrum figure and Hamming product window according to after FFT conversion obtain actual fundamental frequency, are specially:
By put the discrete spectral line sequence number in the nearest left and right sides from theoretical peak, be designated as respectively k
_{1}and k
_{2}, discrete spectral line k
_{1}and k
_{2}corresponding amplitude is respectively y
_{1}and y
_{2}, according to
obtain the discrete spectral line sequence number k of actual firstharmonic
_{0}thereby, obtain fundamental frequency; Wherein, W () is the Hamming product window after FFT conversion.
4. the signal harmonic analytical approach based on Hamming product window as claimed in claim 1, is characterized in that:
Described amplitude, the phase place that obtains firstharmonic and each harmonic according to the signal spectrum figure after FFT conversion and actual fundamental frequency, is specially:
[(the i*f of signal spectrum figure after FFT conversion
_{0}a), (i*f
_{0}+ a)] look for the discrete spectral line k that i subharmonic is corresponding in frequency range
_{i1}and k
_{i2}, i is overtone order; The harmonic wave maximum deviation that a allows according to electric system arranges; Based on discrete spectral line k
_{i1}and k
_{i2}corresponding amplitude y
_{i1}and y
_{i2}, adopt amplitude correction formula and phase place correction formula to obtain amplitude and the phase angle of i subharmonic.
5. the signal harmonic analytic system based on Hamming product window, is characterized in that, comprising:
Analogtodigital conversion module, is used for that power system signal x (t) is carried out to analog to digital conversion and obtains discretize signal x (n), and described power system signal x (t) is voltage signal or current signal;
Windowing and FFT change module, are used for adopting Hamming product window to discretize signal x (n) windowing, and signal after windowing and Hamming product window are carried out to FFT conversion, obtain signal spectrum figure and Hamming product window after FFT conversion; Described Hamming product window is the acquisition of multiplying each other of p Hamming window, and p value is set according to actual conditions;
Actual fundamental frequency obtains module, is used for obtaining actual fundamental frequency according to signal spectrum figure and Hamming product window after FFT conversion;
Harmonic amplitude and harmonic phase obtain module, are used for according to amplitude, the phase place of the signal spectrum figure after FFT conversion and actual fundamental frequency acquisition firstharmonic and each harmonic.
6. the signal harmonic analytic system based on Hamming product window as claimed in claim 5, is characterized in that:
Described harmonic amplitude and harmonic phase obtain module and further comprise:
The discrete spectral line that harmonic wave is corresponding is looked for module, is used for the [(i*f of the signal spectrum figure after FFT conversion
_{0}a), (i*f
_{0}+ a)] look for the discrete spectral line k that i subharmonic is corresponding in frequency range
_{i1}and k
_{i2}, i is overtone order; The harmonic wave maximum deviation that a allows according to electric system arranges;
Amplitude and phase solution module, be used for based on discrete spectral line k
_{i1}and k
_{i2}corresponding amplitude y
_{i1}and y
_{i2}, adopt amplitude correction formula and phase place correction formula to obtain amplitude and the phase angle of i subharmonic.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201410321280.6A CN104062500A (en)  20140704  20140704  Signal harmonic analysis method and system based on Hamming product window 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201410321280.6A CN104062500A (en)  20140704  20140704  Signal harmonic analysis method and system based on Hamming product window 
Publications (1)
Publication Number  Publication Date 

CN104062500A true CN104062500A (en)  20140924 
Family
ID=51550304
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201410321280.6A Pending CN104062500A (en)  20140704  20140704  Signal harmonic analysis method and system based on Hamming product window 
Country Status (1)
Country  Link 

CN (1)  CN104062500A (en) 
Cited By (10)
Publication number  Priority date  Publication date  Assignee  Title 

CN104375111A (en) *  20141116  20150225  甘肃省机械科学研究院  Rapid highprecision refining correction method for intensive frequency spectrum 
CN104897961A (en) *  20150617  20150909  中南民族大学  Three spectral line interpolation FFT harmonic wave analysis method and system based on multiplication window function 
CN105137175A (en) *  20150806  20151209  湖南大学  Windowing shorttime Fourier transform threepoint interpolation dynamic frequency measurement method 
CN105486921A (en) *  20160122  20160413  武汉大学  Kaiser thirdorder mutual convolution window triplespectrumline interpolation harmonic wave and interharmonic wave detection method 
CN105717359A (en) *  20160219  20160629  云南电网有限责任公司电力科学研究院  Harmonic analysis algorithm 
CN106771586A (en) *  20161124  20170531  云南电网有限责任公司电力科学研究院  The loop signal analysis method and device of a kind of DC control protection board 
CN108051682A (en) *  20171222  20180518  西南交通大学  A kind of verification method of singlephase rectifier system impedance model 
CN108680787A (en) *  20180523  20181019  成都玖锦科技有限公司  Real time spectral analysis method based on FPGA 
CN108776263A (en) *  20180502  20181109  三峡大学  Harmonic detecting method based on the peaceful selfconvolution window of the highorder Chinese and improvement interpolation algorithm 
CN113640579A (en) *  20211013  20211112  四川大学  Harmonic measurement method based on double spectral line transformation, electronic device and storage medium 
Citations (7)
Publication number  Priority date  Publication date  Assignee  Title 

CN101261618A (en) *  20080414  20080910  湖南大学  Triangular selfconvolution window weighting algorithm 
JP2009264753A (en) *  20080422  20091112  Yokogawa Electric Corp  Harmonics measuring apparatus 
CN102495285A (en) *  20111124  20120613  中国矿业大学  Method for estimating power harmonic wave parameter by using power gravity center of symmetric window function 
CN102539915A (en) *  20120106  20120704  中国矿业大学  Method for accurately calculating power harmonic wave parameters through adopting time delay Fourier transform frequency measurement method 
CN102818930A (en) *  20120903  20121212  中国矿业大学  Method for quickly calculating power harmonic parameters in highaccuracy mode 
CN103207319A (en) *  20130312  20130717  广东电网公司电力科学研究院  Harmonic wave measurement method of electricity signal of digital substation under nonsynchronous sampling condition 
CN103308766A (en) *  20130515  20130918  湖南大学  Harmonic analysis method based on Kaiser selfconvolution window dualspectrum line interpolation FFT (Fast Fourier Transform) and device thereof 

2014
 20140704 CN CN201410321280.6A patent/CN104062500A/en active Pending
Patent Citations (7)
Publication number  Priority date  Publication date  Assignee  Title 

CN101261618A (en) *  20080414  20080910  湖南大学  Triangular selfconvolution window weighting algorithm 
JP2009264753A (en) *  20080422  20091112  Yokogawa Electric Corp  Harmonics measuring apparatus 
CN102495285A (en) *  20111124  20120613  中国矿业大学  Method for estimating power harmonic wave parameter by using power gravity center of symmetric window function 
CN102539915A (en) *  20120106  20120704  中国矿业大学  Method for accurately calculating power harmonic wave parameters through adopting time delay Fourier transform frequency measurement method 
CN102818930A (en) *  20120903  20121212  中国矿业大学  Method for quickly calculating power harmonic parameters in highaccuracy mode 
CN103207319A (en) *  20130312  20130717  广东电网公司电力科学研究院  Harmonic wave measurement method of electricity signal of digital substation under nonsynchronous sampling condition 
CN103308766A (en) *  20130515  20130918  湖南大学  Harmonic analysis method based on Kaiser selfconvolution window dualspectrum line interpolation FFT (Fast Fourier Transform) and device thereof 
NonPatent Citations (2)
Title 

王玲 等: "一种新的余弦组合窗插值FFT谐波分析算法", 《武汉大学学报（工学版）》, vol. 47, no. 2, 30 April 2014 (20140430), pages 250  254 * 
魏渊: "脉冲雷达的低副瓣数字脉压方法的研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》, no. 7, 15 July 2013 (20130715), pages 136  1104 * 
Cited By (15)
Publication number  Priority date  Publication date  Assignee  Title 

CN104375111B (en) *  20141116  20170329  甘肃省机械科学研究院  The method that quick high accuracy refinement correction is carried out to intensive spectrum 
CN104375111A (en) *  20141116  20150225  甘肃省机械科学研究院  Rapid highprecision refining correction method for intensive frequency spectrum 
CN104897961A (en) *  20150617  20150909  中南民族大学  Three spectral line interpolation FFT harmonic wave analysis method and system based on multiplication window function 
CN105137175B (en) *  20150806  20171107  湖南大学  A kind of adding window Short Time Fourier Transform three point interpolation Dynamic Frequency Measurement method 
CN105137175A (en) *  20150806  20151209  湖南大学  Windowing shorttime Fourier transform threepoint interpolation dynamic frequency measurement method 
CN105486921A (en) *  20160122  20160413  武汉大学  Kaiser thirdorder mutual convolution window triplespectrumline interpolation harmonic wave and interharmonic wave detection method 
CN105717359A (en) *  20160219  20160629  云南电网有限责任公司电力科学研究院  Harmonic analysis algorithm 
CN106771586A (en) *  20161124  20170531  云南电网有限责任公司电力科学研究院  The loop signal analysis method and device of a kind of DC control protection board 
CN106771586B (en) *  20161124  20190726  云南电网有限责任公司电力科学研究院  A kind of the loop signal analysis method and device of DC control protection board 
CN108051682A (en) *  20171222  20180518  西南交通大学  A kind of verification method of singlephase rectifier system impedance model 
CN108051682B (en) *  20171222  20191220  西南交通大学  Verification method of singlephase rectifier system impedance model 
CN108776263A (en) *  20180502  20181109  三峡大学  Harmonic detecting method based on the peaceful selfconvolution window of the highorder Chinese and improvement interpolation algorithm 
CN108776263B (en) *  20180502  20200728  三峡大学  Harmonic detection method based on highorder Hanning selfconvolution window and improved interpolation algorithm 
CN108680787A (en) *  20180523  20181019  成都玖锦科技有限公司  Real time spectral analysis method based on FPGA 
CN113640579A (en) *  20211013  20211112  四川大学  Harmonic measurement method based on double spectral line transformation, electronic device and storage medium 
Similar Documents
Publication  Publication Date  Title 

CN104062500A (en)  Signal harmonic analysis method and system based on Hamming product window  
CN104062528A (en)  Signal harmonic analysis method and system based on Hanning product window  
CN102435844B (en)  Sinusoidal signal phasor calculating method being independent of frequency  
CN202339381U (en)  Harmonic electric energy metering system based on Nuttall selfconvolution window weighed FFT (Fast Fourier Transform)  
CN101852826B (en)  Harmonic analysis method for power system and device thereof  
CN102539915B (en)  Method for accurately calculating power harmonic wave parameters through adopting time delay Fourier transform frequency measurement method  
CN105137180B (en)  Highprecision harmonic analysis method based on six four spectral line interpolations of Cosine Window  
CN103353550A (en)  Method for measuring signal frequency and harmonic parameters of electric power system  
CN101701984B (en)  Fundamental wave and harmonic wave detecting method based on threecoefficient Nuttall windowed interpolation FFT  
CN103308804A (en)  Method for extracting timefrequency parameters of power quality disturbance signals on basis of fast KS (KaiserS) transformation  
CN101261292A (en)  Base wave and harmonic detection method based on fiver item RifeVincent(1)window double spectral line interpolation FFT  
CN103197141A (en)  Method of measuring electrical power system signal frequency and harmonic wave parameters  
CN103308766A (en)  Harmonic analysis method based on Kaiser selfconvolution window dualspectrum line interpolation FFT (Fast Fourier Transform) and device thereof  
CN102331526B (en)  Method for acquiring parameters of electric power harmonic waves by using Hanniing window function continuous frequency spectrum interpolation  
CN103207319A (en)  Harmonic wave measurement method of electricity signal of digital substation under nonsynchronous sampling condition  
CN102495285B (en)  Method for estimating power harmonic wave parameter by using power gravity center of symmetric window function  
CN104391178A (en)  Time shift phase difference steady harmonic signal correction method based on Nuttall window  
CN103399204A (en)  RifeVincent (II) window interpolation FFT (Fast Fourier Transform)based harmonic and interharmonic detection method  
CN101216512A (en)  Nonsine periodic signal real time high precision detection method  
CN106324340B (en)  A kind of method of synchronized phasor and frequency measurement dynamic property  
CN102818930A (en)  Method for quickly calculating power harmonic parameters in highaccuracy mode  
CN103091545A (en)  Sinusoidal signal phasor halfwave computing method irrelevant to frequency  
CN106018956A (en)  Power system frequency calculation method of windowing spectral line interpolation  
CN105486921A (en)  Kaiser thirdorder mutual convolution window triplespectrumline interpolation harmonic wave and interharmonic wave detection method  
CN103543331B (en)  A kind of method calculating electric signal harmonic wave and mAcetyl chlorophosphonazo 
Legal Events
Date  Code  Title  Description 

C06  Publication  
PB01  Publication  
C10  Entry into substantive examination  
SE01  Entry into force of request for substantive examination  
WD01  Invention patent application deemed withdrawn after publication 
Application publication date: 20140924 

WD01  Invention patent application deemed withdrawn after publication 