scholarly journals Multiple sigma-point Kalman smoothers for high-dimensional state-space models

Author(s):  
Jordi Vila-Valls ◽  
Pau Closas ◽  
Angel F. Garcia-Fernandez ◽  
Carles Fernandez-Prades
2021 ◽  
Vol 31 (5) ◽  
Author(s):  
Jacob Vorstrup Goldman ◽  
Sumeetpal S. Singh

AbstractWe propose a novel blocked version of the continuous-time bouncy particle sampler of Bouchard-Côté et al. (J Am Stat Assoc 113(522):855–867, 2018) which is applicable to any differentiable probability density. This alternative implementation is motivated by blocked Gibbs sampling for state-space models (Singh et al. in Biometrika 104(4):953–969, 2017) and leads to significant improvement in terms of effective sample size per second, and furthermore, allows for significant parallelization of the resulting algorithm. The new algorithms are particularly efficient for latent state inference in high-dimensional state-space models, where blocking in both space and time is necessary to avoid degeneracy of MCMC. The efficiency of our blocked bouncy particle sampler, in comparison with both the standard implementation of the bouncy particle sampler and the particle Gibbs algorithm of Andrieu et al. (J R Stat Soc Ser B Stat Methodol 72(3):269–342, 2010), is illustrated numerically for both simulated data and a challenging real-world financial dataset.


2010 ◽  
Vol 30 (2) ◽  
pp. 192-215 ◽  
Author(s):  
Alexander Shkolnik ◽  
Michael Levashov ◽  
Ian R. Manchester ◽  
Russ Tedrake

A motion planning algorithm is described for bounding over rough terrain with the LittleDog robot. Unlike walking gaits, bounding is highly dynamic and cannot be planned with quasi-steady approximations. LittleDog is modeled as a planar five-link system, with a 16-dimensional state space; computing a plan over rough terrain in this high-dimensional state space that respects the kinodynamic constraints due to underactuation and motor limits is extremely challenging. Rapidly Exploring Random Trees (RRTs) are known for fast kinematic path planning in high-dimensional configuration spaces in the presence of obstacles, but search efficiency degrades rapidly with the addition of challenging dynamics. A computationally tractable planner for bounding was developed by modifying the RRT algorithm by using: (1) motion primitives to reduce the dimensionality of the problem; (2) Reachability Guidance, which dynamically changes the sampling distribution and distance metric to address differential constraints and discontinuous motion primitive dynamics; and (3) sampling with a Voronoi bias in a lower-dimensional “task space” for bounding. Short trajectories were demonstrated to work on the robot, however open-loop bounding is inherently unstable. A feedback controller based on transverse linearization was implemented, and shown in simulation to stabilize perturbations in the presence of noise and time delays.


Author(s):  
Oryiema Robert ◽  
David Angwenyi ◽  
Kevin Midenyo

There are several functional forms for non-linear dynamical filters. Extended Kalman filters are algorithms that are used to estimate more accurate values of unknown quantities of internal dynamical systems from a sequence of noisy observation measured over a period of time. This filtering process becomes computationally expensive when subjected to high dimensional data which consequently has a negative impact on the filter performance in real time. This is because integration of the equation of evolution of covariances is extremely costly, especially when the dimension of the problem is huge which is the case in numerical weather prediction.This study has developed a new filter, the First order Extended Ensemble Filter (FoEEF), with a new extended innovation process to improve on the measurement and be able to estimate the state value of high dimensional data. We propose to estimate the covariances empirically, which lends the filter amenable to large dimensional models. The new filter is derived from stochastic state-space models and its performance is tested using Lorenz 63 system of ordinary differential equations and Matlab software.The performance of the newly developed filter is then compared with the performances of three other filters, that is, Bootstrap particle Filter (BPF), First order Extended Kalman Bucy Filter (FoEKBF) and Second order Extended Kalman Bucy Filter (SoEKBF).The performance of the FoEEF improves with the increase in ensemble size. Even with as low number of ensembles as 40, the FoEEF performs as good as the FoEKBF and SoEKBF. This shows, that the proposed filter can register a good performance when used in high-dimensional state-space models.


Acoustics ◽  
2021 ◽  
Vol 3 (3) ◽  
pp. 581-594
Author(s):  
Art J. R. Pelling ◽  
Ennes Sarradj

State-space models have been successfully employed for model order reduction and control purposes in acoustics in the past. However, due to the cubic complexity of the singular value decomposition, which makes up the core of many subspace system identification (SSID) methods, the construction of large scale state-space models from high-dimensional measurement data has been problematic in the past. Recent advances of numerical linear algebra have brought forth computationally efficient randomized rank-revealing matrix factorizations and it has been shown that these factorizations can be used to enhance SSID methods such as the Eigensystem Realization Algorithm (ERA). In this paper, we demonstrate the applicability of the so-called generalized ERA to acoustical systems and high-dimensional input data by means of an example. Furthermore, we introduce a new efficient method of forced response computation that relies on a state-space model in quasi-diagonal form. Numerical experiments reveal that our proposed method is more efficient than previous state-space methods and can even outperform frequency domain convolutions in certain scenarios.


Sign in / Sign up

Export Citation Format

Share Document