work/mixpp/libPF.h

Go to the documentation of this file.

#ifndef PF_H
#define PF_H

#include <itpp/itbase.h>
#include "libBM.h"
#include "libDC.h"

using namespace itpp;

enum RESAMPLING_METHOD { MULTINOMIAL = 0, DETERMINISTIC = 1, RESIDUAL = 2, SYSTEMATIC = 3 };

class PF : public BM { 
protected:
        int n; //number of particles
        vec w; //particle weights
        
public:
        ivec resample(RESAMPLING_METHOD method = SYSTEMATIC);
        //TODO get them on the web
};

class TrivialPF : public PF {
        Array<vec> ptcls;
        
        bool is_proposal;
        mpdf *prop;
        mpdf *par;
        mpdf *obs;
        
        public:
        TrivialPF(mpdf &par, mpdf &obs, mpdf &prop, int n0);
        TrivialPF(mpdf &par, mpdf &obs, int n0);
        void bayes(const vec &dt, bool evalll);
};

class MPF : public TrivialPF {
        Array<BM> Bms;
        public:
        MPF(BM &B, mpdf &prop, mpdf &obs, mpdf &par);
        void bayes(vec &dt);    
};

#endif // KF_H

---