CN110967556A - Real-time harmonic detection method based on feedback neural network - Google Patents
Real-time harmonic detection method based on feedback neural network Download PDFInfo
- Publication number
- CN110967556A CN110967556A CN201911087258.9A CN201911087258A CN110967556A CN 110967556 A CN110967556 A CN 110967556A CN 201911087258 A CN201911087258 A CN 201911087258A CN 110967556 A CN110967556 A CN 110967556A
- Authority
- CN
- China
- Prior art keywords
- harmonic
- neural network
- matrix
- weight
- feedback
- 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.)
- Withdrawn
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 28
- 238000013528 artificial neural network Methods 0.000 title claims abstract description 27
- 239000011159 matrix material Substances 0.000 claims abstract description 46
- 230000005284 excitation Effects 0.000 claims abstract description 8
- 238000005070 sampling Methods 0.000 claims description 11
- 238000000034 method Methods 0.000 claims description 9
- 238000005457 optimization Methods 0.000 claims description 8
- 230000000737 periodic effect Effects 0.000 claims description 4
- 230000017105 transposition Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000004422 calculation algorithm Methods 0.000 abstract description 6
- 238000013461 design Methods 0.000 abstract description 3
- 230000007547 defect Effects 0.000 abstract description 2
- 238000004458 analytical method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 210000002569 neuron Anatomy 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011897 real-time detection Methods 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computing Systems (AREA)
- Computational Linguistics (AREA)
- Biomedical Technology (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Testing Electric Properties And Detecting Electric Faults (AREA)
Abstract
The invention discloses a real-time harmonic detection method based on a feedback neural network, which comprises the steps of constructing a sine basis function feedback neural network structure, designing an excitation function design and a network weight matrix suitable for harmonic detection, constructing a quadratic programming problem, obtaining the optimal weight of the sine basis function weight neural network by operating the feedback neural network, calculating the accurate amplitude and phase of fundamental waves and each subharmonic according to the estimated value of the weight matrix, and reconstructing the fundamental waves and each subharmonic. The real-time harmonic detection method overcomes the defects of complex network structure, multiple iteration times, low algorithm precision and the like of the traditional harmonic detection method, and can obtain a network containing harmonic amplitude and phase information only through a simple feedback network; the detection efficiency is high, the precision is high, and the network structure is simple.
Description
Technical Field
The invention relates to the field of power systems, in particular to a real-time harmonic detection method based on a feedback neural network.
Background
A large number of nonlinear devices exist in the power system to rectify and invert the voltage and current of a power grid to generate time-varying harmonic waves, so that the problem of power harmonic pollution is increasingly serious, the quality of electric energy is seriously influenced, and meanwhile, the safety, stability and high-efficiency operation of the power system are threatened; moreover, due to the randomness generated by the harmonic waves of the power system, the harmonic waves are seriously influenced by the nonlinear complexity of the power grid, and the difficulty in real-time detection of the harmonic waves of the power grid is high, so that the real-time accurate detection of the power harmonic waves at the present stage has important significance.
The traditional harmonic measurement method based on analog filter is eliminated because the ability to resist harmonic distortion rate is large and have voltage with additional phase shift is too weak; the Fourier transform method has the defects of frequency spectrum leakage, barrier effect, incapability of detecting non-stationary harmonic waves and the like; the wavelet detection method has good advantages in detecting transient signals or singularities of signals of a power grid, but signal decomposition performed by the method can cause interleaving between high-pass filter groups and low-pass filter groups, a mixing aliasing phenomenon occurs, and the problem that window energy is not concentrated exists. Other harmonic detection methods based on an intelligent optimization algorithm exist, but the optimization algorithm evolution is carried out after time-frequency transformation is carried out, so that the complexity is high.
The basic idea of the neural network algorithm is to adopt a physically realizable system to simulate the structure and the functional system of human brain nerve cells, the algorithm has good nonlinear expression capability, parallel processing capability, strong robustness and self-organizing self-learning capability, and is widely applied to the fields of signal processing and pattern recognition, the harmonic detection problem is equivalent to a signal detection problem, the feedback neural network has a simple structure and is convenient to realize, and the invention tries to solve the harmonic detection problem by using the feedback neural network method.
Disclosure of Invention
The invention aims to provide a real-time harmonic detection method based on a feedback neural network, which enhances the real-time performance of harmonic detection and achieves the aim of detecting each harmonic.
In order to solve the technical problems, the technical scheme adopted by the invention is as follows: a real-time harmonic detection method based on a feedback neural network comprises the following steps:
1) the periodic signal containing each harmonic in the power electronic system is expressed as:
wherein, w0At fundamental angular frequency, w 02 pi f, f is the fundamental frequency, k is the harmonic order, akAndamplitude and phase of the kth harmonic, respectively; m is the highest harmonic number; discretizing the above expression and expanding the discretized expression into a matrix form, then tiEach sample value is represented asWherein,Tsthe sampling period is generally not more than 100 milliseconds;(ii) a T represents a matrix transposition operation; w is the grid angular frequency;
2) designing an excitation function g (-) of an input and output layer of the feedback neural network: g (x) ═ (x + asin (pi x)); wherein a is more than 0, x is more than infinity and pi is a circumferential rate; sampling value signal x (t)i) As an input signal of the excitation function, an output signal g (x (t) is obtainedi));
3) Establishing a sine basis weight function matrix as follows:
wherein R represents a real number domain;
4) the following quadratic programming optimization problem is constructed:
wherein: represents a 2-norm; d ═ diag [ d (0), d (1), …, d (n)]Diag is a diagonal matrix with main diagonal elements of d (0), d (1), … and d (n), and any other elements are 0, α is a normal number, α ∈ (0,1)]P is the total number of times of feedback of the neural network; w is a weight matrix, and W is a weight matrix,
5) solving the quadratic programming optimization problem to obtain an estimated value of a weight matrix WAnd calculating the accurate amplitude and phase of the fundamental wave and each harmonic according to the estimated value of the weight matrix W, and reconstructing the fundamental wave and each harmonic.
The specific calculation process of the estimated value of the weight matrix W comprises the following steps: order to
a) Derivative W and make its value equal to the zero vector:
b) let the r-th iteration fjWhere j is 0, …, n is fixed, let the first derivative be equal to 0, and a unique solution for W is obtained, then:wherein the superscript-1 represents the matrix inversion operation,is a main diagonal element of | x (j) & ltY2A diagonal matrix of (a); j ═ t1,t2,…,tn;
c) Setting an iterative formula of the r +1 th feedback network weight: w (r +1) ═ 1- μ g (W (r)) + μ W (0); wherein, mu belongs to (0,1), and W (0) is the initial value of the weight;
d) and repeating the steps a) to c) until the value of the cost function J (W) is not reduced any more, and obtaining a weight matrix, thereby obtaining an estimated value of the weight matrix W.
Compared with the prior art, the invention has the beneficial effects that: the invention can enhance the real-time performance of harmonic detection, avoids the complex harmonic detection network structure in the traditional harmonic detection method, obtains the network weight containing harmonic amplitude and phase information through network self-operation, and further extracts the amplitude and phase of fundamental wave and each subharmonic from the weight matrix, thereby achieving the purpose of detecting each subharmonic; the invention has the advantages of high detection precision, good measurement accuracy, strong real-time property and good anti-interference performance.
Drawings
FIG. 1 shows the present invention in analyzing signals asThe fundamental frequency is 50Hz, harmonic waves and phases are randomly generated, and a oscillogram is obtained by analysis under the condition of white Gaussian noise of 10 decibels;
FIG. 2 shows the present invention in analyzing signals asThe fundamental frequency is 50Hz, harmonic waves and phases are randomly generated, and an input and output integral comparison graph is obtained by analysis under the condition of 10 dB Gaussian white noise;
FIG. 3 shows the present invention in analyzing signals asThe fundamental frequency is 50Hz, harmonic waves and phases are randomly generated, and an input and output amplitude error graph is obtained by analyzing under the condition of 10 dB Gaussian white noise.
Detailed Description
The invention provides a real-time harmonic detection method based on a feedback neural network, which comprises the following steps:
(1) according to the existing method, a periodic signal model containing each harmonic in a power electronic system is expanded by adopting a trigonometric function and a difference product formula, so that the periodic signal containing each harmonic in the power electronic system can be expressed as
Wherein, wiAngular frequency of the ith harmonic, w 02 pi f, f being the fundamental frequency, AiAndamplitude and phase of the ith harmonic, respectively(ii) a N is the highest harmonic order and m is the highest harmonic order. Further expressing equation (1) as a trigonometric function transformation equation
Wherein: w is a0Is the fundamental angular frequency, j is the harmonic order, TsIs the sampling period. Discretizing and expanding the formula (2) into a matrix form, tthiA sampling value x (t)i) This can be expressed as follows:
(2) Setting an excitation function g (·), g (x) ═ x + asin (pi x)) of the input and output layers of the neural network; wherein a is more than 0, and x is less than + ∞; pi is the circumference ratio; sampling value signal x (t)i) I-1, 2, …, n as input signal for the excitation function.
(3) Designing a sine basis function neural network weight and constructing a basis function weight matrix;
slave matrixThe amplitude and phase information of the harmonic waves are contained in the matrixAnd the matrix is independent of any signal input sample value; continuing to look at the matrix may find that its elements are composed of Composition if applicableTo accurately obtain the values of the elements of the matrix, the amplitude and the corresponding phase of the ith harmonic can be accurately obtained according to the correlation property of the trigonometric function.
Accordingly, we design the weight matrix of the neural network asIn the form of weight matricesAnd (4) finishing.
Thus, once the weight matrix is obtained, the fundamental amplitude can be calculated asPhase isSimilarly, the amplitude and phase of the obtained kth harmonic are respectivelyThereby reconstructing fundamental wave and each harmonic wave.
Because of CiIs a matrix whose element values are related only to the samples. CiThe function of the method is similar to that of performing function mapping on input sampling points to obtain function outputs, and then the matrix elements can be used as sine triangular basis functions.
Then, according to the input variation of the sampling point, the following sine basis function matrix can be constructed
Wherein: r represents a real number domain, 2m is the number of hidden layer neurons, n is the number of samples,g(x(ti) Is the output signal after the input signal sample has been subjected to the excitation function, i ═ 1,2, …, n;
the design of the neural network is completed as above.
(4) Designing a quadratic programming optimization problem to obtain an estimated value of a weight matrix w;
specifically, the following quadratic programming optimization problem is constructed
Wherein:is a sampling vector of the harmonic input signal, n is a sampling number, the superscript T represents a transposition operation,representing a 2-norm. d ═ diag [ d (0), d (1), …, d (n)]Diag is the main diagonal elements d (0), d (1), …, d (n), and any other elements are diagonal arrays of 0, where And P is the total number of times of the neural network feedback. The purpose of the introduced sigma is to recursively minimize wTAnd w, which does not participate in the algorithm iteration, but gradually reduces the order of the norm by "2" at each feedback.
(5) Derivative W and make its value equal to zero vector
Let, i.e. the kth iteration fjJ is 0, …, n is fixed, let the first derivative equal to 0, and obtain its unique solution, then there is
In the formula: the superscript-1 represents the matrix inversion operation,is formed by rendering a main diagonal element as | x (j)2J is 0, …, n.
(6) Iterative formula for setting weight of (k +1) th feedback network
W(k+1)=(1-μ)g(W(k))+μW(0) (7)
Where μ e (0,1), w (0) is the initial value of the weight, and is generally set to 0.05+ j0.05 for the center tap, and all the other tap values are 0.
The feedback neural network is operated until the value of the cost function J (W) is not reduced any more, namely the feedback network is considered to be converged, and therefore the network weight matrix is obtained.
Further, calculating the accurate amplitude and phase of the fundamental wave and each harmonic according to the estimated value of the weight matrix W, and reconstructing the fundamental wave and each harmonic; in particular, the fundamental amplitude isPhase isSimilarly, the amplitude and phase of the obtained kth harmonic are respectively The fundamental wave and each harmonic can be reconstructed accordingly.
The examples show that: suppose the analysis signal isThe fundamental frequency is 50Hz, the harmonic and the phase are randomly generated, and the noise is white Gaussian noise with 10 dB
FIG. 1 is a graph of waveforms obtained during analysis according to the present invention. Fig. 2 is a comparison of the input and output obtained by the method of the invention. Fig. 3 and 4 are graphs of amplitude and phase errors of input and output subharmonics, respectively, obtained by the method of the present invention. The result shows that the invention can obtain very good harmonic detection effect, has small detection error and high accuracy and has very good application value.
Claims (5)
1. A real-time harmonic detection method based on a feedback neural network is characterized by comprising the following steps:
1) the periodic signal containing each harmonic in the power electronic system is expressed as:
wherein, w0At fundamental angular frequency, w02 pi f, f is the fundamental frequency, k is the harmonic order, akAndamplitude and phase of the kth harmonic, respectively; m is the highest harmonic number; discretizing the above expression and expanding the discretized expression into a matrix form, then tiEach sample value is represented asWherein, i is 1,2, …, n,Tsis a sampling period; (ii) a T represents a matrix transposition operation; w is the grid angular frequency;
2) designing an excitation function g (-) of an input and output layer of the feedback neural network: g (x) ═ (x + asin (pi x)); wherein a is more than 0, x is more than infinity and pi is a circumferential rate; sampling value signal x (t)i) As an input signal of the excitation function, an output signal g is obtained(x(ti));
3) Establishing a sine basis weight function matrix as follows:
wherein R represents a real number domain;
4) the following quadratic programming optimization problem is constructed:
wherein: represents a 2-norm; d ═ diag [ d (0), d (1), …, d (n)]Diag is a diagonal matrix with main diagonal elements of d (0), d (1), … and d (n), and any other elements are 0, l is 1,2, …,2m, α is a normal number, α ∈ (0,1)]P is the total number of times of feedback of the neural network; w is a weight matrix, and W is a weight matrix,
5) solving the quadratic programming optimization problem to obtain an estimated value of a weight matrix WAnd calculating the accurate amplitude and phase of the fundamental wave and each harmonic according to the estimated value of the weight matrix W, and reconstructing the fundamental wave and each harmonic.
2. The method according to claim 1, wherein the specific calculation of the estimated value of the weight matrix W comprises: order to
a) Derivative W and make its value equal to the zero vector:
b) let the r-th iteration fjWhere j is 0, …, n is fixed, let the first derivative be equal to 0, and a unique solution for W is obtained, then:wherein the superscript-1 represents the matrix inversion operation,is a main diagonal element of | x (j) & ltY2A diagonal matrix of (a); j ═ t1,t2,…,tn;
c) Setting an iterative formula of the r +1 th feedback network weight: w (r +1) ═ 1- μ g (W (r)) + μ W (0); wherein, mu belongs to (0,1), and W (0) is the initial value of the weight;
d) and repeating the steps a) to c) until the value of the cost function J (W) is not reduced any more, and obtaining a weight matrix, thereby obtaining an estimated value of the weight matrix W.
5. The feedback neural network-based real-time harmonic detection method of claim 1, wherein T iss≤100ms。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911087258.9A CN110967556A (en) | 2019-11-08 | 2019-11-08 | Real-time harmonic detection method based on feedback neural network |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911087258.9A CN110967556A (en) | 2019-11-08 | 2019-11-08 | Real-time harmonic detection method based on feedback neural network |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110967556A true CN110967556A (en) | 2020-04-07 |
Family
ID=70030492
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911087258.9A Withdrawn CN110967556A (en) | 2019-11-08 | 2019-11-08 | Real-time harmonic detection method based on feedback neural network |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110967556A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112787491A (en) * | 2020-12-28 | 2021-05-11 | 中南大学 | Input current harmonic suppression method of three-stage AC/DC power supply |
CN114252700A (en) * | 2021-10-26 | 2022-03-29 | 深圳市锐风电子科技有限公司 | Power harmonic detection method based on sine and cosine algorithm |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102353839A (en) * | 2011-07-18 | 2012-02-15 | 华北电力大学(保定) | Electric power system harmonics analysis method based on multilayered feedforward neural network |
CN103383413A (en) * | 2013-07-09 | 2013-11-06 | 温州大学 | Real-time harmonic detection method based on direct weight determination method |
CN104833852A (en) * | 2015-05-11 | 2015-08-12 | 重庆大学 | Power system harmonic signal estimation and measurement method based on artificial neural network |
CN105305446A (en) * | 2015-10-22 | 2016-02-03 | 南京亚派科技股份有限公司 | Harmonic current tracking method based on intelligent control |
CN109581054A (en) * | 2018-11-23 | 2019-04-05 | 温州晶彩光电有限公司 | A kind of real-time harmonic rapid detection method of bank base conversion power supply system peculiar to vessel |
-
2019
- 2019-11-08 CN CN201911087258.9A patent/CN110967556A/en not_active Withdrawn
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102353839A (en) * | 2011-07-18 | 2012-02-15 | 华北电力大学(保定) | Electric power system harmonics analysis method based on multilayered feedforward neural network |
CN103383413A (en) * | 2013-07-09 | 2013-11-06 | 温州大学 | Real-time harmonic detection method based on direct weight determination method |
CN104833852A (en) * | 2015-05-11 | 2015-08-12 | 重庆大学 | Power system harmonic signal estimation and measurement method based on artificial neural network |
CN105305446A (en) * | 2015-10-22 | 2016-02-03 | 南京亚派科技股份有限公司 | Harmonic current tracking method based on intelligent control |
CN109581054A (en) * | 2018-11-23 | 2019-04-05 | 温州晶彩光电有限公司 | A kind of real-time harmonic rapid detection method of bank base conversion power supply system peculiar to vessel |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112787491A (en) * | 2020-12-28 | 2021-05-11 | 中南大学 | Input current harmonic suppression method of three-stage AC/DC power supply |
CN114252700A (en) * | 2021-10-26 | 2022-03-29 | 深圳市锐风电子科技有限公司 | Power harmonic detection method based on sine and cosine algorithm |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112395992B (en) | Power harmonic signal denoising method based on improved wavelet threshold | |
CN107102255B (en) | Single ADC acquisition channel dynamic characteristic test method | |
CN106597260A (en) | Simulation circuit fault diagnosis method based on continuous wavelet analysis and ELM network | |
CN106845010A (en) | Based on the low-frequency oscillation dominant pattern discrimination method for improving SVD noise reductions and Prony | |
CN106646121B (en) | A kind of discrimination method of distribution network failure wavefront | |
del Rıo et al. | Comparison between the matrix pencil method and the Fourier transform technique for high-resolution spectral estimation | |
Belega et al. | Amplitude and phase estimation of real-valued sine wave via frequency-domain linear least-squares algorithms | |
CN110967556A (en) | Real-time harmonic detection method based on feedback neural network | |
CN103941072A (en) | Power signal catastrophe parameter measurement method based on real strong tracking filter | |
CN106788060A (en) | Based on the generator parameter identification method for improving SVD Prony algorithms | |
CN111289796A (en) | Method for detecting subsynchronous oscillation of high-proportion renewable energy power system | |
Wang et al. | Two points interpolated DFT algorithm for accurate estimation of damping factor and frequency | |
CN103383413A (en) | Real-time harmonic detection method based on direct weight determination method | |
CN112816779B (en) | Harmonic real signal parameter estimation method for analytic signal generation | |
CN112557751A (en) | Harmonic parameter estimation method based on DFT iteration method | |
Yi et al. | Instantaneous harmonic decomposition technique for three‐phase current based on multiple reference coordinates | |
Xianmin | A new method with high confidence for validation of computer simulation models of flight systems | |
Ruan et al. | Improved Prony method for high-frequency-resolution harmonic and interharmonic analysis | |
CN108917917A (en) | A kind of breaker mechanical vibration signal denoising method based on sparse decomposition | |
CN109581054A (en) | A kind of real-time harmonic rapid detection method of bank base conversion power supply system peculiar to vessel | |
CN114646803A (en) | Harmonic and inter-harmonic frequency estimation method based on parallel factor compressed sensing | |
Liu et al. | A method for blind source separation of multichannel electromagnetic radiation in the field | |
CN111398735A (en) | Transformer substation grounding grid fault detection method based on information entropy | |
CN114258495B (en) | Finite complex signal measurement system and high-precision decomposition method | |
CN115015633B (en) | Frequency estimation method for harmonic wave and inter-harmonic wave in power system |
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 | ||
WW01 | Invention patent application withdrawn after publication | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20200407 |