CN108982966B - Harmonic phase angle analysis method based on linear correction algorithm - Google Patents

Harmonic phase angle analysis method based on linear correction algorithm Download PDF

Info

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
Application number
CN201811015823.6A
Other languages
Chinese (zh)
Other versions
CN108982966A (en
Inventor
傅中君
王建宇
欧云
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiangsu University of Technology
Original Assignee
Jiangsu University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jiangsu University of Technology filed Critical Jiangsu University of Technology
Priority to CN201811015823.6A priority Critical patent/CN108982966B/en
Publication of CN108982966A publication Critical patent/CN108982966A/en
Application granted granted Critical
Publication of CN108982966B publication Critical patent/CN108982966B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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 quasi-synchronous DFT equations
Figure 100004_DEST_PATH_IMAGE002
Analysis ofW+1Obtaining fundamental wave information from data
Figure 100004_DEST_PATH_IMAGE004
And
Figure 100004_DEST_PATH_IMAGE006
(ii) a (3) From the sampling pointi=1Using quasi-synchronous DFT formula analysisW+1Obtaining fundamental wave information from data
Figure 100004_DEST_PATH_IMAGE008
And
Figure 100004_DEST_PATH_IMAGE010
(ii) a (4) Using formulas
Figure 100004_DEST_PATH_IMAGE012
Calculating the frequency drift of a signal
Figure 100004_DEST_PATH_IMAGE014
(ii) a (5) From the sampling pointi=0Using quasi-synchronous DFT formula analysisW+1Obtaining the information of each harmonic wave by the data
Figure 100004_DEST_PATH_IMAGE016
And
Figure 100004_DEST_PATH_IMAGE018
(ii) a (6) Using formulas
Figure 100004_DEST_PATH_IMAGE020
Calculating amplitude angles of the harmonics; (7) using formulas
Figure 100004_DEST_PATH_IMAGE022
And 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

Harmonic phase angle analysis method based on linear correction algorithm
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 high-precision 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 quasi-synchronous sampling technology and the DFT technology can improve the accuracy of harmonic analysis, and the formula is as follows:
Figure 100002_DEST_PATH_IMAGE002
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
Figure 100002_DEST_PATH_IMAGE004
is a primary weighting coefficient;
Figure 100002_DEST_PATH_IMAGE006
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 quasi-synchronous DFT is applied to harmonic analysis, long-range leakage caused by truncation effect and short-range 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 quasi-synchronous DFT. As can be seen from the figure, the harmonic angle of the quasi-synchronous 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 high-precision harmonic phase angle analysis method based on a linear correction algorithm so as to effectively improve the analysis error of a quasi-synchronous DFT harmonic analysis technology and obtain a high-precision 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 (N-1), a complex simpson integration method W = N (N-1)/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 quasi-synchronous DFT formula
Figure 100002_DEST_PATH_IMAGE008
Analysis ofW+1Obtaining fundamental wave information from data
Figure 100002_DEST_PATH_IMAGE010
And
Figure 100002_DEST_PATH_IMAGE012
(3) from the sampling pointi=1Using quasi-synchronous DFT formula analysisW+1Obtaining fundamental wave information from data
Figure 100002_DEST_PATH_IMAGE014
And
Figure 100002_DEST_PATH_IMAGE016
(4) applying the formula:
Figure 100002_DEST_PATH_IMAGE018
calculating the frequency drift of the signal
Figure 100002_DEST_PATH_IMAGE020
(5) From the sampling pointi=0Begin applying quasi-synchronous DFT formula analysisW+1Obtaining the information of each harmonic wave by the data
Figure 100002_DEST_PATH_IMAGE022
And
Figure 100002_DEST_PATH_IMAGE024
(6) using formulas
Figure 100002_DEST_PATH_IMAGE026
Calculating amplitude angles of the harmonics;
(7) using formulas
Figure DEST_PATH_IMAGE028
And linearly correcting harmonic phase angles of the harmonics.
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(N-1)(ii) a If the Simpson integration method is adopted, thenW=n(N-1)/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 iteration
Figure 142132DEST_PATH_IMAGE004
By integration methods, ideal period sampling pointsNAnd number of iterationsnDecision, specific derivation procedures see literature [ wearable in the middle of the world ] several problems in quasi-synchronous sampling applications [ J]Electrical measuring and instrumentation, 1988, (2): 2-7.
Figure 856010DEST_PATH_IMAGE006
Is the sum of all weighting coefficients.
Drift of signal frequency
Figure DEST_PATH_IMAGE030
Based 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 signal
Figure 166906DEST_PATH_IMAGE030
Can also be used for correcting the frequency of fundamental wave and higher harmonic wavef 1 And frequencies of higher harmonicsf k
Figure DEST_PATH_IMAGE032
)。
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-8Stage (fig. 2).
(2) The method provided by the invention fundamentally solves the problem of low analysis precision of quasi-synchronous DFT harmonic phase angles, does not need to perform complicated inversion and correction, and is simple in algorithm.
(3) Compared with quasi-synchronous DFT, the harmonic analysis technology of the invention only needs to add one sampling point to solve the problem of large error of quasi-synchronous 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-8And (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 step-by-step calculation, and the process of the multiple iterations is combined into iteration coefficients in one iteration
Figure 296536DEST_PATH_IMAGE004
The 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 quasi-synchronous 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 (N-1)(ii) a If the Simpson integration method is adopted, thenW=n(N-1)/2
Second, from the sampling pointi=0Start of applying quasi-synchronous DFT formula
Figure 907646DEST_PATH_IMAGE008
Analysis ofW+1Obtaining fundamental wave information from data
Figure 464529DEST_PATH_IMAGE010
And
Figure 262721DEST_PATH_IMAGE012
again, from the sampling pointi=1Using quasi-synchronous DFT formula analysisW+1Obtaining fundamental wave information from data
Figure 258358DEST_PATH_IMAGE014
And
Figure 927237DEST_PATH_IMAGE016
again, the formula is applied:
Figure 655022DEST_PATH_IMAGE018
calculating the frequency drift of the signal
Figure 206089DEST_PATH_IMAGE020
Again, from the sampling pointi=0Begin applying quasi-synchronous DFT formula analysisW+1Obtaining the information of each harmonic from the data
Figure 677521DEST_PATH_IMAGE024
Then, apply the formula
Figure DEST_PATH_IMAGE034
Calculating amplitude angles of the harmonics;
finally, the formula is applied
Figure DEST_PATH_IMAGE036
And linearly correcting harmonic phase angles of the harmonics.
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(N-1)/2
(2) From the sampling pointi=0Start of applying quasi-synchronous DFT formula
Figure DEST_PATH_IMAGE002
Analysis ofW+1Obtaining fundamental wave information from data
Figure DEST_PATH_IMAGE004
And
Figure DEST_PATH_IMAGE006
(3) from the sampling pointi=1Using quasi-synchronous DFT formula analysisW+1Obtaining fundamental wave information from data
Figure DEST_PATH_IMAGE008
And
Figure DEST_PATH_IMAGE010
(4) applying the formula:
Figure DEST_PATH_IMAGE012
calculating the frequency drift of the signal
Figure DEST_PATH_IMAGE014
(5) From the sampling pointi=0Begin applying quasi-synchronous DFT formula analysisW+1Obtaining the information of each harmonic wave by the data
Figure DEST_PATH_IMAGE016
And
Figure DEST_PATH_IMAGE018
(6) using formulas
Figure DEST_PATH_IMAGE020
Calculating amplitude angles of the harmonics;
(7) using formulas
Figure DEST_PATH_IMAGE022
Linearly correcting harmonic angles of each harmonic;
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;
Figure DEST_PATH_IMAGE024
is a primary weighting coefficient;
Figure DEST_PATH_IMAGE026
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.
CN201811015823.6A 2015-05-19 2015-05-19 Harmonic phase angle analysis method based on linear correction algorithm Active CN108982966B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811015823.6A CN108982966B (en) 2015-05-19 2015-05-19 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) 2015-05-19 2015-05-19 A kind of humorous phase angle analysis method
CN201811015823.6A CN108982966B (en) 2015-05-19 2015-05-19 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) 2015-05-19 2015-05-19 A kind of humorous phase angle analysis method

Publications (2)

Publication Number Publication Date
CN108982966A CN108982966A (en) 2018-12-11
CN108982966B true CN108982966B (en) 2020-10-09

Family

ID=53849342

Family Applications (4)

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

Family Applications After (3)

Application Number Title Priority Date Filing Date
CN201811014781.4A Active CN109030942B (en) 2015-05-19 2015-05-19 Harmonic phase angle analysis method
CN201811014772.5A Active CN109142866B (en) 2015-05-19 2015-05-19 Harmonic phase angle analysis method based on linear correction algorithm
CN201510258020.3A Active CN104849551B (en) 2015-05-19 2015-05-19 A kind of humorous phase angle analysis method

Country Status (1)

Country Link
CN (4) CN108982966B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109683016A (en) * 2019-01-18 2019-04-26 江苏理工学院 A kind of harmonic amplitude analysis method
CN110083891A (en) * 2019-04-10 2019-08-02 江苏理工学院 A kind of calculation method of the plesiochronous linear correction factor of DFT amplitude

Family Cites Families (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS56148181A (en) * 1980-04-21 1981-11-17 Toshiba Corp Method of controlling single phase multistage rectifier
US5400366A (en) * 1992-07-09 1995-03-21 Fujitsu Limited Quasi-synchronous detection and demodulation circuit and frequency discriminator used for the same
EP1367794B1 (en) * 1998-01-30 2007-03-14 Matsushita Electric Industrial Co., Ltd. Modulation method and radio communication system
JP2006276006A (en) * 2005-03-01 2006-10-12 Nagoya Institute Of Technology Harmonic analysis method in power system
CN101261293A (en) * 2007-03-08 2008-09-10 国网南京自动化研究院 Electric power steady-state signal tracking measurement based on self-adapting filter
CN102135567A (en) * 2011-02-25 2011-07-27 深圳市业通达实业有限公司 Real-time frequency tracking and harmonic measuring method for AC sampling of power system
CN102156250A (en) * 2011-03-17 2011-08-17 华北电力大学(保定) Dielectric loss factor measurement method based on equivalent model
CN102331526B (en) * 2011-06-10 2013-07-17 中国矿业大学 Method for acquiring parameters of electric power harmonic waves by using Hanniing window function continuous frequency spectrum interpolation
CN102393488B (en) * 2011-08-24 2013-12-11 江苏理工学院 Harmonic analysis method
CN103604989A (en) * 2011-10-21 2014-02-26 蒋春花 Harmonic wave analysis method of power-quality harmonic wave analyzer
CN103592512A (en) * 2011-10-21 2014-02-19 蒋春花 Harmonic analysis method of power quality harmonic analyzer
CN102508026B (en) * 2011-10-21 2014-08-06 江苏理工学院 Harmonic wave analysis method for electric energy quality harmonic wave analyzer
CN102435860B (en) * 2011-10-21 2013-11-06 江苏理工学院 Working method of medium loss current tester
CN103439566B (en) * 2011-10-21 2014-12-24 常州顺创电气科技有限公司 Operating method of MOA resistive current tester with relatively high precision
CN103969507B (en) * 2011-10-21 2016-10-19 江苏理工学院 A kind of electric energy quality harmonic analyzes method
JP5705102B2 (en) * 2011-12-21 2015-04-22 三菱電機株式会社 Insulation deterioration diagnosis device
CN102590615A (en) * 2012-02-16 2012-07-18 安徽理工大学 Method for correcting synchronous phasor under condition of frequency shift of voltage of power grid, and current signal
DE102012204963A1 (en) * 2012-03-28 2013-10-02 Siemens Aktiengesellschaft Method and apparatus for reducing current harmonics
CN103576002B (en) * 2013-11-11 2016-01-20 华北电力大学(保定) A kind of computing method of capacitive insulator arrangement dielectric loss angle
CN103760425A (en) * 2014-01-22 2014-04-30 湖南大学 Method and device for rapidly measuring dielectric loss angle based on time domain quasi-synchronization
CN104181391A (en) * 2014-08-21 2014-12-03 艾德克斯电子(南京)有限公司 Harmonic detection method of digital power meter

Also Published As

Publication number Publication date
CN109030942B (en) 2020-12-04
CN108982966A (en) 2018-12-11
CN104849551A (en) 2015-08-19
CN109142866A (en) 2019-01-04
CN109142866B (en) 2020-10-09
CN104849551B (en) 2018-11-23
CN109030942A (en) 2018-12-18

Similar Documents

Publication Publication Date Title
CN103454494B (en) A kind of high-precision 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 inter-harmonic waves in electric power system
CN103869162A (en) Dynamic signal phasor measurement method based on time domain quasi-synchronization
CN102435860B (en) Working method of medium loss current tester
CN105137180B (en) High-precision 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 m-Acetyl 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 real-time high-precision
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 power-quality 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