Changeset 357 for applications/pmsm/TR2245/pmsm_wishart.cpp
- Timestamp:
- 06/08/09 02:15:30 (15 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
applications/pmsm/TR2245/pmsm_wishart.cpp
r332 r357 17 17 #include <stat/libFN.h> 18 18 19 #include <stat/loggers_ui.h> 20 #include <stat/libEF_ui.h> 21 22 #include "../pmsm_ui.h" 19 #include "user_info.h" 20 #include "../pmsmDS.h" 23 21 24 22 using namespace bdm; … … 28 26 if ( argc>1 ) {fname = argv[1]; } 29 27 else { fname = "pmsm_wishart.cfg"; } 30 UIFile F ( fname ); 28 UI_File F ( fname ); 29 30 double h = 1e-6; 31 int Nsimstep = 125; 32 33 // Kalman filter 34 double k; 35 double l; 36 F.lookupValue ( "k",k); 37 F.lookupValue ( "l",l); 31 38 32 39 int Ndat; 33 40 int Npart; 34 double h = 1e-6;35 int Nsimstep = 125;41 F.lookupValue ( "ndat", Ndat ); 42 F.lookupValue ( "Npart", Npart ); 36 43 37 44 vec Qdiag; 38 45 vec Rdiag; 46 UI::get( Qdiag, F, "dQ" ); //( "1e-6 1e-6 0.001 0.0001" ); //zdenek: 0.01 0.01 0.0001 0.0001 47 UI::get( Rdiag, F, "dR" );// ( "1e-8 1e-8" ); //var(diff(xth)) = "0.034 0.034" 48 49 pmsmDS* DS = UI::build<pmsmDS>(F,"system"); 50 // mpdf* evolQ = UI::build<mpdf>( F, "Qrw" ); 39 51 40 pmsmDS* DS; 41 42 double k; 43 double l; 44 // mpdf* evolQ ; 45 try { 46 // Kalman filter 47 F.lookupValue ( "ndat", Ndat ); 48 F.lookupValue ( "Npart",Npart ); 49 50 F.lookupValue ( "k", k); 51 F.lookupValue ( "l",l); 52 53 // UIbuild ( F.lookup ( "Qrw" ),evolQ ); 54 Qdiag= getvec ( F.lookup ( "dQ" ) ); //( "1e-6 1e-6 0.001 0.0001" ); //zdenek: 0.01 0.01 0.0001 0.0001 55 Rdiag=getvec ( F.lookup ( "dR" ) );// ( "1e-8 1e-8" ); //var(diff(xth)) = "0.034 0.034" 56 57 UIbuild(F.lookup("system"),DS); 58 } 59 catch UICATCH; 60 // internal model 61 62 IMpmsm fxu; 63 // Rs Ls dt Fmag(Ypm) kp p J Bf(Mz) 64 fxu.set_parameters ( 0.28, 0.003465, Nsimstep*h, 0.1989, 1.5 ,4.0, 0.04, 0.0 ); 52 // internal model 53 IMpmsm fxu; 54 // Rs Ls dt Fmag(Ypm) kp p J Bf(Mz) 55 fxu.set_parameters ( 0.28, 0.003465, Nsimstep*h, 0.1989, 1.5 ,4.0, 0.04, 0.0 ); 65 56 // observation model 66 57 OMpmsm hxu; … … 91 82 M.set_rv ( concat ( rQ,rx ) ); 92 83 93 dirfilelog *L ; UIbuild ( F.lookup ( "logger" ), L);// ( "exp/mpf_test",100 );84 dirfilelog *L = UI::build<dirfilelog> ( F, "logger" );// ( "exp/mpf_test",100 ); 94 85 95 86 KFE.set_options ( "logbounds" );