Changeset 352 for applications/pmsm/Barcelona/Mz_mx.m

Timestamp:

06/02/09 10:04:12 (15 years ago)

Author:

smidl

Message:

Barcelona

Files:

• applications/pmsm/Barcelona/Mz_mx.m (modified) (1 diff)

applications/pmsm/Barcelona/Mz_mx.m

@@ -1 @@
-M=pmsm_estim_mex('Mz_mx.cfg');
+if 0
+    M=pmsm_estim_mex('Mz_mx.cfg');
+    save Mz M
+else
+    load Mz
+end
+
+time = 0.000125*[0:size(M.crbzth,1)-1];
 
 figure(1)
 subplot(2,1,1);
-plot(M.x(:,3));
+plot(time,M.xom);
+grid on
+title('Simulated rotor speed profile');
 
 subplot(2,1,2);
 hold off
-plot(M.crbz(:,4));
+plot(time,M.crbzth);
 hold on
-plot(M.crb(:,4),'r');
+plot(time,M.crbth,'r--');
+grid on
+legend('EKF TL','EKF')
+title('Cramer-Rao bounds on rotor position')
 
 figure(2)
 hold off
-plot(M.crbz(:,3));
+plot(time,M.crbzom);
 hold on
-plot(M.crb(:,3),'r');
+plot(time,M.crbom,'r');
 
-ind=100:size(M.crb,1)-100;
-[min(M.crb(ind,3)), max(M.crb(ind,3))]
-[min(M.crbz(ind,3)), max(M.crbz(ind,3))]
+ind=100:size(M.crbth,1)-100;
+[min(M.crbom(ind)), max(M.crbom(ind))]
+[min(M.crbzom(ind)), max(M.crbzom(ind))]
+