47 research outputs found

    Particle filtering in high-dimensional chaotic systems

    Full text link
    We present an efficient particle filtering algorithm for multiscale systems, that is adapted for simple atmospheric dynamics models which are inherently chaotic. Particle filters represent the posterior conditional distribution of the state variables by a collection of particles, which evolves and adapts recursively as new information becomes available. The difference between the estimated state and the true state of the system constitutes the error in specifying or forecasting the state, which is amplified in chaotic systems that have a number of positive Lyapunov exponents. The purpose of the present paper is to show that the homogenization method developed in Imkeller et al. (2011), which is applicable to high dimensional multi-scale filtering problems, along with important sampling and control methods can be used as a basic and flexible tool for the construction of the proposal density inherent in particle filtering. Finally, we apply the general homogenized particle filtering algorithm developed here to the Lorenz'96 atmospheric model that mimics mid-latitude atmospheric dynamics with microscopic convective processes.Comment: 28 pages, 12 figure

    Dimensional reduction in nonlinear filtering: A homogenization approach

    Full text link
    We propose a homogenized filter for multiscale signals, which allows us to reduce the dimension of the system. We prove that the nonlinear filter converges to our homogenized filter with rate ε\sqrt{\varepsilon}. This is achieved by a suitable asymptotic expansion of the dual of the Zakai equation, and by probabilistically representing the correction terms with the help of BDSDEs.Comment: Published in at http://dx.doi.org/10.1214/12-AAP901 the Annals of Applied Probability (http://www.imstat.org/aap/) by the Institute of Mathematical Statistics (http://www.imstat.org

    A stochastic Lagrangian particle model and nonlinear filtering for three dimensional Euler flow with jumps

    Get PDF
    The article of record as published may be found at http://doi.org/10.31390/cosa.5.3.07In this paper we will introduce a stochastic Lagrangian particle model with jumps for the three dimensional Euler flow and study the associated nonlinear filtering problem. We apply results from backward integro-differential equation problem to prove uniqueness of solution to the Zakai equations

    Stochastic Transport in Upper Ocean Dynamics

    Get PDF
    This open access proceedings volume brings selected, peer-reviewed contributions presented at the Stochastic Transport in Upper Ocean Dynamics (STUOD) 2021 Workshop, held virtually and in person at the Imperial College London, UK, September 20–23, 2021. The STUOD project is supported by an ERC Synergy Grant, and led by Imperial College London, the National Institute for Research in Computer Science and Automatic Control (INRIA) and the French Research Institute for Exploitation of the Sea (IFREMER). The project aims to deliver new capabilities for assessing variability and uncertainty in upper ocean dynamics. It will provide decision makers a means of quantifying the effects of local patterns of sea level rise, heat uptake, carbon storage and change of oxygen content and pH in the ocean. Its multimodal monitoring will enhance the scientific understanding of marine debris transport, tracking of oil spills and accumulation of plastic in the sea. All topics of these proceedings are essential to the scientific foundations of oceanography which has a vital role in climate science. Studies convened in this volume focus on a range of fundamental areas, including: Observations at a high resolution of upper ocean properties such as temperature, salinity, topography, wind, waves and velocity; Large scale numerical simulations; Data-based stochastic equations for upper ocean dynamics that quantify simulation error; Stochastic data assimilation to reduce uncertainty. These fundamental subjects in modern science and technology are urgently required in order to meet the challenges of climate change faced today by human society. This proceedings volume represents a lasting legacy of crucial scientific expertise to help meet this ongoing challenge, for the benefit of academics and professionals in pure and applied mathematics, computational science, data analysis, data assimilation and oceanography

    Particle Filtering for Stochastic Navier–Stokes Observed with Linear Additive Noise

    Get PDF
    We consider a nonlinear filtering problem whereby the signal obeys the stochastic Navier–Stokes equations and is observed through a linear mapping with additive noise. The setup is relevant to data assimilation for numerical weather prediction and climate modeling, where similar models are used for unknown ocean or wind velocities. We present a particle filtering methodology that uses likelihood-informed importance proposals, adaptive tempering, and a small number of appropriate Markov chain Monte Carlo steps. We provide a detailed design for each of these steps and show in our numerical examples that they are all crucial in terms of achieving good performance and efficiency

    Parameter estimation, model reduction and quantum filtering

    Get PDF
    This thesis explores the topics of parameter estimation and model reduction in the context of quantum filtering. The last is a mathematically rigorous formulation of continuous quantum measurement, in which a stream of auxiliary quantum systems is used to infer the state of a target quantum system. Fundamental quantum uncertainties appear as noise which corrupts the probe observations and therefore must be filtered in order to extract information about the target system. This is analogous to the classical filtering problem in which techniques of inference are used to process noisy observations of a system in order to estimate its state. Given the clear similarities between the two filtering problems, I devote the beginning of this thesis to a review of classical and quantum probability theory, stochastic calculus and filtering. This allows for a mathematically rigorous and technically adroit presentation of the quantum filtering problem and solution. Given this foundation, I next consider the related problem of quantum parameter estimation, in which one seeks to infer the strength of a parameter that drives the evolution of a probe quantum system. By embedding this problem in the state estimation problem solved by the quantum filter, I present the optimal Bayesian estimator for a parameter when given continuous measurements of the probe system to which it couples. For cases when the probe takes on a finite number of values, I review a set of sufficient conditions for asymptotic convergence of the estimator. For a continuous-valued parameter, I present a computational method called quantum particle filtering for practical estimation of the parameter. Using these methods, I then study the particular problem of atomic magnetometry and review an experimental method for potentially reducing the uncertainty in the estimate of the magnetic field beyond the standard quantum limit. The technique involves double-passing a probe laser field through the atomic system, giving rise to effective non-linearities which enhance the effect of Larmor precession allowing for improved magnetic field estimation. I then turn to the topic of model reduction, which is the search for a reduced computational model of a dynamical system. This is a particularly important task for quantum mechanical systems, whose state grows exponentially in the number of subsystems. In the quantum filtering setting, I study the use of model reduction in developing a feedback controller for continuous-time quantum error correction. By studying the propagation of errors in a noisy quantum memory, I present a computation model which scales polynomially, rather than exponentially, in the number of physical qubits of the system. Although inexact, a feedback controller using this model performs almost indistinguishably from one using the full model. I finally review an exact but polynomial model of collective qubit systems undergoing arbitrary symmetric dynamics which allows for the efficient simulation of spontaneous-emission and related open quantum system phenomenon
    corecore