| 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 | 'options','logll,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) | 
|---|