CN109950903A - A kind of dynamic state estimator method counted and noise statistics are unknown - Google Patents
A kind of dynamic state estimator method counted and noise statistics are unknown Download PDFInfo
- Publication number
- CN109950903A CN109950903A CN201910307752.5A CN201910307752A CN109950903A CN 109950903 A CN109950903 A CN 109950903A CN 201910307752 A CN201910307752 A CN 201910307752A CN 109950903 A CN109950903 A CN 109950903A
- Authority
- CN
- China
- Prior art keywords
- state
- moment
- noise
- generator
- measurement
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 239000011159 matrix material Substances 0.000 claims abstract description 66
- 238000001914 filtration Methods 0.000 claims abstract description 23
- 230000003044 adaptive effect Effects 0.000 claims abstract description 10
- 238000005259 measurement Methods 0.000 claims description 45
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000012937 correction Methods 0.000 claims description 6
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 claims description 3
- 230000000644 propagated effect Effects 0.000 claims description 3
- 229910017435 S2 In Inorganic materials 0.000 claims 1
- 238000005530 etching Methods 0.000 claims 1
- 238000004422 calculation algorithm Methods 0.000 abstract description 9
- 230000015556 catabolic process Effects 0.000 abstract description 4
- 238000006731 degradation reaction Methods 0.000 abstract description 4
- 230000007613 environmental effect Effects 0.000 abstract description 2
- 230000001052 transient effect Effects 0.000 description 6
- 230000003068 static effect Effects 0.000 description 4
- 230000001360 synchronised effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000013016 damping Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000017105 transposition Effects 0.000 description 2
- 230000000052 comparative effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a kind of dynamic state estimator methods that meter and noise statistics are unknown, for realizing the accurate estimation of electric system generator dynamic state variables.This method is in traditional Unscented kalman filtering algorithm (unscented Kalman filter, UKF under frame), the system noise based on innovative information sequence is established respectively and measures noise statistics estimator, dynamically can adjust with update the system noise and measure the covariance matrix of noise satisfaction according to environmental change.Mentioned method (adaptive unscented Kalman filter, AUKF) can effectively solve traditional UKF method because noise matrix setting it is improper caused by state estimation degradation problem, promoted dynamic state estimator device precision of state estimation.The algorithm is because of meter and Practical Project background, and simple and convenient easy to implement, engineering application value with higher.
Description
Technical field
The present invention relates to a kind of dynamic state estimator methods that meter and noise statistics are unknown, belong to Power System Analysis
And monitoring technical field.
Background technique
Accurate state estimation is of great significance for the analysis of electric system with stability contorting.State estimation is generally divided
For two classes, one kind is static estimation, and another kind of is dynamic state estimator.Static state estimation utilizes a certain moment section redundancy
Measurement information realizes the estimation of the system moment state variable.Although static state estimated accuracy is higher, it has ignored electric power
The dynamic characteristic of system.Therefore, the real-time online that static state estimation can not be applied to POWER SYSTEM STATE is estimated.In order to adapt to
The demand of electric system on-line monitoring, the dynamic state estimator method for having estimation and forecast function obtain researcher in recent years
Extensive concern.
Currently, dynamic state estimator device has gradually been applied to the state estimation of electric system generator.The state of use is estimated
Meter method specifically includes that Kalman filtering, Extended Kalman filter, volume Kalman filtering, Unscented kalman filtering etc..It is above-mentioned
These methods improve dynamic state estimator precision to a certain extent.However, it is worth noting that these methods are assumed mostly
System noise and the covariance matrix for measuring noise satisfaction are constant;And in practical power systems, system and measurement noise
Statistical property is difficult accurately to obtain in advance, and noise covariance matrix setting is closely related with the performance of dynamic state estimator device.
So, by degradation mode estimator performance, reducing state estimation essence if system and measurement noise covariance matrix setting are improper
Degree even results in state estimator failure.
Summary of the invention
Goal of the invention: present invention seek to address that improper draw is arranged for system and measurement noise in existing dynamic state estimator device
The degradation problem risen promotes generator dynamic state estimator precision, provides solid number for the safe and stable operation of power grid
It is believed that breath.
The invention discloses a kind of dynamic state estimator methods that meter and noise statistics are unknown, send out for electric system
Motor dynamics state estimation, includes the following steps:
S1: according to generator quadravalence dynamical equation, Generator Status estimation equation is established, and constructs and obtains generator dynamic
The state equation and measurement equation of state estimation:
In formula, f () indicates that Generator Status equation, h () indicate that measurement equation, x, u, z respectively correspond expression state
Variable controls variable and measures vector;Subscript k-1 and k indicate the moment, and w indicates system noise, and v is to measure noise, both for
White Gaussian noise, the mean value met are 0, and corresponding covariance matrix is respectively Q and R, w and v is mutually indepedent and and state
Variable is unrelated;
S2: based on state estimation known to the k-1 momentUsing Unscented transform technology generate the k moment several
Sigma state samples;
S3: it is walked using the status predication of adaptive Unscented kalman filtering, calculates the Generator Status predicted value at k moment
With status predication error co-variance matrix
S4: being based on generator measurement equation, calculates the generator measurement predictor at k momentError association side is predicted with measuring
Poor matrix
S5: the cross-covariance matrix P between status predication value and measurement predictor is calculatedxz,k;
S6: based on adaptive Unscented kalman filtering filtering step, the filtering gain at k moment is obtained;
S7: the measurement information value z at k moment is utilizedk, to status predication valueIt is modified update, state estimation is calculated
Error co-variance matrix
S8: it is based on innovative information valueDynamic corrections update k moment system noise covariance matrix QkWith measurement noise
Covariance matrix Rk;
S9: according to S2-S8 according to measurement information value zkTo electric system generator state dynamic estimation, until k+1 > N, N
State estimation iteration stopping when for the maximum estimated moment, output state estimated result.
Further, 2n+1 Sigma state samples of k moment are generated in S2
In formula, n indicates Generator Status predictor dimension,For state estimation error association side known to the k-1 moment
Poor matrix,It indicates to carry out matrix Cholesky decomposition operation, κ indicates scale parameter.
Further, the Generator Status predicted value at the k moment in S3With status predication error co-variance matrixTable
Show as follows:
In formula,Indicate Sigma state samplesCorrespondence point value after the propagation of Generator Status equation, Qk-1It indicates
The covariance matrix that k-1 moment system noise meets, WlFor the weighted value corresponding to Sigma state samples:
Further, in S4, k moment generator measurement predictorWith measurement predicting covariance matrixIt indicates
Are as follows:
In formula,Indicate correspondence point value of the status predication point after measurement equation is propagated, Rk-1It indicates that the k-1 moment measures to make an uproar
The covariance matrix that sound meets.
Further, in S5, cross-covariance matrix Pxz,kIt indicates are as follows:
Further, in S6, the filtering gain K at k momentkIt indicates are as follows:
Further, in S7, the measurement information value z at k momentk, to status predication valueIt is modified update, and calculates shape
State evaluated error covariance matrixIts calculation formula is:
Further, in S8, dynamic corrections update k moment system noise covariance matrix QkWith measurement noise covariance square
Battle array Rk, calculation formula is as follows:
dk-1=(1-b)/(1-bk) (16)
In formula, b indicates forgetting factor, dk-1For k-1 moment system noise covariance matrix estimator adjustment parameter, diag
() is building diagonal matrix symbol, and exp is indicated using natural constant as the exponential function at bottom.
The utility model has the advantages that a kind of meter of the invention and the unknown dynamic state estimator method of noise statistics, in traditional nothing
Under the frame of mark Kalman filtering algorithm, the system noise based on innovative information sequence is established respectively and measures noise statistics
Estimator, the covariance matrix that can meet according to the adjustment of environmental change dynamic and update the system noise and measurement noise, overcomes
Traditional UKF method because noise matrix setting it is improper caused by state estimation degradation problem, promote dynamic state estimator device
Precision of state estimation.The algorithm has fully considered the actual application background of generator dynamic state estimator, and simple and convenient is easy to
Implement, engineering application value with higher.
Detailed description of the invention
Fig. 1 is flow chart of the method for the present invention;
Fig. 2 is 39 node system structure chart of IEEE10 machine;
Fig. 3 is to utilize tradition UKF method and the method for the present invention to the dynamic estimation result pair of generator's power and angle and angular speed
Than;
Fig. 4 is to utilize tradition UKF method and the method for the present invention to the dynamic estimation Comparative result of generator transient internal voltage;
Fig. 5 is to be compared using tradition UKF method and the method for the present invention to Generator Status estimated result error.
Specific embodiment
The present invention is further explained with reference to the accompanying drawings and examples.
As shown in Figure 1, estimating with the method for the present invention embodiment test macro dynamic variable, including walk as follows
It is rapid:
Step 1: Generator Status estimates model foundation:
The state equation and measurement equation of generator dynamic state estimator, general type can indicate are as follows:
In formula, f () indicates that Generator Status equation, h () indicate that measurement equation, x, u, z respectively correspond expression state
Variable controls variable and measures vector;Subscript k-1 and k indicate the moment, and w indicates system noise, and v is to measure noise, both for
White Gaussian noise, the mean value met are 0, and corresponding covariance matrix is respectively Q and R, w and v is mutually indepedent and and state
Variable is unrelated.
Step 2: the initial parameter value of adaptive Unscented kalman filtering method designed by the setting present invention includes setting
Initial time state variable valueState estimation error co-variance matrixControl variable value u0It is taken as steady-state operation value, scale
Parameter κ;The initial covariance matrix Q that setting system noise and measurement noise are met0, R0And maximum estimated moment N;
Step 3: based on state estimation known to the k-1 momentK moment 2n+1 are generated using Unscented transform technology
Sigma state samplesWherein n indicates Generator Status predictor dimension, and calculation formula is as follows:
In formula,For state estimation error co-variance matrix known to the k-1 moment,It indicates to carry out matrix
Cholesky decomposition operation, κ ∈ [10-4, 1] and indicate scale parameter.
Step 4: it is walked using the status predication of adaptive Unscented kalman filtering, calculates the status predication value at k momentWith
Status predication error co-variance matrixSolution formula is as follows
In formula,Indicate Sigma state samplesCorrespondence point value after the propagation of Generator Status equation, subscript T table
Show the transposition operation of matrix, Qk-1Indicate the covariance matrix that k-1 moment system noise meets, WlTo be adopted corresponding to Sigma state
The weighted value of sampling point, calculation formula are as follows
Step 5: function is measured based on generator, calculates k moment generator measurement predictorIt is assisted with prediction error is measured
Variance matrixCalculation formula is
In formula,Indicate correspondence point value of the status predication point after measurement equation is propagated, Rk-1Indicate that the k-1 moment measures mistake
Poor covariance matrix.
Step 6: calculating status predication and measures the cross-covariance matrix P between predictionxz,k, calculation formula is
Step 7: based on designed adaptive Unscented kalman filtering filtering step, the filtering gain K at k moment is calculatedk, meter
Calculating formula is
Subscript () in formula-1The inversion operation of representing matrix.
Step 8: the measurement information value z at k moment is utilizedk, to status predication valueIt is modified update, and calculates state and estimates
Count error co-variance matrixIts calculation formula is
The transposition operation of subscript T representing matrix in formula.
Step 9: it is based on innovative information valueDynamic corrections update k moment system noise covariance matrix QkAnd measurement
Noise covariance matrix Rk, calculation formula is as follows
dk-1=(1-b)/(1-bk) (16)
B indicates forgetting factor, and 0 < b < 1 in formula, and when the variation degree of system mode is bigger, the value of b is bigger;
dk-1For k-1 moment system noise covariance matrix estimator adjustment parameter, diag () is building diagonal matrix symbol, exp table
Show the exponential function using natural constant the bottom of as.
Step 10: it is dynamic to electric system generator state according to measurement information time series that step is calculated according to (2)-(9)
State estimation, until state estimation iteration stopping when k+1 > N, output state estimated result.
Embodiment:
(a) model foundation
According to generator quadravalence dynamical equation, the Generator Status estimation equation of building is as follows:
In formula: δ indicates generator's power and angle, rad;ω and ω0Respectively angular rate and synchronous rotational speed, pu;e′qWith e 'dPoint
Not Biao Shi generator q axis and d axis transient internal voltage;H indicates generator inertia constant, TmAnd TeRespectively indicate generator mechanical
Power and electromagnetic power, wherein Te=Pe/ω;KDIndicate damping factor, EfdFor stator excitation voltage;T′d0With T 'q0Indicate power generation
Open circuit time constant of the machine machine under d-q coordinate system;xdWith x 'dRespectively indicate generator d axis synchronous reactance and transient state reactance, xq
With x 'qRespectively generator q axis synchronous reactance and transient state reactance;idAnd iqRespectively indicate the stator current of generator d axis and q axis.
When carrying out dynamic estimation to electric system generator dynamic variable, state estimation vector is x=(δ, ω, e 'q,e′d)T;
Choose the electric current i of generator mechanical power, stator excitation voltage and stator R axis and I axisR,iIFor dominant vector, i.e. u=(Tm,
Efd,iR,iI)T;Choose the voltage e of the absolute generator rotor angle of generator, generator angular speed and generator unit stator R axis and I axisR,eIAs
Measuring value, i.e. measurement vector are
Z=(δ, ω, eR,eI)T
Wherein the absolute generator rotor angle of generator and angular speed can directly measure acquisition by PMU measurement equipment, be under this situation
System meets controllability.
(b) embodiment is analyzed
In order to verify meter and noise statistics proposed by the invention unknown dynamic state estimator method validity and reality
With property, the present embodiment chooses 39 node system of IEEE10 machine as test macro, and system structure is shown in Fig. 2.It is verified in algorithm
When, the state variable of the generator G9 using in system is as estimation object, and wherein generator uses quadravalence model.Generator G9's is used
Property time parameter be 34.5, damping factor 2, d axis and q axis open circuit transient time-constant be respectively 4.79 and 1.96.For simulation
Generator transient characteristic, it is assumed that in 20 cycle, three-phase metallic short circuit failure, failure occur for 21 branch of node 16- node
Continued for 6 periods (sampling period 0.02s) disappears afterwards.
With the simulation PMU data acquisition of BPA software, obtains generator and run true value.Metric data value is folded by true value
Random noise is added to be formed.500 cycles (1 cycle is 0.02s) measuring value carries out algorithm before the present invention takes when carrying out emulation experiment
Verifying, i.e. N are 500.The initial value of state variable chooses the quiescent value of last moment when estimation.It is assumed that system noise and measurement noise
Statistical property is unknown, i.e., its value is set as follows: system noise initially defences arranged in matrix jointly as Q0=10-5I4×4, measure noise
Initial covariance matrix is set as R0=10-4I4×4, and the true value of the two is respectively Q0=10-6I4×4And R0=10-6I4×4。
In order to compare and analyze to the estimated result between algorithms of different, the present invention uses average absolute evaluated error
MAE is as performance comparison between index progress algorithm.
MAE (k) is the average value of the sum of absolute error of each state variable estimated result of k moment generator, x in formulai,kFor
The true value (BPA data) of i-th of quantity of state of k moment,Estimated value, N are corresponded to for itsFor total state variable number.
Be utilized respectively traditional Unscented kalman filtering algorithm (UKF) (related parameter values and the method for the present invention needed for it
Initial parameter values are identical) and adaptive Unscented kalman filtering method (AUKF) proposed by the present invention tested.
The state estimation result of generator G9 is respectively such as Fig. 3, shown in Fig. 4.Fig. 5 further provides UKF and AUKF method shape
The mean absolute error of state estimated result compares.It can be seen that from the state estimation result in figure in initial system noise association side
In the case of poor matrix and measurement noise covariance matrix setting improperly, AUKF method designed by the present invention still can be accurate
Generator state variables are tracked, there is higher precision of state estimation compared with UKF.
State estimation result confirms the adaptive Unscented kalman filtering method that is mentioned of the present invention due to being capable of base in figure
In innovative information sequence dynamic corrections and more new system and noise covariance matrix is measured, is overcome because of system noise covariance matrix
Improper caused estimator performance decline is set, generator dynamic state estimator precision is effectively promoted.I.e. the method for the present invention can
The demand of electric system generator dynamic state estimator is better met, realizes the accurate measurements of Electrical Power System Dynamic state.
Claims (8)
1. a kind of dynamic state estimator method counted and noise statistics are unknown, estimates for electric system generator dynamical state
Meter, which comprises the steps of:
S1: according to generator quadravalence dynamical equation, Generator Status estimation equation is established, and constructs and obtains generator dynamical state
The state equation and measurement equation of estimation:
In formula, f () indicates that Generator Status equation, h () indicate that measurement equation, x, u, z respectively correspond expression state variable,
It controls variable and measures vector;Subscript k-1 and k indicate the moment, and w indicates system noise, and v is to measure noise, both for Gauss
White noise, the mean value met are 0, and corresponding covariance matrix is respectively Q and R, w and v is mutually indepedent and and state variable
It is unrelated;
S2: based on state estimation known to the k-1 momentSeveral Sigma states of k moment are generated using Unscented transform technology
Sampled point;
S3: it is walked using the status predication of adaptive Unscented kalman filtering, calculates the Generator Status predicted value at k momentAnd shape
State predicting covariance matrix
S4: being based on generator measurement equation, calculates the generator measurement predictor at k momentWith measurement predicting covariance square
Battle array
S5: the cross-covariance matrix P between status predication value and measurement predictor is calculatedxz,k;
S6: based on adaptive Unscented kalman filtering filtering step, the filtering gain at k moment is obtained;
S7: the measurement information value z at k moment is utilizedk, to status predication valueIt is modified update and obtains state estimation, calculate
To state estimation error co-variance matrix
S8: it is based on innovative information valueDynamic corrections update k moment system noise covariance matrix QkWith measurement noise association side
Poor matrix Rk;
S9: according to S2-S8 according to measurement information value zkTo electric system generator state dynamic estimation, until k+1 > N, N are most
State estimation iteration stopping when the big estimation moment, output state estimated result.
2. a kind of dynamic state estimator method counted and noise statistics are unknown according to claim 1, feature exist
In: 2n+1 Sigma state samples of k moment are generated in S2
In formula, n indicates Generator Status predictor dimension,For state estimation error covariance square known to the k-1 moment
Battle array,It indicates to carry out matrix Cholesky decomposition operation, κ indicates scale parameter.
3. a kind of dynamic state estimator method counted and noise statistics are unknown according to claim 2, feature exist
In: the Generator Status predicted value at the k moment in S3With status predication error co-variance matrixIt is expressed as follows:
In formula,Indicate Sigma state samplesCorrespondence point value after the propagation of Generator Status equation, Qk-1When indicating k-1
The covariance matrix that etching system noise meets, WlFor the weighted value corresponding to Sigma state samples:
4. a kind of dynamic state estimator method counted and noise statistics are unknown according to claim 3, feature exist
In: in S4, k moment generator measurement predictorWith measurement predicting covariance matrixIt indicates are as follows:
In formula,Indicate correspondence point value of the status predication point after measurement equation is propagated, Rk-1It is full to indicate that the k-1 moment measures noise
The covariance matrix of foot.
5. a kind of dynamic state estimator method counted and noise statistics are unknown according to claim 4, feature exist
In: in S5, cross-covariance matrix Pxz,kIt indicates are as follows:
6. a kind of dynamic state estimator method counted and noise statistics are unknown according to claim 5, feature exist
In: in S6, the filtering gain K at k momentkIt indicates are as follows:
7. a kind of dynamic state estimator method counted and noise statistics are unknown according to claim 6, feature exist
In: in S7, the measurement information value z at k momentk, to status predication valueIt is modified update, obtains state estimation, and calculate
State estimation error co-variance matrixIts calculation formula is:
8. a kind of dynamic state estimator method counted and noise statistics are unknown according to claim 7, feature exist
In: in S8, dynamic corrections update k moment system noise covariance matrix QkWith measurement noise covariance matrix Rk, calculation formula is such as
Under:
dk-1=(1-b)/(1-bk) (16)
In formula, b indicates forgetting factor, dk-1For k-1 moment system noise covariance matrix estimator adjustment parameter, diag ()
To construct diagonal matrix symbol, exp is indicated using natural constant as the exponential function at bottom.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910307752.5A CN109950903A (en) | 2019-04-17 | 2019-04-17 | A kind of dynamic state estimator method counted and noise statistics are unknown |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910307752.5A CN109950903A (en) | 2019-04-17 | 2019-04-17 | A kind of dynamic state estimator method counted and noise statistics are unknown |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109950903A true CN109950903A (en) | 2019-06-28 |
Family
ID=67015350
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910307752.5A Pending CN109950903A (en) | 2019-04-17 | 2019-04-17 | A kind of dynamic state estimator method counted and noise statistics are unknown |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109950903A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112260590A (en) * | 2020-09-30 | 2021-01-22 | 西安理工大学 | Control method based on IAUKF in permanent magnet linear synchronous motor |
CN113178898A (en) * | 2021-06-08 | 2021-07-27 | 杭州电力设备制造有限公司 | Method, device, equipment and medium for estimating output state of distributed power supply |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107478990A (en) * | 2017-09-11 | 2017-12-15 | 河海大学 | A kind of generator electromechanical transient process method for dynamic estimation |
CN107590317A (en) * | 2017-08-17 | 2018-01-16 | 河海大学 | A kind of generator method for dynamic estimation of meter and model parameter uncertainty |
CN108155648A (en) * | 2018-01-09 | 2018-06-12 | 河海大学 | Method for estimating state based on the infinite Extended Kalman filter of adaptive H |
-
2019
- 2019-04-17 CN CN201910307752.5A patent/CN109950903A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107590317A (en) * | 2017-08-17 | 2018-01-16 | 河海大学 | A kind of generator method for dynamic estimation of meter and model parameter uncertainty |
CN107478990A (en) * | 2017-09-11 | 2017-12-15 | 河海大学 | A kind of generator electromechanical transient process method for dynamic estimation |
CN108155648A (en) * | 2018-01-09 | 2018-06-12 | 河海大学 | Method for estimating state based on the infinite Extended Kalman filter of adaptive H |
Non-Patent Citations (1)
Title |
---|
YI WANG等: "Robust Dynamic State Estimation of Power Systems with Model Uncertainties Based on Adaptive Unscented H∞ Filter", 《 IET GENERATION TRANSMISSION & DISTRIBUTION》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112260590A (en) * | 2020-09-30 | 2021-01-22 | 西安理工大学 | Control method based on IAUKF in permanent magnet linear synchronous motor |
CN113178898A (en) * | 2021-06-08 | 2021-07-27 | 杭州电力设备制造有限公司 | Method, device, equipment and medium for estimating output state of distributed power supply |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110032812A (en) | A kind of dynamic state estimator method based on adaptive volume Kalman filtering | |
CN107590317B (en) | Generator dynamic estimation method considering model parameter uncertainty | |
CN108155648B (en) | State estimation method based on adaptive H-infinity extended Kalman filtering | |
Ghahremani et al. | Dynamic state estimation in power system by applying the extended Kalman filter with unknown inputs to phasor measurements | |
CN107478990B (en) | A kind of generator electromechanical transient process method for dynamic estimation | |
CN110112770A (en) | A kind of generator dynamic state estimator method based on adaptive H ∞ volume Kalman filtering | |
CN106452247B (en) | Method and device for identifying rotational inertia of permanent magnet synchronous motor | |
CN108574291A (en) | One kind being based on Ensemble Kalman Filter generator dynamic state estimator method | |
CN110222309A (en) | A kind of generator method for dynamic estimation based on robust volume Kalman filtering | |
CN109218073B (en) | Dynamic state estimation method considering network attack and parameter uncertainty | |
CN109950903A (en) | A kind of dynamic state estimator method counted and noise statistics are unknown | |
CN113629729B (en) | Wind power system area inertia estimation method based on frequency measurement point selection | |
CN110768600A (en) | PMSM speed sensorless rotor detection method | |
CN110008638A (en) | A kind of dynamic state estimator method based on adaptive EnKF technology | |
CN109239596A (en) | A kind of dynamic state estimator method based on EKF-IRLS filtering | |
CN112511056A (en) | Robust generator dynamic state estimation method based on phasor measurement | |
CN109918862A (en) | A kind of generator method for dynamic estimation filtered based on robust without mark H infinity | |
CN113224991A (en) | Method, system, terminal and readable storage medium for identifying inductance of synchronous reluctance motor based on unscented Kalman filtering | |
Yin et al. | A speed estimation method for induction motors based on strong tracking extended Kalman filter | |
CN112446010A (en) | Adaptive weak-sensitive rank Kalman filtering method and application thereof | |
CN117013902A (en) | Motor inductance parameter calculation method, device and system, motor and power equipment | |
Zhao et al. | Online estimation of rotor temperature in induction motors based on parameter identification | |
CN109270455A (en) | Induction machine state monitoring method based on hyposensitiveness Ensemble Kalman Filter | |
CN106788064B (en) | Induction motor stator resistance parameter identification method based on EMD-ELM | |
CN110069870A (en) | A kind of generator dynamic state estimator method based on adaptively without mark H ∞ filtering |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190628 |
|
RJ01 | Rejection of invention patent application after publication |