Changeset 221

Show
Ignore:
Timestamp:
01/05/09 19:53:24 (15 years ago)
Author:
smidl
Message:

matlab directory cleanup

Files:
2 added
4 modified
9 moved

Legend:

Unmodified
Added
Removed

  • matlab/pmsm/mpf_test_disp.m

    r220 r221  
    2121        subplot(4,1,i) 
    2222        plot(xthM(:,i)') 
     23    if i<4 
     24        set(gca,'XTick',[]); 
     25    else 
     26        xlabel('sample [t]'); 
     27    end 
     28    ylabel(['Q(' num2str(i) ',' num2str(i) ')']); 
    2329    hold on 
    2430        plot(Qtr(:,i)','--') 

  • pmsm/mpf_test.cpp

    r220 r221  
    1 /* 
     1/*! 
    22  \file 
    3   \brief Models for synchronous electric drive using IT++ and BDM 
     3  \brief TR 2525 file for testing Toy Problem of mpf for Covariance Estimation 
    44  \author Vaclav Smidl. 
    55 
    … …  
    1010  ----------------------------------- 
    1111*/ 
     12 
    1213 
    1314#include <itpp/itbase.h> 
    … …  
    4950 
    5051        vec mu0= "0.0 0.0 0.0 0.0"; 
    51         vec Qdiag ( "1e-6 1e-6 0.001 0.0001" ); //zdenek: 0.01 0.01 0.0001 0.0001 
     52        vec Qdiag ( "1e-5 1e-5 0.0001 0.0001" ); //zdenek: 0.01 0.01 0.0001 0.0001 
    5253        vec Rdiag ( "1e-8 1e-8" ); //var(diff(xth)) = "0.034 0.034" 
    5354        chmat Q ( Qdiag ); 
    … …  
    6566        MPF<EKF_unQ> M ( rx,rQ,evolQ,evolQ,Npart,KFEp ); 
    6667        // initialize 
    67         evolQ.set_parameters ( 10.0, Qdiag, 0.99 ); //sigma = 1/10 mu 
     68        evolQ.set_parameters ( 10.0, Qdiag, 1.0 ); //sigma = 1/10 mu 
    6869        evolQ.condition (Qdiag ); //Zdenek default 
    6970        epdf& pfinit=evolQ._epdf(); 
    … …  
    110111                if (tK==3000)  Qdiag(1)*=10;  
    111112                if (tK==4000) Qdiag(1)/=10;  
    112                 if (tK==5000)  Qdiag(2)*=10;  
    113                 if (tK==6000) Qdiag(2)/=10;  
    114                 if (tK==7000)  Qdiag(3)*=10;  
    115                 if (tK==8000) Qdiag(3)/=10;  
     113                if (tK==5000)  Qdiag(2)*=100;  
     114                if (tK==6000) Qdiag(2)/=100;  
     115                if (tK==7000)  Qdiag(3)*=100;  
     116                if (tK==8000) Qdiag(3)/=100;  
    116117                 
    117118                //estimator 

  • pmsm/pmsm_unkQpf.cpp

    r170 r221  
    1 /* 
    2   \file 
    3   \brief Models for synchronous electric drive using IT++ and BDM 
     1/*! 
     2  \file  
     3  \brief A test for Kalman with unknown Q 
    44  \author Vaclav Smidl. 
    55 

  • pmsm/sim_var.cpp

    r170 r221  
    11/*! 
    22  \file 
    3   \brief Models for synchronous electric drive using IT++ and BDM 
     3  \brief Simulation of disturbances in PMSM model, PWM and integration disturbances are distinguished 
    44  \author Vaclav Smidl. 
    55 

  • pmsm/sim_var_arx.cpp

    r170 r221  
    11/*! 
    22  \file 
    3   \brief Models for synchronous electric drive using IT++ and BDM 
     3  \brief Disturbances of PMSM model are fitted by an ARX model 
    44  \author Vaclav Smidl. 
    55