CN108982966B  Harmonic phase angle analysis method based on linear correction algorithm  Google Patents
Harmonic phase angle analysis method based on linear correction algorithm Download PDFInfo
 Publication number
 CN108982966B CN108982966B CN201811015823.6A CN201811015823A CN108982966B CN 108982966 B CN108982966 B CN 108982966B CN 201811015823 A CN201811015823 A CN 201811015823A CN 108982966 B CN108982966 B CN 108982966B
 Authority
 CN
 China
 Prior art keywords
 sampling
 analysis
 harmonic
 data
 sampling point
 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
Images
Abstract
The invention relates to a harmonic phase angle analysis method based on a linear correction algorithm, which comprises the following steps: (1) samplingW+2Sampling point data (WDetermined by the integration method); (2) from the sampling pointi= 0Using quasisynchronous DFT equationsAnalysis ofW+1Obtaining fundamental wave information from dataAnd(ii) a (3) From the sampling pointi=1Using quasisynchronous DFT formula analysisW+1Obtaining fundamental wave information from dataAnd(ii) a (4) Using formulasCalculating the frequency drift of a signal(ii) a (5) From the sampling pointi=0Using quasisynchronous DFT formula analysisW+1Obtaining the information of each harmonic wave by the dataAnd(ii) a (6) Using formulasCalculating amplitude angles of the harmonics; (7) using formulasAnd linearly correcting harmonic phase angles of the harmonics. The method is beneficial to more accurately obtaining the information of amplitude and phase angles, frequencies and the like of each harmonic wave in the fields of harmonic wave analysis such as power quality monitoring, electronic product production inspection, electrical equipment monitoring and the like.
Description
The application has the following application numbers: 201510258020.3 entitled "a harmonic phase angle analysis method", filed as follows: divisional application of the invention patent application on 19/05/2015.
Technical Field
The invention relates to a highprecision harmonic phase angle analysis method.
Background
The harmonic analysis technology is widely applied to the fields of power quality monitoring, electronic product production inspection, electrical equipment monitoring and the like, and is an important technical means for power grid monitoring, quality inspection and equipment monitoring. The most widely used techniques for harmonic analysis are Discrete Fourier Transform (DFT) and Fast Fourier Transform (FFT) at present. The harmonic analysis technology combining the quasisynchronous sampling technology and the DFT technology can improve the accuracy of harmonic analysis, and the formula is as follows:
in the formula:kfor the order of harmonics to be obtained (e.g. fundamental waves)k=1Harmonic of order 3k=3) (ii) a sin and cos are sine and cosine functions, respectively; whilea _{ k }Andb _{ k }are respectively askReal and imaginary parts of the subharmonics;nis the iteration number;Wdetermined by the integration method, when the complex trapezoidal integration method is adopted,W=nN；is a primary weighting coefficient;is the sum of all weighting coefficients;f(i)is the ith sampling value of the analysis waveform;Nis the number of samples in a cycle.
In engineering applications, harmonic analysis always performs finite point sampling and synchronization sampling which is difficult to achieve strictly. Thus, when the quasisynchronous DFT is applied to harmonic analysis, longrange leakage caused by truncation effect and shortrange leakage caused by barrier effect exist, so that the accuracy of the analysis result is not high, even the analysis result is not credible.
FIG. 1 presents an error plot of harmonic analysis for any given example using quasisynchronous DFT. As can be seen from the figure, the harmonic angle of the quasisynchronous DFT algorithm has extremely large errors except for 50Hz, and is basically not credible.
Disclosure of Invention
The invention aims to solve the technical problem of providing a highprecision harmonic phase angle analysis method based on a linear correction algorithm so as to effectively improve the analysis error of a quasisynchronous DFT harmonic analysis technology and obtain a highprecision harmonic analysis result, thereby improving the effectiveness of quality and state judgment of instruments and equipment in the fields of power quality monitoring, electronic product production inspection, electrical equipment monitoring and the like based on a harmonic analysis theory.
The technical scheme for realizing the aim of the invention is to provide a harmonic phase angle analysis method based on a linear correction algorithm, which comprises the following steps:
(1) w +2 sampling point data { f (i) =0,1, …, W +1} are sampled at equal intervals (W is determined by a selected integration method, the invention does not specify a certain integration method, common integration methods include a complex trapezoidal integration method W = nN, a complex rectangular integration method W = N (N1), a complex simpson integration method W = N (N1)/2 and the like, and a suitable integration method can be selected according to the practical situation of the application of the invention.
(2) From the sampling pointi=0Start of applying quasisynchronous DFT formula
(3) from the sampling pointi=1Using quasisynchronous DFT formula analysisW+1Obtaining fundamental wave information from dataAnd；
(5) From the sampling pointi=0Begin applying quasisynchronous DFT formula analysisW+1Obtaining the information of each harmonic wave by the dataAnd；
The equally spaced sampling is based on the period of the ideal signal for harmonic analysisTSum frequencyfSampling in one cycleNAt a point, i.e. at a sampling frequency off _{ s } =NfAnd is andN≥64。
the samplingW+2The sampling point data is selected according to the selected integration method, if the trapezoidal integration method is adopted, the sampling point data is selected correspondinglyW=nN(ii) a If the complex rectangular integral method is adopted, thenW=n(N1)(ii) a If the Simpson integration method is adopted, thenW=n(N1)/2. Then according to the sampling frequencyf _{ s } =NfObtaining a sequence of sample point data{f(i),i=0,1,…,W+ 1}，nFor the number of iterations, in generaln≥3(ii) a And finally, carrying out harmonic analysis on the data sequence.
Coefficient of first iterationBy integration methods, ideal period sampling pointsNAnd number of iterationsnDecision, specific derivation procedures see literature [ wearable in the middle of the world ] several problems in quasisynchronous sampling applications [ J]Electrical measuring and instrumentation, 1988, (2): 27.
Drift of signal frequencyBased on the phase angle difference of the fundamental wave of adjacent sampling points and the number of sampling points in an ideal periodNObtained by a fixed relationship of the frequency of the signalCan also be used for correcting the frequency of fundamental wave and higher harmonic wavef _{ 1 }And frequencies of higher harmonicsf _{ k }（）。
The invention has the positive effects that: (1) and (5) analyzing the result of the harmonic phase angle with high precision. As for the analysis example given in FIG. 1, the accuracy of the analysis obtained by the present invention is improved to 10^{8}Stage (fig. 2).
(2) The method provided by the invention fundamentally solves the problem of low analysis precision of quasisynchronous DFT harmonic phase angles, does not need to perform complicated inversion and correction, and is simple in algorithm.
(3) Compared with quasisynchronous DFT, the harmonic analysis technology of the invention only needs to add one sampling point to solve the problem of large error of quasisynchronous DFT analysis, and is easy to realize.
(4) The invention is technically feasible for improving the existing instrument and equipment, and the analysis result can be improved to 10 without increasing any hardware expense^{8}And (4) stages.
(5) The method is also suitable for the harmonic analysis process of carrying out multiple iterations instead of one iteration, and only one iteration needs to be decomposed into multiple iterations to realize the harmonic analysis process. One iteration is essentially the same as multiple iterations, except that in the calculation, the multiple iterations are subjected to stepbystep calculation, and the process of the multiple iterations is combined into iteration coefficients in one iterationThe calculation is completed in one time, so the method is also suitable for a plurality of iterative processes.
Drawings
FIG. 1 is a harmonic phase analysis error plot of a quasisynchronous DFT.
FIG. 2 is a harmonic phase angle analysis error map of the present invention.
Detailed Description
(example 1)
The harmonic phase angle analysis method based on the linear correction algorithm in the embodiment comprises the following steps:
first, sampling at equal intervalsW+2Sampling points to obtain a discrete sequence of the signal to be analyzed{f(i),i=0,1,…, Wq+1}。WBy an integration method, number of iterationsnAnd number of samples taken in an ideal periodNAnd (4) jointly determining.
The equal interval sampling means that: based on the frequency of the desired signal for harmonic analysis (e.g. frequency of power frequency signal)fFor 50Hz with a period of 20 mS) determining the sampling frequencyf _{ S } =NfAt the sampling frequencyf _{ S }Is uniformly sampled in one periodNAnd (4) point. Generally, periodic sampling pointsN=64Or more, better harmonic analysis result can be obtained, and the iteration timesn=3~5A more ideal harmonic analysis result can be obtained.
The integration method includes various methods such as a trapezoidal integration method, a rectangular integration method, and a simpson method, and can be selected according to actual conditions. If the method of complex trapezoidal integration is adopted, thenW=nN(ii) a If the complex rectangular integral method is adopted, thenW=n (N1)(ii) a If the Simpson integration method is adopted, thenW=n(N1)/2。
Second, from the sampling pointi=0Start of applying quasisynchronous DFT formula
again, from the sampling pointi=1Using quasisynchronous DFT formula analysisW+1Obtaining fundamental wave information from dataAnd；
Again, from the sampling pointi=0Begin applying quasisynchronous DFT formula analysisW+1Obtaining the information of each harmonic from the data；
It will be appreciated by persons skilled in the art that the above embodiments are only intended to illustrate the present invention, and not to limit the present invention, and that the present invention may be further modified, and that within the spirit and scope of the present invention, changes and modifications to the above described embodiments will fall within the scope of the appended claims.
Claims (2)
1. A harmonic phase angle analysis method based on a linear correction algorithm is characterized by comprising the following steps: the method comprises the following steps:
(1) equidistant samplingW+2Data of each sampling point{f(i),i=0,1,…,W+1}(ii) a The samplingW+2The data of each sampling point adopts a complex Simpson integration method, thenW=n(N1)/2；
(2) From the sampling pointi=0Start of applying quasisynchronous DFT formula
(3) from the sampling pointi=1Using quasisynchronous DFT formula analysisW+1Obtaining fundamental wave information from dataAnd；
(5) From the sampling pointi=0Begin applying quasisynchronous DFT formula analysisW+1Obtaining the information of each harmonic wave by the dataAnd；
in the formula:kthe number of harmonics to be obtained; sin and cos are sine and cosine functions, respectively; whilea _{ k }Andb _{ k }are respectively askReal and imaginary parts of the subharmonics;nis the iteration number;Wdetermined by an integration method;is a primary weighting coefficient;is the sum of all weighting coefficients;f(i)is the ith sampling value of the analysis waveform;Nthe number of sampling times in a period;
in the step (1), the equidistant sampling is based on the period of the ideal signal for harmonic analysisTSum frequencyfSampling in one cycleNAt a point, i.e. at a sampling frequency off _{ s } =NfAnd is andN≥64。
2. the harmonic phase angle analysis method based on the linear correction algorithm according to claim 1, characterized in that: in the step (1), the samplingW+2The data of each sampling point is selected according to the selected integration method and then according to the sampling frequencyf _{ s } =NfObtaining a sequence of sample point data{f(i),i=0,1,…,W+1}；nIn order to be able to perform the number of iterations,n≥3(ii) a And finally, carrying out harmonic analysis on the data sequence.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201811015823.6A CN108982966B (en)  20150519  20150519  Harmonic phase angle analysis method based on linear correction algorithm 
Applications Claiming Priority (2)
Application Number  Priority Date  Filing Date  Title 

CN201510258020.3A CN104849551B (en)  20150519  20150519  A kind of humorous phase angle analysis method 
CN201811015823.6A CN108982966B (en)  20150519  20150519  Harmonic phase angle analysis method based on linear correction algorithm 
Related Parent Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201510258020.3A Division CN104849551B (en)  20150519  20150519  A kind of humorous phase angle analysis method 
Publications (2)
Publication Number  Publication Date 

CN108982966A CN108982966A (en)  20181211 
CN108982966B true CN108982966B (en)  20201009 
Family
ID=53849342
Family Applications (4)
Application Number  Title  Priority Date  Filing Date 

CN201811015823.6A Active CN108982966B (en)  20150519  20150519  Harmonic phase angle analysis method based on linear correction algorithm 
CN201811014781.4A Active CN109030942B (en)  20150519  20150519  Harmonic phase angle analysis method 
CN201811014772.5A Active CN109142866B (en)  20150519  20150519  Harmonic phase angle analysis method based on linear correction algorithm 
CN201510258020.3A Active CN104849551B (en)  20150519  20150519  A kind of humorous phase angle analysis method 
Family Applications After (3)
Application Number  Title  Priority Date  Filing Date 

CN201811014781.4A Active CN109030942B (en)  20150519  20150519  Harmonic phase angle analysis method 
CN201811014772.5A Active CN109142866B (en)  20150519  20150519  Harmonic phase angle analysis method based on linear correction algorithm 
CN201510258020.3A Active CN104849551B (en)  20150519  20150519  A kind of humorous phase angle analysis method 
Country Status (1)
Country  Link 

CN (4)  CN108982966B (en) 
Families Citing this family (2)
Publication number  Priority date  Publication date  Assignee  Title 

CN109683016A (en) *  20190118  20190426  江苏理工学院  A kind of harmonic amplitude analysis method 
CN110083891A (en) *  20190410  20190802  江苏理工学院  A kind of calculation method of the plesiochronous linear correction factor of DFT amplitude 
Family Cites Families (21)
Publication number  Priority date  Publication date  Assignee  Title 

JPS56148181A (en) *  19800421  19811117  Toshiba Corp  Method of controlling single phase multistage rectifier 
US5400366A (en) *  19920709  19950321  Fujitsu Limited  Quasisynchronous detection and demodulation circuit and frequency discriminator used for the same 
EP1367794B1 (en) *  19980130  20070314  Matsushita Electric Industrial Co., Ltd.  Modulation method and radio communication system 
JP2006276006A (en) *  20050301  20061012  Nagoya Institute Of Technology  Harmonic analysis method in power system 
CN101261293A (en) *  20070308  20080910  国网南京自动化研究院  Electric power steadystate signal tracking measurement based on selfadapting filter 
CN102135567A (en) *  20110225  20110727  深圳市业通达实业有限公司  Realtime frequency tracking and harmonic measuring method for AC sampling of power system 
CN102156250A (en) *  20110317  20110817  华北电力大学（保定）  Dielectric loss factor measurement method based on equivalent model 
CN102331526B (en) *  20110610  20130717  中国矿业大学  Method for acquiring parameters of electric power harmonic waves by using Hanniing window function continuous frequency spectrum interpolation 
CN102393488B (en) *  20110824  20131211  江苏理工学院  Harmonic analysis method 
CN103604989A (en) *  20111021  20140226  蒋春花  Harmonic wave analysis method of powerquality harmonic wave analyzer 
CN103592512A (en) *  20111021  20140219  蒋春花  Harmonic analysis method of power quality harmonic analyzer 
CN102508026B (en) *  20111021  20140806  江苏理工学院  Harmonic wave analysis method for electric energy quality harmonic wave analyzer 
CN102435860B (en) *  20111021  20131106  江苏理工学院  Working method of medium loss current tester 
CN103439566B (en) *  20111021  20141224  常州顺创电气科技有限公司  Operating method of MOA resistive current tester with relatively high precision 
CN103969507B (en) *  20111021  20161019  江苏理工学院  A kind of electric energy quality harmonic analyzes method 
JP5705102B2 (en) *  20111221  20150422  三菱電機株式会社  Insulation deterioration diagnosis device 
CN102590615A (en) *  20120216  20120718  安徽理工大学  Method for correcting synchronous phasor under condition of frequency shift of voltage of power grid, and current signal 
DE102012204963A1 (en) *  20120328  20131002  Siemens Aktiengesellschaft  Method and apparatus for reducing current harmonics 
CN103576002B (en) *  20131111  20160120  华北电力大学（保定）  A kind of computing method of capacitive insulator arrangement dielectric loss angle 
CN103760425A (en) *  20140122  20140430  湖南大学  Method and device for rapidly measuring dielectric loss angle based on time domain quasisynchronization 
CN104181391A (en) *  20140821  20141203  艾德克斯电子（南京）有限公司  Harmonic detection method of digital power meter 

2015
 20150519 CN CN201811015823.6A patent/CN108982966B/en active Active
 20150519 CN CN201811014781.4A patent/CN109030942B/en active Active
 20150519 CN CN201811014772.5A patent/CN109142866B/en active Active
 20150519 CN CN201510258020.3A patent/CN104849551B/en active Active
Also Published As
Publication number  Publication date 

CN109030942B (en)  20201204 
CN108982966A (en)  20181211 
CN104849551A (en)  20150819 
CN109142866A (en)  20190104 
CN109142866B (en)  20201009 
CN104849551B (en)  20181123 
CN109030942A (en)  20181218 
Similar Documents
Publication  Publication Date  Title 

CN103454494B (en)  A kind of highprecision harmonic analysis method  
CN102508026B (en)  Harmonic wave analysis method for electric energy quality harmonic wave analyzer  
CN110837001B (en)  Method and device for analyzing harmonic waves and interharmonic waves in electric power system  
CN103869162A (en)  Dynamic signal phasor measurement method based on time domain quasisynchronization  
CN102435860B (en)  Working method of medium loss current tester  
CN105137180B (en)  Highprecision harmonic analysis method based on six four spectral line interpolations of Cosine Window  
CN109239463B (en)  Dielectric loss measurement method based on linear correction algorithm  
CN103543331B (en)  A kind of method calculating electric signal harmonic wave and mAcetyl chlorophosphonazo  
CN108982966B (en)  Harmonic phase angle analysis method based on linear correction algorithm  
CN108896944A (en)  A kind of synchronous measuring apparatus laboratory investment instrument and its synchronous phasor measuring method  
CN103983849B (en)  A kind of Electric Power Harmonic Analysis method of realtime highprecision  
CN104931777B (en)  A kind of signal frequency measuring method based on two DFT plural number spectral lines  
CN105372492B (en)  Signal frequency measuring method based on three DFT plural number spectral lines  
CN109490603B (en)  Method for measuring resistive current fundamental wave  
CN104678170A (en)  Power harmonic analysis method based on harmonic analyzer and harmonic analyzer  
CN103983852A (en)  Harmonic analysis method of power quality harmonic analyzer  
CN110083891A (en)  A kind of calculation method of the plesiochronous linear correction factor of DFT amplitude  
CN103592512A (en)  Harmonic analysis method of power quality harmonic analyzer  
CN103604989A (en)  Harmonic wave analysis method of powerquality harmonic wave analyzer  
CN105372493B (en)  Signal amplitude and Method for Phase Difference Measurement based on three DFT plural number spectral lines  
CN103969507B (en)  A kind of electric energy quality harmonic analyzes method  
CN108957118A (en)  A kind of reactive power calculating method  
CN109683016A (en)  A kind of harmonic amplitude analysis method  
CN116359605B (en)  Harmonic signal analysis method based on secondary weighting  
CN117169590A (en)  Power harmonic analysis method and device based on software variable sampling rate 
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 