Changeset 279 for pmsm/old/pmsm_sim.cpp

Show
Ignore:
Timestamp:
02/24/09 14:12:19 (15 years ago)
Author:
smidl
Message:

Transition of pmsm and libKF

Location:
pmsm/old
Files:
1 added
1 moved

Legend:

Unmodified
Added
Removed
  • pmsm/old/pmsm_sim.cpp

    r271 r279  
    2525public: 
    2626        //! Default constructor 
    27         EKF_unQ ( RV rx, RV ry,RV ru,RV rQ ) :EKFCh ( rx,ry,ru ),BMcond ( rQ ) {}; 
    2827        void condition ( const vec &Q0 ) { 
    2928                Q.setD ( Q0,0 ); 
     
    5352        chmat Q ( Qdiag ); 
    5453        chmat R ( Rdiag ); 
    55         EKFCh KFE ( rx,ry,ru ); 
     54        EKFCh KFE ; 
    5655        KFE.set_parameters ( &fxu,&hxu,Q,R ); 
    5756        KFE.set_est ( mu0, chmat ( 1*ones ( 4 ) ) ); 
    5857 
    5958        RV rQ ( "{Q }","2" ); 
    60         EKF_unQ KFEp ( rx,ry,ru,rQ ); 
     59        EKF_unQ KFEp; 
    6160        KFEp.set_parameters ( &fxu,&hxu,Q,R ); 
    6261        KFEp.set_est ( mu0, chmat ( zeros ( 4 ) ) ); 
    6362 
    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 ); 
    6665        // initialize 
    67         evolQ.set_parameters ( 10.0 ); //sigma = 1/10 mu 
     66        evolQ.set_parameters ( 10.0, "0.01 0.01" ); //sigma = 1/10 mu 
    6867        evolQ.condition ( "0.01 0.01" ); //Zdenek default 
    69         epdf& pfinit=evolQ.posterior(); 
     68        epdf& pfinit=evolQ._epdf(); 
    7069        M.set_est ( pfinit ); 
    71         evolQ.set_parameters ( 10.0 ); 
     70        evolQ.set_parameters ( 10.0,  "0.01 0.01"); 
    7271 
    7372        //