Changeset 279 for pmsm/old/mpf_u_weight.cpp
- Timestamp:
- 02/24/09 14:12:19 (15 years ago)
- Location:
- pmsm/old
- Files:
-
- 1 added
- 1 moved
Legend:
- Unmodified
- Added
- Removed
-
pmsm/old/mpf_u_weight.cpp
r278 r279 28 28 class EKFCh_cond : public EKFCh , public BMcond { 29 29 public: 30 //! Default constructor 31 EKFCh_cond ( RV rx, RV ry,RV ru,RV rC ) :EKFCh ( rx,ry,ru ),BMcond ( rC ) {}; 30 //! condition on value of pfxu 32 31 void condition ( const vec &val ) { 33 32 pfxu->condition ( val ); … … 92 91 chmat Q ( Qdiag ); 93 92 chmat R ( Rdiag ); 94 EKFCh KFE ( rx,ry,ru );93 EKFCh KFE; 95 94 KFE.set_parameters ( &fxu0,&hxu,Q,R ); 96 95 KFE.set_est ( mu0, chmat ( ones ( 4 ) ) ); 97 96 98 97 RV rW ( "{w }" ); 99 EKFCh_cond KFEp ( rx,ry,ru,rW );98 EKFCh_cond KFEp; 100 99 KFEp.set_parameters ( &fxu,&hxu,Q,R ); 101 100 KFEp.set_est ( mu0, chmat ( ones ( 4 ) ) ); 102 101 103 mgamma_fix evolW ( rW,rW );104 MPF<EKFCh_cond> M ( rx,rW,evolW,evolW,Npart,KFEp );102 mgamma_fix evolW ; 103 MPF<EKFCh_cond> M (&evolW,&evolW,Npart,KFEp ); 105 104 // initialize 106 105 vec W0="0.5"; 107 106 evolW.set_parameters ( 10.0, W0, 1.0 ); //sigma = 1/10 mu 108 107 evolW.condition ( W0 ); 109 epdf& pfinit=evolW. posterior();108 epdf& pfinit=evolW._epdf(); 110 109 M.set_est ( pfinit ); 111 110 evolW.set_parameters ( 100.0, W0, 0.99 ); //sigma = 1/10 mu