Changeset 265 for tests/UI

Timestamp:

02/09/09 23:14:58 (15 years ago)

Author:

smidl

Message:

UI in matlab

Location:

tests/UI

Files:

• UIArxDS_test.cfg (modified) (2 diffs)
• UIArxDS_test.cpp (modified) (2 diffs)

tests/UI/UIArxDS_test.cfg

@@ -18 +18 @@
         r = [0.1, 0.0,
              0.0, 1.0];
+        opt="L_theta";
 };
 
@@ -29 +30 @@
 //estimation
 estimator = {
-    type = "ARX";
-        
+   type = "ARXest";
+        ychns = [1];
+        rgrid =  [ 1,  1,  1,  2];
+        delays = [-1, -2, -3, -1];
+
+        //optional fields
+        dV0 = [1e-3, 1e-5, 1e-5, 1e-5, 1e-5]; //default: 1e-3 for y, 1e-5 for rgr
+        nu0 = 8.;      //default: rgrlen + 2
+        frg = .990;    // forgetting, default frg=1.0
 };
 

tests/UI/UIArxDS_test.cpp

@@ -1 +1 @@
 /*!
-\file 
+\file
 \brief Test of UI builders for ARX
 
@@ -9 +9 @@
 
 #include <stat/libDS_ui.h>
+#include <estim/arx_ui.h>
 
 using namespace bdm;
 int main() {
-        UIFile F("UIArxDS_test.cfg");
+        UIFile F ( "UIArxDS_test.cfg" );
 
-        logger* L; 
-        UIbuild(F.lookup("logger"),L);
+        logger* L;
+        UIbuild ( F.lookup ( "logger" ),L );
         ArxDS * DS;
-        UIbuild(F.lookup("system"),DS);
+        UIbuild ( F.lookup ( "system" ),DS );
         int Ndat;
-        F.lookupValue("experiment.ndat",Ndat);
+        F.lookupValue ( "experiment.ndat",Ndat );
         // SET SIMULATOR
+        BM* E;
+        UIbuild ( F.lookup ( "estimator" ),E );
 
-        DS->log_add(*L);
+        DS->log_add ( *L );
+        int L_est= L->add ( E->_epdf()._rv(), "est" ); // estimate
+        int L_lb = L->add ( E->_epdf()._rv(), "lb" ); // lower bound
+        int L_ub = L->add ( E->_epdf()._rv(), "ub" ); // upper bound
         L->init();
-        
+
+        vec dt=zeros ( DS->_drv().count() );   //data variable
+        datalink_e2e dl ( E->_drv(),DS->_drv() ); //datalink between a datasource and estimator
+
         for ( int tK=1;tK<Ndat;tK++ ) {
-                DS->step();
-                DS->logit(*L);
+                DS->step();                                                     // simulator step
+                DS->getdata ( dt );                                     // read data
+                E->bayes ( dl.get_val ( dt ) );         // update estimates
+
+                DS->logit ( *L );
+                L->logit ( L_est, E->_epdf().mean() );
+                L->logit ( L_lb,  E->_epdf().mean()-2*sqrt ( E->_epdf().variance() ) );
+                L->logit ( L_ub,  E->_epdf().mean() +2*sqrt ( E->_epdf().variance() ) );
+
                 L->step();
         }
-        
-        L->finalize(); 
-        
+
+        L->finalize();
+
         delete L;
         delete DS;
+        delete E;
         return 0;
 }