pmsm_disp.m YLim YLim __pad2__ __pad2__ XLim XLim __pad3__ __pad3__ itload ('../pmsm.it') figure(1) itload '../pmsm.it' plot (xthM(1,:)') set(gca plot xthM(1,:)' hold on hold on plot ([0, 10000], 0.1 *[1, 1],'--') legend('\hat plot 0.1 * [1, 1] '--' 2 2 figure () figure subplot (2, 1, 1) plot(xth(3 subplot 2 1 1 xthE(3,: xthE(3,: plot: ('Simulated\omega','Estimated\omega using EKF ', 'Estimated\omega using EKF-MPF 'legend) plot: 'Simulated\omega' legend 'Estimated\omega using EKF ' 'Estimated\omega using EKF-MPF ' legend xlabel ('time index')%set(gca xlabel 'time index' subplot (2, 1, 2) hold off plot(xth(4 subplot 2 1 2 xthE(4,: xthE(4,: plot: ('Simulated\upsilon','Estimated\upsilon using EKF ', 'Estimated\upsilon using EKF-MPF 'legend) plot: 'Simulated\upsilon' legend 'Estimated\upsilon using EKF ' 'Estimated\upsilon using EKF-MPF ' legend itload('../pmsm.it') figure(1); plot(xthM(1,:)') set(gca,'YLim',[0,0.4]) holdon plot([0,10000],0.1*[1,1],'--') legend('\hat{q}_tusingEKF-PF','simulatedvalue') legend('expectedvalueofq_tusingEKF-PF','simulatedvalue') xlabel('timeindex') %set(gcf,'PaperPositionMode','auto'); %print-deps2ekf figure(2); subplot(2,1,1) plot(xth(3,:)','--') holdon plot(xthE(3,:)',':') plot(xthM(5,:)') legend('Simulated\omega','Estimated\omegausingEKF','Estimated\omegausingEKF-MPF'); xlabel('timeindex') %set(gca,'XLim',[2300,3000]); subplot(2,1,2) holdoff plot(xth(4,:)','--') holdon plot(xthE(4,:)',':') plot(xthM(6,:)') legend('Simulated\upsilon','Estimated\upsilonusingEKF','Estimated\upsilonusingEKF-MPF'); xlabel('timeindex') %set(gca,'XLim',[2300,3000]);