Show
Ignore:
Timestamp:
11/04/09 22:55:17 (15 years ago)
Author:
smidl
Message:

New transformation between StateSpace? and ARX

Files:
1 modified

Legend:

Unmodified
Added
Removed
  • library/tests/LQG_test.cpp

    r699 r703  
    3939} 
    4040 
     41TEST(to_state_arx_test) { 
     42        mlnorm<fsqmat> ml; 
     43        mat A="1.1, 2.3"; 
     44        ml.set_parameters(A, vec_1(1.3), eye(1)); 
     45        RV yr=RV("y",1); 
     46        RV ur=RV("u",1); 
     47        ml.set_rv(yr); 
     48        yr.t_plus(-1); 
     49        ml.set_rvc(concat(yr, ur)); 
     50         
     51        shared_ptr<StateFromARX> Stsp=new StateFromARX; 
     52        Stsp->connect_mlnorm(ml); 
     53         
     54        /* results from  
     55        [A,B,C,D]=tf2ss([2.3 0],[1 -1.1]) 
     56        */ 
     57        cout << "---" << endl; 
     58} 
     59 
    4160TEST(arx_LQG_test){ 
    4261        mlnorm<fsqmat> ml; 
    4362        mat A="1.81, -.8189, .00468, .00438"; 
    44         ml.set_parameters(A, vec_1(0), eye(0.0001)); 
     63        ml.set_parameters(A, vec_1(0.0), 0.00001*eye(1)); 
    4564        RV yr=RV("y",1); 
    4665        RV ur=RV("u",1); 
     
    4867        rgr.add(yr.copy_t(-2)); 
    4968        rgr.add(yr.copy_t(-2)); 
     69        rgr.add(ur.copy_t(-2)); 
    5070        rgr.add(ur.copy_t(-1)); 
    51         rgr.add(ur); 
    5271         
    5372        ml.set_rv(yr); 
     
    6483         
    6584        L.redesign(); 
     85        cout << L.to_string()<<endl; 
    6686}