CN108982966A - Humorous phase angle analysis method based on linear correction algorithm - Google Patents
Humorous phase angle analysis method based on linear correction algorithm Download PDFInfo
- Publication number
- CN108982966A CN108982966A CN201811015823.6A CN201811015823A CN108982966A CN 108982966 A CN108982966 A CN 108982966A CN 201811015823 A CN201811015823 A CN 201811015823A CN 108982966 A CN108982966 A CN 108982966A
- Authority
- CN
- China
- Prior art keywords
- analysis
- formula
- phase angle
- data
- harmonic
- 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
Abstract
The humorous phase angle analysis method based on linear correction algorithm that the present invention relates to a kind of, comprises the following steps: (1) sampling W+2 sample point data (W is determined by integration method)；(2) plesiochronous DFT formula is applied from sampled point i=0W+1 data of analysis obtain fundamental informationWith(3) fundamental information is obtained from sampled point i=1 using plesiochronous W+1 data of DFT formula analysisWith(4) formula is appliedCalculate the frequency drift μ of signal；(5) each harmonic information is obtained from sampled point i=0 using plesiochronous W+1 data of DFT formula analysisWith(6) formula is appliedCalculate the width phase angle of each harmonic；(7) formula is appliedThe humorous phase angle of linear amendment each harmonic.Method of the invention facilitates electric energy quality monitoring, electronic product production testing, electric appliances monitoring etc. using information such as the width phase angle of the more accurate acquisition each harmonic in the field of frequency analysis and frequencies.
Description
The application be application No. is: 201510258020.3, invention and created name be " a kind of humorous phase angle analysis method ", Shen
It please day are as follows: the divisional application of application for a patent for invention on May 19th, 2015.
Technical field
The present invention relates to a kind of high-precision humorous phase angle analysis methods.
Background technique
Frequency analysis technology is answered in various fields such as electric energy quality monitoring, electronic product production testing, electric appliances monitorings
It is the important technical for carrying out power system monitor, quality inspection, monitoring of tools with extensive.Frequency analysis is most widely used at present
Technology be discrete Fourier transform (DFT) and Fast Fourier Transform (FFT) (FFT).Quasi-synchronous sampling technique is mutually tied with DFT technique
The frequency analysis technology of conjunction can be improved the precision of frequency analysis, formula are as follows:
In formula: k is the number (such as fundamental wave k=1,3 subharmonic k=3) for needing the harmonic wave obtained；Sin and cos are positive respectively
String and cosine function；And akAnd bkThe respectively real and imaginary parts of k subharmonic；N is the number of iterations；W is determined by integration method, is adopted
When with muiltiple-trapezoid integration method, W=nN；γiFor a weighting coefficient；For the sum of all weighting coefficients；f(i)
For the ith sample value of analysis waveform；N is sampling number in the period.
In engineer application, frequency analysis always carry out the sampling of finite point be difficult to stricti jurise synchronize adopt
Sample.In this way, when the plesiochronous DFT of application carries out frequency analysis, will exist the leakage of the long range as caused by truncation effect and
The short range leakage as caused by fence effect, so that analysis result precision is not high or even insincere.
Fig. 1 gives the Error Graph for carrying out frequency analysis for any given example using plesiochronous DFT.It can be with from figure
It was found that remaining misses by a mile when the humorous phase angle of plesiochronous DFT algorithm is in addition to 50Hz, it is substantially insincere.
Summary of the invention
The technical problem to be solved in the present invention is to provide a kind of high-precision humorous phase angle analyses based on linear correction algorithm
Method, to be efficiently modified the analytical error of plesiochronous DFT frequency analysis technology, obtain high-precision frequency analysis as a result, to
The fields instruments such as electric energy quality monitoring, electronic product production testing, the electric appliances monitoring based on frequency analysis theory are improved to set
The validity of standby quality and state judgement.
It realizes that the technical solution of the object of the invention is to provide a kind of humorous phase angle analysis method based on linear correction algorithm, wraps
Include following steps:
(1) (W is by selected integration method for equidistant W+2 sample point data of sampling { f (i), i=0,1 ..., W+1 }
It determines, the not specified a certain integration method of the present invention, common integration method has muiltiple-trapezoid integration method W=nN, complexification
Rectangular integration method W=n (N-1), iterative Simpson integration method W=n (N-1)/2 etc., the reality that can be applied according to the present invention
Border situation selects suitable integration method.It is generally more satisfactory with muiltiple-trapezoid integration method effect.N is in an ideal period
Sampling number.
(2) plesiochronous DFT formula is applied since sampled point i=0
W+1 data of analysis obtain fundamental informationWith
(3) fundamental information is obtained from sampled point i=1 using plesiochronous W+1 data of DFT formula analysisWith
(4) formula is applied:Calculate the frequency drift μ of signal；
(5) each harmonic information is obtained using plesiochronous W+1 data of DFT formula analysis since sampled point i=0With
(6) formula is appliedCalculate the width phase angle of each harmonic；
(7) formula is appliedThe humorous phase angle of linear amendment each harmonic.
Equal interval sampling is sampled in one cycle according to the cycle T and frequency f of the ideal signal for carrying out frequency analysis
N point, i.e. sample frequency are fs=Nf, and N >=64.
Described W+2 sample point data of sampling is accordingly selected according to selected integration method, according to multiple
Change trapezoidal integration method, then W=nN；According to complexification rectangular integration method, then W=n (N-1)；According to iterative Simpson product
Divide method, then W=n (N-1)/2.Then according to sample frequency fs=Nf, acquisition sample point data sequence f (i), i=0,
1 ..., W+1 }, n is the number of iterations, general n >=3；Frequency analysis finally is carried out to the data sequence.
An iteration coefficient gammaiIt is determined by integration method, ideal period sampled point N and the number of iterations n, specific derivation process
Referring to document [wear some problem [J] the electrical measurement and instrument in the application of quasi-synchro sampling in elder generation, 1988, (2): 2-7.].
For the sum of all weighting coefficients.
The drift μ of signal frequency is the fixation according to sampling number N in neighbouring sample point fundamental wave phase angle difference and ideal period
Relationship and obtain, the drift μ of signal frequency can also be used for amendment fundamental wave and higher hamonic wave frequency f1With the frequency of higher hamonic wave
Rate fk(fk=k μ fs/N)。
The present invention has the effect of positive: (1) high-precision humorous phase angle analyzes result.The analysis reality such as given for Fig. 1
Example, the analysis precision that the present invention obtains are increased to 10-8Grade (Fig. 2).
(2) method of the present invention fundamentally solves the problems, such as that the humorous phase angle analysis precision of plesiochronous DFT is low, and nothing
Complicated inverting and amendment need to be carried out, algorithm is simple.
(3) relative to plesiochronous DFT, frequency analysis technology of the present invention only needs to increase a sampled point and just solves
Plesiochronous DFT analytical error big problem, it is easy to accomplish.
(4) existing instrument and equipment is improved using the present invention, be technically feasible, and do not need to increase any hard
Part expense can be such that analysis result can be improved to 10-8Grade.
(5) this method is similarly also applied for carrying out successive ignition rather than the frequency analysis process of an iteration, at this time only
Needing an iteration to resolve into successive ignition realization can.As an iteration with successive ignition is substantially, only
It is that when calculating, successive ignition carries out decoupled method, and an iteration is that the process of successive ignition is merged into iteration coefficient γi
In once calculate complete, so the present disclosure applies equally to successive ignition processes.
Detailed description of the invention
Fig. 1 is the humorous phase angle analytical error figure of plesiochronous DFT.
Fig. 2 is humorous phase angle analytical error figure of the invention.
Specific embodiment
(embodiment 1)
A kind of humorous phase angle analysis method based on linear correction algorithm of the present embodiment, comprising the following steps:
Firstly, W+2 sampled point of equal interval sampling, with obtain analyzed signal discrete series f (i), i=0,1 ...,
Wq+1}.The value of W is codetermined by sampling number N in integration method, the number of iterations n and ideal period.
Equal interval sampling refers to: according to the frequency for the ideal signal for carrying out frequency analysis, (such as power frequency component frequency f is
50Hz, period 20mS) determine sample frequency fS=Nf, in sample frequency fSUnder the action of equably sample N in one cycle
Point.Generally, periodic sampling point N=64 or more than can obtain preferable frequency analysis as a result, and the number of iterations n=3~5 just
Comparatively ideal frequency analysis result can be obtained.
Integration method has muiltiple-trapezoid integration method, complexification rectangular integration method, Simpson's method etc. a variety of, can basis
Actual conditions are selected.According to muiltiple-trapezoid integration method, then W=nN；According to complexification rectangular integration method, then W=n
(N-1)；According to iterative Simpson integration method, then W=n (N-1)/2.
Secondly, applying plesiochronous DFT formula since sampled point i=0
W+1 data of analysis obtain fundamental informationWith
Again, fundamental information is obtained from sampled point i=1 using plesiochronous W+1 data of DFT formula analysisWith
Again, using formula:Calculate the frequency drift μ of signal；
Again, each harmonic information is obtained using plesiochronous W+1 data of DFT formula analysis since sampled point i=0With
Then, using formulaCalculate the width phase angle of each harmonic；
Finally, using formulaThe humorous phase angle of linear amendment each harmonic.
Those skilled in the art it should be appreciated that more than embodiment be intended merely to illustrate the present invention, and
It is not intended as limitation of the invention, the present invention can also be changing into more modes, as long as in connotation model of the invention
In enclosing, variation, the modification of embodiment described above will all be fallen within the scope of claims of the present invention.
Claims (2)
1. a kind of humorous phase angle analysis method based on linear correction algorithm, characterized by comprising: the following steps are included:
(1) equidistant W+2 sample point data of sampling { f (i), i=0,1 ..., W+1 }；Described W+2 sampling number of sampling
According to using iterative Simpson integration method, then W=n (N-1)/2；
(2) plesiochronous DFT formula is applied since sampled point i=0
W+1 data of analysis obtain fundamental informationWith
(3) fundamental information is obtained from sampled point i=1 using plesiochronous W+1 data of DFT formula analysisWith
(4) formula is applied:Calculate the frequency drift μ of signal；
(5) each harmonic information is obtained using plesiochronous W+1 data of DFT formula analysis since sampled point i=0With
(6) formula is appliedCalculate the width phase angle of each harmonic；
(7) formula is appliedThe humorous phase angle of linear amendment each harmonic；
In formula: k is the number for needing the harmonic wave obtained；Sin and cos is respectively sine and cosine functions；And akAnd bkRespectively k
The real and imaginary parts of subharmonic；N is the number of iterations；W is determined by integration method；γiFor a weighting coefficient；For
The sum of all weighting coefficients；F (i) is the ith sample value of analysis waveform；N is sampling number in the period；
In the step (1), equidistant sampling is according to the cycle T and frequency f of the ideal signal for carrying out frequency analysis, at one
N point is sampled in period, i.e. sample frequency is fs=Nf, and N >=64.
2. the humorous phase angle analysis method according to claim 1 based on linear correction algorithm, it is characterised in that: the step
(1) in, described W+2 sample point data of sampling is accordingly selected according to selected integration method, and then basis is adopted
Sample frequency fs=Nf is obtained sample point data sequence { f (i), i=0,1 ..., W+1 }；N is the number of iterations, n >=3；Finally to this
Data sequence carries out frequency analysis.
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 |
|---|---|---|---|
| CN201811015823.6A CN108982966B (en) | 2015-05-19 | 2015-05-19 | Harmonic phase angle analysis method based on linear correction algorithm |
| CN201510258020.3A CN104849551B (en) | 2015-05-19 | 2015-05-19 | A kind of humorous phase angle analysis method |
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 true CN108982966A (en) | 2018-12-11 |
| CN108982966B CN108982966B (en) | 2020-10-09 |
Family
ID=53849342
Family Applications (4)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| 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 |
| 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 |
Family Applications Before (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| 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 (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811014781.4A Active CN109030942B (en) | 2015-05-19 | 2015-05-19 | Harmonic phase angle analysis method |
Country Status (1)
| Country | Link |
|---|---|
| CN (4) | CN109142866B (en) |
Families Citing this family (2)
| 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 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006276006A (en) * | 2005-03-01 | 2006-10-12 | Nagoya Institute Of Technology | Harmonic analysis method in power system |
| CN102331526A (en) * | 2011-06-10 | 2012-01-25 | 中国矿业大学 | Method for acquiring parameters of electric power harmonic waves by using Hanniing window function continuous frequency spectrum interpolation |
| CN102393488A (en) * | 2011-08-24 | 2012-03-28 | 江苏技术师范学院 | Harmonic analysis method |
| CN102508026A (en) * | 2011-10-21 | 2012-06-20 | 江苏技术师范学院 | Harmonic wave analysis method for electric energy quality harmonic wave analyzer |
| 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 |
| CN103576002A (en) * | 2013-11-11 | 2014-02-12 | 华北电力大学（保定） | Method for calculating dielectric loss angle of compatible insulating device |
| CN104181391A (en) * | 2014-08-21 | 2014-12-03 | 艾德克斯电子（南京）有限公司 | Harmonic detection method of digital power meter |
| JP5705102B2 (en) * | 2011-12-21 | 2015-04-22 | 三菱電機株式会社 | Insulation deterioration diagnosis device |
Family Cites Families (13)
| 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 |
| DE69916151T2 (en) * | 1998-01-30 | 2005-04-21 | Matsushita Electric Ind Co Ltd | Data transmission with a composition of modulation types |
| 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 |
| CN102435860B (en) * | 2011-10-21 | 2013-11-06 | 江苏理工学院 | Working method of medium loss current tester |
| CN103604989A (en) * | 2011-10-21 | 2014-02-26 | 蒋春花 | Harmonic wave analysis method of power-quality harmonic wave analyzer |
| CN103439566B (en) * | 2011-10-21 | 2014-12-24 | 常州顺创电气科技有限公司 | Operating method of MOA resistive current tester with relatively high precision |
| CN103592512A (en) * | 2011-10-21 | 2014-02-19 | 蒋春花 | Harmonic analysis method of power quality harmonic analyzer |
| CN103969507B (en) * | 2011-10-21 | 2016-10-19 | 江苏理工学院 | A kind of electric energy quality harmonic analyzes method |
| DE102012204963A1 (en) * | 2012-03-28 | 2013-10-02 | Siemens Aktiengesellschaft | Method and apparatus for reducing current harmonics |
| CN103760425A (en) * | 2014-01-22 | 2014-04-30 | 湖南大学 | Method and device for rapidly measuring dielectric loss angle based on time domain quasi-synchronization |
-
2015
- 2015-05-19 CN CN201811014772.5A patent/CN109142866B/en active Active
- 2015-05-19 CN CN201510258020.3A patent/CN104849551B/en active Active
- 2015-05-19 CN CN201811015823.6A patent/CN108982966B/en active Active
- 2015-05-19 CN CN201811014781.4A patent/CN109030942B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006276006A (en) * | 2005-03-01 | 2006-10-12 | Nagoya Institute Of Technology | Harmonic analysis method in power system |
| CN102331526A (en) * | 2011-06-10 | 2012-01-25 | 中国矿业大学 | Method for acquiring parameters of electric power harmonic waves by using Hanniing window function continuous frequency spectrum interpolation |
| CN102393488A (en) * | 2011-08-24 | 2012-03-28 | 江苏技术师范学院 | Harmonic analysis method |
| CN103454494A (en) * | 2011-08-24 | 2013-12-18 | 常州顺创电气科技有限公司 | High-precision harmonic analysis method |
| CN102508026A (en) * | 2011-10-21 | 2012-06-20 | 江苏技术师范学院 | Harmonic wave analysis method for electric energy quality harmonic wave analyzer |
| 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 |
| CN103576002A (en) * | 2013-11-11 | 2014-02-12 | 华北电力大学（保定） | Method for calculating dielectric loss angle of compatible insulating device |
| CN104181391A (en) * | 2014-08-21 | 2014-12-03 | 艾德克斯电子（南京）有限公司 | Harmonic detection method of digital power meter |
Non-Patent Citations (2)
| Title |
|---|
| XIANZHONG DAI ET.AL.: "Quasi-synchronous sampling algorithm and its applications 3．high accurate measurement of frequency，frequency deviation and phase angle difference in power systems", 《PROCEEDINGS OF THE 1993 IEEE INSTRUMENTATION》 * |
| 傅中君等: "基于准同步DFT的非整数谐波分析算法", 《仪器仪表学报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108982966B (en) | 2020-10-09 |
| CN109030942A (en) | 2018-12-18 |
| CN109142866B (en) | 2020-10-09 |
| CN109142866A (en) | 2019-01-04 |
| CN109030942B (en) | 2020-12-04 |
| CN104849551A (en) | 2015-08-19 |
| CN104849551B (en) | 2018-11-23 |
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 | |
| CN102435860B (en) | Working method of medium loss current tester | |
| CN109752584B (en) | Method for measuring effective value of periodic signal | |
| CN103983849B (en) | A kind of Electric Power Harmonic Analysis method of real-time high-precision | |
| CN104849569B (en) | A kind of dielectric loss measurement method | |
| CN103543331B (en) | A kind of method calculating electric signal harmonic wave and m-Acetyl chlorophosphonazo | |
| CN104849551B (en) | A kind of humorous phase angle analysis method | |
| CN103983852B (en) | Harmonic analysis method of power quality harmonic analyzer | |
| Salcic et al. | An improved Taylor method for frequency measurement in power systems | |
| CN108957097A (en) | The measurement method of the resistive current first harmonics of metal oxide arrester | |
| CN110083891A (en) | A kind of calculation method of the plesiochronous linear correction factor of DFT amplitude | |
| CN103969507B (en) | A kind of electric energy quality harmonic analyzes method | |
| CN105277780A (en) | Voltage zero-crossing point obtaining method and device | |
| CN103604989A (en) | Harmonic wave analysis method of power-quality harmonic wave analyzer | |
| CN103592512A (en) | Harmonic analysis method of power quality harmonic analyzer | |
| CN103267896B (en) | The measuring method of initial phase angle of periodic signal | |
| CN109490625B (en) | Harmonic signal analysis method based on sliding window and semi-definite programming | |
| Qian et al. | Spectral interpolation for frequency measurement at off-nominal frequencies | |
| CN109683016A (en) | A kind of harmonic amplitude analysis method | |
| CN110879311A (en) | Method for calculating power of power grid | |
| CN104237592B (en) | A kind of distributed Compare System of electrical measurement error and comparison, source analysis method | |
| CN106033098A (en) | AC power measuring method and device based on data collection | |
| CN117743817A (en) | Multiple time-frequency synchronous extrusion method based on short-time fractional Fourier transform | |
| Palkó et al. | Enhanced Spectral Estimation Using FFT in Case of Data Loss |
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 |