Kalman Filter I've got it! Sequential Bayes Filtering is a - - PowerPoint PPT Presentation

kalman filter
SMART_READER_LITE
LIVE PREVIEW

Kalman Filter I've got it! Sequential Bayes Filtering is a - - PowerPoint PPT Presentation

Kalman Filter I've got it! Sequential Bayes Filtering is a general approach to state estimation that gets used all over the place. But, implementations like histogram filters or Kalman filters are computationally complex. Is it


slide-1
SLIDE 1

Kalman Filter

  • Sequential Bayes Filtering is a general

approach to state estimation that gets used all over the place.

  • But, implementations like histogram filters or

Kalman filters are computationally complex.

  • Is it always this way? Is Bayes filtering ever

simple?

I've got it!

slide-2
SLIDE 2

Why Kalman filtering?

Kalman filters are particularly well suited for tracking moving objects In general: state estimation

Image: Thrun et al., Probabilistic Robotics Image: UBC, Kevin Murphy Matlab toolbox

slide-3
SLIDE 3

Review of sequential Bayes filtering

Process update: Measurement update:

slide-4
SLIDE 4

Review of sequential Bayes filtering

Image: Thrun et al., Probabilistic Robotics

slide-5
SLIDE 5

Transition function

Recall the state transition function: – this probability can be expressed as a table: Current state Next state Can also be expressed as a function:

  • r:
slide-6
SLIDE 6

Linear system

A linear system is any system where:

(technically, this is a linear Gaussian system)

Also written:

slide-7
SLIDE 7

Linear system: Example

k b m

How get it into this form? Equation of motion:

slide-8
SLIDE 8

Linear system: Example

k b m

How get it into this form? Equation of motion: Integrate forward one timestep:

slide-9
SLIDE 9

Linear system

A linear system is any system where:

(technically, this is a linear Gaussian system)

Also written: Also, assume that the observation function is linear Gaussian:

slide-10
SLIDE 10

Kalman Idea

update initial position x y x y prediction x y measurement x y

Image: Thrun et al., CS233B course notes

slide-11
SLIDE 11

Kalman Idea

Image: Thrun et al., CS233B course notes

prior Measurement evidence posterior

Image: Thrun et al., CS233B course notes

slide-12
SLIDE 12

Gaussians

Univariate Gaussian: Multivariate Gaussian:

slide-13
SLIDE 13

Playing w/ Gaussians

Suppose: Calculate:

x y x y

slide-14
SLIDE 14

In fact

Suppose: Then:

slide-15
SLIDE 15

Illustration

Image: Thrun et al., CS233B course notes

slide-16
SLIDE 16

And

Suppose: Then: Marginal distribution

slide-17
SLIDE 17

Does this remind us of anything?

slide-18
SLIDE 18

Does this remind us of anything?

Process update (discrete): Process update (continuous):

slide-19
SLIDE 19

Does this remind us of anything?

Process update (discrete): Process update (continuous): prior Transition dynamics

slide-20
SLIDE 20

Does this remind us of anything?

Process update (discrete): Process update (continuous): prior Transition dynamics

slide-21
SLIDE 21

Observation update

Observation update: Where:

slide-22
SLIDE 22

Observation update

Observation update: Where:

slide-23
SLIDE 23

Observation update

Observation update: Where:

slide-24
SLIDE 24

Observation update

But we need:

slide-25
SLIDE 25

Another Gaussian identity...

Suppose: Calculate:

slide-26
SLIDE 26

Observation update

But we need:

slide-27
SLIDE 27

To summarize the Kalman filter

Prior: Process update: Measurement update: System:

slide-28
SLIDE 28

Suppose there is an action term...

Prior: Process update: Measurement update: System:

slide-29
SLIDE 29

To summarize the Kalman filter

Prior: Process update: Measurement update: This factor is often called the “Kalman gain”

slide-30
SLIDE 30

Things to note about the Kalman filter

Process update: Measurement update: – covariance update is independent of observation – Kalman is only optimal for linear-Gaussian systems – the distribution “stays” Gaussian through this update – the error term can be thought of as the different between the

  • bservation and the prediction
slide-31
SLIDE 31

Kalman in 1D

Image: Thrun et al., CS233B course notes

prior Measurement evidence posterior

Image: Thrun et al., CS233B course notes

Process update: Measurement update: System:

slide-32
SLIDE 32

Kalman Idea

Image: Thrun et al., CS233B course notes

initial position prediction measurement

˙ x

x

˙ x

x

˙ x

x

next prediciton

˙ x

x

update

˙ x

x

slide-33
SLIDE 33

Example: estimate velocity

Image: Thrun et al., CS233B course notes

past measurements prediction

slide-34
SLIDE 34

Example: filling a tank

Level of tank Fill rate Process: Observation:

slide-35
SLIDE 35

Example: estimate velocity

slide-36
SLIDE 36

But, my system is NON-LINEAR!

What should I do?

slide-37
SLIDE 37

But, my system is NON-LINEAR!

What should I do? Well, there are some options...

slide-38
SLIDE 38

But, my system is NON-LINEAR!

What should I do? Well, there are some options... But none of them are great.

slide-39
SLIDE 39

But, my system is NON-LINEAR!

What should I do? Well, there are some options... But none of them are great. Here's one: the Extended Kalman Filter

slide-40
SLIDE 40

Extended Kalman filter

Take a Taylor expansion: Where: Where:

slide-41
SLIDE 41

Extended Kalman filter

Take a Taylor expansion: Where: Where: Then use the same equations...

slide-42
SLIDE 42

To summarize the EKF

Prior: Process update: Measurement update:

slide-43
SLIDE 43

Extended Kalman filter

Image: Thrun et al., CS233B course notes

slide-44
SLIDE 44

Extended Kalman filter

Image: Thrun et al., CS233B course notes