1 / 23

Lirong Xia

Approximate inference: Particle filter. Lirong Xia. Tue, April 1, 2014. Reminder. Friday April 4 Project 3 due Project 4 (HMM) out. Recap: Reasoning over Time. Markov models. p( X 1 ) p( X|X -1 ). p(E |X ). Hidden Markov models.

Download Presentation

Lirong Xia

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Approximate inference: Particle filter Lirong Xia Tue, April 1, 2014

  2. Reminder • Friday April 4 • Project 3 due • Project 4 (HMM) out

  3. Recap: Reasoning over Time • Markov models p(X1) p(X|X-1) p(E|X) • Hidden Markov models • Filtering: Given time t and evidences e1,…,et, compute p(Xt|e1:t)

  4. Filtering algorithm • Notation • B(Xt-1)=p(Xt-1|e1:t-1) • B’(Xt)=p(Xt|e1:t-1) • Each time step, we start with p(Xt-1 | previous evidence): • Elapse of time B’(Xt)=Σxt-1p(Xt|xt-1)B(xt-1) • Observe B(Xt) ∝p(et|Xt)B’(Xt) • Renormalize B(Xt)

  5. Prior Sampling (w/o evidences) • Given a Bayesian network (graph and CPTs) • Generate values of random variables sequentially • According to an order compatible with the graph • In each step, generate the value of the top unvalued variable given previously generated values • This sampling procedure is consistent Samples: +c, -s, +r, +w -c, +s, -r, +w

  6. Rejection Sampling • We want p(C|+s) • Do prior sampling, but ignore (reject) samples which don’t have S=+s • It is also consistent for conditional probabilities (i.e., correct in the limit) • Problem • If p(+s) is too small then many samples are wasted +c, -s, +r, +w +c, +s, +r, +w -c, +s, +r, -w +c, -s, +r, +w -c, -s, -r, +w

  7. Likelihood Weighting • Step 1: sampling unvalued variables z (given evidences e) • Step 2: fix the evidences • Step 3: the generated samples have weight • Likelihood weighting is consistent z: C, R e: S, W

  8. Today • Particle filtering • Dynamic Bayesian networks • Viterbi algorithm

  9. Particle Filtering • Sometimes |X| is too big to use exact inference • |X| may be too big to even store B(X) • E.g. X is continuous • |X|2 may be too big to do updates • Solution: approximate inference • Track samples of X, not all values • Samples are called particles • Time per step is linear in the number of samples • But: number needed may be large • In memory: list of particles • This is how robot localization works in practice

  10. Representation: Particles • p(X) is now represented by a list of N particles (samples) • Generally, N << |X| • Storing map from X to counts would defeat the point • p(x) approximated by number of particles with value x • Many x will have p(x)=0! • More particles, more accuracy • For now, all particles have a weight of 1 Particles: (3,3) (2,3) (3,3) (3,2) (3,3) (3,2) (2,1) (3,3) (3,3) (2,1)

  11. Particle Filtering: Elapse Time • Each particle is moved by sampling its next position from the transition model x’= sample(p(X’|x)) • This is like prior sampling – samples’ frequencies reflect the transition probabilities • Here, most samples move clockwise, but some move in another direction or stay in place • This captures the passage of time • If we have enough samples, close to the exact values before and after (consistent)

  12. Particle Filtering: Observe • Slightly trickier: • Don’t do rejection sampling (why not?) • Likelihood weighting • Note that, as before, the probabilities don’t sum to one, since most have been downweighted • in fact they sum to an approximation of p(e)

  13. Particle Filtering: Resample Old Particles: (3,3) w=0.1 (2,1) w=0.9 (2,1) w=0.9 (3,1) w=0.4 (3,2) w=0.3 (2,2) w=0.4 (1,1) w=0.4 (3,1) w=0.4 (2,1) w=0.9 (3,2) w=0.3 New Particles: (2,1) w=1 (2,1) w=1 (2,1) w=1 (3,2) w=1 (2,2) w=1 (2,1) w=1 (1,1) w=1 (3,1) w=1 (2,1) w=1 (1,1) w=1 • Rather than tracking weighted samples, we resample • N times, we choose from our weighted sample distribution (i.e. draw with replacement) • This is analogous to renormalizing the distribution • Now the update is complete for this time step, continue with the next one

  14. Forward algorithm vs. particle filtering Forward algorithm Particle filtering • Elapse of time B’(Xt)=Σxt-1p(Xt|xt-1)B(xt-1) • Observe B(Xt) ∝p(et|Xt)B’(Xt) • Renormalize B(xt) sum up to 1 • Elapse of time x--->x’ • Observe w(x’)=p(et|x) • Resample resample N particles

  15. Robot Localization • In robot localization: • We know the map, but not the robot’s position • Observations may be vectors of range finder readings • State space and readings are typically continuous (works basically like a very fine grid) and so we cannot store B(X) • Particle filtering is a main technique

  16. Dynamic Bayes Nets (DBNs) • We want to track multiple variables over time, using multiple sources of evidence • Idea: repeat a fixed Bayes net structure at each time • Variables from time t can condition on those from t-1 • DBNs with evidence at leaves are HMMs

  17. Exact Inference in DBNs • Variable elimination applies to dynamic Bayes nets • Procedure: “unroll” the network for T time steps, then eliminate variables until p(Xt|e1:t) is computed • Online belief updates: eliminate all variables from the previous time step; store factors for current time only

  18. DBN Particle Filters • A particle is a complete sample for a time step • Initialize: generate prior samples for the t=1 Bayes net • Example particle: G1a=(3,3) G1b=(5,3) • Elapse time: sample a successor for each particle • Example successor: G2a=(2,3) G2b=(6,3) • Observe: weight each entire sample by the likelihood of the evidence conditioned on the sample • Likelihood: p(E1a|G1a)*p(E1b|G1b) • Resample: select prior samples (tuples of values) in proportion to their likelihood

  19. SLAM • SLAM = Simultaneous Localization And Mapping • We do not know the map or our location • Our belief state is over maps and positions! • Main techniques: • Kalman filtering (Gaussian HMMs) • particle methods • http://www.cs.duke.edu/~parr/dpslam/D-Wing.html

  20. HMMs: MLE Queries • HMMs defined by: • States X • Observations E • Initial distribution: p(X1) • Transitions: p(X|X-1) • Emissions: p(E|X) • Query: most likely explanation:

  21. State Path • Graph of states and transitions over time • Each arc represents some transition xt-1→xt • Each arc has weight p(xt|xt-1)p(et|xt) • Each path is a sequence of states • The product of weights on a path is the seq’s probability • Forward algorithm • computing the sum of all paths • Viterbi algorithm • computing the best paths X1X2…XN

  22. Viterbi Algorithm

  23. Example

More Related