An official website of the United States government
To perform power system monitoring and control using synchrophasor measurements, various dynamic state estimators have been proposed in the literature, including the extended Kalman filter (EKF) and the unscented Kalman filter (UKF). However, they are unable to handle system model parameter errors and any type of outliers, precluding them from being adopted for power system real-time applications. In this paper, we develop a robust iterated extended Kalman filter based on the generalized maximum likelihood approach (termed GM-IEKF) for dynamic state estimation. The proposed GM-IEKF can effectively suppress observation and innovation outliers, which may be induced by model parameter gross errors and cyber attacks. We assess its robustness by carrying out extensive simulations on the IEEE 39-bus test system. From the results, we find that the GM-IEKF is able to cope with at least 25% outliers, including in position of leverage.
more »
« less
This content will become publicly available on January 1, 2027
Continuous Venous Oxygen Saturation Estimation: A Robust Population-Informed-Personalized Gaussian Sum-Extended Kalman Filtering Approach
Abstract This article investigates the robustification of the population-informed-personalized Gaussian sum-extended Kalman filter (PI-P-GSEKF) developed in our prior work and its application to continuous venous oxygen saturation (SvO2) estimation. The PI-P-GSEKF was developed to enable state estimation in systems with extremely large variability. It includes a bank of extended Kalman filters (EKFs), whose operating points (i.e., nominal parameter vectors) are selected via generative sampling followed by Markov Chain Monte Carlo (MCMC) sampling with one-time partial state measurement. The state is estimated as the weighted sum of state estimates from all of the EKFs, with the weight for each EKF calculated based on the likelihood of its prediction at every measurement instant. Despite its adequate performance in general, its state estimate can suffer from high-sensitivity operating points whose inaccuracy with respect to the ground truth operating point results in large EKF errors and adversely impacts the PI-P-GSEKF. We explored two ideas to robustify the PI-P-GSEKF against this challenge: (i) penalizing high-sensitivity operating points in MCMC sampling (called robust MCMC sampling) and (ii) calculating the weights for the EKFs based on the likelihood of their predictions in a measurement horizon (called robust Gaussian summing). We examined the efficacy of these ideas in the context of continuous venous oxygen saturation (SvO2) estimation from arterial oxygen saturation measurement, which is important in critical care and cardiopulmonary medicine but is highly invasive and challenging. The results suggested that both ideas could reduce SvO2 estimation error compared with the standard PI-P-GSEKF ((i): 4%; (ii): 13%; (i) + (ii): 16%, on average). However, how to set the length of the sampling interval for weight calculation remains an open challenge.
more »
« less
- Award ID(s):
- 2322533
- PAR ID:
- 10671878
- Publisher / Repository:
- ASME
- Date Published:
- Journal Name:
- ASME Letters in Dynamic Systems and Control
- Volume:
- 6
- Issue:
- 1
- ISSN:
- 2689-6117
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We prove that among all 1-periodic configurations $$\Gamma $$ of points on the real line $$\mathbb{R}$$ the quantities $$\min _{x \in \mathbb{R}} \sum _{\gamma \in \Gamma } e^{- \pi \alpha (x - \gamma )^{2}}$$ and $$\max _{x \in \mathbb{R}} \sum _{\gamma \in \Gamma } e^{- \pi \alpha (x - \gamma )^{2}}$$ are maximized and minimized, respectively, if and only if the points are equispaced and whenever the number of points $$n$$ per period is sufficiently large (depending on $$\alpha $$). This solves the polarization problem for periodic configurations with a Gaussian weight on $$\mathbb{R}$$ for large $$n$$. The first result is shown using Fourier series. The second result follows from the work of Cohn and Kumar on universal optimality and holds for all $$n$$ (independent of $$\alpha $$).more » « less
-
ABSTRACT Oxygen store management underlies dive performance and is dependent on the slow heart rate and peripheral vasoconstriction of the dive response to control tissue blood flow and oxygen uptake. Prior research has revealed two major patterns of muscle myoglobin saturation profiles during dives of emperor penguins. In Type A profiles, myoglobin desaturated rapidly, consistent with minimal muscle blood flow and low tissue oxygen uptake. Type B profiles, with fluctuating and slower declines in myoglobin saturation, were consistent with variable tissue blood flow patterns and tissue oxygen uptake during dives. We examined arterial and venous blood oxygen profiles to evaluate blood oxygen extraction and found two primary patterns of venous hemoglobin desaturation that complemented corresponding myoglobin saturation profiles. Type A venous profiles had a hemoglobin saturation that (a) increased/plateaued for most of a dive's duration, (b) only declined during the latter stages of ascent, and (c) often became arterialized [arterio-venous (a-v) shunting]. In Type B venous profiles, variable but progressive hemoglobin desaturation profiles were interrupted by inflections in the profile that were consistent with fluctuating tissue blood flow and oxygen uptake. End-of-dive saturation of arterial and Type A venous hemoglobin saturation profiles were not significantly different, but did differ from those of Type B venous profiles. These findings provide further support that the dive response of emperor penguins is a spectrum of cardiac and vascular components (including a-v shunting) that are dependent on the nature and demands of a given dive and even of a given segment of a dive.more » « less
-
Zhang, Qichun (Ed.)We develop a general framework for state estimation in systems modeled with noise-polluted continuous time dynamics and discrete time noisy measurements. Our approach is based on maximum likelihood estimation and employs the calculus of variations to derive optimality conditions for continuous time functions. We make no prior assumptions on the form of the mapping from measurements to state-estimate or on the distributions of the noise terms, making the framework more general than Kalman filtering/smoothing where this mapping is assumed to be linear and the noises Gaussian. The optimal solution that arises is interpreted as a continuous time spline, the structure and temporal dependency of which is determined by the system dynamics and the distributions of the process and measurement noise. Similar to Kalman smoothing, the optimal spline yields increased data accuracy at instants when measurements are taken, in addition to providing continuous time estimates outside the measurement instances. We demonstrate the utility and generality of our approach via illustrative examples that render both linear and nonlinear data filters depending on the particular system. Application of the proposed approach to a Monte Carlo simulation exhibits significant performance improvement in comparison to a common existing method.more » « less
-
A newly enhanced, recursive, and robust Bayesian state estimation algorithm for linear and nonlinear systems, referred to as the Multivariate Cauchy Estimator (MCE), is presented. The algorithm enables robust state estimation performance for applications with more volatile system noises than the Gaussian distribution suggests. This is achieved by over-bounding realistic process and measurement noises with additive, heavy-tailed Cauchy random variables. The characteristic function (CF) of the un-normalized conditional probability density function (ucpdf) is propagated as a growing sum of terms in the MCE. Here, the CF is simplified by replacing the original with a representation of linear parameter vectors that operate on bases composed of indicator functions. This insight can lead to eliminating over 99% of terms that previously comprised this CF. Using graphical processing units, the MCE can exploit its parallel mathematical structure and achieve a fast execution rate. A target tracking example shows the robustness of the MCE over the Kalman filter in both heavy-tailed and Gaussian noise.more » « less
