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_unkQpf.cpp

    r278 r279  
    2424public: 
    2525        //! Default constructor 
    26         EKF_unQ ( RV rx, RV ry,RV ru,RV rQ ) :EKFCh ( rx,ry,ru ),BMcond ( rQ ) {}; 
    2726        void condition ( const vec &Q0 ) { 
    2827                Q.setD ( Q0,0 ); 
     
    5352 
    5453        RV rQ ( "{Q}","2" ); 
    55         EKF_unQ KFE ( rx,ry,ru,rQ ); 
     54        EKF_unQ KFE ; 
    5655        KFE.set_parameters ( &fxu,&hxu,Q,R ); 
    5756        KFE.set_est ( mu0, chmat ( 1000*ones ( 4 ) ) ); 
    5857 
    59         mgamma evolQ ( rQ,rQ ); 
     58        mgamma evolQ ; 
    6059        //evolQ.set_parameters ( 10000.0 ); //sigma = 1/10 mu 
    6160 
    62         MPF<EKF_unQ> M ( rx,rQ,evolQ,evolQ,100,KFE ); 
     61        MPF<EKF_unQ> M ( &evolQ,&evolQ,100,KFE ); 
    6362 
    6463        const epdf& KFEep = KFE.posterior(); 
    6564        const epdf& Mep = M.posterior(); 
    6665        // initialize 
    67         evolQ.set_parameters ( 1.0 ); //sigma = 1/10 mu 
     66        evolQ.set_parameters ( 1.0, "0.5 0.5" ); //sigma = 1/10 mu 
    6867        evolQ.condition ( "0.5 0.5" ); 
    69         const epdf& pfinit=evolQ.posterior(); 
     68        const epdf& pfinit=evolQ._epdf(); 
    7069        M.set_est ( pfinit ); 
    71         evolQ.set_parameters ( 1000.0 ); //sigma = 1/10 mu 
     70        evolQ.set_parameters ( 1000.0 , "0.5 0.5"); //sigma = 1/10 mu 
    7271 
    7372        //simulator values