Changeset 357 for applications/pmsm/TR2245/wishart.cpp
- Timestamp:
- 06/08/09 02:15:30 (15 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
applications/pmsm/TR2245/wishart.cpp
r317 r357 17 17 #include <stat/libFN.h> 18 18 19 #include <stat/loggers_ui.h>20 #include <stat/libEF_ui.h>21 22 19 #include "../pmsm.h" 23 20 #include "simulator.h" 24 21 #include "../sim_profiles.h" 22 #include "user_info.h" 23 #include "stat/loggers.h" 25 24 26 25 using namespace bdm; … … 30 29 if ( argc>1 ) {fname = argv[1]; } 31 30 else { fname = "unitsteps.cfg"; } 32 UI File F ( fname );31 UI_File F ( fname ); 33 32 33 double h = 1e-6; 34 34 int Ndat; 35 35 int Npart; 36 double h = 1e-6; 36 F.lookupValue ( "ndat", Ndat ); 37 F.lookupValue ( "Npart",Npart ); 37 38 int Nsimstep = 125; 38 39 40 // Kalman filter 39 41 vec Qdiag; 42 UI::get( Qdiag, F, "dQ" ); //( "1e-6 1e-6 0.001 0.0001" ); //zdenek: 0.01 0.01 0.0001 0.0001 43 40 44 vec Rdiag; 45 UI::get( Rdiag, F, "dR" );// ( "1e-8 1e-8" ); //var(diff(xth)) = "0.034 0.034" 41 46 42 // mpdf* evolQ ; 43 try { 44 // Kalman filter 45 F.lookupValue ( "ndat", Ndat ); 46 F.lookupValue ( "Npart",Npart ); 47 48 // UIbuild ( F.lookup ( "Qrw" ),evolQ ); 49 Qdiag= getvec ( F.lookup ( "dQ" ) ); //( "1e-6 1e-6 0.001 0.0001" ); //zdenek: 0.01 0.01 0.0001 0.0001 50 Rdiag=getvec ( F.lookup ( "dR" ) );// ( "1e-8 1e-8" ); //var(diff(xth)) = "0.034 0.034" 51 } 52 catch UICATCH; 53 // internal model 54 55 IMpmsm fxu; 56 // Rs Ls dt Fmag(Ypm) kp p J Bf(Mz) 57 fxu.set_parameters ( 0.28, 0.003465, Nsimstep*h, 0.1989, 1.5 ,4.0, 0.04, 0.0 ); 47 // internal model 48 IMpmsm fxu; 49 // Rs Ls dt Fmag(Ypm) kp p J Bf(Mz) 50 fxu.set_parameters ( 0.28, 0.003465, Nsimstep*h, 0.1989, 1.5 ,4.0, 0.04, 0.0 ); 58 51 // observation model 59 52 OMpmsm hxu; … … 84 77 M.set_rv ( concat ( rQ,rx ) ); 85 78 86 dirfilelog *L ; UIbuild ( F.lookup ( "logger" ), L);// ( "exp/mpf_test",100 );79 dirfilelog *L = UI::build<dirfilelog> ( F, "logger" );// ( "exp/mpf_test",100 ); 87 80 int l_X = L->add ( rx, "xt" ); 88 81 int l_D = L->add ( concat ( ry,ru ), "" ); … … 103 96 vec xtm=zeros ( 4 ); 104 97 double Ww=0.0; 105 vec vecW=getvec ( F.lookup ( "profile" ) ); 98 vec vecW; 99 UI::get( vecW, F ,"profile" ); 106 100 107 101 mat tQ=diag(Qdiag);