Changeset 283 for doprava/model.h

Timestamp:

02/24/09 14:14:01 (15 years ago)

Author:

smidl

Message:

get rid of BMcond + adaptation in doprava/

Files:

• doprava/model.h (modified) (6 diffs)

doprava/model.h

@@ -5 +5 @@
 
 using namespace bdm;
-
-//Tady se naplni "popis" jednotlivych nahodnych velicin na kterych se pracuje
-// Moznosti je velmi mnoho (viz doc/html/index.html): napriklad
-RV RVstav ( "{stav }", "4");  // Vyrabim stav velikosti 4
-RV RVut ( "{ut }", "2");     // Vstup velikosti 2
-RV RVpozor ( "{I1 O1 }");     //Vystup je intenzita a obsazenost
 
 //! Model stredni hodnoty vyvoje stavu pro k1 
@@ -19 +13 @@
 public:
         //! Constructor
-        IMk1() :diffbifn (RVstav.count(), RVstav, RVut ) {}; 
+        IMk1() :diffbifn () {dimy=4; dimx=3; dimu=1;}; 
         //! set CONSTANT parameters
         void set_parameters ( double alp10,  double alp20) {alp1=alp10; alp2=alp20;}
@@ -25 +19 @@
         vec eval ( const vec &x0, const vec &u0 ) {
                 // napln stav nulami
-                vec xk=zeros ( RVstav.count() );
+                vec xk=zeros ( dimy );
 
                 xk ( 0 ) = 0.2* x0(1) - 0.1* x0(2)+ u0(0);  // vycucane z prstu
@@ -39 +33 @@
 
                 if (full) { // priznak full se nastavi na zacatku => je treba naplnit celou matici
-                        A = eye(RVstav.count());
+                        A = eye(dimy);
                         A(0,1) = 0.2;
                         A(0,2) = -0.1;
@@ -54 +48 @@
 public:
         //! Constructor
-        OMk1() :diffbifn (RVpozor.count(), RVstav, RVut ) {}; 
+        OMk1() :diffbifn ( ) {dimy=2;dimx=3;dimu=2;}; //<======= TODO 
         // Model pozorovani  je trivialni jen se zkopiruji stavy
         vec eval(const vec &x0, const vec &u0 ){
-                vec dt(RVpozor.count());
+                vec dt(dimy);
                 // Pozoruji pouze prvni dva stavy
                 dt(0) = x0(0);
@@ -68 +62 @@
 
                 if (full) { // priznak full se nastavi na zacatku => je treba naplnit celou matici
-                        A = zeros(RVpozor.count(),RVstav.count());
+                        A = zeros(dimy,dimx);
                         A(0,0)=1.0;
                         A(1,1)=1.0;