%   REACDAT
%    "motor_induccion"
%   Este script se ha de ejecutar junto con el modelo motor.mdl

%   $Revision: 0.1 
%   M. I. Capel 14-04-22

% -------------------------------
% Resistencia ohmica del 
% bobinado del rotor
% -------------------------------
Rr = .39;
%[Rr]= Ohmios
% Resistencia ohmica del estator
Rs = .19;
% Inductancia del estator
Lls= .21e-3;
%[Lls]= Henrio = voltio/(Amperio/s.)
% Inductancia del rotor
Llr = .6e-3;
%[Llr]= Henrio = voltio/(Amperio/s.)
% Inductancia de magnetizacion 
Lm = 4e-3;
%[Lm]= Henrio = voltio/(Amperio/s.)
%Frecuencia basica de las fases del motor de induccion
fb=100;
%[fb]= 1/s.
% Numero de polos que producen el campo magnetico giratorio
p= 4;
%Momento de inercia 
J= 0.0226;
%[J]= Newton x metro x (s.^2/rad)
%--------------------------------------
%Inductacia del rotor
Lr= Llr+Lm;
%constante de tiempo del rotor
Tr= Lr/Rr;
%[Tr]= s.
%---------------------------------------
% Calculo de la impedancia y de la velocidad angular
%Velocidad correspondiente a la frecuencia basica
wb=2*pi*fb;
%Impedancia del estator
Xls=wb*Lls;
%Impedancia del rotor
Xlr=wb*Llr;
%Impedancia de magnetizacion
Xm=wb*Lm;
%Impedancia total equivalente
Xmstar=1/(1/Xls+1/Xm+1/Xlr);
%Constate para integrar we
dospi= 2*pi;
%Constante de carga 1000ms.
tau= 1;
YTT=timeseries2timetable(rotor_speed);
XTT=timeseries2timetable(currents);
ZTT=timeseries2timetable(voltages);
YT1=[seconds(YTT.Time),YTT.Data];
XT1=[seconds(XTT.Time),XTT.Data];
ZT1=[seconds(ZTT.Time),ZTT.Data];
WT1=[seconds(ZTT.Time),XTT.Data,ZTT.Data];
[T,w]=tonndata(YT1,false,false);
[W,w]=tonndata(WT1,false,false);


