Nonlinear Filtering. With Polynomial Chaos. Raktim Bhattacharya. Aerospace Engineering, Texas A&M University uq.tamu.edu
|
|
- Norah McDaniel
- 5 years ago
- Views:
Transcription
1 Nonlinear Filtering With Polynomial Chaos Raktim Bhattacharya Aerospace Engineering, Texas A&M University uq.tamu.edu
2 Nonlinear Filtering with PC Problem Setup. Dynamics: ẋ = f(x, ) Sensor Model: ỹ = h(x) + ν where ν is measurement noise with E [ν] = 0, and E [ νν T ] = R. Measurements available at times t k, t k+1, Parametric uncertainty := ( x0 ρ ) Dynamics transformed to ẋ pc = F pc (x pc ) Initial Condition at time t k x pci (t k ) = x0 ϕ i ( )p(k, )d D Workshop on Uncertainty Analysis & Estimation (ACC 2014) 2 /19
3 Filtering with Higher Order Moment Updates 1. Propagation tk=1 x pc (t k+1 ) = x pc (t k ) + F pc (x pc (τ))dτ t k 2. Get Prior Moments from x pc (t k+1 ) recall x pc := vec (X) and X := [x 0 x 1 x N ]. M 1 i M 2 ij M 3 ijk = p = X i0 = X ip X jq ϕ p ϕ q p q X ip X jq X kr ϕ p ϕ q ϕ r and so on... Inner products are with respect to p(t k, ) q r Workshop on Uncertainty Analysis & Estimation (ACC 2014) 3 /19
4 Filtering with Higher Order Moment Updates 3. Update. Incorporate measurements ỹ := ỹ(t k+1 ) and prior moments to get posterior estimates Consider prior state estimate to be ˆx := E [x] = M 1 Let v := ỹ ŷ = h(x) + ν h(ˆx ) Use linear gain K to update moments as K = P xv (P vv ) 1, P xv ij M 1+ M 2+ M 3+ = M 1 + Kv = M 2 + KP vv K T = E [ x i v T j ], P vv ij = E [ v i v T ] j = M 3 + 3K 2 P xvv 3KP xxv K 3 P vvv Workshop on Uncertainty Analysis & Estimation (ACC 2014) 4 /19
5 Filtering with Higher Order Moment Updates 4. Estimation of Posterior PDF max p k+1 ( ) log(p k+1 ( )) d, p k+1 ( ) D subject to p k+1 ( ) d = M 1+ D D Q 3 ( )p k+1 ( ) d = M 3+ Approximate with p k+1 ( ) = D p k+1 ( )d = 1. Q 2 ( )p k+1 ( ) d = M 2+ D Q 4 ( )p k+1 ( ) d = M 4+ D M α i N (µ i, Σ i ), i M α i = 1, p k+1 ( ) 0 α i 0. i Workshop on Uncertainty Analysis & Estimation (ACC 2014) 5 /19
6 Filtering with Higher Order Moment Updates Example Classical duffing oscillator Two state system, x = [x 1, x 2 ] T, Dynamics ẋ 1 = x 2, ẋ 2 = x x 2 x 3 1. Uncertainty in x 0 N ([1, 1], diag(1, 1)) Simulation x 0 = [2, 2] T Scalar measurement model ỹ = x T x + ν, E [ν] = 0 and E [ νν T ] = Workshop on Uncertainty Analysis & Estimation (ACC 2014) 6 /19
7 Results. (a) EKF based estimator with 0.2s update. (b) EKF based estimator with 0.3s update. Figure : Performance of EKF estimators. If data rate is high, no need for nonlinear estimators!! Workshop on Uncertainty Analysis & Estimation (ACC 2014) 7 /19
8 EKF works with faster updates Otherwise we need nonlinear non Gaussian algorithms. (a) gpc based estimator with first two moments updated every 0.3s. (b) gpc based estimator with first three moments updated every 0.5s. Figure : Performance of gpc estimators. If data rate is low, we need nonlinear estimators, with higher order updates! Workshop on Uncertainty Analysis & Estimation (ACC 2014) 8 /19
9 Publications. 1. P. Dutta, R. Bhattacharya, Nonlinear Estimation of Hypersonic Flight Using Polynomial Chaos, AIAA GNC, P. Dutta, R. Bhattacharya, Nonlinear Estimation with Polynomial Chaos and Higher Order Moment Updates, IEEE ACC P. Dutta, R. Bhattacharya, Nonlinear Estimation of Hypersonic State Trajectories in Bayesian Framework with Polynomial Chaos, Journal of Guidance, Control, and Dynamics, vol.33 no.6 ( ), Workshop on Uncertainty Analysis & Estimation (ACC 2014) 9 /19
10 Nonlinear Filtering With Frobenius-Perron Operator Raktim Bhattacharya Aerospace Engineering, Texas A&M University uq.tamu.edu
11 Frobenius-Perron Operator. Given dynamics ẋ = F (t, x), x is augmented state variable captures state and system parameters (including those from KL expansion) p(t 0, x) as the initial state density function. Evolution of density p(t, x) := P t p(t 0, x). P t is defined by p t + (pf ) = 0 Equations ẋ = F (t, x) ṗ = p F Workshop on Uncertainty Analysis & Estimation (ACC 2014) 11 /19
12 Assumptions. Measurements are available at times t 1,, t k 1, t k, t k+1, x k and y k are state and measurement at t k Measurement model y = h(x) + ν E [ν] = 0, E [ νν ] T = R p k ( ) := p(t k, ) p k ( ) is prior at t k p + k ( ) is posterior at t k D x is domain of state augmented Workshop on Uncertainty Analysis & Estimation (ACC 2014) 12 /19
13 A Particle Filter Based Algorithm. 1. Initialize Domain D x is sampled according to p 0 (x) x 0,i samples of r.v. x 0 p 0,i := p 0 (x 0,i ) Recursively apply steps 2,, 6 for k = 1, Workshop on Uncertainty Analysis & Estimation (ACC 2014) 13 /19
14 A Particle Filter Based Algorithm (contd.). 2. Propagate ( xk,i p k,i ) = ( ) tk ( xk 1,i + p k 1,i t k 1 F p F ) dt p k,i because it is prior state PDF Workshop on Uncertainty Analysis & Estimation (ACC 2014) 14 /19
15 A Particle Filter Based Algorithm (contd.) 3. Determine likelihood function p(ỹ k x k = x k,i ) for each grid point i, using Gaussian measurement noise and sensor model y = h(x) + ν. R is the determinant of the covariance matrix of measurement noise l(ỹ k x k = x k,i ) = 1 (2π)m R e 0.5(ỹ k h(x k,i)) T R 1 (ỹ k h(x k,i )), 4. Update: Get Posterior p + k,i := p k(x k = x k,i ỹ k ) = l(ỹ k x k = x k,i )p k (x k = x k,i ) N l(ỹ k x k = x k,j )p k (x k = x k,j ) j=1 Workshop on Uncertainty Analysis & Estimation (ACC 2014) 15 /19
16 A Particle Filter Based Algorithm (contd.) 5. Get State Estimate a Maximum-Likelihood Estimate: Maximize the probability that x k,i = ˆx k ˆx k = mode p + k (x k,i). b Minimum-Variance Estimate: The estimate is the mean of p + k (x k,i) ˆx k = arg min x N N x x k,i 2 p + k (x k,i) = x k,i p + k (x k,i) i=1 i=1 c Minimum-Error Estimate: Minimize maximum x x k,i All same for Gaussian p(t k, x) ˆx = median p + k (x k,i) Workshop on Uncertainty Analysis & Estimation (ACC 2014) 16 /19
17 A Particle Filter Based Algorithm (contd.) 6. Resample Detect degeneracy from p + k (x k,i) p + k (x k,i) < ϵ x k,i is degenerate. Use existing methods for resampling from the new distribution p + k (x k,i). Importance sampling moderate size problems Resampling simple random, multinomial, stratified, systematic Qualitatively, since histogram techniques are not used in determining density functions, this method is less sensitive to the issue of degeneracy. Workshop on Uncertainty Analysis & Estimation (ACC 2014) 17 /19
18 Example. 3 DOF Vinh s Equation Models motion of spacecraft during planetary entry ḣ = V sin(γ) V = ρr 0 V 2 gr 0 sin(γ) 2B c vc 2 γ = ρr0 C L V + gr0 cos(γ) 2B c C D vc 2 ( V R 0 + h 1 V ) R 0 radius of Mars ρ atmospheric density v c escape velocity C L C lift over drag D B c ballistic coefficient h height V velocity γ flight path angle Measurement Model [ ỹ = q = 1 ] 2 ρv 2, Q = kρ 1 2 V 3.15, γ [ ] E [ν] = 0 3 1, E νν T = I 3 scaled Gaussian initial condition uncertainty µ 0 = [54 km, 2.4 km/s, 9 ] T Σ = diag[5.4 km, 240 km/s, 0.9 ] Workshop on Uncertainty Analysis & Estimation (ACC 2014) 18 /19
19 Example (contd.). Compared with generic particle filter and Bootstrap filter All 3 perform equally well FP requires much less number of samples Particle Filter: samples Bootstrap Filter: samples Frobenius-Perron Filter: 7000 samples Generic Particle filter Bootstrap filter FP operator based filter s s s Table : Computational time for each filter Details 1. P. Dutta and R. Bhattacharya, Hypersonic State Estimation Using Frobenius-Perron Operator, AIAA Journal of Guidance, Control, and Dynamics, Volume 34, Number 2, Workshop on Uncertainty Analysis & Estimation (ACC 2014) 19 /19
. Frobenius-Perron Operator ACC Workshop on Uncertainty Analysis & Estimation. Raktim Bhattacharya
.. Frobenius-Perron Operator 2014 ACC Workshop on Uncertainty Analysis & Estimation Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. uq.tamu.edu
More informationNEW ALGORITHMS FOR UNCERTAINTY QUANTIFICATION AND NONLINEAR ESTIMATION OF STOCHASTIC DYNAMICAL SYSTEMS. A Dissertation PARIKSHIT DUTTA
NEW ALGORITHMS FOR UNCERTAINTY QUANTIFICATION AND NONLINEAR ESTIMATION OF STOCHASTIC DYNAMICAL SYSTEMS A Dissertation by PARIKSHIT DUTTA Submitted to the Office of Graduate Studies of Texas A&M University
More informationTSRT14: Sensor Fusion Lecture 8
TSRT14: Sensor Fusion Lecture 8 Particle filter theory Marginalized particle filter Gustaf Hendeby gustaf.hendeby@liu.se TSRT14 Lecture 8 Gustaf Hendeby Spring 2018 1 / 25 Le 8: particle filter theory,
More information1 Kalman Filter Introduction
1 Kalman Filter Introduction You should first read Chapter 1 of Stochastic models, estimation, and control: Volume 1 by Peter S. Maybec (available here). 1.1 Explanation of Equations (1-3) and (1-4) Equation
More informationNonlinear Estimation Techniques for Impact Point Prediction of Ballistic Targets
Nonlinear Estimation Techniques for Impact Point Prediction of Ballistic Targets J. Clayton Kerce a, George C. Brown a, and David F. Hardiman b a Georgia Tech Research Institute, Georgia Institute of Technology,
More informationLecture 7: Optimal Smoothing
Department of Biomedical Engineering and Computational Science Aalto University March 17, 2011 Contents 1 What is Optimal Smoothing? 2 Bayesian Optimal Smoothing Equations 3 Rauch-Tung-Striebel Smoother
More informationData assimilation with and without a model
Data assimilation with and without a model Tim Sauer George Mason University Parameter estimation and UQ U. Pittsburgh Mar. 5, 2017 Partially supported by NSF Most of this work is due to: Tyrus Berry,
More informationData assimilation with and without a model
Data assimilation with and without a model Tyrus Berry George Mason University NJIT Feb. 28, 2017 Postdoc supported by NSF This work is in collaboration with: Tim Sauer, GMU Franz Hamilton, Postdoc, NCSU
More informationA Comparison of the EKF, SPKF, and the Bayes Filter for Landmark-Based Localization
A Comparison of the EKF, SPKF, and the Bayes Filter for Landmark-Based Localization and Timothy D. Barfoot CRV 2 Outline Background Objective Experimental Setup Results Discussion Conclusion 2 Outline
More informationL09. PARTICLE FILTERING. NA568 Mobile Robotics: Methods & Algorithms
L09. PARTICLE FILTERING NA568 Mobile Robotics: Methods & Algorithms Particle Filters Different approach to state estimation Instead of parametric description of state (and uncertainty), use a set of state
More informationLocal Positioning with Parallelepiped Moving Grid
Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 1/?? TA M P E R E U N I V E R S I T Y O F T E C H N O L O G Y M a t h e m a t i c s Local Positioning with Parallelepiped
More informationLecture 6: Bayesian Inference in SDE Models
Lecture 6: Bayesian Inference in SDE Models Bayesian Filtering and Smoothing Point of View Simo Särkkä Aalto University Simo Särkkä (Aalto) Lecture 6: Bayesian Inference in SDEs 1 / 45 Contents 1 SDEs
More informationFrom Bayes to Extended Kalman Filter
From Bayes to Extended Kalman Filter Michal Reinštein Czech Technical University in Prague Faculty of Electrical Engineering, Department of Cybernetics Center for Machine Perception http://cmp.felk.cvut.cz/
More informationUncertainty Quantification for Stochastic Nonlinear Systems using Perron-Frobenius Operator and Karhunen-Loève Expansion
Uncertainty Quantification for Stochastic Nonlinear Systems using Perron-Frobenius Operator and Karhunen-Loève Expansion Parikshit Dutta, Abhishek Halder and Raktim Bhattacharya Abstract In this paper,
More informationA Crash Course on Kalman Filtering
A Crash Course on Kalman Filtering Dan Simon Cleveland State University Fall 2014 1 / 64 Outline Linear Systems Probability State Means and Covariances Least Squares Estimation The Kalman Filter Unknown
More informationAdaptive ensemble Kalman filtering of nonlinear systems
Adaptive ensemble Kalman filtering of nonlinear systems Tyrus Berry George Mason University June 12, 213 : Problem Setup We consider a system of the form: x k+1 = f (x k ) + ω k+1 ω N (, Q) y k+1 = h(x
More informationOrganization. I MCMC discussion. I project talks. I Lecture.
Organization I MCMC discussion I project talks. I Lecture. Content I Uncertainty Propagation Overview I Forward-Backward with an Ensemble I Model Reduction (Intro) Uncertainty Propagation in Causal Systems
More informationTarget tracking and classification for missile using interacting multiple model (IMM)
Target tracking and classification for missile using interacting multiple model (IMM Kyungwoo Yoo and Joohwan Chun KAIST School of Electrical Engineering Yuseong-gu, Daejeon, Republic of Korea Email: babooovv@kaist.ac.kr
More informationBayesian Estimation For Tracking Of Spiraling Reentry Vehicles
Marquette University e-publications@marquette Master's Theses (29 -) Dissertations, Theses, and Professional Projects Bayesian Estimation For Tracking Of Spiraling Reentry Vehicles Juan Esteban Tapiero
More informationIntroduction to Mobile Robotics Bayes Filter Particle Filter and Monte Carlo Localization
Introduction to Mobile Robotics Bayes Filter Particle Filter and Monte Carlo Localization Wolfram Burgard, Cyrill Stachniss, Maren Bennewitz, Kai Arras 1 Motivation Recall: Discrete filter Discretize the
More informationL06. LINEAR KALMAN FILTERS. NA568 Mobile Robotics: Methods & Algorithms
L06. LINEAR KALMAN FILTERS NA568 Mobile Robotics: Methods & Algorithms 2 PS2 is out! Landmark-based Localization: EKF, UKF, PF Today s Lecture Minimum Mean Square Error (MMSE) Linear Kalman Filter Gaussian
More informationIntroduction to Unscented Kalman Filter
Introduction to Unscented Kalman Filter 1 Introdution In many scientific fields, we use certain models to describe the dynamics of system, such as mobile robot, vision tracking and so on. The word dynamics
More informationConstrained State Estimation Using the Unscented Kalman Filter
16th Mediterranean Conference on Control and Automation Congress Centre, Ajaccio, France June 25-27, 28 Constrained State Estimation Using the Unscented Kalman Filter Rambabu Kandepu, Lars Imsland and
More informationA new unscented Kalman filter with higher order moment-matching
A new unscented Kalman filter with higher order moment-matching KSENIA PONOMAREVA, PARESH DATE AND ZIDONG WANG Department of Mathematical Sciences, Brunel University, Uxbridge, UB8 3PH, UK. Abstract This
More informationKalman Filter Computer Vision (Kris Kitani) Carnegie Mellon University
Kalman Filter 16-385 Computer Vision (Kris Kitani) Carnegie Mellon University Examples up to now have been discrete (binary) random variables Kalman filtering can be seen as a special case of a temporal
More informationIN particle filter (PF) applications, knowledge of the computational
Complexity Analysis of the Marginalized Particle Filter Rickard Karlsson, Thomas Schön and Fredrik Gustafsson, Member IEEE Abstract In this paper the computational complexity of the marginalized particle
More informationECE276A: Sensing & Estimation in Robotics Lecture 10: Gaussian Mixture and Particle Filtering
ECE276A: Sensing & Estimation in Robotics Lecture 10: Gaussian Mixture and Particle Filtering Lecturer: Nikolay Atanasov: natanasov@ucsd.edu Teaching Assistants: Siwei Guo: s9guo@eng.ucsd.edu Anwesan Pal:
More informationStochastic Collocation Methods for Polynomial Chaos: Analysis and Applications
Stochastic Collocation Methods for Polynomial Chaos: Analysis and Applications Dongbin Xiu Department of Mathematics, Purdue University Support: AFOSR FA955-8-1-353 (Computational Math) SF CAREER DMS-64535
More information6.4 Kalman Filter Equations
6.4 Kalman Filter Equations 6.4.1 Recap: Auxiliary variables Recall the definition of the auxiliary random variables x p k) and x m k): Init: x m 0) := x0) S1: x p k) := Ak 1)x m k 1) +uk 1) +vk 1) S2:
More informationLagrangian Data Assimilation and Its Application to Geophysical Fluid Flows
Lagrangian Data Assimilation and Its Application to Geophysical Fluid Flows Laura Slivinski June, 3 Laura Slivinski (Brown University) Lagrangian Data Assimilation June, 3 / 3 Data Assimilation Setup:
More informationNew Advances in Uncertainty Analysis and Estimation
New Advances in Uncertainty Analysis and Estimation Overview: Both sensor observation data and mathematical models are used to assist in the understanding of physical dynamic systems. However, observational
More informationPATTERN RECOGNITION AND MACHINE LEARNING CHAPTER 13: SEQUENTIAL DATA
PATTERN RECOGNITION AND MACHINE LEARNING CHAPTER 13: SEQUENTIAL DATA Contents in latter part Linear Dynamical Systems What is different from HMM? Kalman filter Its strength and limitation Particle Filter
More informationPROPAGATION OF PDF BY SOLVING THE KOLMOGOROV EQUATION
PROPAGATION OF PDF BY SOLVING THE KOLMOGOROV EQUATION Puneet Singla Department of Mechanical & Aerospace Engineering University at Buffalo, Buffalo, NY-1460 Workshop: New Advances in Uncertainty Analysis
More informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond January 18, 2017 Contents 1 Batch and Recursive Estimation 2 Towards Bayesian Filtering 3 Kalman Filter and Bayesian Filtering and Smoothing
More information4 Derivations of the Discrete-Time Kalman Filter
Technion Israel Institute of Technology, Department of Electrical Engineering Estimation and Identification in Dynamical Systems (048825) Lecture Notes, Fall 2009, Prof N Shimkin 4 Derivations of the Discrete-Time
More informationParticle Filters. Outline
Particle Filters M. Sami Fadali Professor of EE University of Nevada Outline Monte Carlo integration. Particle filter. Importance sampling. Degeneracy Resampling Example. 1 2 Monte Carlo Integration Numerical
More informationExtended Kalman Filter Tutorial
Extended Kalman Filter Tutorial Gabriel A. Terejanu Department of Computer Science and Engineering University at Buffalo, Buffalo, NY 14260 terejanu@buffalo.edu 1 Dynamic process Consider the following
More informationCERTAIN THOUGHTS ON UNCERTAINTY ANALYSIS FOR DYNAMICAL SYSTEMS
CERTAIN THOUGHTS ON UNCERTAINTY ANALYSIS FOR DYNAMICAL SYSTEMS Puneet Singla Assistant Professor Department of Mechanical & Aerospace Engineering University at Buffalo, Buffalo, NY-1426 Probabilistic Analysis
More informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond Department of Biomedical Engineering and Computational Science Aalto University January 26, 2012 Contents 1 Batch and Recursive Estimation
More informationAUTOMOTIVE ENVIRONMENT SENSORS
AUTOMOTIVE ENVIRONMENT SENSORS Lecture 5. Localization BME KÖZLEKEDÉSMÉRNÖKI ÉS JÁRMŰMÉRNÖKI KAR 32708-2/2017/INTFIN SZÁMÚ EMMI ÁLTAL TÁMOGATOTT TANANYAG Related concepts Concepts related to vehicles moving
More information2D Image Processing (Extended) Kalman and particle filter
2D Image Processing (Extended) Kalman and particle filter Prof. Didier Stricker Dr. Gabriele Bleser Kaiserlautern University http://ags.cs.uni-kl.de/ DFKI Deutsches Forschungszentrum für Künstliche Intelligenz
More informationKalman Filters with Uncompensated Biases
Kalman Filters with Uncompensated Biases Renato Zanetti he Charles Stark Draper Laboratory, Houston, exas, 77058 Robert H. Bishop Marquette University, Milwaukee, WI 53201 I. INRODUCION An underlying assumption
More informationEnKF-based particle filters
EnKF-based particle filters Jana de Wiljes, Sebastian Reich, Wilhelm Stannat, Walter Acevedo June 20, 2017 Filtering Problem Signal dx t = f (X t )dt + 2CdW t Observations dy t = h(x t )dt + R 1/2 dv t.
More informationA6523 Modeling, Inference, and Mining Jim Cordes, Cornell University
A6523 Modeling, Inference, and Mining Jim Cordes, Cornell University Lecture 19 Modeling Topics plan: Modeling (linear/non- linear least squares) Bayesian inference Bayesian approaches to spectral esbmabon;
More informationGaussian Mixture Filter in Hybrid Navigation
Digest of TISE Seminar 2007, editor Pertti Koivisto, pages 1-5. Gaussian Mixture Filter in Hybrid Navigation Simo Ali-Löytty Institute of Mathematics Tampere University of Technology simo.ali-loytty@tut.fi
More informationSTATISTICAL ORBIT DETERMINATION
STATISTICAL ORBIT DETERMINATION Satellite Tracking Example of SNC and DMC ASEN 5070 LECTURE 6 4.08.011 1 We will develop a simple state noise compensation (SNC) algorithm. This algorithm adds process noise
More informationThe Kalman filter. Chapter 6
Chapter 6 The Kalman filter In the last chapter, we saw that in Data Assimilation, we ultimately desire knowledge of the full a posteriori p.d.f., that is the conditional p.d.f. of the state given the
More informationPosition Estimation and Modeling of a Flexible Industrial Robot
Position Estimation and Modeling of a Flexible Industrial Robot Rickard Karlsson, Mikael Norrlöf, Division of Automatic Control Department of Electrical Engineering Linköpings universitet, SE-581 83 Linköping,
More informationUAVBook Supplement Full State Direct and Indirect EKF
UAVBook Supplement Full State Direct and Indirect EKF Randal W. Beard March 14, 217 This supplement will explore alternatives to the state estimation scheme presented in the book. In particular, we will
More informationSuggestions for Making Useful the Uncertainty Quantification Results from CFD Applications
Suggestions for Making Useful the Uncertainty Quantification Results from CFD Applications Thomas A. Zang tzandmands@wildblue.net Aug. 8, 2012 CFD Futures Conference: Zang 1 Context The focus of this presentation
More informationGaussian Process Approximations of Stochastic Differential Equations
Gaussian Process Approximations of Stochastic Differential Equations Cédric Archambeau Dan Cawford Manfred Opper John Shawe-Taylor May, 2006 1 Introduction Some of the most complex models routinely run
More informationRaktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Preliminaries
. AERO 632: of Advance Flight Control System. Preliminaries Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Preliminaries Signals & Systems Laplace
More informationAutonomous Navigation for Flying Robots
Computer Vision Group Prof. Daniel Cremers Autonomous Navigation for Flying Robots Lecture 6.2: Kalman Filter Jürgen Sturm Technische Universität München Motivation Bayes filter is a useful tool for state
More informationEKF, UKF. Pieter Abbeel UC Berkeley EECS. Many slides adapted from Thrun, Burgard and Fox, Probabilistic Robotics
EKF, UKF Pieter Abbeel UC Berkeley EECS Many slides adapted from Thrun, Burgard and Fox, Probabilistic Robotics Kalman Filter Kalman Filter = special case of a Bayes filter with dynamics model and sensory
More informationEKF, UKF. Pieter Abbeel UC Berkeley EECS. Many slides adapted from Thrun, Burgard and Fox, Probabilistic Robotics
EKF, UKF Pieter Abbeel UC Berkeley EECS Many slides adapted from Thrun, Burgard and Fox, Probabilistic Robotics Kalman Filter Kalman Filter = special case of a Bayes filter with dynamics model and sensory
More informationEvery real system has uncertainties, which include system parametric uncertainties, unmodeled dynamics
Sensitivity Analysis of Disturbance Accommodating Control with Kalman Filter Estimation Jemin George and John L. Crassidis University at Buffalo, State University of New York, Amherst, NY, 14-44 The design
More informationKalman filtering and friends: Inference in time series models. Herke van Hoof slides mostly by Michael Rubinstein
Kalman filtering and friends: Inference in time series models Herke van Hoof slides mostly by Michael Rubinstein Problem overview Goal Estimate most probable state at time k using measurement up to time
More informationComputer Vision Group Prof. Daniel Cremers. 9. Gaussian Processes - Regression
Group Prof. Daniel Cremers 9. Gaussian Processes - Regression Repetition: Regularized Regression Before, we solved for w using the pseudoinverse. But: we can kernelize this problem as well! First step:
More informationGaussian Processes. Le Song. Machine Learning II: Advanced Topics CSE 8803ML, Spring 2012
Gaussian Processes Le Song Machine Learning II: Advanced Topics CSE 8803ML, Spring 01 Pictorial view of embedding distribution Transform the entire distribution to expected features Feature space Feature
More informationParameter Estimation
1 / 44 Parameter Estimation Saravanan Vijayakumaran sarva@ee.iitb.ac.in Department of Electrical Engineering Indian Institute of Technology Bombay October 25, 2012 Motivation System Model used to Derive
More informationSequential Monte Carlo Methods (for DSGE Models)
Sequential Monte Carlo Methods (for DSGE Models) Frank Schorfheide University of Pennsylvania, PIER, CEPR, and NBER October 23, 2017 Some References These lectures use material from our joint work: Tempered
More informationSequential Monte Carlo and Particle Filtering. Frank Wood Gatsby, November 2007
Sequential Monte Carlo and Particle Filtering Frank Wood Gatsby, November 2007 Importance Sampling Recall: Let s say that we want to compute some expectation (integral) E p [f] = p(x)f(x)dx and we remember
More informationComputer Vision Group Prof. Daniel Cremers. 4. Gaussian Processes - Regression
Group Prof. Daniel Cremers 4. Gaussian Processes - Regression Definition (Rep.) Definition: A Gaussian process is a collection of random variables, any finite number of which have a joint Gaussian distribution.
More informationOn Underweighting Nonlinear Measurements
On Underweighting Nonlinear Measurements Renato Zanetti The Charles Stark Draper Laboratory, Houston, Texas 7758 Kyle J. DeMars and Robert H. Bishop The University of Texas at Austin, Austin, Texas 78712
More informationExtension of Farrenkopf Steady-State Solutions with Estimated Angular Rate
Extension of Farrenopf Steady-State Solutions with Estimated Angular Rate Andrew D. Dianetti and John L. Crassidis University at Buffalo, State University of New Yor, Amherst, NY 46-44 Steady-state solutions
More informationBayesian Statistics and Data Assimilation. Jonathan Stroud. Department of Statistics The George Washington University
Bayesian Statistics and Data Assimilation Jonathan Stroud Department of Statistics The George Washington University 1 Outline Motivation Bayesian Statistics Parameter Estimation in Data Assimilation Combined
More informationA Novel Gaussian Sum Filter Method for Accurate Solution to Nonlinear Filtering Problem
A Novel Gaussian Sum Filter Method for Accurate Solution to Nonlinear Filtering Problem Gabriel Terejanu a Puneet Singla b Tarunraj Singh b Peter D. Scott a Graduate Student Assistant Professor Professor
More informationUsing the Kalman Filter to Estimate the State of a Maneuvering Aircraft
1 Using the Kalman Filter to Estimate the State of a Maneuvering Aircraft K. Meier and A. Desai Abstract Using sensors that only measure the bearing angle and range of an aircraft, a Kalman filter is implemented
More informationUncertainty Evolution In Stochastic Dynamic Models Using Polynomial Chaos
Noname manuscript No. (will be inserted by the editor) Uncertainty Evolution In Stochastic Dynamic Models Using Polynomial Chaos Umamaheswara Konda Puneet Singla Tarunraj Singh Peter Scott Received: date
More informationLeast Squares. Ken Kreutz-Delgado (Nuno Vasconcelos) ECE 175A Winter UCSD
Least Squares Ken Kreutz-Delgado (Nuno Vasconcelos) ECE 75A Winter 0 - UCSD (Unweighted) Least Squares Assume linearity in the unnown, deterministic model parameters Scalar, additive noise model: y f (
More informationInvariant Extended Kalman Filter: Theory and application to a velocity-aided estimation problem
Invariant Extene Kalman Filter: Theory an application to a velocity-aie estimation problem S. Bonnabel (Mines ParisTech) Joint work with P. Martin (Mines ParisTech) E. Salaun (Georgia Institute of Technology)
More informationTime Series Prediction by Kalman Smoother with Cross-Validated Noise Density
Time Series Prediction by Kalman Smoother with Cross-Validated Noise Density Simo Särkkä E-mail: simo.sarkka@hut.fi Aki Vehtari E-mail: aki.vehtari@hut.fi Jouko Lampinen E-mail: jouko.lampinen@hut.fi Abstract
More informationProportional, Integral & Derivative Control Design. Raktim Bhattacharya
AERO 422: Active Controls for Aerospace Vehicles Proportional, ntegral & Derivative Control Design Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University
More informationProbabilistic Fundamentals in Robotics. DAUIN Politecnico di Torino July 2010
Probabilistic Fundamentals in Robotics Gaussian Filters Basilio Bona DAUIN Politecnico di Torino July 2010 Course Outline Basic mathematical framework Probabilistic models of mobile robots Mobile robot
More informationNonlinear Estimation with Perron-Frobenius Operator and Karhunen-Loève Expansion
Nonlinear Estimation with Perron-Frobenius Operator and Karhunen-Loève Expansion Parikshit Dutta a, Abhishek Halder b, Raktim Bhattacharya b. a INRIA Rhône Alpes & Laboratoire Jean Kuntzmann, 655 avenue
More informationExtension of the Sparse Grid Quadrature Filter
Extension of the Sparse Grid Quadrature Filter Yang Cheng Mississippi State University Mississippi State, MS 39762 Email: cheng@ae.msstate.edu Yang Tian Harbin Institute of Technology Harbin, Heilongjiang
More informationA New Nonlinear Filtering Method for Ballistic Target Tracking
th International Conference on Information Fusion Seattle, WA, USA, July 6-9, 9 A New Nonlinear Filtering Method for Ballistic arget racing Chunling Wu Institute of Electronic & Information Engineering
More informationSimultaneous Localization and Mapping (SLAM) Corso di Robotica Prof. Davide Brugali Università degli Studi di Bergamo
Simultaneous Localization and Mapping (SLAM) Corso di Robotica Prof. Davide Brugali Università degli Studi di Bergamo Introduction SLAM asks the following question: Is it possible for an autonomous vehicle
More informationV Control Distance with Dynamics Parameter Uncertainty
V Control Distance with Dynamics Parameter Uncertainty Marcus J. Holzinger Increasing quantities of active spacecraft and debris in the space environment pose substantial risk to the continued safety of
More informationThe Kalman Filter. Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience. Sarah Dance
The Kalman Filter Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience Sarah Dance School of Mathematical and Physical Sciences, University of Reading s.l.dance@reading.ac.uk July
More informationRESEARCH ARTICLE. Online quantization in nonlinear filtering
Journal of Statistical Computation & Simulation Vol. 00, No. 00, Month 200x, 3 RESEARCH ARTICLE Online quantization in nonlinear filtering A. Feuer and G. C. Goodwin Received 00 Month 200x; in final form
More informationROBUST CONSTRAINED ESTIMATION VIA UNSCENTED TRANSFORMATION. Pramod Vachhani a, Shankar Narasimhan b and Raghunathan Rengaswamy a 1
ROUST CONSTRINED ESTIMTION VI UNSCENTED TRNSFORMTION Pramod Vachhani a, Shankar Narasimhan b and Raghunathan Rengaswamy a a Department of Chemical Engineering, Clarkson University, Potsdam, NY -3699, US.
More informationInitial Orbit Determination Using Stereoscopic Imaging and Gaussian Mixture Models
Initial Orbit Determination Using Stereoscopic Imaging and Gaussian Mixture Models Keith A. LeGrand K. DeMars, H.J. Pernicka Department of Mechanical and Aerospace Engineering Missouri University of Science
More informationState Estimation Based on Nested Particle Filters
Cleveland State University EngagedScholarship@CSU ETD Archive 2013 State Estimation Based on Nested Particle Filters Swathi Srinivasan Cleveland State University How does access to this work benefit you?
More informationNonlinear State Estimation! Particle, Sigma-Points Filters!
Nonlinear State Estimation! Particle, Sigma-Points Filters! Robert Stengel! Optimal Control and Estimation, MAE 546! Princeton University, 2017!! Particle filter!! Sigma-Points Unscented Kalman ) filter!!
More informationFactor Analysis and Kalman Filtering (11/2/04)
CS281A/Stat241A: Statistical Learning Theory Factor Analysis and Kalman Filtering (11/2/04) Lecturer: Michael I. Jordan Scribes: Byung-Gon Chun and Sunghoon Kim 1 Factor Analysis Factor analysis is used
More informationMoving Horizon Estimation with Dynamic Programming
Cleveland State University EngagedScholarship@CSU ETD Archive 2013 Moving Horizon Estimation with Dynamic Programming Mohan Kumar Ramalingam Cleveland State University How does access to this work benefit
More informationEvolution Strategies Based Particle Filters for Fault Detection
Evolution Strategies Based Particle Filters for Fault Detection Katsuji Uosaki, Member, IEEE, and Toshiharu Hatanaka, Member, IEEE Abstract Recent massive increase of the computational power has allowed
More informationContent.
Content Fundamentals of Bayesian Techniques (E. Sucar) Bayesian Filters (O. Aycard) Definition & interests Implementations Hidden Markov models Discrete Bayesian Filters or Markov localization Kalman filters
More informationWhy do we care? Examples. Bayes Rule. What room am I in? Handling uncertainty over time: predicting, estimating, recognizing, learning
Handling uncertainty over time: predicting, estimating, recognizing, learning Chris Atkeson 004 Why do we care? Speech recognition makes use of dependence of words and phonemes across time. Knowing where
More informationData Assimilation for Dispersion Models
Data Assimilation for Dispersion Models K. V. Umamaheswara Reddy Dept. of Mechanical and Aerospace Engg. State University of New Yor at Buffalo Buffalo, NY, U.S.A. venatar@buffalo.edu Yang Cheng Dept.
More informationCertain Thoughts on Uncertainty Analysis for Dynamical Systems
Certain Thoughts on Uncertainty Analysis for Dynamical Systems!"#$$%&'(#)*+&!"#$%"&$%'()$%(*+&$,$-$+.(",/(0,&122341,&(5$#$67(8'.&1-.(!"#$%%&'()*+,-.+/01'&2+,304 5,#')67,-642849,:!'-(:'&4;4
More informationA Bank of Maximum A Posteriori (MAP) Estimators for Target Tracking
A Bank of Maximum A Posteriori (MAP) Estimators for Target Tracking Guoquan Huang, Ke Zhou, Nikolas Trawny, and Stergios I. Roumeliotis Abstract Nonlinear estimation problems, such as range-only and bearing-only
More informationGaussian with mean ( µ ) and standard deviation ( σ)
Slide from Pieter Abbeel Gaussian with mean ( µ ) and standard deviation ( σ) 10/6/16 CSE-571: Robotics X ~ N( µ, σ ) Y ~ N( aµ + b, a σ ) Y = ax + b + + + + 1 1 1 1 1 1 1 1 1 1, ~ ) ( ) ( ), ( ~ ), (
More informationBayes Filter Reminder. Kalman Filter Localization. Properties of Gaussians. Gaussians. Prediction. Correction. σ 2. Univariate. 1 2πσ e.
Kalman Filter Localization Bayes Filter Reminder Prediction Correction Gaussians p(x) ~ N(µ,σ 2 ) : Properties of Gaussians Univariate p(x) = 1 1 2πσ e 2 (x µ) 2 σ 2 µ Univariate -σ σ Multivariate µ Multivariate
More informationChris Bishop s PRML Ch. 8: Graphical Models
Chris Bishop s PRML Ch. 8: Graphical Models January 24, 2008 Introduction Visualize the structure of a probabilistic model Design and motivate new models Insights into the model s properties, in particular
More informationROBOTICS 01PEEQW. Basilio Bona DAUIN Politecnico di Torino
ROBOTICS 01PEEQW Basilio Bona DAUIN Politecnico di Torino Probabilistic Fundamentals in Robotics Gaussian Filters Course Outline Basic mathematical framework Probabilistic models of mobile robots Mobile
More informationMini-Course 07 Kalman Particle Filters. Henrique Massard da Fonseca Cesar Cunha Pacheco Wellington Bettencurte Julio Dutra
Mini-Course 07 Kalman Particle Filters Henrique Massard da Fonseca Cesar Cunha Pacheco Wellington Bettencurte Julio Dutra Agenda State Estimation Problems & Kalman Filter Henrique Massard Steady State
More informationI. INTRODUCTION 338 IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS VOL. 33, NO. 1 JANUARY 1997
VI. CONCLUSION We have shown that dilution of precision terms for relative positioning, using double-difference processing of GPS satellite signals, are bounded by the corresponding dilution of precision
More informationTactical Ballistic Missile Tracking using the Interacting Multiple Model Algorithm
Tactical Ballistic Missile Tracking using the Interacting Multiple Model Algorithm Robert L Cooperman Raytheon Co C 3 S Division St Petersburg, FL Robert_L_Cooperman@raytheoncom Abstract The problem of
More information