UNIVERSIDAD POLITECNICA SALESIANA

ECUACIONES DIFERENCIALES

RESOLUCION DE UN CIRCUITO RLC

CUARTO NIVEL G-1

INTEGRANTES:

PAOLA FLOR

ALEXANDRA ROJAS

JOSE LUIS TAFUR

ISRAEL VERGARA

MATLAB contiene dos funciones para calcular soluciones numérica de ecuaciones diferenciales

ordinarias; "ode23" y "ode45".

Veremos como calcular funciones con ode45.

[x,y] = ode45('función', a, b ,inicial) 

Esta instrucción regresa un conjunto de coordenadas "x" y "y" que representan a la función y=f(x), los valores se calculan a través de métodos Runge-Kuta de cuarto y quinto orden.El nombre "función", define una función que representa a una ecuación diferencial ordinaria, ODE45 proporciona los valores de la ecuación diferencial y'=g(x,y).Los valores "a" y "b" especifican los extremos del intervalo en el cual se desea evaluar a la función y=f(x).El valor inicial y = f(a) especifica el valor de la función en el extremo izquierdo del intervalo [a,b].

CIRCUITO RLC

FUNCIÓN:

ECUACION:

PASOS PARA RESOLVER UN CIRCUITO RLC

Primero ingresamos la funcion de esta manera Function B=cirlcr(t,A)

Despues creamos una matriz(2x1) B=zeros(2,1)

Asignamos valores de carga

B(1)=A(2);

Ingresamos la funcion B(2)=23*sin(t)-3*B(1)-A(1)/6.5;

NOS QUEDA DE LA SIGUIENTE MANERA

Despues de esto guardamos el archivo con el nombre por default que nos da Matlab.

El siguiente codigo lo pegamos en la ventana command window:

[t,A]=ode45('cirlcr',[0 10],[0 0]); q=A(:,1); i=A(:,2); plot(t,q) title('q vs t') xlabel('t(s)') ylabel(‘q(C)') figure (2) plot(t,i) title('i vs t') xlabel('t(s)') ylabel('i(A)')

NOS QUEDA DE LA SIGUIENTE MANERA

GRAFICA 1 Q(C) VS T(S)

GRAFICA 2 I(A) VS T(S)