%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%  rotmmf.m - simulates rotor produced air gap traveling mmf
%             wave for one magnetic pole-pair pitch of movement
%             in 10 magnetic (electric) degree increments. Rotor
%             position is indicated by cage.
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear; colordef white;
p=input('How many poles    ');       % No. of poles
s=input('Rotor slip(per unit)   ');  % Rotor slip
nb=input('How many rotor bars    '); % No. of rotor bars 
theta=linspace(0,2*pi,90); n=length(theta);
rs=5; rr=4;                          % Bore & rotor radii
rgr=4.5;                             % Air gap center grid radius
Rs=rs*ones(1,n);Rgr=rgr*ones(1,n);Rr=rr*ones(1,n); cla;
F=zeros(1,n); phi=0; rb=0.96*rr; rotbar=[]; rotang=[]; 
for k=1:nb; rotbar=[rotbar;rb]; rotang=[rotang;0]; end 
for i=1:37; pause(0.001);
   for j=1:n; F(j)=rgr+0.45*sin(p/2*theta(j)-phi); end 
ang=(1-s)*2/p*phi; rotang=[];
   for k=1:nb; rotang=[rotang;ang+(k-1)*2*pi/nb]; end
polar(theta,Rs,'-.'); hold on;
polar(theta,Rgr,'m');polar(theta,Rr,'g-.');
polar(rotang,rotbar,'go'); polar(theta,F,'r.');
text(1.6,rs+0.45,'Stator bore');text(-1.6,rr-0.9,'Rotor');
name=['Rotor produced mmf wave & rotor position';
'                                        ';
'                                        ']; title(name);
hold off; phi=phi+pi/18;
end;