No Arabic abstract
Most data based state and parameter estimation methods require suitable initial values or guesses to achieve convergence to the desired solution, which typically is a global minimum of some cost function. Unfortunately, however, other stable solutions (e.g., local minima) may exist and provide suboptimal or even wrong estimates. Here we demonstrate for a 9-dimensional Lorenz-96 model how to characterize the basin size of the global minimum when applying some particular optimization based estimation algorithm. We compare three different strategies for generating suitable initial guesses and we investigate the dependence of the solution on the given trajectory segment (underlying the measured time series). To address the question of how many state variables have to be measured for optimal performance, different types of multivariate time series are considered consisting of 1, 2, or 3 variables. Based on these time series the local observability of state variables and parameters of the Lorenz-96 model is investigated and confirmed using delay coordinates. This result is in good agreement with the observation that correct state and parameter estimation results are obtained if the optimization algorithm is initialized with initial guesses close to the true solution. In contrast, initialization with other exact solutions of the model equations (different from the true solution used to generate the time series) typically fails, i.e. the optimization procedure ends up in local minima different from the true solution. Initialization using random values in a box around the attractor exhibits success rates depending on the number of observables and the available time series (trajectory segment).
Observability of state variables and parameters of a dynamical system from an observed time series is analyzed and quantified by means of the Jacobian matrix of the delay coordinates map. For each state variable and each parameter to be estimated a measure of uncertainty is introduced depending on the current state and parameter values, which allows us to identify regions in state and parameter space where the specific unknown quantity can (not) be estimated from a given time series. The method is demonstrated using the Ikeda map and the Hindmarsh-Rose model.
Optimization-based state estimation is useful for nonlinear or constrained dynamic systems for which few general methods with established properties are available. The two fundamental forms are moving horizon estimation (MHE) which uses the nearest measurements within a moving time horizon, and its theoretical ideal, full information estimation (FIE) which uses all measurements up to the time of estimation. Despite extensive studies, the stability analyses of FIE and MHE for discrete-time nonlinear systems with bounded process and measurement disturbances, remain an open challenge. This work aims to provide a systematic solution for the challenge. First, we prove that FIE is robustly globally asymptotically stable (RGAS) if the cost function admits a property mimicking the incremental input/output-to-state stability (i-IOSS) of the system and has a sufficient sensitivity to the uncertainty in the initial state. Second, we establish an explicit link from the RGAS of FIE to that of MHE, and use it to show that MHE is RGAS under enhanced conditions if the moving horizon is long enough to suppress the propagation of uncertainties. The theoretical results imply flexible MHE designs with assured robust stability for a broad class of i-IOSS systems. Numerical experiments on linear and nonlinear systems are used to illustrate the designs and support the findings.
We present a method to determine the relative parameter mismatch in a collection of nearly identical chaotic oscillators by measuring large deviations from the synchronized state. We demonstrate our method with an ensemble of slightly different circle maps. We discuss how to apply our method when there is noise, and show an example where the noise intensity is comparable to the mismatch.
We present a general framework for sensitivity optimization in quantum parameter estimation schemes based on continuous (indirect) observation of a dynamical system. As an illustrative example, we analyze the canonical scenario of monitoring the position of a free mass or harmonic oscillator to detect weak classical forces. We show that our framework allows the consideration of sensitivity scheduling as well as estimation strategies for non-stationary signals, leading us to propose corresponding generalizations of the Standard Quantum Limit for force detection.
The development of advanced closed-loop irrigation systems requires accurate soil moisture information. In this work, we address the problem of soil moisture estimation for the agro-hydrological systems in a robust and reliable manner. A nonlinear state-space model is established based on the discretization of the Richards equation to describe the dynamics of agro-hydrological systems. We consider that model parameters are unknown and need to be estimated together with the states simultaneously. We propose a consensus-based estimation mechanism, which comprises two main parts: 1) a distributed extended Kalman filtering algorithm used to estimate several model parameters; and 2) a distributed moving horizon estimation algorithm used to estimate the state variables and one remaining model parameter. Extensive simulations are conducted, and comparisons with existing methods are made to demonstrate the effectiveness and superiority of the proposed approach. In particular, the proposed approach can provide accurate soil moisture estimate even when poor initial guesses of the parameters and the states are used, which can be challenging to be handled using existing algorithms.