Algoritmo para modelar la máquina de inducción (secuencia qd0 en función de flujos)

Algoritmo para modelar la máquina de inducción
(secuencia qd0 en función de ()
Por:
Arturo Chávez González

Referencia del capítulo 4 del libro: “Analysis of Electric Machinery and DriveSystems”, Paul C. Krause, Oleg Wasynczuk and Scott D. Sudohoff.

Secuencia “qd0”.

Dar de alta los datos de la máquina a modelar.

f=60; %Frecuencia de alimentación
w=2*pi*f;%Velocidad angular
v= 220; %Voltaje linea-neutro
rs=0.435; %Resistencia de los devanados del estator
rr=0.816; %Resistencia de los devanados del rotor
xls=0.754; %Reactancia dedispersion del estator
xlr=0.754; %Reactancia de dispersion del rotor
xm=26.13; %Reactancia de magnetizacion
J=0.089; %Inercia del rotor
P=4; %Numero de polos

Calcular:v= sqrt(2)*vl_l/sqrt(3); %Voltaje linea-neutro
w=2*pi*f; %Frecuencia angular
Lls=xls/w;
Llr=xlr/w;
lms=xm/w;
LM=(3/2)*lms;

Declarar variables que no dependen de los estados nidel tiempo.

alf=(2*pi)/3;
RS=diag([rs rs rs]); %Matriz de resistencias del estator
RR=diag([rr rr rs]); %Matriz de resistencias del rotor
Ls= diag( [Lls+LM Lls+LM Lls]);%Matriz deinductancias del estator
Lr= diag( [Llr+LM Llr+LM Llr]); %Matriz de inductancias del rotor
Lsr= diag( [LM LM 0]); %Matriz de inductancias mutuas
Lrs=transpose(Lsr);
R=[RS zeros(3);zeros(3)RR]; %Matriz de resistencias
L=[Ls Lsr;Lrs Lr]; %Matriz de inductancias
Plantear y resolver las ODE´s del modelo con condiciones iniciales cero en un intervalo de tiempo adecuado (de 0 a 1seg. por ejemplo).

Iniciar el ciclo de la ODE con t=0 plantear el vector de estados con sus valores iniciales (cero) y los vectores de las variables restantes que dependen de los estados o deltiempo.

Y= [0 0 0 0 0 0 0];

Ks=(2/3)*[cos(w*t) cos(w*t-alf) cos(w*t+alf);
sin(w*t) sin(w*t-alf) sin(w*t+alf);
1/2 1/2 1/2 ];

om1=[0 w 0;-w 0 0;0 0 0];...