Bayesian Filtering Library  Generated from SVN r
Public Member Functions | Protected Member Functions | Protected Attributes
MixtureBootstrapFilter< StateVar, MeasVar > Class Template Reference

Particular mixture particle filter : Proposal PDF = SystemPDF. More...

#include <mixtureBootstrapFilter.h>

Inheritance diagram for MixtureBootstrapFilter< StateVar, MeasVar >:
MixtureParticleFilter< StateVar, MeasVar > Filter< StateVar, MeasVar >

List of all members.

Public Member Functions

 MixtureBootstrapFilter (Mixture< StateVar > *prior, int resampleperiod=0, double resamplethreshold=0, int resamplescheme=DEFAULT_RS, int maintainMixturePeriod=1)
 Constructor.
 MixtureBootstrapFilter (Mixture< StateVar > *prior, Mixture< StateVar > *post, int resampleperiod=0, double resamplethreshold=0, int resamplescheme=DEFAULT_RS, int maintainMixturePeriod=1)
 Constructor.
virtual ~MixtureBootstrapFilter ()
 Destructor.
virtual void Reset (Mixture< StateVar > *prior)
 Reset Filter.
virtual void Reset (Pdf< StateVar > *prior)
 Reset Filter.
virtual void ProposalSet (ConditionalPdf< StateVar, StateVar > *const cpdf)
 Set the proposal density.
ConditionalPdf< StateVar,
StateVar > * 
ProposalGet ()
 Get a pointer to the proposal density.
virtual Mixture< StateVar > * PostGet ()
 Get Posterior density.
virtual bool Update (SystemModel< StateVar > *const sysmodel, const StateVar &u, MeasurementModel< MeasVar, StateVar > *const measmodel, const MeasVar &z, const StateVar &s)
 Full Update (system with inputs/sensing params)
virtual bool Update (SystemModel< StateVar > *const sysmodel, MeasurementModel< MeasVar, StateVar > *const measmodel, const MeasVar &z, const StateVar &s)
 Full Update (system without inputs, with sensing params)
virtual bool Update (SystemModel< StateVar > *const sysmodel, MeasurementModel< MeasVar, StateVar > *const measmodel, const MeasVar &z)
 Full Update (system without inputs/sensing params)
virtual bool Update (SystemModel< StateVar > *const sysmodel, const StateVar &u, MeasurementModel< MeasVar, StateVar > *const measmodel, const MeasVar &z)
 Full Update (system with inputs, without sensing params)
virtual bool Update (SystemModel< StateVar > *const sysmodel, const StateVar &u)
 System Update (system with inputs)
virtual bool Update (SystemModel< StateVar > *const sysmodel)
 System Update (system without inputs)
virtual bool Update (MeasurementModel< MeasVar, StateVar > *const measmodel, const MeasVar &z, const StateVar &s)
 Measurement Update (system with "sensing params")
virtual bool Update (MeasurementModel< MeasVar, StateVar > *const measmodel, const MeasVar &z)
 Measurement Update (system without "sensing params")
int TimeStepGet () const
 Get current time.

Protected Member Functions

virtual bool UpdateInternal (SystemModel< StateVar > *const sysmodel, const StateVar &u, MeasurementModel< MeasVar, StateVar > *const measmodel, const MeasVar &z, const StateVar &s)
 Actual implementation of updateinternal.
virtual bool ProposalStepInternal (SystemModel< StateVar > *const sysmodel, const StateVar &u, MeasurementModel< MeasVar, StateVar > *const measmodel, const MeasVar &z, const StateVar &s)
 Proposal step.
virtual bool ProposalStepInternalOne (int component, SystemModel< StateVar > *const sysmodel, const StateVar &u, MeasurementModel< MeasVar, StateVar > *const measmodel, const MeasVar &z, const StateVar &s)
 Proposal step for one component.
virtual bool UpdateWeightsInternal (SystemModel< StateVar > *const sysmodel, const StateVar &u, MeasurementModel< MeasVar, StateVar > *const measmodel, const MeasVar &z, const StateVar &s)
 Update Weights.
virtual bool UpdateWeightsInternalOne (int component, SystemModel< StateVar > *const sysmodel, const StateVar &u, MeasurementModel< MeasVar, StateVar > *const measmodel, const MeasVar &z, const StateVar &s)
 Update Weights for one component.
virtual bool DynamicResampleStep ()
 Resample if necessary.
virtual bool DynamicResampleStepOne (int component)
 Resampling for one component.
virtual bool StaticResampleStep ()
 Resample if wanted.
virtual bool Resample ()
 Actual Resampling happens here;.
virtual bool ResampleOne (int component)
 Actual Resampling for one component;.
virtual bool MaintainMixtureStep ()
 Maintain Mixture if wanted.
virtual bool MaintainMixture ()
 Actual mixture maintainance happens here;.

Protected Attributes

ConditionalPdf< StateVar,
StateVar > * 
_proposal
 Pointer to the Proposal Density.
WeightedSample< StateVar > _sample
 While updating use sample<StateVar>
vector< vector< WeightedSample
< StateVar > > > 
_old_samplesVec
 While updating store list of old samples.
vector< vector< WeightedSample
< StateVar > > > 
_new_samplesVec
 While updating store list of new samples.
vector< vector< Sample
< StateVar > > > 
_new_samples_unweightedVec
 While resampling.
vector< WeightedSample
< StateVar > >::iterator 
_os_it
 Iterator for old list of samples.
vector< WeightedSample
< StateVar > >::iterator 
_ns_it
 Iterator for new list of samples.
vector< Probability_newMixtureWeights
 Vector containing the new mixture weights during update step.
vector< Probability_sumWeights
 Vector containing the sum of weights during update step.
int _resamplePeriod
 Number of timestep between resampling from the Posterior Pdf.
double _resampleThreshold
 Threshold used when dynamic resampling.
int _resampleScheme
 Which resample algorithm (see top of particle.h for #defines)
bool _dynamicResampling
 Dynamic resampling or fixed period resampling?
bool _proposal_depends_on_meas
 Proposal depends on last measurement?
bool _created_post
 created own post
int _maintainMixturePeriod
 Number of timestep between mixture maintainance of the Posterior.
Pdf< StateVar > * _prior
 prior Pdf
Pdf< StateVar > * _post
 Pointer to the Posterior Pdf.
int _timestep
 Represents the current timestep of the filter.

Detailed Description

template<typename StateVar, typename MeasVar>
class BFL::MixtureBootstrapFilter< StateVar, MeasVar >

Particular mixture particle filter : Proposal PDF = SystemPDF.

This is one (simple) particular implementation of a particle filter, in which the proposal density is equal to the pdf describing the system model (aka as SystemPdf), and involving a resampling step

The reason why I chose the name bootstrap filter is the fact that this is the name used in the book by Doucet et al.

    @Book{		  doucet_book,
    editor	= {Doucet, Arnaud and de Freytas, Nando and Gordon, Neil},
    title		= {{S}equential {M}onte {C}arlo {M}ethods in {P}ractice},
    publisher	= {Springer--Verlag},
    year		= {2001},
    series	= {Statistics for engineering and information science},
    month		= {january},
    annote	= {see http://www-sigproc.eng.cam.ac.uk/~ad2/book.html}
    }
    

(and I presume this will become a/the standard book about particle filtering). Typical for the bootstrap filter is the fact that the proposal density is chosen to be the SystemPdf of the SystemModel. So there is no proposal density in the constructor here

Todo:
The implementation is very slow for the moment. It would probably be much faster to add a vector<WeightedSample> to the private members of this class.
See also:
Pdf

Definition at line 71 of file mixtureBootstrapFilter.h.


Constructor & Destructor Documentation

MixtureBootstrapFilter ( Mixture< StateVar > *  prior,
int  resampleperiod = 0,
double  resamplethreshold = 0,
int  resamplescheme = DEFAULT_RS,
int  maintainMixturePeriod = 1 
)

Constructor.

Precondition:
you created the necessary models and the prior
Parameters:
priorpointer to the Monte Carlo Pdf prior density
resampleperiodfixed resampling period (if desired)
resamplethresholdthreshold used when dynamic resampling
resampleschemeresampling scheme, see header file for different defines and their meaning
maintainMixturePeriodfixed mixture maintainance period
MixtureBootstrapFilter ( Mixture< StateVar > *  prior,
Mixture< StateVar > *  post,
int  resampleperiod = 0,
double  resamplethreshold = 0,
int  resamplescheme = DEFAULT_RS,
int  maintainMixturePeriod = 1 
)

Constructor.

Precondition:
you created the necessary models and the prior
Parameters:
priorpointer to the Monte Carlo Pdf prior density
postpointer to the Monte Carlo Pdf post density
resampleperiodfixed resampling period (if desired)
resamplethresholdthreshold used when dynamic resampling
resampleschemeresampling scheme, see header file for different defines and their meaning
maintainMixturePeriodfixed mixture maintainance period

Member Function Documentation

bool DynamicResampleStepOne ( int  component) [protected, virtual, inherited]

Resampling for one component.

Parameters:
componentcomponent number of component to update
Bug:
let the user implement her/his own resamplescheme

Definition at line 368 of file mixtureParticleFilter.h.

bool MaintainMixtureStep ( ) [protected, virtual, inherited]

Maintain Mixture if wanted.

Bug:
let the user implement her/his own mixture maintaince scheme

Definition at line 508 of file mixtureParticleFilter.h.

Mixture< SV > * PostGet ( ) [virtual, inherited]

Get Posterior density.

Get the current Posterior density

Returns:
a pointer to the current posterior

Reimplemented from Filter< StateVar, MeasVar >.

Definition at line 502 of file mixtureParticleFilter.h.

ConditionalPdf< SV, SV > * ProposalGet ( ) [inherited]

Get a pointer to the proposal density.

Returns:
a pointer to the proposal density

Definition at line 170 of file mixtureParticleFilter.h.

void ProposalSet ( ConditionalPdf< StateVar, StateVar > *const  cpdf) [virtual, inherited]

Set the proposal density.

Parameters:
cpdfthe new proposal density. The order of the conditional arguments is fixed and should be: x (state), u (input), z (measurement), s (sensor param). Off course all of them are optional

Definition at line 164 of file mixtureParticleFilter.h.

bool ProposalStepInternal ( SystemModel< StateVar > *const  sysmodel,
const StateVar &  u,
MeasurementModel< MeasVar, StateVar > *const  measmodel,
const MeasVar &  z,
const StateVar &  s 
) [protected, virtual, inherited]

Proposal step.

Implementation of proposal step

Parameters:
sysmodelpointer to the used system model
uinput param for proposal density
measmodelpointer to the used measurementmodel
zmeasurement param for proposal density
ssensor param for proposal density
Bug:
Make sampling method variable. See implementation.

Definition at line 177 of file mixtureParticleFilter.h.

bool ProposalStepInternalOne ( int  component,
SystemModel< StateVar > *const  sysmodel,
const StateVar &  u,
MeasurementModel< MeasVar, StateVar > *const  measmodel,
const MeasVar &  z,
const StateVar &  s 
) [protected, virtual, inherited]

Proposal step for one component.

Implementation of proposal step

Parameters:
componentcomponent number of component to update
sysmodelpointer to the used system model
uinput param for proposal density
measmodelpointer to the used measurementmodel
zmeasurement param for proposal density
ssensor param for proposal density
Bug:
Make sampling method variable. See implementation.

Definition at line 193 of file mixtureParticleFilter.h.

bool ResampleOne ( int  component) [protected, virtual, inherited]

Actual Resampling for one component;.

Parameters:
componentcomponent number of component to update

Definition at line 457 of file mixtureParticleFilter.h.

bool StaticResampleStep ( ) [protected, virtual, inherited]

Resample if wanted.

Bug:
let the user implement her/his own resamplescheme

Definition at line 404 of file mixtureParticleFilter.h.

int TimeStepGet ( ) const [inherited]

Get current time.

Get the current time of the filter

Returns:
the current timestep

Definition at line 50 of file filter.h.

bool Update ( SystemModel< SVar > *const  sysmodel,
const SVar &  u,
MeasurementModel< MVar, SVar > *const  measmodel,
const MVar &  z,
const SVar &  s 
) [virtual, inherited]

Full Update (system with inputs/sensing params)

Parameters:
sysmodelpointer to the system model to use for update
uinput to the system
measmodelpointer to the measurement model to use for update
zmeasurement
s"sensing parameter"

Definition at line 56 of file filter.h.

bool Update ( SystemModel< SVar > *const  sysmodel,
MeasurementModel< MVar, SVar > *const  measmodel,
const MVar &  z,
const SVar &  s 
) [virtual, inherited]

Full Update (system without inputs, with sensing params)

Parameters:
sysmodelpointer to the system model to use for update
measmodelpointer to the measurement model to use for update
zmeasurement
s"sensing parameter"

Definition at line 66 of file filter.h.

bool Update ( SystemModel< SVar > *const  sysmodel,
MeasurementModel< MVar, SVar > *const  measmodel,
const MVar &  z 
) [virtual, inherited]

Full Update (system without inputs/sensing params)

Parameters:
sysmodelpointer to the system model to use for update
measmodelpointer to the measurement model to use for update
zmeasurement

Definition at line 76 of file filter.h.

bool Update ( SystemModel< SVar > *const  sysmodel,
const SVar &  u,
MeasurementModel< MVar, SVar > *const  measmodel,
const MVar &  z 
) [virtual, inherited]

Full Update (system with inputs, without sensing params)

Parameters:
sysmodelpointer to the system model to use for update
uinput to the system
measmodelpointer to the measurement model to use for update
zmeasurement

Definition at line 85 of file filter.h.

bool Update ( SystemModel< SVar > *const  sysmodel,
const SVar &  u 
) [virtual, inherited]

System Update (system with inputs)

Parameters:
sysmodelpointer to the system model to use for update
uinput to the system

Definition at line 95 of file filter.h.

bool Update ( SystemModel< SVar > *const  sysmodel) [virtual, inherited]

System Update (system without inputs)

Parameters:
sysmodelpointer to the system model to use for update

Definition at line 103 of file filter.h.

bool Update ( MeasurementModel< MVar, SVar > *const  measmodel,
const MVar &  z,
const SVar &  s 
) [virtual, inherited]

Measurement Update (system with "sensing params")

Parameters:
measmodelpointer to the measurement model to use for update
zmeasurement
s"sensing parameter"

Definition at line 110 of file filter.h.

bool Update ( MeasurementModel< MVar, SVar > *const  measmodel,
const MVar &  z 
) [virtual, inherited]

Measurement Update (system without "sensing params")

Parameters:
measmodelpointer to the measurement model to use for update
zmeasurement

Definition at line 119 of file filter.h.

bool UpdateWeightsInternal ( SystemModel< StateVar > *const  sysmodel,
const StateVar &  u,
MeasurementModel< MeasVar, StateVar > *const  measmodel,
const MeasVar &  z,
const StateVar &  s 
) [protected, virtual, inherited]

Update Weights.

Parameters:
sysmodelpointer to the used system model
uinput param for proposal density
measmodelpointer to the used measurementmodel
zmeasurement param for proposal density
ssensor param for proposal density

Definition at line 248 of file mixtureParticleFilter.h.

bool UpdateWeightsInternalOne ( int  component,
SystemModel< StateVar > *const  sysmodel,
const StateVar &  u,
MeasurementModel< MeasVar, StateVar > *const  measmodel,
const MeasVar &  z,
const StateVar &  s 
) [protected, virtual, inherited]

Update Weights for one component.

Parameters:
componentcomponent number of component to update
sysmodelpointer to the used system model
uinput param for proposal density
measmodelpointer to the used measurementmodel
zmeasurement param for proposal density
ssensor param for proposal density

Definition at line 275 of file mixtureParticleFilter.h.


Member Data Documentation

int _maintainMixturePeriod [protected, inherited]

Number of timestep between mixture maintainance of the Posterior.

By choosing this period, one can execute period maintainance of mixtures. You can implement maintainMixture such that target splitting, appearing and disappearing is handled

Definition at line 213 of file mixtureParticleFilter.h.

Pdf<StateVar>* _post [protected, inherited]

Pointer to the Posterior Pdf.

The Posterior Pdf represents the subjective belief of the person applying the filter AFTER processing inputs and measurements. A filter does not maintain the beliefs at all timesteps t, since this leads to non-constant (or ever growing if you prefer) memory requirements. However, it is possible, to copy the Posterior density at all timesteps in your application by means of the PostGet() member function

See also:
PostGet()

Definition at line 95 of file filter.h.

ConditionalPdf<StateVar,StateVar>* _proposal [protected, inherited]

Pointer to the Proposal Density.

Every particle filter (or more correct: every Sequential Importance Sampling method) uses a proposal density to do the forward sampling step

Definition at line 168 of file mixtureParticleFilter.h.

int _resamplePeriod [protected, inherited]

Number of timestep between resampling from the Posterior Pdf.

By choosing this period, one can avoid numerical instability (aka Degeneration of the particle filter

Definition at line 191 of file mixtureParticleFilter.h.

int _timestep [protected, inherited]

Represents the current timestep of the filter.

Todo:
Check wether this really belongs here

Definition at line 100 of file filter.h.


The documentation for this class was generated from the following files: