1.
Realice una investigación acerca de los fundamentos de las Ecuaciones diferenciales
ordinarias y sus aplicaciones del Mathlab
Ecuaciones Diferenciales Ordinarias (ODES)
Matlab posee multiples rutinas para resolver ODES, entre ellas se
encuentran ode45 y ode23tb. Ode45 es un método de solución de ecuaciones
diferenciales ordinarias y su sintaxis es:
[t,x]=ode45(@función,tinit,xinit)
Donde @función es la función que contiene la ecuación diferencial a resolver,
tinit es el intervalo de tiempo sobre el que se resuelve la ecuación y xinit es un
vector que actúa como condiciones iniciales para el problema. Las salidas t y x
corresponden, respectivamente, a los valores del tiempo y la solución de la
función.
Es importante recordar que la solución de una ecuación diferencial es una
función del tiempo, es decir, una trayectoria tal que para cada valor de t, existe
un valor de la función.
Con el fin de configurar más salidas, es posible usar options. Options es
configurada por medio del comando odeset, que genera salidas tales como la
jacobiana, gráficos y estadísticas. Un ejemplo de esto es:
options=odeset('outputfcn','odeplot')
[t,x]=ode45(@función,tinit,xinit, options)
Options hará que se arroje un gráfico de la solución al problema. Las opciones
son muy amplias y la explicación sobre ellas en la página oficial de Matlab es
detallada.
A continuación, se muestra un ejemplo del uso de Ode45 para resolver un
problema:
En un script escribimos la función:
function dxdt=deq(t,x)
dxdt=3-x
end
En un segundo script:
tinit=[0 20]
xinit=[-1.5];
options=odeset('outputfcn','odeplot')
[t,x]=ode45(@deq,tinit,xinit,options)
La solución a los problemas de ecuaciones diferenciales varía con respecto a
las condiciones iniciales y a los parámetros que posea la función. Por esto, a
continuación, se muestra una forma simple de incluir un cambio en los
parámetros en la solución. Para las condiciones iniciales el proceso es similar,
por lo que no se explica.
En un primer script definimos la función:
Sistemas de Ecuaciones diferenciales Ordinarias
Anteriormente vimos el uso de Ode45 para resolver una ecuación diferencial de
primer orden. Sin embargo, esta rutina permite dar solución a ecuaciones hasta
de orden n. Para esto, se debe transformar la ecuación en un sistema de
ecuaciones ordinarias de primer orden y luego implementarlo en Matlab. A
continuación, se muestra un ejemplo de cómo hacer esto:
Dada la ecuación:
Se sustituye:
Note que son tres ecuaciones porque la variable y se encontraba derivada tres
veces. Para implementar este sistema en Matlabs se crean tres ecuaciones con
los términos del sistema dentro de una función y agrega una más que
representa un vector de tres filas por una columna, como se muestra a
continuación:
En un nuevo script:
function dydt=sist(t,x)
dydt=zeros(3,1)
dydt=x(1)
dydt=x(2)
dydt=1-4^3 - 2*x(3) + x(2) + 2*x(1)
end
En general, el sistema tiene tantas ecuaciones como se encuentre derivada la
variable de interés y el tamaño del vector de ceros es también determinado por
esta cantidad. Para darle solución, se sigue el mismo procedimiento que con
las ODES. Sin embargo, hay que tener cuidado al defini xinit, pues esta debe
tener tantas condiciones iniciales como ecuaciones tiene el sistema.
En un segundo script:
tinit=0:0.5:10
xinit=[1 0 0]
[t,x]=ode45(@sist,tinit,xinit)
