Local Positioning with Parallelepiped Moving Grid

Similar documents
Gaussian Mixture Filter in Hybrid Navigation

Gaussian Mixture Filter Allowing Negative Weights and its Application to Positioning Using Signal Strength Measurements

Lecture 2: From Linear Regression to Kalman Filter and Beyond

Lecture 2: From Linear Regression to Kalman Filter and Beyond

Lecture 6: Bayesian Inference in SDE Models

Bayesian Methods in Positioning Applications

The Kalman Filter ImPr Talk

Introduction to Mobile Robotics Bayes Filter Particle Filter and Monte Carlo Localization

TSRT14: Sensor Fusion Lecture 8

Lecture 7: Optimal Smoothing

PATTERN RECOGNITION AND MACHINE LEARNING CHAPTER 13: SEQUENTIAL DATA

Tampere University of Technology Tampere Finland

A Comparison of the EKF, SPKF, and the Bayes Filter for Landmark-Based Localization

Posterior Cramer-Rao Lower Bound for Mobile Tracking in Mixed Line-of-Sight/Non Line-of-Sight Conditions

Gaussian Mixture Filters in Hybrid Positioning

Autonomous Mobile Robot Design

1 Kalman Filter Introduction

2D Image Processing (Extended) Kalman and particle filter

CONSTRAINT KALMAN FILTER FOR INDOOR BLUETOOTH LOCALIZATION

Parametric Models. Dr. Shuang LIANG. School of Software Engineering TongJi University Fall, 2012

Particle lter for mobile robot tracking and localisation

Lecture: Gaussian Process Regression. STAT 6474 Instructor: Hongxiao Zhu

Gaussian Process Approximations of Stochastic Differential Equations

Sensor Fusion: Particle Filter

Kalman filtering and friends: Inference in time series models. Herke van Hoof slides mostly by Michael Rubinstein

Kalman Filter Computer Vision (Kris Kitani) Carnegie Mellon University

Introduction to Unscented Kalman Filter

Particle Filters. Outline

Nonlinear Filtering. With Polynomial Chaos. Raktim Bhattacharya. Aerospace Engineering, Texas A&M University uq.tamu.edu

Online tests of Kalman filter consistency

Mobile Robot Localization

AUTOMOTIVE ENVIRONMENT SENSORS

Markov localization uses an explicit, discrete representation for the probability of all position in the state space.

Extension of the Sparse Grid Quadrature Filter

L06. LINEAR KALMAN FILTERS. NA568 Mobile Robotics: Methods & Algorithms

EKF and SLAM. McGill COMP 765 Sept 18 th, 2017

Organization. I MCMC discussion. I project talks. I Lecture.

Partitioned Update Kalman Filter

DRAGON ADVANCED TRAINING COURSE IN ATMOSPHERE REMOTE SENSING. Inversion basics. Erkki Kyrölä Finnish Meteorological Institute

Kalman Filters, Dynamic Bayesian Networks

Time Series Prediction by Kalman Smoother with Cross-Validated Noise Density

CIS 390 Fall 2016 Robotics: Planning and Perception Final Review Questions

(Extended) Kalman Filter

L09. PARTICLE FILTERING. NA568 Mobile Robotics: Methods & Algorithms

Lagrangian Data Assimilation and Manifold Detection for a Point-Vortex Model. David Darmon, AMSC Kayo Ide, AOSC, IPST, CSCAMM, ESSIC

Self Adaptive Particle Filter

Dynamic Multipath Estimation by Sequential Monte Carlo Methods

The Unscented Particle Filter

Mobile Robot Localization

Particle Filters. Pieter Abbeel UC Berkeley EECS. Many slides adapted from Thrun, Burgard and Fox, Probabilistic Robotics

Sequential Monte Carlo Methods for Bayesian Computation

Probabilistic Graphical Models

Bayes Filter Reminder. Kalman Filter Localization. Properties of Gaussians. Gaussians. Prediction. Correction. σ 2. Univariate. 1 2πσ e.

Sequential Monte Carlo in the machine learning toolbox

A Comparison of Particle Filters for Personal Positioning

Introduction to Mobile Robotics Probabilistic Robotics

Strong Lens Modeling (II): Statistical Methods

SEQUENTIAL MONTE CARLO METHODS WITH APPLICATIONS TO COMMUNICATION CHANNELS. A Thesis SIRISH BODDIKURAPATI

L11. EKF SLAM: PART I. NA568 Mobile Robotics: Methods & Algorithms

Introduction to Probabilistic Graphical Models: Exercises

Robotics. Mobile Robotics. Marc Toussaint U Stuttgart

CS491/691: Introduction to Aerial Robotics

Data Fusion Kalman Filtering Self Localization

System identification and sensor fusion in dynamical systems. Thomas Schön Division of Systems and Control, Uppsala University, Sweden.

Inferring biological dynamics Iterated filtering (IF)

ECE295, Data Assimila0on and Inverse Problems, Spring 2015

Dynamic System Identification using HDMR-Bayesian Technique

Outline. 1 Introduction. 2 Problem and solution. 3 Bayesian tracking model of group debris. 4 Simulation results. 5 Conclusions

Autonomous Navigation for Flying Robots

ECE276A: Sensing & Estimation in Robotics Lecture 10: Gaussian Mixture and Particle Filtering

Probabilistic Graphical Models

Robotics. Lecture 4: Probabilistic Robotics. See course website for up to date information.

Terrain Navigation Using the Ambient Magnetic Field as a Map

Localization of DECT mobile phones based on a new nonlinear filtering technique

An introduction to Sequential Monte Carlo

Stochastic Models, Estimation and Control Peter S. Maybeck Volumes 1, 2 & 3 Tables of Contents

Nonlinear Estimation Techniques for Impact Point Prediction of Ballistic Targets

Introduction to Gaussian Processes

Introduction to Mobile Robotics Bayes Filter Kalman Filter

Rao-Blackwellized Particle Filter for Multiple Target Tracking

Cheng Soon Ong & Christian Walder. Canberra February June 2018

Parallel Particle Filter in Julia

Modeling and state estimation Examples State estimation Probabilities Bayes filter Particle filter. Modeling. CSC752 Autonomous Robotic Systems

Recursive Noise Adaptive Kalman Filtering by Variational Bayesian Approximations

A Note on Auxiliary Particle Filters

Large Scale Bayesian Inference

Par$cle Filters Part I: Theory. Peter Jan van Leeuwen Data- Assimila$on Research Centre DARC University of Reading

Robot Localization and Kalman Filters

Implementation of Particle Filter-based Target Tracking

Lecture 9. Time series prediction

Sequential Monte Carlo Samplers for Applications in High Dimensions

Introduction to Machine Learning

Evaluating the Consistency of Estimation

SLAM Techniques and Algorithms. Jack Collier. Canada. Recherche et développement pour la défense Canada. Defence Research and Development Canada

From Bayes to Extended Kalman Filter

Smoothers: Types and Benchmarks

Lecture Particle Filters

Log Gaussian Cox Processes. Chi Group Meeting February 23, 2016

NONLINEAR BAYESIAN FILTERING FOR STATE AND PARAMETER ESTIMATION

Blind Equalization via Particle Filtering

Transcription:

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 Moving Grid Niilo Sirola and Simo Ali-Löytty Tampere University of Technology, Finland

Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 2/?? Personal positioning for example, a mobile handset sources of positioning information: navigation equipment: GPS, IMU cellular network, WLAN, Bluetooth,... digital compass, step counter, barometer aiding data (altitude, cell sector,... ) digital maps etc.

Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 3/?? Challenges large linearization errors scarce measurements underdetermined non-normal noise structure?

Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 3/?? Challenges large linearization errors scarce measurements underdetermined non-normal noise structure make use of all available information?

Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 3/?? Challenges large linearization errors scarce measurements underdetermined non-normal noise structure make use of all available information exploit the time dependency?

Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 3/?? Challenges large linearization errors scarce measurements underdetermined non-normal noise structure make use of all available information exploit the time dependency allow flexible measurement models?

Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 4/?? Discrete nonlinear filtering problem Distribution of initial state x 0 given. Motion model: Measurement model: x k+1 = f k (x k ) + w k y k = h k (x k ) + v k Optimal solution: recursive Bayesian filter

Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 5/?? Recursive Bayesian filter Prediction step: p(x k Y k 1 ) = p(x k x k 1 ) p(x }{{} k 1 Y k 1 ) }{{} R d motion model previous dx k 1

Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 5/?? Recursive Bayesian filter Prediction step: p(x k Y k 1 ) = p(x k x k 1 ) p(x }{{} k 1 Y k 1 ) }{{} R d Update step: motion model p(x k Y k ) p(y k x k ) p(x }{{} k Y k 1 ) measurement previous dx k 1

Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 5/?? Recursive Bayesian filter Prediction step: p(x k Y k 1 ) = p(x k x k 1 ) p(x }{{} k 1 Y k 1 ) }{{} R d Update step: motion model p(x k Y k ) p(y k x k ) p(x }{{} k Y k 1 ) measurement previous dx k 1 In general, exact solution requires infinite time and memory!

Numerical pdf approximations point mass grid mass true pdf Gaussian Monte Carlo Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 6/??

Numerical pdf approximations point mass grid mass true pdf Gaussian Monte Carlo Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 6/??

Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 7/?? Moving grid algorithm 1. prior approximation

Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 7/?? Moving grid algorithm 1. prior approximation 2. measurement likelihood

Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 7/?? Moving grid algorithm 1. prior approximation 2. measurement likelihood 3. approximate likelihood

Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 7/?? Moving grid algorithm 1. prior approximation 2. measurement likelihood 3. approximate likelihood 4. multiply to get posterior

Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 7/?? Moving grid algorithm 1. prior approximation 2. measurement likelihood 3. approximate likelihood 4. multiply to get posterior 5. propagate with motion model

Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 7/?? Moving grid algorithm 1. prior approximation 2. measurement likelihood 3. approximate likelihood 4. multiply to get posterior 5. propagate with motion model 6. repeat from step 2...

Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 8/?? Example run true track reference x 0 base station

Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 8/?? Example run x 0 true track EKF reference base station

Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 8/?? Example run moving grid x 0 true track EKF reference base station

Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 9/?? Grid accuracy 2500 2000 error from track cell radius error / m 1500 1000 500 error from reference 0 0 20 40 60 80 100 120 time / s The 2D error of the moving grid filter mean estimate compared to the cell radius

Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 10/?? How to compare filters? Ideally: compare the posterior distribution to the ideal one In practise we compute the mean and covariance estimates error from true track error from reference track consistency (error vs. estimated covariance)

Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 11/?? Some numbers mean error 95 % CERP inconsistent time grid 200 m 574 m <0.1 % 100 EKF 105 m 370 m 4 % 1 SMC, 10k 75 m 260 m 1 % 30 SMC, 2M 75 m 255 m <0.1 % 10000

Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 12/?? Conclusions grid filter gives rough but reliable results fair and expressive comparison of nonlinear filters still an open problem how to weigh the overall score? compare to true track or reference track? accuracy vs. reliability accuracy vs. computation time?

Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 12/?? Conclusions grid filter gives rough but reliable results fair and expressive comparison of nonlinear filters still an open problem how to weigh the overall score? compare to true track or reference track? accuracy vs. reliability accuracy vs. computation time? Questions?