| 1 | %load graf1.txt -ascii |
|---|
| 2 | load graf2.txt -ascii |
|---|
| 3 | |
|---|
| 4 | % t, Isx, Isy, Isd, Isdw, Isq, Isqw, M, Fs, poloha, rychlost, rychlost_w, |
|---|
| 5 | % Urm, Urm_max, beta, zbeta, Uc |
|---|
| 6 | |
|---|
| 7 | ind = 1:size(graf1,1); |
|---|
| 8 | |
|---|
| 9 | % t, isx,isy, usx usy fme rychlost |
|---|
| 10 | figure |
|---|
| 11 | plot(graf1(ind,1),graf1(ind,2)) % Isx |
|---|
| 12 | hold on |
|---|
| 13 | plot(graf1(ind,1),graf2(ind,3),'r') % compare with g2 |
|---|
| 14 | legend('Isx','Isy'); |
|---|
| 15 | |
|---|
| 16 | figure |
|---|
| 17 | plot(graf1(ind,1),graf2(ind,4)) % Isx |
|---|
| 18 | hold on |
|---|
| 19 | plot(graf1(ind,1),graf2(ind,5),'r') % compare with g2 |
|---|
| 20 | legend('Usx','Usy'); |
|---|
| 21 | |
|---|
| 22 | figure |
|---|
| 23 | plot(graf1(ind,1),graf1(ind,25)/2/pi) % fme_ekf |
|---|
| 24 | hold on |
|---|
| 25 | plot(graf1(ind,1),graf1(ind,11)/(2*pi),'r') % fme skutecna |
|---|
| 26 | plot(graf1(ind,1),graf2(ind,6)/(2*pi),'g') % fme MPF |
|---|
| 27 | legend('fme ekf','fme', 'fme mpf'); |
|---|
| 28 | |
|---|
| 29 | figure |
|---|
| 30 | plot(graf1(ind,1),graf1(ind,26)) % theta_elektricka_ekf |
|---|
| 31 | hold on |
|---|
| 32 | plot(graf1(ind,1),graf1(ind,10),'r') % Poloha elektricka skutecna |
|---|
| 33 | plot(graf1(ind,1),graf2(ind,7),'g') % theta_elektricka_mpf |
|---|
| 34 | legend('Theta ekf','Theta real', 'Theta MPF'); |
|---|
| 35 | |
|---|
