Changeset 28 for tests

Timestamp:

02/22/08 16:40:12 (16 years ago)

Author:

smidl

Message:

prelozitelna verze

Location:

tests

Files:

• test0.cpp (modified) (2 diffs)
• testKF.cpp (modified) (2 diffs)
• testSmp.cpp (modified) (1 diff)

tests/test0.cpp

@@ -14 +14 @@
         RV th = RV ( "1 2","{a b }","1 1","0 0","0 0" );
         RV r = RV ( "3 4" );
-        int i =0;
 
         cout << th << r << endl;
@@ -42 +41 @@
 
         cout << "Test ldform()"<<endl;
-        mat V = "1 2; 2 3"; 
+        mat V = "1 2; 2 13"; 
         ldmat lV(V);
         ldmat ilV(V);  

tests/testKF.cpp

@@ -22 +22 @@
 
         bool xxx= fin.seek( "d" );
+        if (!xxx){ it_error("testKF.it not found");}
         fin >>Dt;
         fin.seek( "A" ); 
@@ -43 +44 @@
         XtE=zeros( 2,Ndat );
 
-        Kalman<fsqmat> KF( A,B,C,D,fsqmat(R),fsqmat(Q),fsqmat(P0),mu0 );
 //      cout << KF;
-        KalmanFull KF2( A,B,C,D,R,Q,P0,mu0 );
         RV rx("1","{x}","2","0","0");
         RV ru("2","{u}","1","0","0");
         RV ry("3","{y}","1","0","0");
+        //
+        Kalman<ldmat> KF(rx,ry,ru);
+        KF.set_parameters(A,B,C,D,ldmat(R),ldmat(Q));
+        KF.set_est(mu0,ldmat(P0) );
+        //
+        Kalman<fsqmat> KFf(rx,ry,ru);
+        KFf.set_parameters(A,B,C,D,fsqmat(R),fsqmat(Q));
+        KFf.set_est(mu0,fsqmat(P0) );
+        //
+        KalmanFull KF2( A,B,C,D,R,Q,P0,mu0 );
+        //
         bilinfn fxu(rx,ru,A,B);
         bilinfn hxu(rx,ru,C,D);
-        EKF<fsqmat> KFE(&fxu,&hxu,Q,R,mu0,P0);
+        EKF<ldmat> KFE(rx,ry,ru);
+        KFE.set_parameters(&fxu,&hxu,Q,R);
+        KFE.set_est(mu0,P0);
 
-        Xt.set_col( 0,KF.mu );
-        Xt2.set_col( 0,KF.mu );
-        XtE.set_col( 0,KF.mu );
+        Xt.set_col( 0,*((enorm<ldmat>*)(KF._epdf()))->_mu() );
+        Xt2.set_col( 0,KF2.mu );
+        XtE.set_col( 0,*((enorm<ldmat>*)(KFE._epdf()))->_mu() );
         for ( int t=1;t<Ndat;t++ ) {
+                KFf.bayes( Dt.get_col( t ));
                 KF.bayes( Dt.get_col( t ));
                 KF2.bayes( Dt.get_col( t ));
-                KFE.bayes( Dt.get_col( t ));
-                Xt.set_col(t,KF.mu);
+//              KFE.bayes( Dt.get_col( t ));
+                Xt.set_col(t,*((enorm<ldmat>*)(KF._epdf()))->_mu());
                 Xt2.set_col(t,KF2.mu);
-                XtE.set_col(t,KFE.mu);
-//              Kmu = KF.mu;
-//              cout <<  "t:" <<t<< "  " << dt<<"  "<<Kmu <<endl;
+                XtE.set_col(t,*((enorm<ldmat>*)(KFE._epdf()))->_mu());
         }
 

tests/testSmp.cpp

@@ -17 +17 @@
         mat V0("1.2 0.3; 0.3 5");
         ldmat R = ldmat(V0);
-        enorm<ldmat> eN(rv,mu0,R);
+        enorm<ldmat> eN(rv);
+        eN.set_parameters(mu0,R);
         
         mat Smp = eN.sample(N);