Changeset 279 for pmsm/old/pmsm_unkQpf.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_unkQpf.cpp
r278 r279 24 24 public: 25 25 //! Default constructor 26 EKF_unQ ( RV rx, RV ry,RV ru,RV rQ ) :EKFCh ( rx,ry,ru ),BMcond ( rQ ) {};27 26 void condition ( const vec &Q0 ) { 28 27 Q.setD ( Q0,0 ); … … 53 52 54 53 RV rQ ( "{Q}","2" ); 55 EKF_unQ KFE ( rx,ry,ru,rQ );54 EKF_unQ KFE ; 56 55 KFE.set_parameters ( &fxu,&hxu,Q,R ); 57 56 KFE.set_est ( mu0, chmat ( 1000*ones ( 4 ) ) ); 58 57 59 mgamma evolQ ( rQ,rQ );58 mgamma evolQ ; 60 59 //evolQ.set_parameters ( 10000.0 ); //sigma = 1/10 mu 61 60 62 MPF<EKF_unQ> M ( rx,rQ,evolQ,evolQ,100,KFE );61 MPF<EKF_unQ> M ( &evolQ,&evolQ,100,KFE ); 63 62 64 63 const epdf& KFEep = KFE.posterior(); 65 64 const epdf& Mep = M.posterior(); 66 65 // initialize 67 evolQ.set_parameters ( 1.0 ); //sigma = 1/10 mu66 evolQ.set_parameters ( 1.0, "0.5 0.5" ); //sigma = 1/10 mu 68 67 evolQ.condition ( "0.5 0.5" ); 69 const epdf& pfinit=evolQ. posterior();68 const epdf& pfinit=evolQ._epdf(); 70 69 M.set_est ( pfinit ); 71 evolQ.set_parameters ( 1000.0 ); //sigma = 1/10 mu70 evolQ.set_parameters ( 1000.0 , "0.5 0.5"); //sigma = 1/10 mu 72 71 73 72 //simulator values