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 | |
---|