Ecuaciones diferenciales en busca de estabilidad

Alberdi Celaya, Elisabete

Matematikan lizentziatua eta doktoregaia

La estabilidad es un tesoro preciado. La vida es una búsqueda constante de estabilidad: búsqueda de un trabajo fijo, de una relación fuerte y larga, etc. Aunque sean variables, se suele denominar estable a aquel que se encuentra dentro de parámetros que no presentan cambios bruscos. Sinónimo de aburrimiento para algunos, un tesoro deseado en general.
Ecuaciones diferenciales en busca de estabilidad
01/11/2009 | Alberdi Celaya, Elisabete | Matematikan licenciado eta doctorando

(Foto: ©Fotolia)
También en matemáticas hay disciplinas en las que la estabilidad es fundamental. Los que queremos fijarnos tienen cabida en los programas de las carreras científicas y técnicas, se utilizan en la formulación de modelos matemáticos y para algunos se busca un resultado analítico, mientras que otros se resuelven con métodos numéricos --aproximaciones - y con el ordenador como acompañante. A algunos alumnos les gusta, pero a la mayoría les causan quebraderos de cabeza. ¿Qué son? Son ecuaciones diferenciales. Las ecuaciones diferenciales, al igual que el ser humano en la vida, buscan el mayor campo de estabilidad posible.

Utilidad de las ecuaciones diferenciales

Un modelo matemático es un dispositivo que describe un sistema o suceso de la vida. La formulación de un modelo matemático comienza con la identificación de las variables que influyen en el sistema, es decir, producen un cambio en el sistema. A continuación se establecen hipótesis razonables sobre el sistema, identificando las leyes empíricas aplicables. Algunas de estas hipótesis indican la medida de la variación de algunas de las variables previamente definidas. El enunciado matemático de estas hipótesis será una ecuación o un sistema de ecuaciones en el que aparezcan las derivadas, que es el sistema de ecuaciones diferenciales.

Su resolución concreta nos permitirá conocer el comportamiento del sistema. Pero, ¿de qué hablamos cuando hablamos de modelo o modelo matemático?

En la cinética de las reacciones químicas interesa su evolución a lo largo del tiempo. Al ser las velocidades derivadas del tiempo de cualquier variable, la cinética de las reacciones se modeliza mediante ecuaciones diferenciales. Ejemplo de ello es el caso de que ambas sustancias generan un tercero. Las variables son concentraciones de sustancias, mientras que las leyes empíricas son la ley de acción de masas y la ley de conservación de masas. La primera nos dice que el producto de las concentraciones de los reactivos es proporcional al producto de las concentraciones de los productos, y las segundas, la suma de las masas de los reactivos es igual a la de los productos.

En la refrigeración de los cuerpos también se utilizan ecuaciones diferenciales. Por ejemplo, si sacamos un pastel del horno a 150ºC y queremos saber su temperatura en cualquier momento. Si estamos investigando el homicidio o el asesinato de una persona y queremos calcular a qué hora murió, utilizaremos la ley de refrigeración de Newton. Esta ley establece que el cambio de temperatura de la superficie corporal es proporcional a la diferencia entre la temperatura corporal y la temperatura ambiente. Así, si la temperatura del cuerpo en el instante t es y ( t ) y la temperatura en el medio T, según la ley de Newton, se cumplirá la siguiente ecuación diferencial:

siendo k 0 la proporcionalidad constante.

¿Queremos calcular los intereses que nos dará una cantidad de dinero que ingresamos en el banco? Deberemos utilizar ecuaciones diferenciales. Supongamos que introducimos en el banco la cantidad D 0 y nos pagan un tipo de interés r. Llamaremos D ( t ) a la cantidad que tendremos dentro de t años. La variación de la cantidad será la suma de la variación por acumulación de intereses y la variación de ingresos que realizamos en el banco, p de fórmula.

Hay, por tanto, que modelar en la vida cotidiana utilizando ecuaciones diferenciales.

Dificultades de resolución

Distribución de los autovalores en un plano complejo. A la derecha de la línea verde, la ecuación diferencial es inestable y a la izquierda estable. La zona rosa es la zona de estabilidad de un método. Las estrellas rosas se encuentran en la zona de estabilidad, mientras que las verdes, multiplicadas por un número, las trasladamos a la zona (amarillas). Graf. : Elisabete Alberdi.

Una vez obtenida la ecuación diferencial se debe resolver. En ocasiones, no es posible su liberación analítica, y mediante métodos numéricos se obtiene un resultado aproximado. Las ecuaciones diferenciales incluyen unos autovalores complejos que denominaremos l. Cuando la parte real de estos autovalores es positiva, la ecuación diferencial es inestable y el resultado aproximado que obtendríamos utilizando un método numérico para su liberación no tiene por qué parecerse a la solución real, ya que un pequeño cambio tiene un fuerte impacto en ellos. Por ello, los métodos numéricos se aplican a ecuaciones diferenciales estables, autobalioides de parte real negativa. Por otra parte, el propio método también tiene una zona de estabilidad. Y si se quiere conseguir un buen resultado, los autovalores deben situarse allí. Por tanto, los autovalores de las ecuaciones diferenciales estables no presentes en esta zona se trasladan a la zona de estabilidad, multiplicando los autovalores por un número h medido (0,1) que denominaremos medida del paso.

La medida del paso es fundamental para resolver la ecuación y, por supuesto, no puede ser cualquier número. En términos de fiabilidad, el método de resolución debe ser lo suficientemente pequeño como para mantenerse dentro de la zona de estabilidad y lo suficientemente grande como para dar el menor número de pasos posibles de cara al trabajo. ¿Cómo conseguir ese equilibrio?

Equilibrio del área de estabilidad

Para saber qué pasos se pueden dar en el método numérico hay que fijarse en el error. Al utilizar métodos numéricos se genera un error local en cada paso. Pero hay otro error que aparece en la longitud, que es la reproducción de los errores de cada paso, y que está formado por las potencias de un número llamado factor de amplificación. Para evitar el aumento del error, el factor de amplificación debe ser menor que la unidad. La zona que cumple esta condición es la conocida como zona de estabilidad del método. Pues bien, el autovalor l debe estar en la zona de estabilidad tras multiplicarse por h para que el factor de amplificación sea menor que la unidad. Pero no sólo eso, una vez que el producto l .h se encuentra en la zona de estabilidad, se ve si el número h debe ser menor, ya que cualquier número que sea menor que h que hemos encontrado hace posible que después de multiplicar el autovalor se lleve a la zona. De los números que permiten llevar el autovalor al campo, el número h que mantiene el error local dentro de una tolerancia es el que se elige como medida del paso. Esta medida de paso permite tener controlados los dos errores (local y longitudinal). Cuanto mayor sea este número, mayor será el número de pasos que podremos dar y más rápido será el resultado, un objetivo que el gran campo de estabilidad permite alcanzar.

Cuando la zona de estabilidad es elevada, se podrán utilizar medidas mayores del paso h para llevar los autovalores a la zona, es decir, se podrán dar pasos mayores a través del método. En las pequeñas áreas de estabilidad se deben dar pequeños pasos cuando los autovalores son grandes. Ed. : ©Fotolia.
Pero la medida de la zona de estabilidad también está condicionada. Depende del orden. El orden es el parámetro que se utiliza para medir la precisión del método y es más preciso cuanto más alto es el orden. Normalmente, a medida que aumenta el orden, la zona de estabilidad disminuye. Esto nos obliga a dar pasos menores, aunque el resultado sea más preciso.

El propio método también condiciona la medida de la zona de estabilidad. Los métodos implícitos suelen tener un área de estabilidad mayor que los explícitos, pero también presentan desventajas: la principal es que se requieren operaciones más o menos complejas que en los métodos explícitos. Por ello, cuando los autovalores no son muy grandes, se prefiere el método explícito, ya que las operaciones a realizar serán más sostenibles.

El sueño sería encontrar la máxima estabilidad en órdenes altos y métodos explícitos, pero tampoco en este campo de la ciencia hay rareza. Con el fin de conseguir que las ecuaciones diferenciales se suelten fácilmente por métodos numéricos, son numerosos los trabajos realizados para aumentar el campo de estabilidad. Los primeros son los métodos de Adams Bashforth y Moulton, en los que se utiliza tanto la información del último paso como la información de otros pasos más rápidos para construir el siguiente paso. De esta forma se consiguió subir el orden del método numérico, y el de orden 1 tiene un campo mayor que el de Euler explícito (también de orden 1). Una de las propuestas más importantes fue la denominada BDF (fórmula backward differentiation) realizada por Gear hacia 1971, un método implícito que utilizaba información en varios pasos más ligeros. El método BDF consiguió que las zonas de estabilidad fueran también importantes en órdenes elevados. Últimamente predominan los métodos que utilizan la derivada 2 o los llamados puntos de superfuturo, ya que tienen grandes áreas de estabilidad, aunque incrementan el trabajo a realizar en cada paso.

A la izquierda, Adams Bashforth-Moulton PECE (Predict Evaluate Correctate Evaluate), campo de estabilidad del método explícito para diferentes órdenes. En este caso, el área de estabilidad es el interior de las líneas. A la derecha, la zona de estabilidad de los métodos implícitos denominados FMI para diferentes órdenes, siendo la zona exterior de las líneas la zona de estabilidad. Para indicar el orden se ha utilizado la letra k. En ambos casos, a medida que aumenta el orden, la zona de estabilidad disminuye. Graf. : Elisabete Alberdi.

Muchos pequeños pasos o pocos grandes pasos, ahí está la clave. En esta competencia, la primera opción, con muchos pequeños pasos, tiene la ventaja de la fiabilidad y la segunda la optimización del trabajo. Unificar ambas opciones en una sola ocasión sería excelente, ya que se conseguiría "fiabilidad sin trabajo excesivo". La clave para romper este equilibrio radica en métodos de gran campo de estabilidad, que no ponen límites tan rígidos a la medida del paso. Por eso es posible dar pasos no tan fiables y grandes pasos fiables.

Ecuaciones diferenciales y métodos numéricos en historia
Historia de las ecuaciones diferenciales XVII. Comenzó a finales del siglo XX con obras de Isaac Newton (1642-1727) y Gottfried Wilhelm Leibniz (1646-1716). Gracias al desarrollo de la mecánica, la física y el electromagnetismo, las ecuaciones diferenciales se convirtieron en una herramienta muy útil para describir los fenómenos naturales. Pero XIX. Hasta el siglo XIX la cuestión que los científicos manejaron fue buscar soluciones explícitas a estas ecuaciones diferenciales. En esa época se encuentran Leonhard Paul Euler (1707-1783), Daniel Bernoulli (1700-1782), Jean le Rond d'Alembert (1717-1783), Joseph Louis Lagrange (1736-1813) y Pierre Simon de Laplace (1749-1827).
Dos pasos del método de Euler explícito e implícito. En el explícito se da el paso moviéndose por la recta tangente del punto anterior (t1,y1) para obtener el punto (t2,y2). Implícitamente, se busca el punto t=t2 en directo (t2,y2), de manera que la recta tangente (t1,y1) pase por el punto. Graf. : Elisabete Alberdi.
En los casos más simples, el resultado se obtenía integrando la ecuación diferencial, pero no siempre era fácil o no era posible encontrar la solución analítica. Por ello, en 1760, Euler fue quien propuso la primera aproximación numérica para liberar ecuaciones diferenciales. La idea en la que se basaba este método numérico era muy sencilla: el intervalo en el que se debía encontrar el resultado se divide en subgrupos de tamaño h, siendo h el tamaño del paso. Para dar el paso 1 se conoce el valor inicial del primer tramo y se calcula el valor final del tramo moviendo por la recta tangente del punto inicial, es decir: y n+1 = y n + hf ( t n , y n ) .
Zonas de estabilidad de los métodos de Euler pintadas de rosa. Se observa que el área de estabilidad del implícito es mayor que la del explícito. Graf. : Elisabete Alberdi.
Para dar el paso 2 se conoce el valor inicial del tramo 2 y se repite el procedimiento anterior. Así en todos los pasos. Este método es explícito ya que utiliza los datos del paso anterior para encontrar el siguiente valor. Sin embargo, dado que el campo de estabilidad de este método es pequeño, los autovalores que no estaban en la zona debían multiplicarse por un número muy pequeño h para que entraran en la zona. Por ello, un cambio en el método es el conocido como método implícito de Euler, en el que el punto final del rango se calcula de forma que la recta tangente pase por su punto inicial intermedio: y n+1 = y n + hf ( t n+1 , y n+1 ). Este método es implícito porque en el cálculo de un desconocido aparece el mismo desconocido.
Resultado de la ecuación diferencial de valor inicial y' = -50(y-cost), y(0) = 0: en azul, el resultado analítico, y los resultados obtenidos mediante la aplicación de diferentes pasos al método implícito o explícito de Euler en los colores indicados. El método implícito de Euler, gracias a su mayor área de estabilidad, requiere pocos pasos para que el resultado aproximado sea similar al resultado analítico. GraF. : Elisabete Alberdi.
En ambos casos, el resultado se construye paso a paso y la posibilidad de dar pasos mayores permite obtener resultados más rápidos. La ventaja del método implícito de Euler es que tiene una mayor zona de estabilidad, ya que se pueden utilizar números mayores para llevar los autovalores a la zona.
Alberdi Celaya, Elisabete
Servicios
258
2009
Resultados
022
Matemáticas
Libre
Formación
Babesleak
Eusko Jaurlaritzako Industria, Merkataritza eta Turismo Saila