Changeset 279 for pmsm/old/pmsm_sim.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/pmsm_sim.cpp
r271 r279 25 25 public: 26 26 //! Default constructor 27 EKF_unQ ( RV rx, RV ry,RV ru,RV rQ ) :EKFCh ( rx,ry,ru ),BMcond ( rQ ) {};28 27 void condition ( const vec &Q0 ) { 29 28 Q.setD ( Q0,0 ); … … 53 52 chmat Q ( Qdiag ); 54 53 chmat R ( Rdiag ); 55 EKFCh KFE ( rx,ry,ru );54 EKFCh KFE ; 56 55 KFE.set_parameters ( &fxu,&hxu,Q,R ); 57 56 KFE.set_est ( mu0, chmat ( 1*ones ( 4 ) ) ); 58 57 59 58 RV rQ ( "{Q }","2" ); 60 EKF_unQ KFEp ( rx,ry,ru,rQ );59 EKF_unQ KFEp; 61 60 KFEp.set_parameters ( &fxu,&hxu,Q,R ); 62 61 KFEp.set_est ( mu0, chmat ( zeros ( 4 ) ) ); 63 62 64 mgamma evolQ ( rQ,rQ );65 MPF<EKF_unQ> M ( rx,rQ,evolQ,evolQ,Npart,KFEp );63 mgamma evolQ ; 64 MPF<EKF_unQ> M (&evolQ,&evolQ,Npart,KFEp ); 66 65 // initialize 67 evolQ.set_parameters ( 10.0 ); //sigma = 1/10 mu66 evolQ.set_parameters ( 10.0, "0.01 0.01" ); //sigma = 1/10 mu 68 67 evolQ.condition ( "0.01 0.01" ); //Zdenek default 69 epdf& pfinit=evolQ. posterior();68 epdf& pfinit=evolQ._epdf(); 70 69 M.set_est ( pfinit ); 71 evolQ.set_parameters ( 10.0 70 evolQ.set_parameters ( 10.0, "0.01 0.01"); 72 71 73 72 //