skip to main content


Title: Peak Estimation for Uncertain and Switched Systems
Peak estimation bounds extreme values of a function of state along trajectories of a dynamical system. This paper focuses on extending peak estimation to continuous and discrete settings with time-independent and time-dependent uncertainty. Techniques from optimal control are used to incorporate uncertainty into an existing occupation measure-based peak estimation framework, which includes special consideration for handling switching-type (polytopic) uncertainties. The resulting infinite-dimensional linear programs can be solved approximately with Linear Matrix Inequalities arising from the moment-SOS hierarchy.  more » « less
Award ID(s):
1808381 1646121 2038493
NSF-PAR ID:
10349423
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
60th IEEE Conf. Decision and Control
Page Range / eLocation ID:
3222 to 3228
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. This work bounds extreme values of state functions for a class of input-affine continuous-time systems that are affected by polyhedral-bounded uncertainty. Instances of these systems may arise in data-driven peak estimation, in which the state function must be bounded for all systems that are consistent with a set of state-derivative data records corrupted under L-infinity bounded noise. Existing occupation measure-based methods form a convergent sequence of outer approximations to the true peak value, given an initial set, by solving a hierarchy of semidefinite programs in increasing size. These techniques scale combinatorially in the number of state variables and uncertain parameters. We present tractable algorithms for peak estimation that scale linearly in the number of faces of the uncertainty-bounding polytope rather than combinatorially in the number of uncertain parameters by leveraging convex duality and a theorem of alternatives (facial decomposition). The sequence of decomposed semidefinite programs will converge to the true peak value under mild assumptions (convergence and smoothness of dynamics). 
    more » « less
  2. In this paper, we present a method for updating robotic state belief through contact with uncertain surfaces and apply this update to a Kalman filter for more accurate state estimation. Examining how guard surface uncertainty affects the time spent in each mode, we derive a novel guard saltation matrix- which maps perturbations prior to hybrid events to perturbations after - accounting for additional variation in the resulting state. Additionally, we propose the use of parameterized reset functions - capturing how unknown parameters change how states are mapped from one mode to the next - the Jacobian of which accounts for additional uncertainty in the resulting state. The accuracy of these mappings is shown by simulating sampled distributions through uncertain transition events and comparing the resulting covariances. Finally, we integrate these additional terms into the “uncertainty aware Salted Kalman Filter”, uaSKF, and show a peak reduction in average estimation error by 24–60% on a variety of test conditions and systems. 
    more » « less
  3. Abstract

    The near‐surface air temperature lapse rate is the predominant source of spatial temperature variability in mountains and controls snowfall and snowmelt regimes, glacier mass balance, and species distributions. Lapse rates are often estimated from observational data, however there is little guidance on best practices for estimating lapse rates. We use observational and synthetic datasets to evaluate the error and uncertainty in lapse rate estimates stemming from sample size, dataset noise, covariate collinearity, domain selection, and estimation methods. We find that lapse rates estimated from small sample sizes (<5) or datasets with high noise or collinearity can have errors of several °C km−1. Uncertainty in lapse rates due to non‐elevation related large‐scale temperature variability was reduced by correcting for spatial temperature gradients and restricting domains based on spatial clusters of stations. We generally found simple linear regression to be more robust than multiple linear regression for lapse rate estimation. Finally, lapse rates had lower error and uncertainty when estimated from a sample of topoclimatically self‐similar stations. Motivated by these results, we outline a set of best practices for lapse rate estimation that include using quality controlled temperature observations from as many locations as possible within the study domain, accounting for and minimizing non‐elevational sources of climatic gradients, and calculating lapse rates using simple linear regression across topoclimatically self‐similar samples of stations which are roughly 80% of the station population size.

     
    more » « less
  4. Abstract Rapid earthquake magnitude estimation from real-time space-based geodetic observation streams provides an opportunity to mitigate the impact of large and potentially damaging earthquakes by issuing low-latency warnings prior to any significant and destructive shaking. Geodetic contributions to earthquake characterization and rapid magnitude estimation have evolved in the last 20 yr, from post-processed seismic waveforms to, more recently, improved capacity of regional geodetic networks enabled real-time Global Navigation Satellite System seismology using precise point positioning (PPP) displacement estimates. In addition, empirical scaling laws relating earthquake magnitude to peak ground displacement (PGD) at a given hypocentral distance have proven effective in rapid earthquake magnitude estimation, with an emphasis on performance in earthquakes larger than ∼Mw 6.5 in which near-field seismometers generally saturate. Although the primary geodetic contributions to date in earthquake early warning have focused on the use of 3D position estimates and displacements, concurrent efforts in time-differenced carrier phase (TDCP)-derived velocity estimates also have demonstrated that this methodology has utility, including similarly derived empirical scaling relationships. This study builds upon previous efforts in quantifying the ambient noise of three-component ground-displacement and ground-velocity estimates. We relate these noise thresholds to expected signals based on published scaling laws. Finally, we compare the performance of PPP-derived PGD to TDCP-derived peak ground velocity (PGV), given several rich event datasets. Our results indicate that TDCP-PGV is more likely than PPP-PGD to detect intermediate magnitude (∼Mw 5.0–6.0) earthquakes, albeit with greater magnitude estimate uncertainty and across smaller epicentral distances. We conclude that the computationally lightweight TDCP-derived PGV magnitude estimation is complementary to PPP-derived PGD magnitude estimates, which could be produced at the network edge at high rates and with increased sensitivity to ground motion than current PPP estimates. 
    more » « less
  5. Abstract

    In light of the significant damage observed after earthquakes in Japan and New Zealand, enhanced performing seismic force‐resisting systems and energy dissipation devices are increasingly being utilized in buildings. Numerical models are needed to estimate the seismic response of these systems for seismic design or assessment. While there have been studies on modeling uncertainty, selecting the model features most important to response can remain ambiguous, especially if the structure employs less well‐established lateral force‐resisting systems and components. Herein, a global sensitivity analysis was used to address modeling uncertainty in specimens with elastic spines and force‐limiting connections (FLCs) physically tested at full‐scale at the E‐Defense shake table in Japan. Modeling uncertainty was addressed for both model class and model parameter uncertainty by varying primary models to develop several secondary models according to pre‐established uncertainty groups. Numerical estimates of peak story drift ratio and floor acceleration were compared to the results from the experimental testing program using confidence intervals and root‐mean‐square error. Metrics such as the coefficient of variation, variance, linear Pearson correlation coefficient, and Sobol index were used to gain intuition about each model feature's contribution to the dispersion in estimates of the engineering demands. Peak floor acceleration was found to be more sensitive to modeling uncertainty compared to story drift ratio. Assumptions for the spine‐to‐frame connection significantly impacted estimates of peak floor accelerations, which could influence future design methods for spines and FLC in enhanced lateral‐force resisting systems.

     
    more » « less