%E316p6_10.m
%rjb 3/24/03
%code for prob 6.10 
clear
format compact
%Parameters
R1 = 0.103; R2 = 0.225; X1 = 1.10; X2 = 1.13; Xm = 59.4;
Pfw = 265; Pcore = 220;
V = 460; f = 60; p =4; 

%%%%%%%%%%%%%%% Prob 6-10 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%loop over slip
for slip =1:3
    if slip == 1
        s =.01;
        disp('%%%%%%%%%%%%%%%%%%%%%% s = .01  %%%%%%%%%%%%%%%%%%%%%%%')
    elseif slip == 2
        s =.02;
        disp('%%%%%%%%%%%%%%%%%%%%%% s = .02  %%%%%%%%%%%%%%%%%%%%%%%')
    elseif slip == 3
         s=.03;
        disp('%%%%%%%%%%%%%%%%%%%%%% s = .03  %%%%%%%%%%%%%%%%%%%%%%%')
    end

%Calculations
    Zt = R1+j*X1 + ((j*Xm*((R2/s) + j*X2))/((R2/s) +j*(Xm+X2)));
    pf = cos(angle(Zt))
    I1 = V/(sqrt(3)*Zt);
    I2 = I1 * j*Xm / ((R2/s)+j*(X2 +Xm));
    Pdev = 3  *(abs(I2)^2)* R2 * ((1-s)/s);
    Pout = Pdev - Pcore - Pfw
    Psta = 3  *(abs(I1)^2)* R1;
    Prot = 3  *(abs(I2)^2)* R2;
    Pin = sqrt(3)*V*abs(I1)*pf;
    Pin2 = Pdev + Psta + Prot;
    Eff = Pout/Pin 
    ns = 120*f/p; ws = ns*pi/30;
    wr = (1-s)*ws; nr = (1-s)*ns
    Tout = Pout/wr
end