%Three phase induction motor

%Three phase induction motor
r1=0.200;
x1=0.250;
r2=0.10;
x2=0.250;
xm=25.0;
v_phase=208/sqrt(3);
nsync=1800;
wsync=2*pi*nsync/60;
vth=v_phase*(xm/sqrt(r1^2 +(x1+xm)^2));
zth=((j*xm)*(r1+j*x1))/(r1+j*(x1+xm));
rth=real(zth);
xth=imag(zth);
s=(0:1:50)/50;
s(1)=0.001;
nm=(1-s)*nsync;
for i = 1:51
td(i)=(3*vth^2*r2/s(i))/(wsync*((rth+r2/s(i))^2 +(xth+x2)^2));
end
figure(1);
plot(nm, td);
grid on;
xlabel(‘n(r/m)’)
ylabel(‘Td(Voltage Ratio (p.u)’)