Revision 47, 0.8 kB
(checked in by smidl, 17 years ago)
|
|
-
Property svn:eol-style set to
native
-
Property svn:keywords set to
Rev Author Date
|
Line | |
---|
1 | itload('../pmsm.it') |
---|
2 | figure(1); |
---|
3 | plot(xthM(1,:)') |
---|
4 | set(gca,'YLim',[0,0.4]) |
---|
5 | hold on |
---|
6 | plot([0,10000],0.1*[1,1],'--') |
---|
7 | legend('\hat{q}_t using EKF-PF','simulated value') |
---|
8 | legend('expected value of q_t using EKF-PF','simulated value') |
---|
9 | xlabel('time index') |
---|
10 | %set(gcf,'PaperPositionMode','auto'); |
---|
11 | %print -deps2 ekf |
---|
12 | |
---|
13 | figure(2); |
---|
14 | subplot(2,1,1) |
---|
15 | plot(xth(3,:)','--') |
---|
16 | hold on |
---|
17 | plot(xthE(3,:)',':') |
---|
18 | plot(xthM(5,:)') |
---|
19 | legend('Simulated \omega','Estimated \omega using EKF ', 'Estimated \omega using EKF-MPF '); |
---|
20 | xlabel('time index') |
---|
21 | %set(gca,'XLim',[2300,3000]); |
---|
22 | |
---|
23 | subplot(2,1,2) |
---|
24 | hold off |
---|
25 | plot(xth(4,:)','--') |
---|
26 | hold on |
---|
27 | plot(xthE(4,:)',':') |
---|
28 | plot(xthM(6,:)') |
---|
29 | legend('Simulated \upsilon','Estimated \upsilon using EKF ', 'Estimated \upsilon using EKF-MPF '); |
---|
30 | xlabel('time index') |
---|
31 | %set(gca,'XLim',[2300,3000]); |
---|
32 | |
---|