%%
%--------------------------------------------------------------------------
%--------------------------------------------------------------------------
%
%
% Modeling and simulation of a 4 phase variable reluctance motor
%
%
%--------------------------------------------------------------------------
%--------------------------------------------------------------------------
clc;
clear all;
close all;
dirsig=-1;%clockwise rotation
sim('VRM')
figure (1)
subplot(2,2,1);plot(t(1:1000),I1(1:1000),'b','linewidth',2);grid on;ylabel('I1');xlabel('Time(sec)');axis([0 0.1 0 8])
subplot(2,2,2);plot(t(1:1000),I2(1:1000),'b','linewidth',2);grid on;ylabel('I2');xlabel('Time(sec)');axis([0 0.1 0 8])
subplot(2,2,3);plot(t(1:1000),I3(1:1000),'b','linewidth',2);grid on;ylabel('I3');xlabel('Time(sec)');axis([0 0.1 0 8])
subplot(2,2,4);plot(t(1:1000),I4(1:1000),'b','linewidth',2);grid on;ylabel('I4');xlabel('Time(sec)');axis([0 0.1 0 8])
figure (2)
plot(t,n_rpm,'b','linewidth',2);grid on;axis([0 12 -8000 0])
dirsig=1;%clockwise rotation
sim('VRM')
figure (3)
subplot(2,2,1);plot(t(1:1000),I1(1:1000),'b','linewidth',2);grid on;ylabel('I1');xlabel('Time(sec)');axis([0 0.1 0 8])
subplot(2,2,2);plot(t(1:1000),I2(1:1000),'b','linewidth',2);grid on;ylabel('I2');xlabel('Time(sec)');axis([0 0.1 0 8])
subplot(2,2,3);plot(t(1:1000),I3(1:1000),'b','linewidth',2);grid on;ylabel('I3');xlabel('Time(sec)');axis([0 0.1 0 8])
subplot(2,2,4);plot(t(1:1000),I4(1:1000),'b','linewidth',2);grid on;ylabel('I4');xlabel('Time(sec)');axis([0 0.1 0 8])
figure (4)
plot(t,n_rpm,'b','linewidth',2);grid on;axis([0 12 0 8000])
%--------------------------------------------------------------------------
%--------------------------------------------------------------------------