Modelado maquina de cd en matlab

Proyecto Maquinas CD
Objetivo: De las ecuaciones de las maquinas de CD, calcular la corriente de
campo (Ia), la corriente de armadura (If), el par electromecánico (Te), la velocidad angular (Wm) yel voltaje inducido (eind). A través de un método numérico; así como graficar cada una de las variables.

Introducción
Maquinas de CD Son maquinas que convierten energía mecánica a eléctrica yeléctrica a mecánica. Las ecuaciones necesarias para el análisis de las maquinas de cd son las siguientes:

1.-

2.3.4.De 1 despejamos la derivada y reacomodando en forma de una matriz, quedando lasecuaciones de la sig. manera:

a)

Sustituyendo valores en la matriz

Ampers

Inicialmente se tiene If = Ia= w=0  Al calcular nuestras nuevas Ia e If; podemos observar que nuestra corriente dearranque es muy alta y la corriente de campo es aproximadamente 0 Sustituimos las corrientes calculadas en la ecuación del par = 8.4639e-005 rpm

 Nuestro par es aproximadamente de 0 Calculando lavelocidad angular tenemos:

=  La velocidad en la primera iteración es igual al par por lo tanto es aproximadamente de 0

Aplicando el Runge Kutta de 4to orden a las ecuaciones tenemos:
%datosmotor de CD Vf Va Rf Ra Lf La J K b If Ia W eind h t(1) = 230; = 230; = 10; = 0.4; = 120; = 0.02; = 10; = 0.0368; = 1; = 0; = 0; = 0; = 0; = .01; = 0; %voltaje de campo %voltaje de armadura%resistencia de campo %resistencia de armadura %inductancia de campo %inductancia de armadura %momento inercia %constante proporcionalidad %coeficiente de amortiguación %corriente de campo %corriente dearmadura %velocidad angular %voltaje inducido %paso de iteración

%Declarando Matriz V = [Va/La;Vf/Lf]; B = [Ra/La,(K/La)*W;0,Rf/Lf]; I = [Ia;If]; y = I; for i=1:1:100/h; %aplicando el método runge Kutta alas ecuaciones de la corriente K1 = V-[Ra/La,(K/La)*W;0,Rf/Lf]*y; a1 = y+(h/2)*K1; %en la primera iteración se calcula las Ia' e If' con losK2 = V-[Ra/La,(K/La)*W;0,Rf/Lf]*a1; %-valores iníciales...