Show
Ignore:
Timestamp:
08/27/10 17:05:55 (14 years ago)
Author:
smidl
Message:

pmsm stuff fro bierman

Files:
1 modified

Legend:

Unmodified
Added
Removed
  • applications/pmsm/simulator_zdenek/ekf_example/ekf_obj.h

    r384 r1168  
    6565        //! Default constructor 
    6666        EKFfixed ():BM(),E(),Ry(2,2){ 
    67         int i; 
    68  for(i=0;i<16;i++){Q[i]=0;} 
    69  for(i=0;i<4;i++){R[i]=0;} 
     67                int i; 
     68                for(i=0;i<16;i++){Q[i]=0;} 
     69                for(i=0;i<4;i++){R[i]=0;} 
    7070 
    71  for(i=0;i<4;i++){x_est[i]=0;} 
    72  for(i=0;i<4;i++){x_pred[i]=0;} 
    73  for(i=0;i<16;i++){P_pred[i]=0;} 
    74  for(i=0;i<16;i++){P_est[i]=0;} 
    75  P_est[0]=0x7FFF; 
    76  P_est[5]=0x7FFF; 
    77  P_est[10]=0x7FFF; 
    78  P_est[15]=0x7FFF; 
    79  for(i=0;i<2;i++){Y_mes[i]=0;} 
    80  for(i=0;i<2;i++){ukalm[i]=0;} 
    81  for(i=0;i<8;i++){Kalm[i]=0;} 
     71                for(i=0;i<4;i++){x_est[i]=0;} 
     72                for(i=0;i<4;i++){x_pred[i]=0;} 
     73                for(i=0;i<16;i++){P_pred[i]=0;} 
     74                for(i=0;i<16;i++){P_est[i]=0;} 
     75                P_est[0]=0x7FFF; 
     76                P_est[5]=0x7FFF; 
     77                P_est[10]=0x7FFF; 
     78                P_est[15]=0x7FFF; 
     79                for(i=0;i<2;i++){Y_mes[i]=0;} 
     80                for(i=0;i<2;i++){ukalm[i]=0;} 
     81                for(i=0;i<8;i++){Kalm[i]=0;} 
    8282 
    83  for(i=0;i<16;i++){PSI[i]=0;} 
    84 }; 
     83                for(i=0;i<16;i++){PSI[i]=0;} 
     84 
     85                set_dim(4); 
     86                E._mu()=zeros(4); 
     87                E._R()=zeros(4,4); 
     88                init_ekf(0.000125); 
     89        }; 
    8590        //! Here dt = [yt;ut] of appropriate dimensions 
    86         void bayes ( const vec &dt ); 
     91        void bayes ( const vec &yt, const vec &ut ); 
    8792        //!dummy! 
    88         epdf& posterior(){return E;}; 
    89         void condition ( const vec &Q0 ) { 
    90                  
    91                 Q[0]=prevod(minQ(Q0(0)),15);       // 0.05 
    92                 Q[5]=prevod(minQ(Q0(1)),15); 
    93                 Q[10]=prevod(minQ(Q0(2)),15);      // 1e-3 
    94                 Q[15]=prevod(minQ(Q0(3)),15);      // 1e-3 
    95  
    96         } 
     93        const epdf& posterior() const {return E;}; 
     94         
    9795}; 
    9896 
     97UIREGISTER(EKFfixed); 
    9998 
    10099#endif // KF_H