CN112698092A - Rapid broadband measuring device and method based on ESPRIT algorithm - Google Patents

Rapid broadband measuring device and method based on ESPRIT algorithm Download PDF

Info

Publication number
CN112698092A
CN112698092A CN202011461130.7A CN202011461130A CN112698092A CN 112698092 A CN112698092 A CN 112698092A CN 202011461130 A CN202011461130 A CN 202011461130A CN 112698092 A CN112698092 A CN 112698092A
Authority
CN
China
Prior art keywords
signal
data
module
matrix
sampling
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202011461130.7A
Other languages
Chinese (zh)
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.)
State Grid Corp of China SGCC
State Grid Liaoning Electric Power Co Ltd
Original Assignee
State Grid Liaoning Electric Power Co Ltd
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 State Grid Liaoning Electric Power Co Ltd filed Critical State Grid Liaoning Electric Power Co Ltd
Priority to CN202011461130.7A priority Critical patent/CN112698092A/en
Publication of CN112698092A publication Critical patent/CN112698092A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/02Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage

Abstract

The invention discloses a rapid broadband measuring device and method based on an ESPRIT algorithm. The rapid sampling and processing module adopts a rapid broadband measurement method based on ESPRIT. And the signal frequency component is estimated through the kurtosis of the singular value of the signal, so that the performance under the condition of a low signal-to-noise ratio is improved while the fast-changing broadband signal has short response time and high resolution. The practical broadband measurement technology is promoted, and a new monitoring means is provided for the power electronic power grid.

Description

Rapid broadband measuring device and method based on ESPRIT algorithm
Technical Field
The invention relates to a broadband measuring device, in particular to a rapid broadband measuring device and a rapid broadband measuring method based on an ESPRIT algorithm.
Background
With the development and utilization of large-scale renewable energy sources and the development of smart power grids, China currently builds an ultra-large-scale complex interconnected power system. The high permeability power electronic device makes the dynamic characteristics of the power system signal become more and more complex, and the frequency band range presents the characteristic of wide frequency
Power electronics can also cause frequency components up to 300Hz to the system, which also causes new grid stability problems. And a large amount of field wave recording data analysis shows that the frequency of the 0-300Hz component is dense, and the frequency and the amplitude can change rapidly along with time (the amplitude can change by 50 percent within 1 s). In addition, more and more power electronic equipment is put into operation, and the operation of the equipment improves the flexibility and the reliability of the operation of the power grid, but also generates a large amount of harmonic pollution to influence the power quality of the power grid. Power electronics also generate a large number of higher harmonics (up to 50 harmonics), which can cause severe harmonic pollution in power systems. Therefore, it is necessary to perform real-time synchronous measurement on the broadband signal of 0-2500Hz to provide data for research on broadband signal source, propagation path, and security control.
In the aspect of a measurement method of a broadband signal, Fast Fourier Transform (FFT) and an improved algorithm thereof are widely used for electrical quantity spectrum analysis of a power system due to the advantages of small calculated amount, easiness in hardware implementation and the like. However, the frequency resolution of the FFT is proportional to the time window length, i.e., a longer time window is required to achieve higher frequency resolution. This results in that it cannot take into account both the measurement accuracy and the tracking speed when measuring the fast changing frequency component of 0-300 Hz. And the method has a good measuring effect on integer subharmonics with large frequency intervals and relatively stable frequency intervals. In addition, wavelet transform and Prony algorithm are also common spectrum analysis methods, but these 2 methods have poor stability and cannot maintain good performance under the condition of low signal-to-noise ratio. The rotation invariant subspace (ESPRIT) algorithm is a spatial spectrum estimation method, which utilizes the orthogonal characteristics of the signal subspace and the noise subspace of a sampling signal, can ensure a shorter time window compared with the method while ensuring a higher frequency resolution, can better track the signal fast dynamic process and reduce the averaging effect of the time window, and has a good application prospect. However, this method is susceptible to noise when dividing the signal subspace and the noise subspace, which affects the performance of the algorithm.
Disclosure of Invention
The invention aims to provide a rapid broadband measuring device and a measuring method based on an ESPRIT algorithm. The method can accurately measure the rapidly-changed broadband signal and rapidly track the rapidly-changed broadband signal. The device has the advantages of short response time and high resolution, and meanwhile, the measurement accuracy of the device under the condition of low signal-to-noise ratio is improved.
In order to solve the existing technical problems. The technical scheme adopted by the invention is as follows:
the fast broadband measuring device based on the ESPRIT algorithm comprises a power supply module, a CPU module, a fast sampling and processing module, a broadband oscillation detection module, a high-speed bus backboard, a man-machine interface module and a data transmission module, wherein the power supply module is used for providing working voltage required by the device; the high-speed bus back board is an access board for all the plug-ins and is used for transmitting data among the plug-ins; the CPU module is inserted and connected on a high-speed bus backboard, converts a small signal sent by the alternating current input module into a digital signal through an analog-to-digital conversion chip AD to finish current and voltage signal sampling, calculates fundamental wave, harmonic wave, inter-harmonic wave, synchronous phasor and power of voltage and current according to the digital sampling signal of the current and voltage, and sends the voltage and the current to the liquid crystal module through the high-speed bus backboard for displaying; the man-machine interface module exchanges data with other modules through a high-speed bus backboard; the rapid sampling and processing module adopts an ESPRIT-based rapid broadband measurement algorithm; the broadband oscillation detection module collects data including a dominant component, an equivalent value and a transient energy flow. The data transmission module is divided into two links of an in-station link and a master station link. And transmitting measurement data, alarm events and file data to the in-station monitoring system, wherein various data are modeled according to a DL/T860 standard. Phasor, inter-harmonic waves and alarm events are transmitted to the master station in real time according to the GB/T26865.2 standard, and a recording wave data file is uploaded. And each dominant component of the inter-harmonic waves and the oscillation power is carried by real-time data frame time division multiplexing, so that compatible transmission with synchronous phasor data is realized.
Further, the man-machine interface module comprises a liquid crystal and a keyboard.
The ESPRIT algorithm of the rapid broadband measuring device based on the ESPRIT algorithm comprises the following steps: firstly, a dynamic model of a broadband multi-frequency signal in the power system is composed of k frequency components and noise components at time t, namely
Figure BDA0002831831900000031
In the formula: f. ofi,xmi(t) and θi(t) are each si(t) has a frequency, amplitude and phase angle of xm1(t)≥xm2(t)≥…≥xmk(t) of (d). Arranging the sampled data of the voltage or current signals into a sampled data matrix X, and arranging the sampled data of the voltage or current signals into a matrix X consisting of the first L-1 elements of the n +1 th column vector of X because each sampled value can be regarded as the superposition of the instantaneous value of each frequency component and the noise componentn1=[x(n),x(n+1),…,x(n+L-2)]TCan be expressed as:
Figure BDA0002831831900000032
in the formula: ts is the sampling interval.
The formula (2) is simplified into a matrix form
Xn1=A1S+N1 (3)
In the formula: a. the1Representing the phase angle difference of each frequency component of the reference element when the sampling elements meet the sampling characteristic matrix; s is a column vector consisting of each frequency component; n is a radical of1Is a noise vector.
For X, the same principle appliesnThe matrix X composed of the last L-1 elementsn2Comprises the following steps:
Xn2=A2S+N2 (4)
in the formula: a. the2Is a sampling characteristic matrix; n is a radical of2Is a noise vector.
Thirdly, the corresponding amplitude x can be obtained by the least square method through the frequency of each component of the signalmiAngle of sum
Figure BDA0002831831900000041
Matrix parameters are defined as shown in equations (5) and (6):
Figure BDA0002831831900000042
Xs=[x(0),x(1),…,x(L+M-2)]T (6)
fourthly, according to the least square method, a matrix S containing phasor information of each frequency component of the signal is
S=(WTW)-1WTXs=[s1,s2,…,sk]T (7)
The amplitude and phase angle corresponding to each component of the signal are:
Figure BDA0002831831900000043
and fifthly, in order to improve the measurement performance of the ESPRIT algorithm under the condition of containing noise, a signal frequency component estimation method based on kurtosis is provided.
Performing SVD processing on the L multiplied by M sampling data matrix X:
X=UPV (9)
in the formula: p is singular value matrix, P ═ diag (P)1,p2,…,pL)。
Sixthly, regarding the singular value of X from small to large as a discrete signal, the kurtosis K of the mth singular valuemThe expression is as follows:
Figure BDA0002831831900000044
in the formula: p ═ pm-n+1,pm-n,…,pmThe m-n +1 th to m-th singular values of the sampled data matrix X are arranged from small to large;
Figure BDA0002831831900000051
is the mean value of all elements in p; σ is the standard deviation of p, and the expression is shown below.
Figure BDA0002831831900000052
The following equations (10) and (11) can be obtained:
Figure BDA0002831831900000053
when the elements in the group p are noise singular values, the kurtosis is due to the fact that no significant impact part of the noise singular values appears overall
Figure BDA0002831831900000056
Will not generate too large fluctuation and is always at the limit value KaThe following; when m is1+1 singular values into the array P, kurtosis
Figure BDA0002831831900000054
It will increase significantly. Thus by judging kurtosis value KmCan estimate the frequency of the signal multiplied by a fraction k, i.e.
Figure BDA0002831831900000055
The invention has the advantages and beneficial effects that:
the invention relates to a rapid broadband measuring device based on an ESPRIT algorithm, which consists of a power supply module, a CPU module, a rapid sampling and processing module, a broadband oscillation detection module, a high-speed bus back plate, a data transmission module and a man-machine interface module. The fast sampling and processing module adopts an ESPRIT-based fast broadband measurement method, can accurately measure and fast track fast-changing broadband signals, estimates the frequency component of the signals through the kurtosis of singular values of the signals, has short response time and high resolution ratio for the fast-changing broadband signals, and improves the performance under the condition of low signal-to-noise ratio. The practical broadband measurement technology is promoted, and a new monitoring means is provided for the power electronic power grid.
Drawings
Fig. 1 is a schematic structural diagram of a fast broadband measurement device based on the ESPRIT algorithm according to the present invention.
Detailed Description
As shown in fig. 1, the fast broadband measurement device based on the ESPRIT algorithm of the present invention includes a power module 1, a CPU module 2, a fast sampling and processing module 3, a broadband oscillation detection module 4, a high-speed bus backplane 5, a human-computer interface module 6, and a data transmission module 7, wherein the power module is used for providing a working voltage required by the device; the high-speed bus back board is an access board for all the plug-ins and is used for transmitting data among the plug-ins; the system is provided with a man-machine interface, so that a rapid sampling and processing module can be combined with a rapid broadband measurement algorithm of ESPRIT, and accurate measurement and rapid tracking of broadband data signals are realized; the CPU module is inserted and connected on a high-speed bus backboard, converts a small signal sent by the alternating current input module into a digital signal through an analog-to-digital conversion chip AD, completes current and voltage signal sampling, calculates parameters such as voltage, fundamental wave, harmonic wave, inter-harmonic wave, synchronous phasor and power of current according to the digital sampling signal of the current and voltage, and sends the parameters to the liquid crystal module for display through the high-speed bus backboard; the man-machine interface module exchanges data with other modules through a high-speed bus backboard; the rapid sampling and processing module adopts an ESPRIT-based rapid broadband measurement algorithm, and has shorter response time and higher resolution ratio for rapidly-changed broadband signals; the broadband oscillation detection module collects data including a dominant component, an equivalent value and a transient energy flow. Detecting the dominant component according to the frequency resolution of 1Hz to realize the accurate estimation of frequency, amplitude and phase; the equivalent value is the comprehensive evaluation quantity of subsynchronous/supersynchronous current, voltage total energy and distortion power thereof, and is used as out-of-limit alarm quantity to be monitored in real time; the data transmission module is divided into two links of an in-station link and a master station link. And transmitting measurement data, alarm events and file data to the in-station monitoring system, wherein various data are modeled according to a DL/T860 standard. Phasor, inter-harmonic and alarm events are transmitted in real time with the main station according to the GB/T26865.2 standard, and a recording wave data file is uploaded. And each dominant component of (inter) harmonic waves and oscillation power is carried by real-time data frame time division multiplexing, so that compatible transmission with synchronous phasor data is realized. The man-machine interface module comprises a liquid crystal and a keyboard.
The ESPRIT algorithm of the rapid broadband measuring device based on the ESPRIT algorithm comprises the following steps:
the dynamic model of the broadband multi-frequency signal in the power system is composed of k frequency components and noise components at time t, that is:
Figure BDA0002831831900000071
in the formula: f. ofi,xmi(t) and θi(t) are each si(t) has a frequency, amplitude and phase angle of xm1(t)≥xm2(t)≥…≥xmk(t)。
Arranging the sampled data of the voltage or current signals into a sampled data matrix X, and arranging the sampled data of the voltage or current signals into a matrix X consisting of the first L-1 elements of the n +1 th column vector of X because each sampled value can be regarded as the superposition of the instantaneous value of each frequency component and the noise componentn1=[x(n),x(n+1),…,x(n+L-2)]TCan be expressed as:
Figure BDA0002831831900000072
in the formula: ts is the sampling interval.
The equation (2) is simplified to a matrix form:
Xn1=A1S+N1 (3)
in the formula: a. the1For sampling the characteristic matrix, characterizing the meeting of the sampling elements to the reference element at each frequencyRate component phase angle difference; s is a column vector consisting of each frequency component; n is a radical of1Is a noise vector.
For X, the same principle appliesnThe matrix X composed of the last L-1 elementsn2Comprises the following steps:
Xn2=A2S+N2 (4)
in the formula: a. the2Is a sampling characteristic matrix; n is a radical of2Is a noise vector.
Thirdly, the corresponding amplitude x can be obtained by the least square method through the frequency of each component of the signalmiAngle of sum
Figure BDA0002831831900000074
Matrix parameters are defined as shown in equations (5) and (6):
Figure BDA0002831831900000073
Xs=[x(0),x(1),…,x(L+M-2)]T (6)
fourthly, according to a least square method, a matrix S containing phasor information of each frequency component of the signal is as follows:
S=(WTW)-1WTXs=[s1,s2,…,sk]T (7)
the amplitude and phase angle corresponding to each component of the signal are:
Figure BDA0002831831900000081
and fifthly, in order to improve the measurement performance of the ESPRIT algorithm under the condition of containing noise, a signal frequency component estimation method based on kurtosis is provided.
Performing SVD processing on the L multiplied by M sampling data matrix X:
X=UPV (9)
in the formula: p is singular value matrix, P ═ diag (P)1,p2,…,pL)。
Sixthly, the singularity that X is arranged from small to largeThe value is treated as a discrete signal whose m-th singular value has a kurtosis KmThe expression is as follows:
Figure BDA0002831831900000082
in the formula: p ═ pm-n+1,pm-n,…,pmThe m-n +1 th to m-th singular values of the sampled data matrix X are arranged from small to large;
Figure BDA0002831831900000083
is the mean value of all elements in p; σ is the standard deviation of p, and the expression is shown below.
Figure BDA0002831831900000084
The following equations (10) and (11) can be obtained:
Figure BDA0002831831900000085
when the elements in the group p are noise singular values, the kurtosis is due to the fact that no significant impact part of the noise singular values appears overall
Figure BDA0002831831900000093
Will not generate too large fluctuation and is always at the limit value KaThe following; when m is1+1 singular values into the array P, kurtosis
Figure BDA0002831831900000091
It will increase significantly. Thus by judging kurtosis value KmThe frequency multiplied by the fraction k of the signal can be estimated as:
Figure BDA0002831831900000092

Claims (3)

1. quick wide band measuring device based on ESPRIT algorithm, its characterized in that: the system comprises a power module (1), a CPU module (2), a rapid sampling and processing module (3), a broadband oscillation detection module (4), a high-speed bus backboard (5), a man-machine interface module (6) and a data transmission module (7); the power supply module is used for providing working voltage required by the device; the high-speed bus back board is an access board for all the plug-ins and is used for transmitting data among the plug-ins; the CPU module is inserted and connected on a high-speed bus backboard, converts a small signal sent by the alternating current input module into a digital signal through an analog-to-digital conversion chip AD to finish current and voltage signal sampling, calculates fundamental wave, harmonic wave, inter-harmonic wave, synchronous phasor and power of voltage and current according to the digital sampling signal of the current and voltage, and sends the voltage and the current to the liquid crystal module through the high-speed bus backboard for displaying; the man-machine interface module exchanges data with other modules through a high-speed bus backboard; the rapid sampling and processing module adopts an ESPRIT-based rapid broadband measurement algorithm; the broadband oscillation detection module acquires data comprising a dominant component, an equivalent value and a transient energy flow; the data transmission module is divided into an in-station link and a master station link, and transmits measurement data, alarm events and file data to an in-station monitoring system, and various types of data are modeled according to a DL/T860 standard; phasor, inter-harmonic waves and alarm events are transmitted to the master station in real time according to the GB/T26865.2 standard, and a wave recording data file is uploaded; and each dominant component of the inter-harmonic waves and the oscillation power is carried by real-time data frame time division multiplexing, so that compatible transmission with synchronous phasor data is realized.
2. The ESPRIT algorithm-based fast broadband measurement device according to claim 1, wherein: the man-machine interface module comprises a liquid crystal and a keyboard.
3. The measurement algorithm of the ESPRIT algorithm based fast broadband measurement device according to claim 1, comprising the steps of:
the dynamic model of the broadband multi-frequency signal in the power system is composed of k frequency components and noise components at time t, that is:
Figure FDA0002831831890000021
in the formula: f. ofi,xmi(t) and θi(t) are each si(t) has a frequency, amplitude and phase angle of xm1(t)≥xm2(t)≥…≥xmk(t);
Arranging the sampled data of the voltage or current signals into a sampled data matrix X, and arranging the sampled data of the voltage or current signals into a matrix X consisting of the first L-1 elements of the n +1 th column vector of X because each sampled value can be regarded as the superposition of the instantaneous value of each frequency component and the noise componentn1=[x(n),x(n+1),…,x(n+L-2)]TCan be expressed as:
Figure FDA0002831831890000022
in the formula: ts is a sampling interval;
the equation (2) is simplified to a matrix form:
Xn1=A1S+N1 (3)
in the formula: a. the1Representing the phase angle difference of each frequency component of the reference element when the sampling elements meet the sampling characteristic matrix; s is a column vector consisting of each frequency component; n is a radical of1Is a noise vector;
for X, the same principle appliesnThe matrix X composed of the last L-1 elementsn2Comprises the following steps:
Xn2=A2S+N2 (4)
in the formula: a. the2Is a sampling characteristic matrix; n is a radical of2Is a noise vector;
thirdly, the corresponding amplitude x can be obtained by the least square method through the frequency of each component of the signalmiAngle of sum
Figure FDA0002831831890000024
Matrix parameters are defined as shown in equations (5) and (6):
Figure FDA0002831831890000023
Xs=[x(0),x(1),…,x(L+M-2)]T (6)
fourthly, according to a least square method, a matrix S containing phasor information of each frequency component of the signal is as follows:
S=(WTW)-1WTXs=[s1,s2,…,sk]T (7)
the amplitude and phase angle corresponding to each component of the signal are:
Figure FDA0002831831890000031
in order to improve the measurement performance of the ESPRIT algorithm under the condition of containing noise, a signal frequency component estimation method based on kurtosis is provided;
performing SVD processing on the L multiplied by M sampling data matrix X:
X=UPV (9)
in the formula: p is singular value matrix, P ═ diag (P)1,p2,…,pL)。
Sixthly, regarding the singular value of X from small to large as a discrete signal, the kurtosis K of the mth singular valuemThe expression is as follows:
Figure FDA0002831831890000032
in the formula: p ═ pm-n+1,pm-n,…,pmThe m-n +1 th to m-th singular values of the sampled data matrix X are arranged from small to large;
Figure FDA0002831831890000033
is the mean value of all elements in p; sigma being pStandard deviation, the expression is as follows:
Figure FDA0002831831890000034
the following equations (10) and (11) can be obtained:
Figure FDA0002831831890000035
when the elements in the group p are noise singular values, the kurtosis is due to the fact that no significant impact part of the noise singular values appears overall
Figure FDA0002831831890000041
Will not generate too large fluctuation and is always at the limit value KaThe following; when m is1+1 singular values into the array P, kurtosis
Figure FDA0002831831890000042
Will increase significantly; thus by judging kurtosis value KmThe frequency multiplied by the fraction k of the signal can be estimated as:
Figure FDA0002831831890000043
CN202011461130.7A 2020-12-11 2020-12-11 Rapid broadband measuring device and method based on ESPRIT algorithm Pending CN112698092A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011461130.7A CN112698092A (en) 2020-12-11 2020-12-11 Rapid broadband measuring device and method based on ESPRIT algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011461130.7A CN112698092A (en) 2020-12-11 2020-12-11 Rapid broadband measuring device and method based on ESPRIT algorithm

Publications (1)

Publication Number Publication Date
CN112698092A true CN112698092A (en) 2021-04-23

Family

ID=75507321

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011461130.7A Pending CN112698092A (en) 2020-12-11 2020-12-11 Rapid broadband measuring device and method based on ESPRIT algorithm

Country Status (1)

Country Link
CN (1) CN112698092A (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101252382A (en) * 2008-03-17 2008-08-27 成都国恒空间技术工程有限公司 Wide frequency range signal polarizing and DOA estimating method and apparatus
US20120200284A1 (en) * 2011-02-09 2012-08-09 Smart Energy Instruments, Inc. Power measurement device
CN103716262A (en) * 2012-10-09 2014-04-09 王晓安 Channel estimation by time-domain parameter extraction
CN105301557A (en) * 2015-11-06 2016-02-03 中国石油大学(华东) Direction-of-arrival estimate configuration method
CN108469545A (en) * 2017-12-13 2018-08-31 南京南瑞继保电气有限公司 A kind of substation's wideband measuring device
CN110231514A (en) * 2019-06-26 2019-09-13 北京四方继保自动化股份有限公司 A kind of synchronous phasor measuring method suitable for wideband measurement
CN111190042A (en) * 2019-12-30 2020-05-22 中国电力科学研究院有限公司 Intelligent sensing terminal and method for measuring full frequency band of power grid signal
CN111521898A (en) * 2020-05-06 2020-08-11 中国电力科学研究院有限公司 Power grid broadband signal intelligent sensing terminal and implementation method

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101252382A (en) * 2008-03-17 2008-08-27 成都国恒空间技术工程有限公司 Wide frequency range signal polarizing and DOA estimating method and apparatus
US20120200284A1 (en) * 2011-02-09 2012-08-09 Smart Energy Instruments, Inc. Power measurement device
CN103716262A (en) * 2012-10-09 2014-04-09 王晓安 Channel estimation by time-domain parameter extraction
CN105301557A (en) * 2015-11-06 2016-02-03 中国石油大学(华东) Direction-of-arrival estimate configuration method
CN108469545A (en) * 2017-12-13 2018-08-31 南京南瑞继保电气有限公司 A kind of substation's wideband measuring device
CN110231514A (en) * 2019-06-26 2019-09-13 北京四方继保自动化股份有限公司 A kind of synchronous phasor measuring method suitable for wideband measurement
CN111190042A (en) * 2019-12-30 2020-05-22 中国电力科学研究院有限公司 Intelligent sensing terminal and method for measuring full frequency band of power grid signal
CN111521898A (en) * 2020-05-06 2020-08-11 中国电力科学研究院有限公司 Power grid broadband signal intelligent sensing terminal and implementation method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘灏等: "基于K-ESPRIT的快速宽频测量方法", 《电力系统自动化》 *
吴艳平等: "多功能宽频测量装置的设计与实现", 《电力系统自动化》 *

Similar Documents

Publication Publication Date Title
CN102279327B (en) On-line monitoring and state evaluation system for photovoltaic grid-connected power generation
CN102565574B (en) Ship electric energy quality monitoring system
CN102495281B (en) Method for measuring phasor frequency of power system
CN108490257B (en) Short-time window inter-harmonic measurement method based on spectral curve fitting
CN105004939A (en) Composite electric energy quality disturbance signal quantitative analysis method
CN107064630B (en) A kind of power system frequency measurement method and device
CN104793053A (en) DFT (discrete Fourier transform) based synchronous phaser phase angle measurement method
CN104217112A (en) Multi-type signal-based power system low-frequency oscillation analysis method
CN106324340A (en) Method of measuring dynamic performance based on phasor and frequency
CN105353330A (en) On-line verifying system and algorithm for electronic current transformer based on virtual instrument technology
Stanisavljević et al. Voltage dips detection in a system with grid-tie inverter
CN112018784B (en) Subsynchronous resonance tracing method based on synchronous phasor measurement data
CN112698092A (en) Rapid broadband measuring device and method based on ESPRIT algorithm
CN110320400B (en) Voltage flicker envelope parameter extraction method for quasi-synchronous sampling and improved energy operator
CN103560509B (en) Voltage sag detection device based on wavelet analysis and control method of the device
CN103630744B (en) PMU phase-angle detection method and system
Rane et al. Real-Time Electromechanical Mode Identification through Energy-sorted Matrix Pencil Method
Jiang et al. A power cable fault location method combining with wavelet analysis and curve fitting
CN102680783B (en) Dislocation analysis weighting method capable of improving interharmonic calculation accuracy
Liu et al. An approach to power system harmonic analysis based on triple-line interpolation discrete Fourier transform
Xia et al. Research on field testing method of digital energy meters based on digital reference meter with high accuracy
Chen et al. Single-phase short-circuit voltage sag source location method based on wavelet transform and traveling wave ranging
Zhong et al. Measuring Method of Harmonics and Interharmonics Based on Ratio Correction and Matrix Inversion
Wang et al. ESPRIT-based fast broadband measurement method
Tomse et al. Detection and Monitoring System of Disturbances in Electrical Networks

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
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20211115

Address after: 125099 No. 1, Longwan South Street, Longgang District, Huludao City, Liaoning Province

Applicant after: HULUDAO POWER SUPPLY COMPANY OF STATE GRID LIAONING ELECTRIC POWER Co.,Ltd.

Applicant after: STATE GRID LIAONING ELECTRIC POWER SUPPLY Co.,Ltd.

Applicant after: STATE GRID CORPORATION OF CHINA

Address before: 125000 Longwan South Street Nanyan Road, Longgang District, Huludao City, Liaoning Province

Applicant before: HULUDAO POWER SUPPLY COMPANY OF STATE GRID LIAONING ELECTRIC POWER Co.,Ltd.

RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20210423