| 1 | % script for estimation of Mz | 
|---|
| 2 | addpath ~/work/git/mixpp/library/mex/ | 
|---|
| 3 | addpath ~/work/git/mixpp/applications/pmsm/ | 
|---|
| 4 |  | 
|---|
| 5 | dQ=[0.017, 0.017, 0.00000004, 1e-14, 1e-3]; | 
|---|
| 6 | dR=[0.025, 0.025]; | 
|---|
| 7 |  | 
|---|
| 8 | % experiment | 
|---|
| 9 | S.experiment.ndat = 32000; | 
|---|
| 10 |  | 
|---|
| 11 | Orv.type  = 'rv'; | 
|---|
| 12 | Orv.names = { 'X' }; | 
|---|
| 13 | Orv.sizes = [ 4 ]; | 
|---|
| 14 |  | 
|---|
| 15 | Oekf.type = 'rv'; | 
|---|
| 16 |  | 
|---|
| 17 | S.system= struct('type', 'ITppFileDS',... | 
|---|
| 18 |     'filename','tec0013ALL.it',... | 
|---|
| 19 |     'varname','D',... | 
|---|
| 20 |     'rv', Orv); | 
|---|
| 21 |  | 
|---|
| 22 | S.estimator = {struct('type','EKF',... | 
|---|
| 23 |          'IM', struct('type','pmsmIM',... | 
|---|
| 24 |                 'params', struct('type', 'external', 'filename', '../cfg/zcu.cfg','path' ,'pmsm107'),... | 
|---|
| 25 |         'variant', 'Mz'),... | 
|---|
| 26 |          'OM', struct('type','pmsmOM'),... | 
|---|
| 27 |          'dQ',dQ,... | 
|---|
| 28 |          'dR',dR,... | 
|---|
| 29 |          'dP0',[1., 1., 1., 10., 1e-10],... | 
|---|
| 30 |     'drv', Orv,... | 
|---|
| 31 |     'rv', struct('type','rv', 'names', {{'EKFMz'}},'sizes',[5]),... | 
|---|
| 32 |          'log_level','logevidence,logbounds')}; | 
|---|
| 33 |          | 
|---|
| 34 | %S.logger= struct(      'type','mexlog','maxlen',32000); | 
|---|
| 35 | S.logger= struct(       'type','dirfilelog','maxlen',32000,'dirname','t2'); | 
|---|
| 36 |  | 
|---|
| 37 | %pmsm_estim_mex(S) | 
|---|