1 | load graf1.txt -ascii
|
---|
2 |
|
---|
3 | % t, Isx, Isy, Isd, Isdw, Isq, Isqw, M, Fs, poloha, rychlost, rychlost_w,
|
---|
4 | % Urm, Urm_max, beta, zbeta, Uc
|
---|
5 |
|
---|
6 | figure
|
---|
7 | plot(graf1(:,1),graf1(:,2)) % Isx
|
---|
8 | % hold on
|
---|
9 | % plot(graf1(:,1),graf1(:,3),'r')
|
---|
10 | legend('Isx');
|
---|
11 |
|
---|
12 | figure
|
---|
13 | plot(graf1(:,1),graf1(:,4)) % Isd
|
---|
14 | hold on
|
---|
15 | plot(graf1(:,1),graf1(:,5),'r')
|
---|
16 | legend('Isd','Isdw');
|
---|
17 |
|
---|
18 | figure
|
---|
19 | plot(graf1(:,1),graf1(:,6)) % Isq
|
---|
20 | hold on
|
---|
21 | plot(graf1(:,1),graf1(:,7),'r')
|
---|
22 | legend('Isq','Isqw');
|
---|
23 |
|
---|
24 | figure
|
---|
25 | plot(graf1(:,1),graf1(:,8)) % M
|
---|
26 | hold on
|
---|
27 | plot(graf1(:,1),graf1(:,16),'r') % M_model
|
---|
28 | legend('Moment','Moment model');
|
---|
29 |
|
---|
30 | figure
|
---|
31 | plot(graf1(:,1),graf1(:,9)) % Fs
|
---|
32 | hold on
|
---|
33 | plot(graf1(:,1),graf1(:,15),'r') % Fs_model
|
---|
34 | legend('Fs','Fs model');
|
---|
35 |
|
---|
36 | figure
|
---|
37 | plot(graf1(:,1),graf1(:,10)) % Poloha
|
---|
38 | title('Poloha theta elektricka');
|
---|
39 |
|
---|
40 | figure
|
---|
41 | plot(graf1(:,1),graf1(:,11)/(2*pi)) % fme
|
---|
42 | hold on
|
---|
43 | plot(graf1(:,1),graf1(:,12)/(2*pi),'r')
|
---|
44 | legend('fme','fmew');
|
---|
45 |
|
---|
46 | figure
|
---|
47 | plot(graf1(:,1),graf1(:,13)) % Urm
|
---|
48 | hold on
|
---|
49 | plot(graf1(:,1),graf1(:,14),'r')
|
---|
50 | legend('Urm','Urm max');
|
---|
51 |
|
---|
52 | figure
|
---|
53 | plot(graf1(:,1),graf1(:,17)/pi*180) % zatezny uhel (beta)
|
---|
54 | legend('zatezny uhel (beta)');
|
---|
55 |
|
---|
56 | figure
|
---|
57 | plot(graf1(:,1),graf1(:,18)/pi*180) % zatezny uhel (zbeta) - vypocet z napeti
|
---|
58 | legend('zatezny uhel (zbeta)');
|
---|
59 |
|
---|
60 | figure
|
---|
61 | plot(graf1(:,1),graf1(:,19)) % Uc
|
---|
62 | legend('Uc');
|
---|
63 |
|
---|
64 | figure
|
---|
65 | plot(graf1(:,1),graf1(:,20)) % Um*2/Ucn
|
---|
66 | legend('Um*2/Ucn');
|
---|
67 |
|
---|
68 | figure
|
---|
69 | plot(graf1(:,1),graf1(:,21)) % Upila
|
---|
70 | hold on
|
---|
71 | plot(graf1(:,1),graf1(:,22),'r') % Ua
|
---|
72 | legend('Upila','Ua');
|
---|
73 |
|
---|
74 | figure
|
---|
75 | plot(graf1(:,1),graf1(:,23)) % Ukomp
|
---|
76 | legend('Ukomp');
|
---|
77 |
|
---|
78 | figure
|
---|
79 | plot(graf1(:,1),graf1(:,24)) % ALGORITMUS
|
---|
80 | legend('ALGORITMUS');
|
---|
81 |
|
---|
82 | figure
|
---|
83 | plot(graf1(:,1),graf1(:,25)/2/pi) % fme_ekf
|
---|
84 | hold on
|
---|
85 | plot(graf1(:,1),graf1(:,11)/(2*pi),'r') % fme skutecna
|
---|
86 | legend('fme ekf','fme');
|
---|
87 |
|
---|
88 | figure
|
---|
89 | plot(graf1(:,1),graf1(:,26)) % theta_elektricka_ekf
|
---|
90 | hold on
|
---|
91 | plot(graf1(:,1),graf1(:,10),'r') % Poloha elektricka skutecna
|
---|
92 | legend('Theta ekf','Theta real');
|
---|
93 |
|
---|