ECUACIONES DIFERENCIALES EN MATLAB

ECUACIONES
DIFERENCIALES
EN

MATLAB
Daniel Parcero Sánchez
Rocío Salgueiro Fernández

Ecuaciones Diferenciales en Matlab®
Matlab ofrece varios algoritmos numéricos para resolver una extensa variedad de
ecuaciones diferenciales. Esta demostración enseña la formulación y solución para dos
tipos distintos de ecuaciones diferenciales usando Matlab.
-

El problema del valor inicial.
Problemas devalor limite.

a)El Problema del Valor Inicial
Para esto utilizaremos la ecuación de Van Der Pol, en Matlab VANDERPOLDEMO es
una función que define la ecuación de Van Der Pol.
La ecuación de Van Der Pol tiene una larga historia de utilización tanto en física como
en biología. Por ejemplo, en biología, Fitzhugh y Nagumo ampliaron la ecuación como
un modelo para la acción neuronal. En cuanto a lafísica, la ecuación ha sido utilizada en
sismología para modelizar las dos placas en una falla geológica.
En la dinámica, el oscilador de Van Der Pol es un oscilador no conservativo con
amortiguamiento no lineal Se desarrolla en el tiempo de acuerdo con la siguiente
ecuación diferencial de 2º orden:

Donde x es la coordenada de posición - que es función del tiempo t - y µ es un
parámetro escalar queindica la fuerza de la amortiguación. Se puede demostrar a través
del Teorema de Lienard que existe un ciclo límite para el oscilador de Van Der Pol no
forzado, por lo que es un ejemplo de un sistema de Lienard.

Resultados para el oscilador no forzado
Dos regímenes de interés para las características del oscilador no forzado son:
•

Cuando µ = 0, es decir, no hay ninguna función de amortiguación, laecuación
es:

Esta es una forma del oscilador armonico simple, y por tanto se conserva la
energia mecanica.

•

Cuando µ> 0, el sistema entrará en un ciclo límite, donde la energía sigue siendo
conservada, pero parte de la energía mecanica se perdera en forma de calor.
Cerca del origen x = dx / dt = 0 el sistema es inestable, y lejos del origen el
sistema se amortigua.

El oscilador forzado VanDer Pol
El oscilador de Van Der Pol forzado o dirigido, toma la función original, y añade una
función de corrección A.sen (ωt) para dar una ecuación diferencial de la forma:

Donde A es la amplitud, t el tiempo y ω es su velocidad angular.

type vanderpoldemo
dydt = [y(2); Mu*(1-y(1)^2)*y(2)-y(1)];
La ecuación se escribe como un sistema de dos funciones ODE de primer
orden. Estas son evaluadas paradistintos valores del parámetro Mu.
Para una integración más rápida, elegimos un método de solución basado
en el valor del parámetro Mu.

Para Mu = 1, cualquiera de los métodos de solución ODE Matlab puede resolver la
ecuación de Van Der Pol eficientemente. el método de solución ODE45 usado a
continuación es un ejemplo. La ecuación es resuelta en el dominio [0, 20].
- tspan: es el vector queespecifica el intervalo de integración.
– y0: es el vector de condiciones iniciales.
- Mu: es un parámetro, al que le damos el valor que queramos, 1 en este caso..
- ODE: es un método de solución que emplea Matlab para resolver ecuaciones
diferenciales ordinarias.
– [t,y]: Nos dice el método empleado por Matlab para la resolución del problema, y las
variables de las que depende (ode, tspan y y0), yaexplicadas.
tspan = [0, 20];
y0 = [2; 0];
Mu = 1;
ode = @(t,y) vanderpoldemo(t,y,Mu);
[t,y] = ode45(ode, tspan, y0);
%

Gráfico de la solución

- plot: Comando que ordena a Matlab representar graficamente la función que le damos.
- xlabel: Nos permite ponerle el titulo que queramos al eje x.
- ylabel: Nos permite ponerle el titulo que queramos al eje y.
- title: Nos permite ponerle a la gráfica eltitulo que queramos.

plot(t,y(:,1))
xlabel('t')
ylabel('Solución y')
title('Ecuación Van Der Pol, \mu = 1')

Para magnitudes más grandes de Mu, el problema se vuelve más rígido, es decir, mas
dificil de resolver numericamente. Para una integración rapida son necesarios métodos
numéricos especiales como ODE15S, ODE23S, ODE23T, y ODE23TB, que pueden
resolver problemas rígidos eficientemente....