%CurvaParVelocidad%

disp('A continuación se le va a pedir que ingrese los datos del modelo equivalente del MI, versión estrella equivalente')
R1 = input('Ingrese la resistenica estatórica: ');
R2e = input('Ingrese la resistenica rotórica vista del estator: ');
X = input('Ingrese la reactancia total de fugas "X1+X2e": ');
%Xo = input('Ingrese la reactancia magnetizante "Xo": ');
%Ro = input('Ingrese la resistencia de vacio "Ro": ');
U = input('Ingrese la tensión de linea de la red: ');
p = input('Ingrese la cantidad de pares de polos del MI: ');
f = 50;
g =[.001:.001:0.5];
ns = 60*f/p;
n = ns.*(1-g);
ws = 2*pi*ns/60;
Cm = ((R2e/ws)./g).*((U^2)./((R1+R2e./g).^2+(X^2)));
subplot(2,1,1)
plot(n,Cm) 
title('Curva par velocidad'); 
ylabel('C(N.m.)'); 
xlabel('N(rpm)');
subplot(2,1,2) 
plot(g*100,Cm); 
ylabel('C(N.m.)');
xlabel('g(%)')