|
Revision 47, 0.8 kB
(checked in by smidl, 18 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 | |
|---|